Effect-directed analysis for identification of river basin specific pollutants.
Ph.D. Dissertation
Zuzana Toušová
Masaryk University, Faculty of Science,
Centre RECETOX
Brno, Czech Republic 2020
Supervisor: Dr. Ing. Jaroslav Slobodník
2
BIBLIOGRAPHIC ENTRY
Author: Mgr. Zuzana Toušová Faculty of Science, Masaryk University Centre RECETOX
Title of Dissertation: Effect-directed analysis for identification of river basin specific pollutants
Degree Programme: Environmental Health Sciences
Specialization: Environmental Chemistry and Toxicology
Supervisor: Dr. Ing. Jaroslav Slobodník
Supervisor specialist: prof. RNDr. Luděk Bláha, Ph.D.
Year: 2020
Keywords: contaminants of emerging concern, effect directed analysis, in-vitro bioassays, prioritization, non- target screening.
3
BIBLIOGRAFICKÝ ZÁZNAM
Autor: Mgr. Zuzana Toušová Přírodovědecká fakulta, Masarykova univerzita Centrum RECETOX
Název práce: Užití efektem řízené analýzy pro identifikaci specifických znečišťujících látek v povodích
Studijní program: Životní prostředí a zdraví
Specializace: Environmentální chemie a toxikologie
Školitel: Dr. Ing. Jaroslav Slobodník
Školitel specialista: prof. RNDr. Luděk Bláha, Ph.D.
Rok: 2020
Klíčová slova: nově se objevující znečišťující látky, efektem řízená analýza, biotesty in-vitro, prioritizace, necílená analýza.
4
© Zuzana Toušová, Masaryk University, 2020 5
ACKNOWLEDGEMENTS
I would like to thank my supervisor Dr.. Jaroslav Slobodník for his great guidance, encouragement, support and patience during my postgraduate study. My deep gratitude also goes to prof. RNDr. Luděk Bláha, Ph.D., who carefully co- supervised my work, inspired me, provided me with great scientific training and contributed to the direction and richness of this research. Prof. Bláha was always ready to give advice and discuss immediate problems, which helped me greatly to progress smoothly with both laboratory work and writing papers. I am also very grateful to doc. Mgr. Klára Hilscherová, PhD., who largely contributed to the study design and results’ interpretation in publications II and III. Doc. Hilscherová kindly shared her rich practical experience with me and devoted a lot of her time to help me finalize the publications. My thanks also go to all co-authors of our common publications, for their great work on the research and writing papers. Namely, I would like to thank Dr. Tobias Schulze and doc. Ing. Branislav Vrana, PhD. for their great work on the sampling part of this research. I gratefully acknowledge Dr. Peter Oswald and Dr. Natalia Glowacka., colleagues from the EI, who contributed on the analytical part and helped me with lots of practical issues. I extend my thanks to all members of the EDA EMERGE Project for all those inspirational discussions, training courses, hospitality and fun. My thanks go to all RECETOX colleagues for making such a friendly environment at the workplace and being helpful at all times. I express my deep gratitude to all my family. I would like to thank my parents and grandparents for their huge support throughout my studies. Great thanks go to my husband for his encouragement, support, patience and time with the kids. I am very grateful to my mother and especially to my late mother- in-law for being such good grandmas and all the baby-sitting, without which this thesis would never had taken shape.
The PhD training programme was funded by the European Union under the Marie Curie Actions—Initial Training Networks, FP7-PEOPLE-2011-ITN, EDA-EMERGE project, Grant Agreement No. 290100.
6
ABSTRACT
Pollution of surface waters with contaminants of emerging concern (CECs) has become a pressing environmental issue of global importance. CECs is a group of organic compounds of diverse structure, physicochemical properties and usage patterns comprising pesticides, pharmaceuticals, personal care products, flame retardants and others. CECs can be found in surface and waste waters, they are not monitored or regulated and raise concerns of the research community and public due to their extensive occurrence and possible adverse effects on biota, which are mostly unknown. Innovative approaches to monitor surface waters, with regard to CECs, have been proposed. These include novel sampling techniques, effect assessment with sensitive bioassays, non-target chemical screening and use of advanced software tools for data interpretation. Further, a powerful technique to identify bioactive compounds in complex environmental mixtures, known as effect-directed analysis (EDA), is based on the combination of bioassays, fractionation and chemical analysis. The objective of this thesis was to contribute to the identification of relevant toxicity drivers in surface waters by advancing the EDA utilization. The specific objectives, addressed by particular studies, were to apply a newly developed simplified EDA protocol (study I and IV), identify compounds responsible for algal growth and acetylcholine esterase inhibition in a WWTP effluent extract (study II and V) and characterize contamination of the Bosna River with CECs (study III). Novel sampling device, based on solid phase extraction of large water volumes (LVSPE), was tested and optimized in a sampling campaign spanning across four European river basins (Studies I and IV). In the effect assessment of extracts, the most pronounced effects were estrogenicity, toxicity to algae and fish embryo toxicity, whereas the major portions of the observed effects could not be explained by analyzed compounds. Simplified risk assessment procedure with analyzed compounds enabled identification and prioritization of 21 candidate compounds for future monitoring efforts. Effect directed analysis was applied to identify compounds contained in a wastewater treatment plant (WWTP) effluent extract that were responsible for growth inhibition of green algae (Study II) and acetylcholinesterase inhibition in vitro (Study V). Our results suggest that pesticides and their transformation products, pharmaceuticals (barbiturate derivatives and macrolide antibiotics e.g. azithromycin), industrial compounds or caffeine and its metabolites were the most likely toxicity drivers for green algae in study II. 7
A combined chemical and effect screening of water quality in the River Bosna was carried out, to characterize, in some detail, its pollution with CECs (Study III). The assessment of cumulative pollutant concentrations and hazard profiles enabled to identify the major source of contamination with CECs and hotspots of biological potency. Simplified risk assessment procedure of detected target compounds suggested that 7 compounds, namely diazinon, diclofenac, 17-β estradiol, estrone, benzo[k]fluoranthene, fluoranthene and benzo[k]fluoranthene, might pose serious risks to aquatic biota in the Bosna River. The presented dissertation demonstrates that combination of novel sampling techniques with effect-based methods and chemical screening is a functional strategy to address CECs in surface waters. This thesis shows that EDA is a useful tool to reduce the complexity of environmental mixtures and to identify relevant toxicity drivers. Risk assessment and prioritization of detected chemicals is a crucial step in the process of identifying river basin specific pollutants.
8
ABSTRAKT
Znečištění povrchových vod nově se objevujícími znečišťujícími látkami (Contaminants of emerging concern - CECs) se stalo vážným environmentálním problémem celosvětového rozsahu. CECs jsou velmi rozmanitou skupinou organických látek, co se týče struktury, fyzikálně chemických vlastností a oblastí jejich využití. Mezi CECs řadíme například pesticidy, farmaka, výrobky pro osobní hygienu, zpomalovače hoření a další. CECs se vyskytují v odpadních a povrchových vodách, přičemž nejsou pravidelně monitorovány ani nijak regulovány. Díky častému výskytu a možným nežádoucím účinkům na živé organismy vyvolává skupina CECs znepokojení veřejnosti a velkou pozornost vědecké komunity. Nově navrhované přístupy k monitoringu povrchových vod s ohledem na CECs proto zahrnují speciální vzorkovací techniky, sledování biologických účinků citlivými biotesty, necílený chemický screening a užití pokročilých softwarových nástrojů k interpretaci získaných dat. Dále také efektem řízená analýza (EDA – effect directed analysis), založená na kombinaci biotestů, frakcionace, a chemické analýzy, je účinnou metodou k identifikaci látek zodpovědných za pozorované biologické účinky v komplexních environmentálních směsích. Cílem této disertace bylo přispět k identifikaci látek zodpovědných za toxické účinky sledované v povrchových vodách užitím a dalším vývojem metody EDA. Dílčími cíli pak byly i) využití zjednodušeného protokolu pro efektem řízenou analýzu (studie I a IV); ii) identifikace inhibitorů růstu zelených řas a acetylcholinesterázy v extraktu přečištěné odpadní vody (studie II a V) a iii) charakterizace znečištění řeky Bosny nově se objevujícími znečišťujícími látkami V rámci studie I a IV proběhla vzorkovací kampaň ve čtyřech povodích Evropských řek, přičemž bylo testováno a optimalizováno nové vzorkovací zařízení založené na extrakci velkého objemu vody na pevné fázi. Nejčastěji sledovaným biologickým efektem byla estrogenita, dále toxicita pro zelené řasy a rybí embrya. Převážnou část těchto efektů nebylo možné vysvětlit výskytem cílových analytů. Na základě zjednodušené analýzy rizik cílových analytů bylo identifikováno a prioritizováno 21 látek zajímavých pro budoucí monitoring povrchových vod. K identifikací látek způsobujících toxicitu pro zelené řasy (Studie II) a inhibici acetylcholinesterázy in vitro (studie V) byla využita efektem řízená analýza extraktu přečištěné odpadní vody. Výsledky studie II naznačují, že pro inhibici růstu zelených řas jsou určující pesticidy a produkty jejich přeměny, 9 některá farmaka (deriváty barbiturátů, makrolidy), průmyslové látky a kofein včetně jeho metabolitů. Znečištění řeky Bosny CECs bylo zkoumáno využitím kombinace screeningu biologických účinků a chemických analýz provedené na výluzích pasivních vzorkovačů (Studie III). Hlavní zdroj kontaminace řeky Bosny CECs a ohniska sledované biologické aktivity se podařilo identifikovat porovnáním kumulativních koncentrací cílových analytů a profilů nebezpečnosti. Zjednodušená analýza rizik cílových analytů poukázala na 7 látek (diazinon, diclofenac, 17-β estradiol, estrone, benzo[k]fluoranthene, fluoranthene and benzo[k]fluoranthene), které mohou představovat riziko pro vodní biotu. Tato disertční práce demonstruje, že kombinace nových vzorkovacích technik s metodami k posouzení biologických účinků a chemickými analýzami je dobrou strategií, jak přistupovat ke znečištění povrchových vod CECs. Tato práce také prokazuje užitečnost metody EDA pro snížení komplexity environmentálních směsí a identifikaci látek způsobujících toxicitu. Analýza rizik a prioritizace detekovaných látek je důležitým krokem v procesu hledání specifických znečišťujících látek povodí.
10
LIST OF ORIGINAL PUBLICATIONS AND THE AUTHOR’S CONTRIBUTION
Publication I (see ANNEX I) Tousova, Z., Oswald, P., Slobodnik, J., Blaha, L., Muz, M., Hu, M., Brack, W., Krauss, M., Di Paolo, C., Tarcai, Z., Seiler, T.-B., Hollert, H., Koprivica, S., Ahel, M., Schollée, J.E., Hollender, J., Suter, M.J.-F., Hidasi, A.O., Schirmer, K., Sonavane, M., Ait-Aissa, S., Creusot, N., Brion, F., Froment, J., Almeida, A.C., Thomas, K., Tollefsen, K.E., Tufi, S., Ouyang, X., Leonards, P., Lamoree, M., Torrens, V.O., Kolkman, A., Schriks, M., Spirhanzlova, P., Tindall, A., Schulze, T., 2017. European demonstration program on the effect-based and chemical identification and monitoring of organic pollutants in European surface waters. Sci. Total Environ. 601–602. https://doi.org/10.1016/j.scitotenv.2017.06.032
Zuzana Toušová participated in the planning of the study, conducted designated parts of the study (sampling in Danube RB, sample processing, algal bioassays), evaluated and interpreted relevant results (linking of biological effects and detected compounds, risk assessment and prioritization) and prepared the manuscript.
Publication II (see ANNEX II) Tousova, Z., Froment, J., Oswald, P., Slobodník, J., Hilscherova, K., Thomas, K.V., Tollefsen, K.E., Reid, M., Langford, K., Blaha, L., 2018. Identification of algal growth inhibitors in treated waste water using effect-directed analysis based on non-target screening techniques. J. Hazard. Mater. 358. https://doi.org/10.1016/j.jhazmat.2018.05.031
Zuzana Toušová participated in the planning of the study, conducted designated parts of the study (sampling, sample processing, algal bioassays, fractionation, toxicity data search in literature and databases, ECOSAR modelling), evaluated and interpreted relevant results and prepared the manuscript.
Publication III (see ANNEX III) Toušová, Z., Vrana, B., Smutná, M., Novák, J., Klučárová, V., Grabic, R., Slobodník, J., Giesy, J.P., Hilscherová, K., 2019. Analytical and bioanalytical assessments of organic micropollutants in the Bosna River using a combination of passive sampling, bioassays and multi-residue analysis. Sci. Total Environ. 650. https://doi.org/10.1016/j.scitotenv.2018.08.336 11
Zuzana Toušová participated in the planning of the study, conducted designated parts of the study (sample processing, in vitro bioassays, literature and databases, ECOSAR modelling), evaluated and interpreted relevant results (linking of biological effects and detected compounds, contamination profiling and hazard assessment), and prepared the manuscript.
Publication IV (see ANNEX IV) Schulze, T., Ahel, M., Ahlheim, J., Aït-Aïssa, S., Brion, F., Di Paolo, C., Froment, J., Hidasi, A.O., Hollender, J., Hollert, H., Hu, M., Kloß, A., Koprivica, S., Krauss, M., Muz, M., Oswald, P., Petre, M., Schollée, J.E., Seiler, T.-B., Shao, Y., Slobodnik, J., Sonavane, M., Suter, M.J.-F., Tollefsen, K.E., Tousova, Z., Walz, K.-H., Brack, W., 2017. Assessment of a novel device for onsite integrative large-volume solid phase extraction of water samples to enable a comprehensive chemical and effect-based analysis. Sci. Total Environ. 581–582. https://doi.org/10.1016/j.scitotenv.2016.12.140
Zuzana Toušová participated in the planning of the study, conducted designated parts of the study (sampling in Danube RB, sample processing, algal bioassay) and contributed to writing of the manuscript.
Publication V (see ANNEX V) Ouyang, X., Leonards, P.E.G., Tousova, Z., Slobodnik, J., De Boer, J., Lamoree, M.H., 2016. Rapid Screening of Acetylcholinesterase Inhibitors by Effect- Directed Analysis Using LC × LC Fractionation, a High Throughput in Vitro Assay, and Parallel Identification by Time of Flight Mass Spectrometry. Anal. Chem. 88. https://doi.org/10.1021/acs.analchem.5b04311
Zuzana Toušová carried out the sampling and sample processing, and contributed to writing of the manuscript. 12
LIST OF ABBREVIATIONS
AR androgen receptor AChE acetylcholinesterase BEQ bioanalytical equivalent concentration CA concentration addition CECs contaminants of emerging concern CI contamination index CUPs currently used pesticides [1-(3,4-dichlorophenyl)-3-methylurea], diuron DCPMU transformation product Dex-EQ dexamethasone equivalent DHT dihydrotestosterone DHT-EQ dihydrotestosterone equivalent DMSO dimethylsulfoxid E1 estrone E2 17ß-estradiol E2-EQ 17-ß-estradiol equivalent E3 estriol effect concentration 50 - concentration at which the effect EC50 reaches 50% of the effect in untreated control EDA effect directed analysis EDP European demonstration program EE2 17a-ethinylestradiol EQS environmental quality standard ER estrogen receptor EU European union Flu-EQ flutamide equivalent GC gas chromatography (two dimensional) gas chromatography coupled to mass GC(xGC)-MS spectrometry HI hazard index HQ hazard quotient HR-ToF MS high resolution time of flight mass spectrometry IC inhibitory concentration IA independent action LC liquid chromatography liquid chromatography coupled to high resolution mass LC-HRMS spectrometry LOD limit of detection; 13
LOEC lowest observed effect concentration LOQ limit of quantification LVSPE large volume solid phase extraction MEC95 95th percentile of maximum environmental concentration MeOH methanol MTPexp maximum toxicity potential - experimental
MTPpred maximum toxicity potential - predicted with ECOSAR OCPs organochlorine pesticides OH-Tam -EQ hydroxytamoxifen equivalent OPE organophosphate esters PAH polycyclic aromatic hydrocarbons PBDE polybrominated diphenylethers PCBs polychlorinated biphenyls PCPs personal care products PEC predicted environmental concentration PFAS per- and polyfluoroalkyl substances PFOA perfluorooctanoic acid PFOS perfluorooctanesulfonate PNEC predicted no effect concentration POCIS polar organic chemical integrative sampler PPP plant protection product PRC performance reference compound RB river basin RBSP river basin specific pollutants REF relative enrichment factor REP relative effect potency RP-HPLC reverse phase - high performance liquid chromatography RP-SPE reverse phase-solid phase extraction SPMD Semi permeable membrane device 2,3,7,8-Tetrachlorodibenzo-p-dioxin; TCDD-EQ – TCDD TCDD equivalent TU toxic unit UHPLC ultra-high-performance liquid chromatography US United states WFD Water framework directive WW wastewater WWTP wastewater treatment plant
14
TABLE OF CONTENTS
1 INTRODUCTION ...... 16 1.1 Monitoring of water quality ...... 16 1.2 Contaminants of emerging concern (CECs) ...... 17 Main classes of CECs: ...... 18 1.2.1 Artificial sweeteners ...... 18 1.2.2 Per- and polyfluoroalkyl compounds (PFASs) ...... 19 1.2.3 Pharmaceuticals and steroids ...... 19 1.2.4 Personal care products (PCPs) ...... 20 1.2.5 Brominated and emerging flame retardants ...... 20 1.2.6 Benzotriazoles and benzothiazoles ...... 21 1.2.7 Pesticides ...... 21 1.2.8 Others ...... 21 1.3 Sampling ...... 22 1.3.1 Large volume solid phase extraction (LVSPE) ...... 22 1.3.2 Passive sampling ...... 23 1.4 Chemical analyses and non-target screening ...... 25 1.5 Bioassays ...... 25 1.6 Effect directed analysis (EDA) ...... 26 1.7 Risk assessment (RA) ...... 27 2 Aims of the study ...... 29 3 SUMMARY OF MATERIALS AND METHODS ...... 30 3.1 European demonstration program (EDP) ...... 30 3.1.1 LVSPE sampling ...... 30 3.1.2 Effect assessment ...... 30 3.1.3 Chemical analyses ...... 31 3.1.4 Linking effects and detected compounds ...... 32 3.1.5 Risk assessment and prioritization ...... 32 3.2 Effect directed analysis of WWTP effluent extract ...... 33 15
3.2.1 Fractionation ...... 33 3.2.2 Algal growth inhibition assay ...... 33 3.2.3 Non-target screening ...... 33 3.2.4 Linking effects and detected compounds ...... 33 3.2.5 LCxLC system for identification of AChE inhibitors ...... 34 3.3 Analytical and bioanalytical assessment of the Bosna River ...... 34 3.3.1 Passive sampling ...... 34 3.3.2 Chemical analyses ...... 34 3.3.3 In vitro bioassays ...... 35 3.3.4 Contamination profiling and hazard assessment ...... 35 4 SUMMARY OF RESULTS AND DISCUSSION ...... 36 4.1 European demonstration program (EDP) ...... 36 4.1.1 Sampling ...... 36 4.1.2 Effect assessment ...... 37 4.1.3 Chemical analyses ...... 37 4.1.4 Risk assessment and prioritization ...... 39 4.2 Effect directed analysis of the WWTP effluent extract ...... 41 4.2.1 Toxicity to algae ...... 41 4.2.2 AChE inhibition ...... 43 4.3 Analytical and bioanalytical assessment of the Bosna River ...... 45 5 CONCLUSIONS AND FUTURE PROSPECTS...... 48 6 REFERENCES ...... 50 7 ANNEXES...... 60
16
1 INTRODUCTION
1.1 Monitoring of water quality Environmental quality monitoring of surface waters is fundamental to the sustainable management of water resources and to reducing risks posed by multiple anthropogenic stressors (Geissen et al., 2015). In spite of the marked improvement of water quality in European surface water bodies, the ultimate goal of the European Union (EU) water policy, i.e. reaching good ecological and chemical status, still remains a challenge in most cases (European Environment Agency, 2018). The current monitoring scheme delineated by the Water Framework Directive (WFD), based on a few priority substances and biological indicators, fails to address a wide range of chemicals occurring in complex environmental mixtures, which may elicit adverse biological effects (Directive 2000/60/EC, Directive 2008/1005/EC, Directive 2013/39/EC). Novel holistic approaches to water quality monitoring with focus on the mixture risks and identification of the main toxicity drivers have been proposed (Fig. 1). These empathize the use of tailored sampling techniques, effect-based methods, chemical suspect and non-target screening followed by linking of the chemical and biological data (Altenburger et al., 2019). Integrative monitoring strategies reflect the complex character of water pollution and aim to bring relevant information to the water management bodies about the deleterious biological effects of chemical contaminants at reasonable cost (Brack et al., 2019).
Figure 1: Scheme of the current and proposed future approaches to surface water quality monitoring under the WFD (adapted from European Environment Agency, 2018). EQS – Environmental quality standard, RBSP – River basin specific pollutant INTRODUCTION 17
1.2 Contaminants of emerging concern (CECs) Chemicals known as contaminants of emerging concern (CECs) comprise a diverse group of compounds including personal care products, human and veterinary pharmaceuticals, surfactants and surfactant-derived compounds, X-ray contrast media, pesticides, disinfection by-products, algal toxins, flame retardants, plasticizers, UV-filters, industrial compounds and transformation products (Sima et al., 2014). These compounds typically occur at rather small (down to sub-ng L-1) concentrations. However, some of these are sufficiently potent or have the potential to be accumulated to concentrations such that they can elicit biological effects. Due to the broad range of physicochemical properties and usage patterns, CECs constitute a major challenge for analytical chemists and environmental regulators (Ginebreda et al., 2014). The main sources of CECs entering surface waters are wastewater treatment plant (WWTP) effluents, agricultural and urban runoff (Fairbairn et al., 2018). The occurrence and fate of CECs in surface waters varies widely, depending on their physicochemical properties, regional usage and efficiency of wastewater treatment. In surface waters, CECs attenuate naturally by biotransformation, photolysis, sorption, volatilization, and dispersion, or a combination there of (Pal et al., 2010). While some CECs have short half-lives e.g. certain pharmaceuticals, others are susceptible to poor removal during conventional wastewater treatment processes e.g. artificial sweeteners (Loos et al., 2013) Eventually, some compounds e.g. bisphenol A, became virtually ubiquitous as they occur in remote areas with no direct anthropogenic impact (Weissinger et al., 2018). Exposure of CECs to aquatic biota is essentially constant, even for those with short half- lives, because the feed from wastewater treatment plants is continuous (Vera-Candioti et al., 2008). CECs are an issue of major concern as they may have negative impact on human health and aqueous ecosystems through various modes of toxic action e.g. endocrine disruption, immunotoxicity, oxidative stress or genotoxicity and embryotoxicity (Connon et al., 2012). These specific biological effects were repeatedly observed for surface waters, drinking waters, raw and treated wastewaters and technical systems (Dizer et al., 2002, Loos et al., 2013; Jálová et al., 2013; Osman et al., 2015; Zhao et al., 2011) and some can be associated with particular compounds as the main toxicity drivers, e.g., estrogenicity due to E2 and EE2 (Aerni et al., 2004; Koenig et al., 2017; Miège et al., 2009) or phytotoxicity to herbicides (Tang and Escher, 2014). However, the cause of many effects observed in bioassays remains unexplained by the chemicals detected in the same samples (Burgess et al., 2013; Neale et al., 2015). 18
CECs occur in highly complex and variable environmental mixtures. Combined effects of individual constituents, phenomenon known as ‘cocktail effect’, present a major challenge, since considering one substance at a time may lead to underestimation of the overall environmental toxicity (Posthuma et al., 2019). Research is moving towards new strategies to evaluate the effects of chemical mixtures. Two main approaches are available to study the toxicity of mixtures: “whole mixture” and “component-based”. The first directly assesses the toxicity of a chemical mixture, considers all the chemicals in the mixture and evaluates the realistic exposure scenario. However, the results cannot be extrapolated to different samples and are only applicable to the investigated mixture. The component-based approach determines the toxicity of the mixture on the basis of the effects of its individual components, following two concepts of potential interaction between chemicals: Concentration Addition (CA) and Independent Action (IA). CA assumes similar activity and a common mode of action for all the compounds in the mixture, while IA implies that the compounds act on different molecular targets that lead to a common toxicological effect (Riva et al., 2019). CECs are largely unregulated, although several CECs were added to the EU WFD list of priority substances for chemical monitoring and three compounds, i.e. diclofenac, 17-beta-estradiol and 17-alpha-ethinylestradiol, have been monitored as a part of the WFD Watch list (Directive 2008/1005/EC, Directive 2013/39/EC). However, the focus on a few routinely monitored compounds encourages reduction in their use, while replacement of these chemicals by alternatives that pose similar hazards is an unresolved problem e.g. atrazine vs. terbuthylazine in the EU (Brack et al., 2019).
Main classes of CECs:
1.2.1 Artificial sweeteners Artificial sweeteners such as sucralose, acesulfame, saccharin, aspartame, etc. are widely used in high concentrations as sugar substitutes in food, drinks and pharmaceuticals (Scheurer et al., 2009). Artificial sweeteners are rarely metabolized and primarily enter the environment through treated wastewater. Due to their stability and high levels found in the environment (units of µg/L in river waters), sucralose and acesulfame are recognized as potential tracers of anthropogenic inputs into environmental waters (Richardson and Ternes, 2018). Artificial sweeteners are not expected to bioaccumulate, but their impacts on the behavior of aquatic organisms and toxicity of sweetener transformation products were shown in several studies. INTRODUCTION 19
Environmentally relevant levels (0.05 and 155 μg/L) of sucralose caused changes in oxidative biomarkers on the the gills, muscle, brain, and liver of carp and also oxidative damage in lipids and proteins (Saucedo-Vence et al., 2017). Increased oxidative stress in the liver of Carassius auratus after the exposure to phototransformation products of acesulfame was reported by Ren et al., (2016).
1.2.2 Per- and polyfluoroalkyl compounds (PFASs) PFASs such as perfluorooctanoic acid (PFOA) and perfluorooctanesulfonate (PFOS) are substances containing a fully fluorinated hydrophobic alkylchain with variable number of carbon atoms. Due to their unusual chemical properties, repelling both grease and water, PFASs are used in food packaging, fabrics (e.g., waterproof jackets), carpets, nonstick cooking pans, paints, adhesives, electronics, personal care products, and firefighting foams. PFASs are highly stable in the environment, can be transported over long distances and accumulate in red blood cells (Richardson and Ternes, 2018). PFASs may enter the aquatic environment directly from production facilities, further due to the use of products with their content, from wastewater treatment plants and furthermore from degradation of long-range- transported volatile precursors (Bečanová et al., 2016). The most common PFASs (PFOA and PFOS) have been regulated in the EU and US as they were identified as persistent, bioaccumulative and toxic (OECD, 2002), however, the replacement PFASs remain a concern. While the acute toxicity of PFASs is moderate, toxicological studies reported adverse effects of PFASs spanning from hepatotoxicity, neurotoxicity, immunotoxicity, developmental and reproduction toxicity to endocrine disruption (Rainieri et al., 2017).
1.2.3 Pharmaceuticals and steroids Pharmaceuticals bear diverse classes of organic compounds such as analgesics, antiepileptics, antihyperlipidemics, non-steroidal anti-inflammatory drugs, synthetic hormones, antimicrobials and many more. Pharmaceutical administrations are not always entirely metabolized within the body and can be excreted as unaltered or an active metabolite of the parent compound, which may undergo further transformation or remain intact when they pass through WWTPs (Salimi et al., 2017). Globally, the major source of pharmaceuticals in surface waters is the discharge of urban wastewater, although emissions from the pharmaceutical industry, agriculture and aquaculture can be important on a local scale (Richardson and Ternes, 2018). Pharmaceuticals, their metabolites and transformation products may pose both acute and long-term chronic 20
effects on aquatic biota at low, environmentally relevant concentrations (Patel et al., 2019). Adverse biological effects of pharmaceuticals on biota comprise e.g. embryotoxicity, mutagenicity, gill alterations, growth inhibition and changes in foraging behavior. Among all, the great advance of antibiotic resistance in environmental microbial communities raises major concerns (Patel et al., 2019). Natural endogenous (e.g. 17β-estradiol, estrone, estriol, progesterone, cortisone, testosterone and thyroid hormone) and synthetic steroids (17α- ethinyloestradiol, mestranol) excreted by humans and livestock enter surface waters by direct discharges or WWTP effluents (Thomaidis et al., 2012). Steroids together with non-steroidal chemicals from various classes known as endocrine disruptors (e.g. bisphenol A, phthalates, phenothrin) cause adverse endocrine effects such as fish feminization.
1.2.4 Personal care products (PCPs) Personal care products (PCPs) are a diverse group of chemicals including disinfectants (e.g. triclosan), fragrances (e.g. musks), insect repellants (e.g.DEET), preservatives (e.g. parabens) and UV filters (Benzophenone-3), which are used in consumer products like soaps, lotions, toothpaste, fragrances, sunscreens etc., and they are applied directly on the human body to change appearance, taste or odour (Brausch and Rand, 2011). PCPs enter the aquatic environments through recreational activities such as swimming and also via showering and bathing as well as other technological process (Tijani et al., 2016). PCPs are ubiquitous and even though some may have restricted lifetime in the aquatic environment, their levels in surface waters are very stable as they are released constantly (Ohoro et al., 2019). PCPs were reported to cause a wide range of adverse effects in aquatic organisms such as endocrine disruption (parabens), changes in developmental rates of invertebrates and fish (polycyclic musks), or acute toxicity to algae (disinfectant triclosan), (Brausch and Rand, 2011).
1.2.5 Brominated and emerging flame retardants Flame retardants are chemical additives used in plastics, textiles, electronic circuitry, and other materials, from the class of polybrominated diphenyl ethers (PBDEs), organophosphate esters (OPEs) and chloro- organophosphates (Thomaidis et al., 2012). PBDEs were banned in some countries, due to their widespread presence in the environment and adverse effects on biota (developmental neurotoxicity and others), but they are still being released to the environment from older products (Linares et al., 2015). INTRODUCTION 21
New replacement flame retardants, OPEs and chloro-organophosphates, mainly released from WWTPs, have emerged in surface waters (Salimi et al., 2017). However, these compounds were also reported to cause neurotoxicity, embryotoxicity or disturbance of thyroid signaling at higher than current environmental levels, and therefore need to be carefully monitored (Greaves and Letcher, 2017).
1.2.6 Benzotriazoles and benzothiazoles Benzotriazoles and benzothiazoles are widely used as corrosion inhibitors, UV stabilizers in plastics and polymers or in rubber production (Richardson and Ternes, 2018) . They are highly soluble in water, resistant to biological degradation, and not well removed in wastewater treatment. In surface waters, they are present at low ng/L levels but can reach up to tens of µg/L (Shi et al., 2019). Benzotriazoles have high bio-accumulation potential and they were shown to elicit multiple toxic effects including changes in mitochondrial bioenergetics, embryonic development, and locomotor activity of zebrafish, in vitro anti-androgenic and in vivo estrogenic effects, neurotoxicity and hepatotoxicity in fish or alteration of molting frequency in Daphnia (Liang et al., 2019; Shi et al., 2019; Tangtian et al., 2012).
1.2.7 Pesticides Pesticides are a wide class of chemicals used to limit, inhibit, prevent the growth or repel harmful animals, insects, weeds, invasive plants, and fungi. Agricultural activity is the main source of pesticides in the environment, however, industrial emissions during their production and urban run-off also significantly contribute to the levels of pesticides in surface waters (Pietrzak et al., 2019). The levels of pesticides in surface waters fluctuate in time, depending on the application pattern and weather. Pesticides in surface waters may undergo degradation by hydrolysis, oxidation, biodegradation, or photolysis, while for some compounds the transformation products can occur at higher levels and show higher toxicity than the parent compounds (Houtman, 2010). Pesticides elicit a multitude of toxic effects to aquatic life, threaten human health and water pollution due to pesticides is a priority issue of global concern (Pietrzak et al., 2019).
1.2.8 Others Besides the main, above mentioned groups, CECs comprise other compounds or compound classes which have been increasingly detected in surface waters and therefore should not stay out of the scientific and regulatory 22
focus, e.g. halogenated methanesulfonic acids, ionic liquids, algal toxins, dioxane, disinfection byproducts, siloxanes and bisphenol A.
1.3 Sampling Water is an extremely heterogeneous matrix where the distribution and mixing of waterborne chemicals are affected by the hydrodynamics of the water, the sorption partition coefficients of the chemicals, and the amount of organic matter. The concentrations of contaminants vary spatially and temporally, because episodic events from surface runoff, spills, and other point source contamination can result in isolated or short-lived chemical pulses in the water (Alvarez and Jones-Lepp, 2010). Traditional sampling methods i.e. grab sampling have several limitations for monitoring of CECs in surface waters. Modern chemical analytical instrumentation allows for the analysis of small water volumes with no or only low sample enrichment for most of the typical water micropollutants pollutants, which occur at trace levels, while the analysis of some priority substances with very low EQS values as well as in vivo and in vitro tests may require greater enrichment and larger water volumes (Neale et al., 2015). Grab samples represent the concentration of CECs only at the instant of sampling and even repetitive grab sampling schemes usually fail to capture the concentration pulses during episodic events (Vrana et al., 2005). To overcome these limitations, integrative sampling techniques like passive sampling or solid phase extraction have been developed (Dimpe and Nomngongo, 2016; Jones et al., 2015).
1.3.1 Large volume solid phase extraction (LVSPE) The combination of targeted and nontargeted chemical screening analysis in combination with bioassays has been recommended for the identification of (eco-)toxicologically active compounds and mixtures by several studies (Altenburger et al., 2015; Brack et al., 2015). Large sample volumes are required for such a combined approach, which may present logistic, technical, economic and scientific issues related to the storage and transport. To prevent these problems, a novel automated sampling device based on solid phase extraction (SPE) was developed by Schulze et al. (2017). Solid phase extraction (SPE) is the most powerful sampling and enrichment technique for complex mixtures of known and unknown contaminants, which are captured and stabilized on the sorbent during the sampling. Several well- tested and widely used solid phases, based on C18 or polystyrene- divinylbenzene (co-)polymers, that trap organic compounds with a broad range INTRODUCTION 23
of properties (nonpolar to polar, neutral to charged) are commercially available (Fontanals et al., 2010). The newly introduced on-site large volume solid phase extraction (LVSPE) device combines a pre-filtration cartridge (separation of suspended particulate matter (SPM) from the water phase) with a set of tailor- made columns that allow customizable selection of sorbents to cover a broad range of compounds with different properties (Fig. 2). To force the water through the extraction columns, the LVSPE device uses a pressurized system which is powered by a car battery and controlled by 12 V electronic components with low energy consumption (Schulze et al., 2017).
Figure 2: (a) Picture of the LVSPE50 device; (1): Dosing system (500 mL), (2): pre-filter (3): ball valve, (4): pressure chamber (550 mL), (5):extraction cartridge, (6): controller (Photo by MAXX GmbH); (b) LVSPE device during field sampling at the WWTP effluent in Brno – Modřice, Czech. Rep. in August 2014.
1.3.2 Passive sampling Passive sampling is a sampling technique based on free flow of analyte molecules from the sampled medium to a collecting medium, as a result of a difference in chemical potentials of the analyte between the two media. (Górecki and Namienik, 2002). Passive samplers are non-mechanical devices, that require minimal resources of personnel and equipment for sampling and avoid almost every disadvantage of active sampling and/or of the methods of preparing the samples mentioned above. Materials used in passive sampler construction have constant composition and well-defined diffusion and partition properties a constant uptake capacity (Taylor et al., 2019). Passive sampling provides a sensitive measurement of dissolved concentrations that is integrated over time (Cfree). Due to its proportionality to the chemical activity and chemical potential, Cfree is considered a key parameter in understanding 24
chemical’s exposure of aquatic organisms (Reichenberg and Mayer, 2006). Passive sampling enables integrative collection of contaminants over an extended period of time and captures residues from episodic events, which typically remain undetected when using grab sampling (Alvarez et al., 2004; Vrana et al., 2005). Nowadays, various passive samplers’ designs, multiple receiving phases and diffusion membranes are available. Passive samplers can be used for sampling of a wide variety of compounds, e.g. semipermeable membrane device (SPMD) for hydrophobic substances such as PAHs or PCBs (Huckins et al., 1993) and polar organic chemical integrative sampler (POCIS) for hydrophilic substances such as polar pesticides and pharmaceuticals (Alvarez et al., 2004), (Figure 3). Because of the integrative character of sampling, passive samplers accumulate a sufficient amount of sampled chemicals for detection of small concentrations in water and samples for multiple analyses including bioassays (Jones et al., 2015; Moschet et al., 2014; Vrana et al., 2014). Passive sampling has been successfully used for combined chemical and effect analyses of CECs in surface waters in many studies (Emelogu et al., 2013; Jalova et al., 2013; Jarosova et al., 2012; Novák et al., 2018). Nevertheless, the recognition of the added value of passive sampling by regulators is still rather limited. In Europe, the strict monitoring requirements of the EU WFD hinder the implementation of passive sampling for regulatory purposes, whereas environmental managers in the United States have more freedom to apply approaches in line with the current scientific insights (Booij et al., 2016). The main weakness of passive samplers from the regulators’ perspective is insufficient quality control related to variability in the reported results due to inaccuracies of the partition coefficients of target analytes (Booij et al., 2016). This issue is more pronounced in case of polar compounds (Křesinová et al., 2016).
Figure 3: Stainless steel protective cage with POCIS discs in triplicate (left); POCIS disc (middle) and SPMD membrane stretched in the stainless steel protective cage (right) with visible biofouling after deployment in surface water (photos by Dr. Vrana) INTRODUCTION 25
1.4 Chemical analyses and non-target screening Chemical analyses of CECs in surface waters require highly sensitive analytical methods, which have become available with the advance of gas chromatography (GC) and (ultra)-high performance liquid chromatography [(U)HPLC] coupled either to tandem mass spectrometry (MS/MS), or a wide range of high resolution mass spectrometers (García-Córcoles et al., 2019; Richardson and Ternes, 2018). The selection for GC or LC separation usually depends on the physicochemical properties of the analytes. Thus, the less volatile and polar compounds are normally separated with LC, whereas GC is often used for the volatile compounds (Martín-Pozo et al., 2019). Conventional multi-residue target analyses of complex environmental samples like WWTP effluents provide only a limited picture of a mixture’s complete chemical composition. Many compounds are co-extracted in the extraction process, they remain unnoticed during the analysis, but they may contribute significantly to the observed biological effects. Non-target screening techniques using GC and LC-MS together with high-resolution mass spectrometers (HR-MS) have therefore been increasingly applied in environmental analysis (Gómez et al., 2009; Ibáñez et al., 2008; Schymanski et al., 2015; Zedda and Zwienerbu, 2012). A study by Schymanski et al. (2015) showed that while the non-target screening analytical techniques were substantially harmonized among laboratories, the data-processing and evaluation part still presented a challenge. Currently, powerful chromatographic deconvolution and structure elucidation software, in silico prediction tools for fragmentation, retention time and ionization, as well as MS and MS/MS libraries of environmental contaminants are being developed rapidly and enable tentative identification of unknowns in a reasonable time frame with a reasonable confidence (Krauss et al., 2010; Schymanski et al., 2015). Non-target screening has a great potential for characterization of complex environmental mixtures, identification of unknown toxicity drivers, prioritization of substances for monitoring programs and assessment of environmental quality (Du et al., 2017; Hollender et al., 2019).
1.5 Bioassays Adverse biological effects of CECs occurring in surface waters have been studied extensively in the recent decade (Richardson and Ternes, 2018; Rykowska and Wasiak, 2015). To assess these effects, a large number of bioassays indicative of different endpoints at various levels of biological organization have been developed (Escher et al., 2018). In vivo-bioassays using 26
whole organisms like green alga Raphidocelis subcapitata, crustacean Daphnia magna or fish Danio rerio investigate the effects of CECs on apical endpoints such as survival, growth and reproduction. In vitro assays generally measure effects at the cellular level and are mode-of action specific. They allow rapid and sensitive detection and are designed for high throughput applications in the laboratory. Further, several in vitro approaches allow to assess specific endpoints, relevant to CECs, including neurotoxicity, immunotoxicity, oxidative stress response, xenobiotic metabolism regulation, genotoxicity, dioxin-like and hormonal activities like (anti-)estrogenicity, (anti-)androgenicity, glucocorticoidand thyroid activity, etc. (Connon et al., 2012; Neale et al., 2015). These specific biological effects were repeatedly observed in the samples from surface waters, drinking waters, raw and treated wastewaters and technical systems (Dizer et al., 2002, Loos et al., 2013; Jálová et al., 2013; Osman et al., 2015; Zhao et al., 2011). Some can be associated with particular compounds as the main toxicity drivers, e.g., estrogenicity due to E2 and EE2 (Aerni et al., 2004; Koenig et al., 2017; Miège et al., 2009) or phytotoxicity to herbicides (Tang and Escher, 2014). However, the cause of many effects observed in bioassays remains unexplained by the chemicals detected in the same samples (Burgess et al., 2013; Neale et al., 2015). By applying a panel of cellular and small-scale whole- organism assays it is possible to obtain a more holistic profile of the effects of all chemicals present in a water sample without identifying the causative compounds individually. (Escher et al., 2018)
1.6 Effect directed analysis (EDA) Effect-directed analysis (EDA) is a method, which combines bioassays, fractionation and chemical analysis with the aim to identify compounds responsible for the observed biological effects (Brack, 2003). The procedure is typically applied to organic extracts of the sample and starts with biotesting. In case significant effects are detected, the complexity of the sample is sequentially reduced by fractionation, repeated biotesting and elimination of fractions with low or no biological activity (Fig. 4). Several fractionation steps can be carried out until the isolated toxic fractions are subjected to the identification of toxicity drivers by target chemical analysis and non-target screening. Finally, the identified toxicity drivers need to be confirmed as the cause of the measured effect (Brack et al., 2016). EDA is applied in drug discovery, toxicology, forensics, and environmental sciences, where it proved to INTRODUCTION 27
be a powerful tool to facilitate identification of unknown toxicants (Brack et al., 2016; Burgess et al., 2013). EDA studies are very costly, tedious and time consuming. However, novel sampling, bioanalytical, chemical and software tools make EDA studies increasingly feasible through miniaturized and automated high-throughput formats, hyphenated tools, lowered detection limits, and optimized multi-target and non-target screening (Guijarro et al., 2015; Brack et al., 2013).
Figure 4: General scheme of effect directed analysis. Redrawn from Brack (2003)
1.7 Risk assessment (RA) The assessment of potential risk of CECs to aquatic ecosystems has been addressed by several studies (Ginebreda et al., 2014; Kuzmanovic et al., 2014; Smital et al., 2013). These studies combine the established methods of risk assessment of chemicals with novel approaches to overcome the greatest challenge which is data scarcity (Von der Ohe et al., 2011). The risk assessment process considers the exposure level, often referred to as – for example – Predicted Environmental Concentration (PEC), and potential effect of a given substance, referred to as Predicted No-effect Concentration (PNEC). The PECs can be derived from available measured data (concentrations in the environment) and/or model calculations. The PNEC values are ecological safety thresholds usually determined on the basis of results from single species laboratory tests or, in a few cases, established effect and/or no-effect concentrations from model ecosystem tests (e.g. pesticides). PNEC is regarded as a concentration below which an unacceptable effect will most likely not 28
occur (European Commission, 2003). PEC/PNEC risk ratios above 1 indicate that the substance poses risk to aquatic life in the investigated area. Due to the great number of CECs potentially occurring in the environment, there is a need to prioritize these compounds for the management optimization purposes (Kuzmanovic et al., 2014; Von der Ohe et al., 2011). Existing prioritization schemes, mainly based on monitoring and ecotoxicity data as well as predictive modelling, are varied and often highlight different compounds of concern (Letsinger and Kay, 2019). A methodology to prioritize CECs based on the risk of individual compounds was developed within the NORMAN Network (Dulio and Von der Ohe, 2013). This scheme seems most suitable for CECs as it reflects the ever growing list of the compounds of interest along with the evolving availability of monitoring and toxicity data (Dulio et al., 2018). The prioritization is based on the lowest NORMAN PNECs, which were determined with the use of either experimental data, existing environmental quality standards (EQSs), or in silico predictions. Within this prioritization approach, maximum environmental concentrations 95 (MEC95), which are the 95th percentiles of the measured concentrations at the investigated area, are compared to the lowest NORMAN PNEC values to determine frequency of PNEC exceedance (indicator 1), and the extent of exceedance (indicator 2). Indicators 1 and 2 are given risk scores which are used to derive the final NORMAN risk score. Measured compounds on the priority list are then ranked in order of the resulting NORMAN risk score. This ranking than allows to efficiently prioritize various CECs based on their potential risks to aquatic biota.
AIMS OF THE STUDY 29
2 AIMS OF THE STUDY
The overall aims of the present study were to apply and assess different combinations of effect directed analysis (EDA) tools to identify and prioritize river basin specific pollutants (RBSPs) from the class of contaminants of emerging concern (CECs), particularly endocrine disrupters, compounds toxic to microalgae and neurotoxic compounds. Specific objectives of the study were then as follows:
1. to evaluate applicability of a newly developed simplified EDA protocol for effect-based monitoring of surface waters including the use a novel sampling device, set of bioassays, extensive multi-residue analysis and non-target screening (publication I and IV) 2. to identify compounds responsible for algal growth inhibition and acetylcholinesterase inhibition in a WWTP effluent extract using complex EDA approaches based on non-target screening (publication II) and 2- dimensional liquid chromatography system (publication V) 3. to asses water quality in the Bosna River by use of a combination of passive sampling, a set of in vitro bioassays and multi-residue analysis of several compound classes (publication III)
30
3 SUMMARY OF MATERIALS AND METHODS
3.1 European demonstration program (EDP) – a case study on effect-based water monitoring (publications I and IV) Surface water samples were collected at 18 sampling sites in four European river basins in six countries (Figure 5) from July 2013 to August 2014 by means of grab sampling and large volume solid phase extraction (LVSPE). The grab samples (2L) were designated for in-vivo thyroid assay with Xenopus tadpoles. The LVSPE extracts were subject to a set of bioassays addressing 9 endpoints, multi-residue target analysis addressing 151 target compounds and GC-MS non-target screening. Resulting effect data were linked to the measured concentrations of target compounds, which were also used for basic risk assessment and prioritization.
3.1.1 LVSPE sampling An onsite LVSPE device (UFZ, Leipzig, Germany; Maxx GmbH, Rangendingen, Germany) was used to extract organic micropollutants from 50L of river water. The extraction device comprised three different sorbents in sequence, designed to capture neutral, weakly acidic and weakly basic organic compounds. Laboratory blanks were prepared by percolation of 2L of mineralized LC-grade distilled water through the LVSPE device for 100 cycles (equivalent of 50L).
3.1.2 Effect assessment A set of 4 in vitro and 3 in vivo bioassays was applied to screen for 10 endpoints covering both non-specific and specific toxicity (Tab. 1). The effective concentrations of surface water samples were expressed in relative enrichment factors (REFs) as proposed by Escher et al. (2006). The LVSPE extracts were reconstituted in a solvent and mixed with the test media to reach final REFs from 1 to 100, meaning that the tested range covered the original river concentrations (REF=1) as well as concentrations up to 100-times higher (REF=100). The in vivo thyroid activity assay with Xenopus was performed with whole raw water samples. For active samples in the receptor-mediated assays, effect equivalents of standard agonists or antagonists were calculated.
SUMMARY OF MATERIALS AND METHODS 31
Table 1: Overview of bioanalytical methods
Bioassay Biological model Endpoint Method reference
in vitro U2-OS (GR- Glucocorticoid GR – CALUX® Sonneveld et al. (2005) CALUX®) activity Estrogenic MELN (MCF7- Balaguer et al. (1999) ER-mediated activity ERE-Luciferase- activity Antiestrogenic Neomycin) Creusot et al. (2016) activity Androgenic Wilson et al. (2002) AR-mediated activity MDA-kb2 activity Antiandrogenic Creusot et al. (2014) activity Purified enzyme Enzyme Ellman et al. (1961) AChE inhibition AChE inhibition Galgani and Bocquene (1991) in vivo Algal growth Raphidocelis Growth rate OECD guideline 201 inhibition subcapitata inhibition Rojickova and Dvorakova (1998) Survival OECD guideline 236 Zebrafish embryo Danio rerio - Sum of acute toxicity embryo sublethal OECD guideline 236 endpoints Thyroid activity Xenopus laevis - Thyroid Fini et al. (2007) embryo activity
3.1.3 Chemical analyses Neutral LVSPE extracts were subject to target analysis of 151 compounds and to GC-MS non-target screening analysis using state-of-the-art LC-(HR)MS/MS and GC-MS tools (Tab. 2). The list of target compounds was set up to cover several classes of environmentally and toxicologically relevant CECs and their transformation products.
32
Table 2: Overview of chemical analytical methods
Analytical method Target compound Nr. of Method reference /Instrumentation group targets LC-ESI-QFT-HRMS/MS (Thermo Q-Exactive Polar micropollutants 44 Schymanski et al. (2014) Orbitrap) LC-QgQ-MS/MS Industrial solvents, 4 developed for EDP (Thermo TSQ-Vantage) pesticide LC-ESI-HRMS/MS (Thermo LTQ-FT Glucocorticoids 20 Schriks et al. (2010) Orbitrap) LC-ESI-MS/MS Phenolic (Thermo TSQ Quantum 6 Petrovic et al. (2003) xenoestrogens Ultra AM) LC-ESI-ITFT-HRMS (Thermo LTQ Orbitrap Polar micropollutants 47 Hug et al. (2014) XL) LC-MS/MS Steroids 27 Griffith et al. (2014) (AB Sciex Qtrap 6500) GC-EI-MS Legacy pesticides, WFD (Agilent 5795C) priority compounds, 22 modified ISO 6468:1996 industrial compounds GC-EI-MS Non-target screening NA Slobodnik et al. (2012) (Agilent 5795C)
3.1.4 Linking effects and detected compounds For receptor mediated in vitro assays, mass balance calculations were conducted using relative effect potencies (REPs) of known agonists according to Kinani et al. (2010) to quantify the contribution of the detected target compounds to the observed biological activity. For in vivo assays, the zebrafish embryo acute toxicity assay and algal growth inhibition assay, the link between measured compounds and observed effects was calculated using toxic units (TU) and their sum, similarly to Booij et al. (2014) and Kuzmanovic et al. (2014). Toxicity of detected compounds for fish and
green algae was estimated by the program ECOSAR (v1.11) and the 96h-EC50 value
of the most toxic ECOSAR chemical class (minECOSAR 96h-EC50i) was selected (US-EPA, 2012).
3.1.5 Risk assessment and prioritization Risk assessment and prioritization of detected target compounds was based on the NORMAN lowest predicted no effect concentrations (PNECs) developed within the NORMAN Network (Dulio and Von der Ohe, 2013). Target compounds were prioritized in order of NORMAN risk scores, which take into account the frequency and extent of the NORMAN lowest PNEC exceedance by SUMMARY OF MATERIALS AND METHODS 33
each compound measured at individual sampling sites and their 95th percentile (MEC95), respectively.
3.2 Effect directed analysis of WWTP effluent extract (publication II and V) The LVSPE extract of WWTP effluent in Brno Modřice (Czech Rep.) from publication I elicited toxic effects in the algal growth and AChE inhibition assays. The sample was subject to a full higher-tier EDA study to identify AChE inhibitors and a new LVSPE extract was collected at the site to perform a higher-tier EDA study focused on toxicity to algae.
3.2.1 Fractionation (publication II) The LVSPE extract was fractionated in three fractionation steps. The first fractionation was achieved on site by the LVSPE50 sampling device followed by 5 elution steps in the laboratory (5 fractions; F1-F5). The second fractionation of the bioactive fraction F1, was performed with RP-SPE, followed by stepwise elution with water and MeOH mixtures with increasing eluotropic strength (9 fractions; F1.1 - F1.9). The third fractionation of pooled bioactive fractions F1.2, F1.3 and F1.4, was conducted with RP-HPLC resulting in 31 fractions (F1.4.1-F1.4.31).
3.2.2 Algal growth inhibition assay (publication II) Miniaturized algal growth inhibition assay with unicellular green alga Pseudokirchneriella subcapitata (syn. Raphidocelis subcapitata) was carried out according to a modified OECD method 201 (2011) as described by Rojickova and Dvorakova (1998). Sample extracts in DMSO (0.5% v/v) were tested at a relative enrichment factor (REF) up to 200.
3.2.3 Non-target screening (publication II) Both GC- and LC-MS non target screening techniques were applied for analysis of the bioactive fractions. Specific workflows were followed to tentatively identify proposed chemical structures. The workflows are described in detail in the supplementary material of publication II.
3.2.4 Linking effects and detected compounds (publication II) To link the observed bioactivity and the tentatively identified compounds, several experimental datasets and ECOSAR modelling were used. Firstly, effect concentrations (ECx, ICx and LOEC) from laboratory studies with algae and cyanobacteria were extracted from the US-EPA Ecotox database (US-EPA, 2015). 34
Secondly, experimental toxicity data for photosynthesizing organisms, i.e. algae and cyanobacteria, were collected from open literature and thereafter an internal database with more than 1700 experimental data points was built to derive the lowest effective experimental concentration value for each compound (maximum
toxicity potential - MTPexp). Thirdly, green algae 96h EC50s were predicted by the ECOSAR program, v1.11(US-EPA, 2012) to fill the gaps in experimental data. The
toxicity of the ECOSAR class with the lowest EC50 value was selected to represent
the predicted toxicity of each compound (maximum toxicity potential - MTPpred). The detected compounds were ranked in order of their maximum toxicity potential (MTP), whereas experimental data were always given priority over the predicted ones
3.2.5 LCxLC system for microfractionation, high throughput AChE bioassay and parallel HR-TOF detection (publication V) The primary fractionation was achieved on site by the LVSPE50 sampling device followed by 3 elution steps in the laboratory (3 fractions; F1-F3). Neutral eluate (F1), which elicited AChE inhibition, was subject to microfractionation (384 fractions) using a 2-dimensional LC system. The flow of LCxLC effluent was split into the in vitro AChE bioassay (four 96-well plates, 80%) and into the parallel high-resolution time of flight (HR-TOF) mass spectrometric detection (20%).
3.3 Analytical and bioanalytical assessment of the Bosna River (publication III)
3.3.1 Passive sampling Passive samplers (SPMD and POCIS) were deployed at 10 locations (S1- S10) along the Bosna River, BiH, in 2012. At each site 3 POCIS discs and 3 SPMD replicates were co-deployed in a protective cage. Concentrations of compounds from SPMDs were estimated according to an approach based on the sampler-
water partition coefficients (KSW) of the compounds and their sampling rate (RS) calculated on the basis of dissipation of performance reference compounds (PRCs) as described in Vrana et al. (2014), Rusina et al. (2010) and Booij and Smedes (2010). Concentrations of compounds from POCIS samplers were estimated using
median RS value of 0.2 L d-1 according to Harman et al. (2012).
3.3.2 Chemical analyses Chemical analyses of 168 target compounds (134 in POCIS and 34 in SPMD extracts) in 5 compound classes, namely hydrophobic compounds (PAHs, PCBs, SUMMARY OF MATERIALS AND METHODS 35
OCPs), current use pesticides, estrogens and pharmaceuticals, were conducted by use of state-of-the-art GC-MS(-MS) and HPLC-MS(-MS) techniques.
3.3.3 In vitro bioassays Three in vitro cell-based reporter gene bioassays were used to assess ER-, (anti)AR- and AhR-mediated potencies and cytotoxicity of the SPMD and POCIS extracts (Tab. 3). Cytotoxicity of the sample extracts was tested with combination of three dyes (Alamar Blue, CFDA-AM, neutral red) according to Schirmer et al. (1998).
Table 3: Overview of bioanalytical methods Effect Cell line Standard reference Method reference Wilson et al. Androgenicity MDA-kb2 cells DHT: 3.3 pM–100 nM (2002) Antiandrogenicity MDA-kb2 cells FLU: 110 nM – 100 µM Jálová et al. (2013) Demirpence et al. Estrogenicity MVLN cells E2: 1–500 pM (1993) H4G1.1c2 cells Dioxin-like activity (AhR) TCDD: 1–500 pM Nagy et al. (2002) (CAFLUX assay)
3.3.4 Contamination profiling and hazard assessment Toxicity profiles based on the results of applied bioassays were translated into site-specific contamination profiles according to an approach outlined by Hamers et al. (2010). The ratio between the biological response of downstream sites (S2-S10) and a reference site (S1), contamination index (CI), was regarded as a measure of contamination by each of the tested endocrine effects. Assessments of hazards of detected target compounds were conducted by use of the lowest predicted no effect concentrations (PNEC), values derived by the NORMAN Network (Dulio and Von der Ohe, 2013). Hazard quotients (HQs) were calculated as a ratio of the estimated dissolved concentration of an individual compound at a particular sampling site and the NORMAN lowest PNEC value, whereas compounds with HQs exceeding 1 might pose risk to aquatic life. Overall hazard index (HI) was calculated by summation of all HQs of compounds detected at each sampling site.
36
4 SUMMARY OF RESULTS AND DISCUSSION
4.1 European demonstration program (EDP) – a case study on effect-based water monitoring (publications I and IV)
4.1.1 Sampling The novel LVSPE device itself (Fig. 2) and sample processing steps were optimized within the study and its repeated application in field during the EDP sampling campaign (18 sites shown in Fig. 5) demonstrated the suitability of LVSPE for combined chemical and effect-based screening of water quality. Within 2-4 hours, the LVSPE device enabled extraction of 50L, which proved to be sufficient volume for a set of 6 bioassays (with required enrichment up to REF 100) and chemical analyses. The LVSPE device proved to be useful for sampling of complex environmental mixtures of known water contaminants, as the recoveries of 251 spiked compounds (log D range -3,6 – 9,4 at pH 7, concentration range 1-2400 ng.L-1) were acceptable (60-123%) for most compounds (Fig. 6).
Figure 5: Overview map of the EDP sampling locations (clockwise from top left to bottom left) in the Saale RB (Germany), Danube RB (Czech Rep., Slovakia and Hungary), Sava RB (Croatia) and Emme RB (Switzerland). SUMMARY OF RESULTS AND DISCUSSION 37
The effect assessment showed the LVSPE extract enriched at REF 100 allowed a discrimination of active from non-active samples in case of 8 out of 10 toxic endpoints. The neutral (HR-X) sorbents performed better than the weak acidic (HR-XAW) or weak basic (HR-XCW) ones in terms of recovery and capture of toxic effects.
Fig. 6: Scatterplot of the total recoveries (in %) of compounds (N = 251) spiked in a pristine water sample of Wormsgraben (Harz Mountains, Germany) and extracted with the LVSPE50 device versus the water-octanol partition coefficient at pH 7.0 corrected for the speciation (log D).
4.1.2 Effect assessment The most frequently observed effects were estrogenicity, zebrafish acute embryo toxicity and algal growth inhibition (Tab. 4). The biological effects assessed in the EDP program were environmentally relevant, except for the AChE inhibition, which was measurable only in highly enriched extracts
(EC50s>REF 100). With the receptor mediated assays, effects could be detected even for diluted (estrogenicity) or moderately enriched samples (REF<20, androgenicity, anti-estrogenicity and glucocorticoid activity). In the algal growth inhibition and the FET assay, the EC50s were never lower than REF 10 and at some cases exceeded REF 100. In vivo thyroid activity in transgenic Xenopus, using whole water samples, was detected at one sampling site (Saale RB).
4.1.3 Chemical analyses Target analysis of 151 compounds in neutral LVSPE extracts (HR-X) provided extensive chemical characterization of the sampling sites. Out of the total 151 compounds, 107 compounds were detected at least at one sampling 38
site, while 44 compounds, mainly glucocorticoids, steroids and WFD priority compounds, were below detection limits. Twelve compounds occurred at all sampling sites with the exception of the reference site (Saale RB). The concentrations of detected target compounds ranged from a few ng L-1 to a few µg L-1, which is in line with a similar study conducted by Kuzmanovic et al. (2014) and Smital et al., (2013). The highest median concentrations were detected for 1H-benzotriazole, sucralose, phenylbenzimidazolesulfonic acid, 4- toluenesulfonamide and 5-methyl-1H-benzotriazole.
Table 4: Overview of the results of effect assessment linked to the results of target analyses (results of mass balance calculations and identification of the main effect drivers) Range of measured bioactivity Main effect drivers (percentage of Bioactivity (Nr. of active/examined sites) effect explained) Glucocorticoid activity cortisol (29%) 0.3 - 30.5 (4/18) Dex-EQ [ng L−1] 6α-methylprednisolone (2%) estrone (E1) nonylphenoxyacetic Estrogenicity 0.06-1.85 (15/18) acid, bisphenol A and octylphenol E2-EQ [ng L−1] (0 - 77 %) Antiestrogenic activity 35 - 185 ng/L (5/13) unexplained OH-Tam EQ [ng L−1] Androgenic activity 0.93 - 2.7 ng L-1 (4/18) unexplained DHT-EQs [ng L−1] Antiandrogenic activity n.d. n.a. FLU-EQ [pg L−1] AChE inhibition IC50 > 100 (5/18) n.a. [REF] Algal growth inhibition diuron, atrazine and IC50 = 17 – 259 (14/18) [REF] terbuthylazine (30-166%) nonylphenol, nonylphenoxyacetic Zebrafish embryo LC50 = 12-83 (15/18) acid, octylphenol, diazinon, acute toxicity [REF] fipronil (0-20.6%) Thyroid activity 1/18 unexplained n.d. – never detected, n.a. – not available
SUMMARY OF RESULTS AND DISCUSSION 39
4.1.4 Risk assessment and prioritization NORMAN PNEC values were available for 94 out of 151 target compounds. The MEC95 of 15 compounds exceeded the PNEC value and measured concentrations of another 6 compounds exceeded their PNEC values at least once. Target compounds were ranked according to their NORMAN risk score and listed in Table 5. The EDP priority list contained 21 compounds including 4 pesticides and 2 pesticide transformation products, 7 pharmaceuticals (3 of them macrolide antibiotics), 3 surfactant-derived compounds and 1 PAH, biocide, repellent and plasticizer each. 40
Table 5: The list of priority compounds from the EDP case study based on NORMAN prioritization methodology Frequency of Frequency of MEC95 NORMAN PNEC MEC95 Median occurrence PNEC Rank Compound name Compound group / Usage pattern exceedance risk [ng/L] [ng/L] [ng/L] N=18 exceedance of PNEC score [%] [%] Perfluorooctanesulfonic 1 surfactant, WFD priority 0.13 31 4.2 Effect directed analysis of the WWTP effluent extract (publication II and V) 4.2.1 Toxicity to algae (Publication II) Algal growth inhibition was observed in two primary fractions of the WWTP extract, namely F1 (EC50=18.6 REF) and F3 (EC50=56.2 REF). Consequently, F1 was subject to further fractionation by RP-SPE creating 9 fractions, out of which F1.3 was the most active. F1.3 and its neighbouring fractions F1.2, and F1.4 were pooled (EC50=32.9) and further fractionated on RP-HPLC. Out of the total 31 fractions, F1.4.7 (EC50=82.3 REF), F1.4.8 (EC50=84.3) and F1.4.31 (EC50=105.5) were toxic to algae. Toxicity of treated waste waters to algae was described in earlier studies (Köhler et al., 2006; Magdaleno et al., 2014; Maselli et al., 2013), and it remains a concern regarding the chemical and ecological status of the recipient surface waters. In the GC(xGC)-MS non-target screening, 41 and 22 compounds were detected and tentatively identified in the phytotoxic primary fraction F1 and combined tertiary fractions F.1.4.7 and F1.4.8, respectively. These compounds comprised of organophosphates, pharmaceuticals, personal care products, pesticides and industrial pollutants. Two-dimensional GC separation enabled identification of several compounds with co-eluting peaks in the one- dimensional setup, e.g. ticlopidine and fluconazole, promethryn and clorophene, tramadol and desmethyltramadol, venlafaxine and norvenlafaxine. In the LC-HRMS non-target screening of three tertiary fractions active in the bioassay (F1.4.7, F1.4.8, F1.4.31), the total of 345 compounds, including mostly pharmaceuticals, industrial chemicals and by-products, biocides, herbicides and their transformation products, were detected and tentatively identified. The collection of experimental and predicted ecotoxicity data for the detected compounds resulted in the total of 377 Maximum toxicity potential (MTP) values. Experimental values (MTPexp) were available for 68 compounds (17.8%). ECOSAR predicted maximum toxicity potential (MTVpred) values were used for the remaining 309 compounds, while no MTP could be derived for 6 chemicals. Subsequently, detected compounds were ranked according to their MTPs to identify those with the greatest potential to contribute to the observed toxicity. The list of top 10 compounds (Tab. 6), identified as the most likely toxicity drivers, was dominated by herbicides and their transformation products, which is in line with earlier studies focused on phytotoxicity of 42 surface and waste waters (Bengtson Nash et al., 2005; Beate I Escher et al., 2006; Tang and Escher, 2014; Vermeirssen et al., 2010). Table 6: List of the top 10 compounds detected and tentatively identified by GC- and LC-MS non-target screening of phytotoxic fractions of WWTP effluent extract (F1, F1.4.7, F1.4.8 and F1.4.31) ranked according to their estimated maximum toxicity potential. Status under EU Reference to Name Use pesticide MTP [µg/L] the MTP regulation* selective triazine herbicide, US EPA- 1 Terbutryn not approved 0.02 photosystem II inhibitor ECOTOX chlorotriazine herbicide, Booij et al., 2 Terbuthylazine approved 0.40 photosystem II inhibitor 2014 selective triazine herbicide, US EPA- 3 Prometryn not approved 0.41 photosystem II inhibitor ECOTOX industrial formulation of 2-(2-heptadec-8- lubricant additives, 4 enyl-2-imidazolin-1- n.a. 0.84 ECOSAR lubricants, functional fluids yl)ethanol and greases chloroacetanilide herbicide, Nagai and 5 Acetochlor not approved 1.40 elongase inhibitor Taya, 2015 6- chlorinated degradation parent compound US EPA- 6 2 deisopropylatrazine product herbicide atrazine not approved ECOTOX selective herbicide, 7 Flurochloridone biosynthesis of carotenoids approved 3 ECOSAR inhibitor broad- US EPA- 8 Carbendazim spectrum benzimidazole not approved 3.3 ECOTOX fungicide central nervous system (CNS) stimulant of the US EPA- 9 Caffeine n.a. 5 methylxanthine class of ECOTOX psychoactive drugs parent compound metabolite of herbicide not approved in US EPA- 10 DCPMU 5 diuron most EU member ECOTOX states Information on the EU regulatory status retrieved from the Pesticide Properties Database (PPDB) by Lewis et al., 2015 Several pesticides, detected in the active fractions, were not approved in the Czech Republic for use in either plant protection products (PPPs) or biocides i.e. – prometryn and acetochlor. The triazine herbicide, prometryn, detected by GC(xGC)-MS, was ranked as number 3 in the list of likely toxicity drivers in the studied sample. This is an interesting finding because the use of prometryn in PPPs has been banned in the EU since 2003 (or since 2007 in some countries) and the compound is not registered for use in biocides either. Prometryn was reported to be present in effluents and the concentrations in the post ban years did not decrease (Quednow and Püttmann, 2009) and similar trend was reported also for atrazine (Vonberg et al., 2014). Possible SUMMARY OF RESULTS AND DISCUSSION 43 explanations for the continuous input of banned compounds to the environment could be the use of the old stocks, leaching from reservoirs in contaminated soils or the run-off from roof paints, plasters or other building materials enhanced by these compounds (Quednow and Püttmann, 2009; Venzmer, 2008). The results of this study indicated that active substances entering the aquatic environment from biocidal products should be carefully considered. Among the other top-ranking compounds following the group of herbicides, were industrial chemicals, other pesticides, caffeine and its metabolites, barbiturate drugs, macrolide antibiotics and other pharmaceuticals. 4.2.2 AChE inhibition (publication V) The optimized LCxLC system enhanced peak capacity and near perfect orthogonality to provide the basis for high resolution microfractionation in a relatively short time to support fast and accurate identification of active compounds after bioassay screening. Out of total 384 microfractions, 7 were active in the AChE bioassay (threshold 10% effect) and high resolution microfractionation greatly reduced the number of candidate chemicals in the active fractions (Fig. 7). Several pharmaceuticals and their metabolites were tentatively identified in the bioactive microfractions and among those, three pharmaceuticals with applications in neurology/psychiatry were selected for further investigation. Tiapride, lamotrigine and amisulpride were identified as the drivers of AChE inhibition based on the bioanalytical confirmation and earlier findings (Ferrer and Thurman, 2010; Fontaine and Reuse, 1980; Waldmeier et al., 1995; Wode et al., 2015). 44 Figure 7: Heatmap of AChE inhibition (%, n = 3) of the 384 microfractions of the Brno WWTP effluent extract after LC × LC fractionation. The numbers on the top are the plate column numbers (1−12) and the letter and number combinations to the left of the graph show the plate number together with the row number (A-H) of the 4 plates. SUMMARY OF RESULTS AND DISCUSSION 45 4.3 Analytical and bioanalytical assessment of the Bosna River Of the 168 target compounds, 103 compounds were detected in extracts of samplers from at least one sampling site. Specifically, 71 out of 134 compounds (52.9%) were found in POCIS and 32 out of 34 (94.1%) in SPMD extracts. 65 (38.7%) compounds never exceeded their LOQ. There is a clear trend of decreasing cumulative concentration from S3 downstream to S10 in the POCIS samples, while no such pattern can be seen in case of hydrophobic compounds determined in the SPMD. Most analysed compounds were undetectable at the reference site S1 (spring of Bosna) with only a few compounds detected in concentrations near their LOQs. Within the class of PAHs, compounds with 3 and 4 condensed aromatic rings were detected at the greatest concentrations and occurred at all sampling sites. In case of CUPs (20 out of 52 compounds detected), the greatest concentrations were detected for diazinon and prometryn (53 ng L−1 and 17 ng L−1 at sites S2 and S5, respectively) and these two compounds occurred at all sampling sites except for the reference site (S1). Natural estrogens E1, E2 and E3 were detected at all sites in a range of 2.0×10-2-5.8 ng L−1 except for sites S1 and S5. EE2, a synthetic estrogen contained in hormonal contraceptives, was less than the LOQ of 1.0×10-2 – 3.0×10-2 ng L−1 at all sampling sites. This might be related to the 8- to 9-fold lesser prevalence of hormonal contraceptives in BiH compared to western European countries like Belgium or the UK due to cultural and economic reasons as well as limited availability of hormonal contraceptives (Boussen, 2012). Analyses of 77 pharmaceuticals in 11 subclasses resulted in detection of 47 compounds (61%) at least once. The most frequently detected compounds, disopyramide, carbamazepine and caffeine, occurred at all sampling sites. Bioanalytical assessment of the samples resulted in detection of anti- /androgenic, estrogenic and dioxin-like activities (Tab. 7) and the results were in line with earlier studies (Bain et al., 2014; Creusot et al., 2013; Jálová et al., 2013; König et al., 2017; Miège et al., 2009). A profile of integrated effects of mixtures at various locations in the Bosna RB, based on potencies observed in the three bioassays for extracts of the two types of sampler at each location were developed based on comparison to the reference site (S1). Contamination indices (CI), the ratio between the response of downstream sites (S2-S10) and a reference site (S1), in combination with the overall cumulative concentration and number of detected target compounds for 46 each sampling site are shown in Fig. 8. None of the sampling sites downstream of S1 can be considered as uncontaminated as the reference site because all sites exceeded a CI of 1.0 for at least two endpoints. The CI profiles, as well as the cumulative concentration and the number of detected hydrophobic compounds, differed less between individual sites in extracts of SPMD than in extracts of POCIS. Contamination indices indicate that the most contaminated sites were S2 and S3. The greatest cumulative concentrations and numbers of detected compounds were observed for the extract of the POCIS at S3 (a complete analysis for S2 was not available). This result implies that the major source of pollution to the Bosna River was Sarajevo (S2), the capital with a population of about 300,000 (Milinovic, 2013). The trend of decreasing cumulative concentrations in extracts of POCIS samplers downstream of S3 cannot be clearly seen in CIs and no patterns between CIs and cumulative concentrations in extracts of SPMD were observed, despite extensive, multi- residue analyses. Hazard quotients based on the NORMAN lowest PNEC values exceeded the threshold value 1.0 in case of 7 compounds, namely diazinon, diclofenac, E1, E2, benzo[b]fluoranthene, fluoranthene and benzo[k]fluoranthene. The overall hazard index (HI), resulting from the summation of all HQs at each sampling site, indicates that all sites downstream of the reference site S1 might cause adverse effects to aquatic biota as their HIs exceed 1.0. Table 7. Results of the bioanalytical assessment. Sampler Range of measured Nr. of Effect drivers (percentage Effect type bioactivity (Nr. of compounds with of effect explained) active sites) REP Androgenicity SPMD 4.2 – 10 (2) n.a. n.a. DHT-EQs [pg L−1] POCIS 2.1×102 -1.7×103(2) n.a. n.a. SPMD 3.4×104 -5.1×104 2 benzo[a]pyrene (0.03- Antiandrogenicity (5) 0.05%) FLU-EQ [pg L−1] POCIS 2.8×106 - 10 diazinon (0.04-0.06%) 3.2×106(3) SPMD n.d. Estrogenicity POCIS 2.3 ×102 - 2.5×103 8 natural estrogens - E1, E2 E2-EQ [pg L−1] (8) and E3 (0.84-305%) SPMD 9 benzo[b]fluoranthene, 2.9 - 7.3 (9) benzo[k]fluoranthene and Dioxin-like activity chrysene (6.1-24%) TCDD-EQ [pg L−1] POCIS 7 propiconazole (0.02- 31-2.2×102 (6) 0.04%) n.a. – not available, n.d. – not detected SUMMARY OF RESULTS AND DISCUSSION 47 SPMD 400 35 ] pM 30 [ 300 25 20 200 15 100 10 5 Sumconcentration 0 0 S1 S2 S3 S4 S5 S6 S7 S8 S9 S10 Numbercompoundsofdetected Androgenicity 1 74.1 30.0 1 1 1 1 1 1 1 Anti-androgenicity 1 1 1 25.4 22.8 20.0 1 17.1 1 21.6 Estrogenicity 1 1 1 1 1 1 1 1 1 1 Dioxin-like activity 1 17.3 14.8 15.5 13.2 8.43 7.33 7.93 18.3 11.8 ] POCIS pM 5000 70 [ 60 4000 50 3000 40 2000 30 20 1000 10 Sum Sum concentration 0 0 S1 S2* S3 S4 S5 S6 S7 S8 S9 S10 Number of detected compounds Androgenicity 1 5628 707 1 1 1 1 1 1 1 Anti-androgenicity 1 1 1 1 1 1 31.8 30.6 1 27.2 Estrogenicity 1 42.6 194 83.6 86.4 71.9 122 25.5 17.7 1 Dioxin-like activity 1 2.52 1.11 0.87 0.36 0.95 1 1 1 1 0.1-1 Sum concentration 1 Nr of detected compounds 1-10 10-100 100-1000 Reference >1000 Fig. 8: Contamination profiles of sampling sites S2-S10 combined with the total number and cumulative concentration (pM) of detected compounds in SPMD (top) and POCIS (bottom) extracts. Colours green to red indicate to what extent the bioassay responses exceed the bioassay response of the reference site S1 (on a logarithmic scale). *77 pharmaceuticals and 16 target CUPs were not analysed in the POCIS extract from site S2 48 5 CONCLUSION AND FUTURE PROSPECTS The presented dissertation demonstrated a successful application of a newly developed simplified effect-directed analysis (EDA) protocol. The novel on-site LVSPE sampling device proved to be a valuable instrument for integrative effect-based and chemical monitoring purposes. The set of selected bioassays enabled detection of effects relevant for surface waters on an EU- wide scale whereas the most frequently observed were estrogenicity, fish embryo toxicity and toxicity to algae. The findings of this dissertation confirmed that the non-target screening techniques may bring significant added value to the classical chemical target monitoring anchored in the current EU legislation (WFD). In the case of LC-MS non-target screening, however, the data have only qualitative character, and further research is needed to quantify the real concentrations of the identified compounds. This may be challenging especially in case of transformation products, for which analytical standards are hardly available. In this thesis, several compounds were identified as novel candidates for future surface water monitoring campaigns across Europe, e.g. ibuprofen, terbuthylazine, triphenylphosphate and nonylphenoxyacetic acid, as well as macrolide antibiotics such as azithromycin. The results also showed that many WFD priority compounds are of lesser concern, while several CECs, particularly from the class of pesticides and pharmaceuticals may pose risk to aquatic life. On the other hand, our studies indicated that several pesticides, highly prioritized by our research, were no longer registered for use in plant protection products or biocides, which indicated that their non-agricultural input into aquatic environment via WWTPs is not negligible and the sources of these compounds should be explored further. Most of the biological effects observed and described in the particular studies could not be explained by an extensive list of target compounds. which clearly demonstrated the need for implementation of the effect-based tools in surface water monitoring programs. Further investigation using higher tier EDA studies could possibly identify the yet unknown toxicity drivers of the observed effects. Another important aspect, which was only partially addressed by the research in this dissertation, is mixture effects. Understanding the interactions between compounds in complex environmental mixtures is crucial to proper evaluation of the water quality assessment based on individual components of the mixture. CONCLUSIONS AND FUTURE PROSPECTS 49 In conclusion, the objectives of the present dissertation were successfully addressed by conducted research. The studies demonstrated the advantages and needs to combine novel sampling tools, effects-based methods and chemical analyses including non—target screening, in order to address the myriad of contaminants of emerging concern and to achieve a comprehensive assessment of surface water quality. In future, the combined approaches should be further enhanced by application of modelling of mixtures responses and effect-directed analysis together with non-target identification of chemicals to prioritize the most relevant toxicants and effect drivers. 50 6 REFERENCES Aerni, H.R., Kobler, B., Rutishauser, B. V., Wettstein, F.E., Fischer, R., Giger, W., Hungerbuhler, A., Marazuela, M.D., Peter, A., Schonenberger, R., Vogeli, A.C., Suter, M.J.F., Eggen, R.I.L., 2004. Combined biological and chemical assessment of estrogenic activities in wastewater treatment plant effluents. Anal. Bioanal. Chem. 378, 688–696. https://doi.org/10.1007/s00216-003-2276-4 Altenburger, R., Ait-Aissa, S., Antczak, P., Backhaus, T., Barceló, D., Seiler, T.-B., Brion, F., Busch, W., Chipman, K., de Alda, M.L., de Aragão Umbuzeiro, G., Escher, B.I., Falciani, F., Faust, M., Focks, A., Hilscherova, K., Hollender, J., Hollert, H., Jäger, F., Jahnke, A., Kortenkamp, A., Krauss, M., Lemkine, G.F., Munthe, J., Neumann, S., Schymanski, E.L., Scrimshaw, M., Segner, H., Slobodnik, J., Smedes, F., Kughathas, S., Teodorovic, I., Tindall, A.J., Tollefsen, K.E., Walz, K.-H., Williams, T.D., Van den Brink, P.J., van Gils, J., Vrana, B., Zhang, X., Brack, W., 2015. Future water quality monitoring — Adapting tools to deal with mixtures of pollutants in water resource management. Sci. Total Environ. 512–513, 540–551. https://doi.org/10.1016/j.scitotenv.2014.12.057 Altenburger, R., Brack, W., Burgess, R.M., Busch, W., Escher, B.I., Focks, A., Hewitt, L.M., Jacobsen, B.N., Alda, M.L. De, Aissa, S.A., Backhaus, T., Ginebreda, A., Hilscherová, K., Hollender, J., Hollert, H., Neale, P.A., Schulze, T., 2019. Future water quality monitoring : improving the balance between exposure and toxicity assessments of real ‑ world pollutant mixtures. Environ. Sci. Eur. 1–17. https://doi.org/10.1186/s12302-019-0193-1 Alvarez, D., Jones-Lepp, T., 2010. Sampling and Analysis of Emerging Pollutants, Water Quality Concepts, Sampling, and Analyses. CRC Press. https://doi.org/10.1201/b10157-12 Alvarez, D.A., Petty, J.D., Huckins, J.N., Jones-Lepp, T.L., Getting, D.T., Goddard, J.P., Manahan, S.E., 2004. Development of a passive, in situ, integrative sampler for hydrophilic organic contaminants in aquatic environments. Environ. Toxicol. Chem. 23, 1640–1648. https://doi.org/10.1897/03-603 Bain, P., Williams, M., Kumar, A., 2014. Assessment of multiple hormonal activities in wastewater at different stages of treatment. Environ. Toxicol. Chem. 33, 2297–2307. https://doi.org/10.1002/etc.2676 Balaguer, P., François, F., Comunale, F., Fenet, H., Boussioux, A.M., Pons, M., Nicolas, J.C., Casellas, C., 1999. Reporter cell lines to study the estrogenic effects of xenoestrogens. Sci. Total Environ. 233, 47–56. https://doi.org/10.1016/S0048-9697(99)00178-3 Bečanová, J., Komprdová, K., Vrana, B., Klánová, J., 2016. Annual dynamics of perfluorinated compounds in sediment: A case study in the Morava River in Zlín district, Czech Republic. Chemosphere 151, 225–233. https://doi.org/10.1016/j.chemosphere.2016.02.081 Bengtson Nash, S.M., Schreiber, U., Ralph, P.J., Müller, J.F., 2005. The combined SPE:ToxY-PAM phytotoxicity assay; application and appraisal of a novel biomonitoring tool for the aquatic environment. Biosens. Bioelectron. 20, 1443–1451. https://doi.org/10.1016/j.bios.2004.09.019 Booij, K., Robinson, C.D., Burgess, R.M., Mayer, P., Roberts, C.A., Ahrens, L., Allan, I.J., Brant, J., Jones, L., Kraus, U.R., Larsen, M.M., Lepom, P., Petersen, J., Pröfrock, D., Roose, P., Schäfer, S., Smedes, F., Tixier, C., Vorkamp, K., Whitehouse, P., 2016. Passive Sampling in Regulatory Chemical Monitoring of Nonpolar Organic Compounds in the Aquatic Environment. Environ. Sci. Technol. 50, 3–17. https://doi.org/10.1021/acs.est.5b04050 Booij, K., Smedes, F., 2010. An improved method for estimating in situ sampling rates of nonpolar passive samplers. Environ. Sci. Technol. 44, 6789–6794. Booij, P., Vethaak, D., Leonards, P.E.G., Sjollema, S.B., Kool, J., De Voogt, P., Lamoree, M.H., 2014. Identification of photosynthesis inhibitors of pelagic marine algae using 96-well plate microfractionation for enhanced throughput in effect-directed analysis. Environ. Sci. Technol. 48, 8003–8011. https://doi.org/10.1021/es405428t Brack, W., 2003. Effect-directed analysis: a promising tool for the identification of organic REFERENCES 51 toxicants in complex mixtures? Anal. Bioanal. Chem. 377, 397–407. Brack, W., Aissa, S.A., Backhaus, T., Dulio, V., Escher, B.I., Faust, M., Hilscherova, K., Hollender, J., Hollert, H., Müller, C., Munthe, J., 2019. Effect ‑ based methods are key . The European Collaborative Project SOLUTIONS recommends integrating effect ‑ based methods for diagnosis and monitoring of water quality. Environ. Sci. Eur. 4–9. https://doi.org/10.1186/s12302-019-0192-2 Brack, W., Ait-Aissa, S., Burgess, R.M., Busch, W., Creusot, N., Di Paolo, C., Escher, B.I., Mark Hewitt, L., Hilscherova, K., Hollender, J., Hollert, H., Jonker, W., Kool, J., Lamoree, M., Muschket, M., Neumann, S., Rostkowski, P., Ruttkies, C., Schollee, J., Schymanski, E.L., Schulze, T., Seiler, T.B., Tindall, A.J., De Aragão Umbuzeiro, G., Vrana, B., Krauss, M., 2016. Effect-directed analysis supporting monitoring of aquatic environments - An in-depth overview. Sci. Total Environ. 544, 1073–1118. https://doi.org/10.1016/j.scitotenv.2015.11.102 Brack, W., Altenburger, R., Schüürmann, G., Krauss, M., López Herráez, D., van Gils, J., Slobodnik, J., Munthe, J., Gawlik, B.M., van Wezel, A., Schriks, M., Hollender, J., Tollefsen, K.E., Mekenyan, O., Dimitrov, S., Bunke, D., Cousins, I., Posthuma, L., van den Brink, P.J., López de Alda, M., Barceló, D., Faust, M., Kortenkamp, A., Scrimshaw, M., Ignatova, S., Engelen, G., Massmann, G., Lemkine, G., Teodorovic, I., Walz, K.H., Dulio, V., Jonker, M.T.O., Jäger, F., Chipman, K., Falciani, F., Liska, I., Rooke, D., Zhang, X., Hollert, H., Vrana, B., Hilscherova, K., Kramer, K., Neumann, S., Hammerbacher, R., Backhaus, T., Mack, J., Segner, H., Escher, B., de Aragão Umbuzeiro, G., 2015. The SOLUTIONS project: Challenges and responses for present and future emerging pollutants in land and water resources management. Sci. Total Environ. 503–504, 22–31. https://doi.org/10.1016/j.scitotenv.2014.05.143 Brack, W., Govender, S., Schulze, T., Krauss, M., Hu, M., Muz, M., Hollender, J., Schirmer, K., Schollee, J., Hidasi, A., Slobodnik, J., Rabova, Z., Ait-Aissa, S., Sonavane, M., Carere, M., Lamoree, M., Leonards, P., Tufi, S., Ouyang, X., Schriks, M., Thomas, K., Almeida, A.C. De, Froment, J., Hammers-Wirtz, M., Ahel, M., Koprivica, S., Hollert, H., Seiler, T.-B., Paolo, C. Di, Tindall, A., Spirhanzlova, P., 2013. EDA-EMERGE: an FP7 initial training network to equip the next generation of young scientists with the skills to address the complexity of environmental contamination with emerging pollutants. Environ. Sci. Eur. 25, 1–7. https://doi.org/10.1186/2190-4715-25-18 Brausch, J.M., Rand, G.M., 2011. A review of personal care products in the aquatic environment: Environmental concentrations and toxicity. Chemosphere 82, 1518–1532. https://doi.org/10.1016/j.chemosphere.2010.11.018 Burgess, R.M., Ho, K.T., Brack, W., Lamoree, M., 2013. Effects-directed analysis (EDA) and toxicity identification evaluation (TIE): Complementary but different approaches for diagnosing causes of environmental toxicity. Environ. Toxicol. Chem. 32, 1935–1945. https://doi.org/10.1002/etc.2299 Connon, R.E., Geist, J., Werner, I., 2012. Effect-based tools for monitoring and predicting the ecotoxicological effects of chemicals in the aquatic environment. Sensors (Switzerland) 12, 12741–12771. https://doi.org/10.3390/s120912741 Creusot, N., Aït-Aïssa, S., Tapie, N., Pardon, P., Brion, F., Sanchez, W., Thybaud, E., Porcher, J.M., Budzinski, H., 2014. Identification of synthetic steroids in river water downstream from pharmaceutical manufacture discharges based on a bioanalytical approach and passive sampling. Environ. Sci. Technol. 48, 3649–3657. https://doi.org/10.1021/es405313r Creusot, N., Dévier, M.H., Budzinski, H., Aït-Aïssa, S., 2016. Evaluation of an extraction method for a mixture of endocrine disrupters in sediment using chemical and in vitro biological analyses. Environ. Sci. Pollut. Res. 23, 10349–10360. https://doi.org/10.1007/s11356- 016-6062-1 Creusot, N., Tapie, N., Piccini, B., Balaguer, P., Porcher, J.M., Budzinski, H., Aït-Aïssa, S., 2013. Distribution of steroid- and dioxin-like activities between sediments, POCIS and SPMD in a French river subject to mixed pressures. Environ. Sci. Pollut. Res. 20, 2784–2794. https://doi.org/10.1007/s11356-012-1452-5 Demirpence, E., Duchesne, M.-J., Badia, E., Gagne, D., Pons, M., 1993. MVLN Cells: A bioluminescent MCF-7-derived cell line to study the modulation of estrogenic activity. J. Steroid Biochem. Mol. Biol. 46, 355–364. https://doi.org/10.1016/0960-0760(93)90225-L 52 Dimpe, K.M., Nomngongo, P.N., 2016. Current sample preparation methodologies for analysis of emerging pollutants in different environmental matrices. TrAC - Trends Anal. Chem. 82, 199–207. https://doi.org/10.1016/j.trac.2016.05.023 DIRECTIVE 2000/60/EC OF THE EUROPEAN PARLIAMENT AND OF THE COUNCIL of 23 October 2000 establishing a framework for Community action in the field of water policy., 2000. . Off. J. Eur. Union L327, 1e77. DIRECTIVE 2008/105/EC OF THE EUROPEAN PARLIAMENT AND OF THE COUNCIL of 16 December 2008 on environmental quality standards in the field of water policy, amending and subsequently repealing Council Directives 82/176/EEC, 83/513/EEC, 84/156/EEC, 84/491/EEC, 2008. . Off. J. Eur. Union L 348/84. DIRECTIVE 2013/39/EU OF THE EUROPEAN PARLIAMENT AND OF THE COUNCIL of 12 August 2013 amending Directives 2000/60/EC and 2008/105/EC as regards priority substances in the field of water policy, 2013. . Off. J. Eur. Union L 226/15. Dizer, H., Wittekindt, E., Fischer, B., Hansen, P.D., 2002. The cytotoxic and genotoxic potential of surface water and wastewater effluents as determined by bioluminescence, umu-assays and selected biomarkers. Chemosphere 46, 225–233. https://doi.org/10.1016/S0045- 6535(01)00062-5 Du, B., Lofton, J.M., Peter, K.T., Gipe, A.D., James, C.A., McIntyre, J.K., Scholz, N.L., Baker, J.E., Kolodziej, E.P., 2017. Development of suspect and non-target screening methods for detection of organic contaminants in highway runoff and fish tissue with high-resolution time-of-flight mass spectrometry. Environ. Sci. Process. Impacts 19, 1185–1196. https://doi.org/10.1039/c7em00243b Dulio, V., van Bavel, B., Brorström-Lundén, E., Harmsen, J., Hollender, J., Schlabach, M., Slobodnik, J., Thomas, K., Koschorreck, J., 2018. Emerging pollutants in the EU: 10 years of NORMAN in support of environmental policies and regulations. Environ. Sci. Eur. 30, 5. https://doi.org/10.1186/s12302-018-0135-3 Dulio, V., Von der Ohe, P.C. (Eds.), 2013. NORMAN prioritisation framework for emerging substances. [WWW Document]. URL http://www.norman- network.net/sites/default/files/files/Publications/NORMAN_prioritisation_Manual_15 April2013_final for website-f.pdf (accessed 10.16.15). Ellman, G.L., Courtney, K.D., Andres, V., Featherstone, R.M., 1961. A new and rapid colorimetric determination of acetylcholinesterase activity. Biochem. Pharmacol. 7, 88–95. Emelogu, E.S., Pollard, P., Dymond, P., Robinson, C.D., Webster, L., McKenzie, C., Dobson, J., Bresnan, E., Moffat, C.F., 2013. Occurrence and potential combined toxicity of dissolved organic contaminants in the Forth estuary and Firth of Forth, Scotland assessed using passive samplers and an algal toxicity test. Sci. Total Environ. 461–462, 230–239. https://doi.org/10.1016/j.scitotenv.2013.05.011 Escher, B.I., Aїt-Aїssa, S., Behnisch, P.A., Brack, W., Brion, F., Brouwer, A., Buchinger, S., Crawford, S.E., Du Pasquier, D., Hamers, T., Hettwer, K., Hilscherová, K., Hollert, H., Kase, R., Kienle, C., Tindall, A.J., Tuerk, J., van der Oost, R., Vermeirssen, E., Neale, P.A., 2018. Effect-based trigger values for in vitro and in vivo bioassays performed on surface water extracts supporting the environmental quality standards (EQS) of the European Water Framework Directive. Sci. Total Environ. 628–629, 748–765. https://doi.org/10.1016/j.scitotenv.2018.01.340 Escher, B.I., Pronk, W., Suter, M.J.F., Maurer, M., 2006. Monitoring the removal efficiency of pharmaceuticals and hormones in different treatment processes of source-separated urine with bioassays. Environ. Sci. Technol. 40, 5095–5101. Escher, Beate I, Quayle, P., Muller, R., Schreiber, U., Mueller, J.F., 2006. Passive sampling of herbicides combined with effect analysis in algae using a novel high-throughput phytotoxicity assay (Maxi-Imaging-PAM). J. Environ. Monit. 8, 456–464. https://doi.org/10.1039/b517512g European Commission, 2003. Document on Risk Assessment. Tech. Guid. Doc. Ris Assess. Part II 337. https://doi.org/10.1002/mp.12308 European Environment Agency, 2018. European waters — Assessment of status and pressures 2018. https://doi.org/0.2800/303664 REFERENCES 53 Fairbairn, D.J., Elliott, S.M., Kiesling, R.L., Schoenfuss, H.L., Ferrey, M.L., Westerhoff, B.M., 2018. Contaminants of emerging concern in urban stormwater: Spatiotemporal patterns and removal by iron-enhanced sand filters (IESFs). Water Res. 145, 332–345. https://doi.org/10.1016/j.watres.2018.08.020 Ferrer, I., Thurman, E.M., 2010. Identification of a New Antidepressant and its Glucuronide Metabolite in Water Samples Using Liquid Chromatography/Quadrupole Time-of-Flight Mass Spectrometry. Anal. Chem. 82, 8161–8168. https://doi.org/10.1021/ac1014645 Fini, J.B., Le Mével, S., Turque, N., Palmier, K., Zalko, D., Cravedi, J.P., Demeneix, B.A., 2007. An in vivo multiwell-based fluorescent screen for monitoring vertebrate thyroid hormone disruption. Environ. Sci. Technol. 41, 5908–5914. Fontaine, J., Reuse, J., 1980. The effects of substituted benzamides on frog rectus abdominis. Eur. J. Pharmacol. 68, 55–60. https://doi.org/10.1016/0014-2999(80)90060-6 Fontanals, N., Marcé, R.M., Borrull, F., 2010. Núria Fontanals, Rosa Maria Marcé, Francesc Borrull_2010 6, 199–213. https://doi.org/10.2436/20.7010.01.97 Galgani, F., Bocquene, G., 1991. Semi-automated colorimetric and enzymaic assays for aquatic organisms using microplate readers. Water Res. 25, 147–150. García-Córcoles, M.T., Rodríguez-Gómez, R., de Alarcón-Gómez, B., Çipa, M., Martín-Pozo, L., Kauffmann, J.M., Zafra-Gómez, A., 2019. Chromatographic Methods for the Determination of Emerging Contaminants in Natural Water and Wastewater Samples: A Review. Crit. Rev. Anal. Chem. 49, 160–186. https://doi.org/10.1080/10408347.2018.1496010 Geissen, V., Mol, H., Klumpp, E., Umlauf, G., Nadal, M., van der Ploeg, M., van de Zee, S.E. a. T.M., Ritsema, C.J., 2015. Emerging pollutants in the environment: A challenge for water resource management. Int. Soil Water Conserv. Res. 3, 57–65. https://doi.org/10.1016/j.iswcr.2015.03.002 Ginebreda, A., Kuzmanovic, M., Guasch, H., de Alda, M.L., López-Doval, J.C., Muñoz, I., Ricart, M., Romaní, A.M., Sabater, S., Barceló, D., 2014. Assessment of multi-chemical pollution in aquatic ecosystems using toxic units: Compound prioritization, mixture characterization and relationships with biological descriptors. Sci. Total Environ. 468–469, 715–723. https://doi.org/10.1016/j.scitotenv.2013.08.086 Gómez, M.J., Gómez-Ramos, M.M., Agüera, a., Mezcua, M., Herrera, S., Fernández-Alba, a. R., 2009. A new gas chromatography/mass spectrometry method for the simultaneous analysis of target and non-target organic contaminants in waters. J. Chromatogr. A 1216, 4071–4082. https://doi.org/10.1016/j.chroma.2009.02.085 Górecki, T., Namienik, J., 2002. Passive sampling. TrAC - Trends Anal. Chem. 21, 276–291. https://doi.org/10.1016/S0165-9936(02)00407-7 Greaves, A.K., Letcher, R.J., 2017. A Review of Organophosphate Esters in the Environment from Biological Effects to Distribution and Fate. Bull. Environ. Contam. Toxicol. 98, 2–7. https://doi.org/10.1007/s00128-016-1898-0 Griffith, D.R., Kido Soule, M.C., Matsufuji, H., Eglinton, T.I., Kujawinski, E.B., Gschwend, P.M., 2014. Measuring free, conjugated, and halogenated estrogens in secondary treated wastewater effluent. Environ. Sci. Technol. 48, 2569–2578. https://doi.org/10.1021/es402809u Guijarro, C., Fuchs, K., Bohrn, U., Stütz, E., Wölfl, S., 2015. Simultaneous detection of multiple bioactive pollutants using a multiparametric biochip for water quality monitoring. Biosens. Bioelectron. 72, 71–79. https://doi.org/10.1016/j.bios.2015.04.092 Hamers, T., Leonards, P.E.G., Legler, J., Dick Vethaak, A., Schipper, C.A., 2010. Toxicity profiling: An integrated effect-based tool for site-specific sediment quality assessment. Integr. Environ. Assess. Manag. 6, 761–773. https://doi.org/10.1002/ieam.75 Harman, C., Allan, I.J., Vermeirssen, E.L.M., 2012. Calibration and use of the polar organic chemical integrative sampler-a critical review. Environ. Toxicol. Chem. 31, 2724–2738. https://doi.org/10.1002/etc.2011 Hollender, J., van Bavel, B., Dulio, V., Farmen, E., Furtmann, K., Koschorreck, J., Kunkel, U., Krauss, M., Munthe, J., Schlabach, M., Slobodnik, J., Stroomberg, G., Ternes, T., Thomaidis, N.S., Togola, A., Tornero, V., 2019. High resolution mass spectrometry-based non-target screening can support regulatory environmental monitoring and chemicals management. Environ. Sci. Eur. 31. https://doi.org/10.1186/s12302-019-0225-x 54 Houtman, C.J., 2010. Emerging contaminants in surface waters and their relevance for the production of drinking water in Europe. J. Integr. Environ. Sci. 7, 271–295. https://doi.org/10.1080/1943815X.2010.511648 Huckins, J.N., Manuweera, G.K., Petty, J.D., Mackay, D., Lebo, J.A., 1993. Lipid-containing semipermeable membrane devices for monitoring organic contaminants in water. Environ. Sci. Technol. 27, 2489–2496. Hug, C., Ulrich, N., Schulze, T., Brack, W., Krauss, M., 2014. Identification of novel micropollutants in wastewater by a combination of suspect and nontarget screening. Environ. Pollut. 184, 25–32. https://doi.org/10.1016/j.envpol.2013.07.048 Ibáñez, M., Sancho, J. V., Hernández, F., McMillan, D., Rao, R., 2008. Rapid non-target screening of organic pollutants in water by ultraperformance liquid chromatography coupled to time- of-light mass spectrometry. TrAC - Trends Anal. Chem. 27, 481–489. https://doi.org/10.1016/j.trac.2008.03.007 Jálová, V., Jarošová, B., Bláha, L., Giesy, J.P., Ocelka, T., Grabic, R., Jurčíková, J., Vrana, B., Hilscherová, K., 2013. Estrogen-, androgen- and aryl hydrocarbon receptor mediated activities in passive and composite samples from municipal waste and surface waters. Environ. Int. 59, 372–383. https://doi.org/10.1016/j.envint.2013.06.024 Jalova, V., Jarosova, B., Blaha, L., Giesy, J.P.P., Ocelka, T., Grabic, R., Jurcikova, J., Vrana, B., Hilscherova, K., Jálová, V., Jarošová, B., Bláha, L., Giesy, J.P.P., Ocelka, T., Grabic, R., Jurčíková, J., Vrana, B., Hilscherová, K., 2013. Estrogen-, androgen- and aryl hydrocarbon receptor mediated activities in passive and composite samples from municipal waste and surface waters. Environ. Int. 59, 372–83. https://doi.org/10.1016/j.envint.2013.06.024 Jarosova, B., Blaha, L., Vrana, B., Randak, T., Grabic, R., Giesy, J.P., Hilscherova, K., 2012. Changes in concentrations of hydrophilic organic contaminants and of endocrine-disrupting potential downstream of small communities located adjacent to headwaters. Environ. Int. 45, 22–31. https://doi.org/10.1016/j.envint.2012.04.001 Jones, L., Ronan, J., McHugh, B., McGovern, E., Regan, F., 2015. Emerging priority substances in the aquatic environment: a role for passive sampling in supporting WFD monitoring and compliance. Anal. Methods 7, 7976–7984. https://doi.org/10.1039/C5AY01059D Kinani, S., Bouchonnet, S., Creusot, N., Bourcier, S., Balaguer, P., Porcher, J.-M., Aït-Aïssa, S., 2010. Bioanalytical characterisation of multiple endocrine- and dioxin-like activities in sediments from reference and impacted small rivers. Environ. Pollut. 158, 74–83. https://doi.org/10.1016/j.envpol.2009.07.041 Koenig, M., Escher, B.I., Neale, P.A., Krauss, M., Hilscherova, K., Novak, J., Teodorovic, I., Schulze, T., Seidensticker, S., Kamal Hashmi, M.A., Ahlheim, J., Brack, W., 2017. Impact of untreated wastewater on a major European river evaluated with a combination of in??vitro bioassays and chemical analysis. Environ. Pollut. 220, 1220–1230. https://doi.org/10.1016/j.envpol.2016.11.011 Köhler, A., Hellweg, S., Escher, B.I., Hungerbühler, K., 2006. Organic pollutant removal versus toxicity reduction in industrial wastewater treatment: The example of wastewater from fluorescent whitening agent production. Environ. Sci. Technol. 40, 3395–3401. https://doi.org/10.1021/es060555f König, M., Escher, B.I., Neale, P.A., Krauss, M., Hilscherová, K., Novák, J., Teodorović, I., Schulze, T., Seidensticker, S., Kamal Hashmi, M.A., Ahlheim, J., Brack, W., 2017. Impact of untreated wastewater on a major European river evaluated with a combination of in vitro bioassays and chemical analysis. Environ. Pollut. 220, 1220–1230. https://doi.org/10.1016/j.envpol.2016.11.011 Krauss, M., Singer, H., Hollender, J., 2010. LC-high resolution MS in environmental analysis: From target screening to the identification of unknowns. Anal. Bioanal. Chem. 397, 943–951. https://doi.org/10.1007/s00216-010-3608-9 Křesinová, Z., Petrů, K., Lhotský, O., Rodsand, T., Cajthaml, T., 2016. Passive sampling of pharmaceuticals and personal care products in aquatic environments. Eur. J. Environ. Sci. 6, 43–56. https://doi.org/10.14712/23361964.2016.8 Kuzmanovic, M., Ginebreda, A., Petrovic, M., Barcelo, D., 2014. Risk assessment based prioritization of 200 organic micropollutants in 4 Iberian rivers. Sci. Total Environ. 504, REFERENCES 55 289–299. https://doi.org/10.1016/j.scitotenv.2014.06.056 Letsinger, S., Kay, P., 2019. Comparison of Prioritisation Schemes for Human Pharmaceuticals in the Aquatic Environment. Environ. Sci. Pollut. Res. 26, 3479–3491. https://doi.org/10.1007/s11356-018-3834-9 Liang, X., Adamovsky, O., Souders, C.L., Martyniuk, C.J., 2019. Biological effects of the benzotriazole ultraviolet stabilizers UV-234 and UV-320 in early-staged zebrafish (Danio rerio). Environ. Pollut. 245, 272–281. https://doi.org/10.1016/j.envpol.2018.10.130 Linares, V., Bellés, M., Domingo, J.L., 2015. Human exposure to PBDE and critical evaluation of health hazards. Arch. Toxicol. 89, 335–356. https://doi.org/10.1007/s00204-015-1457-1 Loos, R., Carvalho, R., António, D.C., Comero, S., Locoro, G., Tavazzi, S., Paracchini, B., Ghiani, M., Lettieri, T., Blaha, L., Jarosova, B., Voorspoels, S., Servaes, K., Haglund, P., Fick, J., Lindberg, R.H., Schwesig, D., Gawlik, B.M., 2013. EU-wide monitoring survey on emerging polar organic contaminants in wastewater treatment plant effluents. Water Res. 47, 6475–6487. https://doi.org/10.1016/j.watres.2013.08.024 Magdaleno, A., Juárez, Á.B., Dragani, V., Saenz, M.E., Paz, M., Moretton, J., 2014. Ecotoxicological and Genotoxic Evaluation of Buenos Aires City ( Argentina ) Hospital Wastewater 2014. Martín-Pozo, L., de Alarcón-Gómez, B., Rodríguez-Gómez, R., García-Córcoles, M.T., Çipa, M., Zafra- Gómez, A., 2019. Analytical methods for the determination of emerging contaminants in sewage sludge samples. A review. Talanta 192, 508–533. https://doi.org/10.1016/j.talanta.2018.09.056 Maselli, B.D.S., Luna, L.A.V. De, Palmeira, J.D.O., Babosa, S., Beijo, L.A., Umbuzeiro, G.D.A., Kummrow, F., 2013. Ecotoxicity of raw and treated effluents generated by a veterinary medicine industry. Ambient. e Agua - An Interdiscip. J. Appl. Sci. 8, 795–804. https://doi.org/10.4136/ambi-agua.1121 Miège, C., Gabet, V., Coquery, M., Karolak, S., Jugan, M.L., Oziol, L., Levi, Y., Chevreuil, M., 2009. Evaluation of estrogenic disrupting potency in aquatic environments and urban wastewaters by combining chemical and biological analysis. TrAC - Trends Anal. Chem. 28, 186–195. https://doi.org/10.1016/j.trac.2008.11.007 Milinovic, Z., 2013. PRELIMINARY RESULTS Of the 2013 Census of Population, Households and Dwellings in Bosnia and Herzegovina [WWW Document]. Agency Stat. Bosnia Herzegovina. URL http://www.bhas.ba/obavjestenja/Preliminarni_rezultati_bos.pdf (accessed 10.10.17). Moschet, C., Vermeirssen, E.L.M., Seiz, R., Pfefferli, H., Hollender, J., 2014. Picogram per liter detections of pyrethroids and organophosphates in surface waters using passive sampling. Water Res. 66, 411–422. https://doi.org/10.1016/j.watres.2014.08.032 Nagy, S.R., Sanborn, J.R., Hammock, B.D., Denison, M.S., 2002. Development of a green fluorescent protein-based cell bioassay for the rapid and inexpensive detection and characterization of Ah receptor agonists. Toxicol. Sci. 65, 200–210. https://doi.org/10.1093/toxsci/65.2.200 Neale, P.A., Ait-Aissa, S., Brack, W., Creusot, N., Denison, M.S., Deutschmann, B., Hilscherova, K., Hollert, H., Krauss, M., Novak, J., Schulze, T., Seiler, T.B., Serra, H., Shao, Y., Escher, B.I., 2015. Linking in Vitro Effects and Detected Organic Micropollutants in Surface Water Using Mixture-Toxicity Modeling. Environ. Sci. Technol. 49, 14614–14624. https://doi.org/10.1021/acs.est.5b04083 Novák, J., Vrana, B., Rusina, T., Okonski, K., Grabic, R., Neale, P.A., Escher, B.I., Macová, M., Ait- Aissa, S., Creusot, N., Allan, I., Hilscherová, K., 2018. Effect-based monitoring of the Danube River using mobile passive sampling. Sci. Total Environ. 636, 1608–1619. https://doi.org/10.1016/j.scitotenv.2018.02.201 OECD, 2013. Test No. 236: Fish Embryo Acute Toxicity (FET) Test [WWW Document]. https://doi.org/http://dx.doi.org/10.1787/9789264203709-en OECD, 2011. Test No. 201: Freshwater Alga and Cyanobacteria, Growth Inhibition Test [WWW Document]. https://doi.org/http://dx.doi.org/10.1787/9789264069923-en OECD, 2002. Hazard assessment of perfluorooctane sulfonate (PFOS) and its salts.Co-operation on existing chemicals. Dev. 362. Ohoro, Adeniji, Okoh, Okoh, 2019. Distribution and Chemical Analysis of Pharmaceuticals and Personal Care Products (PPCPs) in the Environmental Systems: A Review. Int. J. Environ. 56 Res. Public Health 16, 3026. https://doi.org/10.3390/ijerph16173026 Osman, A.G.M., AbouelFadl, K.Y., Krüger, A., Kloas, W., 2015. Screening of multiple hormonal activities in water and sediment from the river Nile, Egypt, using in vitro bioassay and gonadal histology. Environ. Monit. Assess. 187, 317. https://doi.org/10.1007/s10661-015- 4553-z Pal, A., Gin, K.Y.-H., Lin, A.Y.-C., Reinhard, M., 2010. Impacts of emerging organic contaminants on freshwater resources: review of recent occurrences, sources, fate and effects. Sci. Total Environ. 408, 6062–9. https://doi.org/10.1016/j.scitotenv.2010.09.026 Patel, M., Kumar, R., Kishor, K., Mlsna, T., Pittman, C.U., Mohan, D., 2019. Pharmaceuticals of emerging concern in aquatic systems: Chemistry, occurrence, effects, and removal methods. Chem. Rev. 119, 3510–3673. https://doi.org/10.1021/acs.chemrev.8b00299 Petrovic, M., Barceló, D., Diaz, A., Ventura, F., 2003. Low nanogram per liter determination of halogenated nonylphenols, nonylphenol carboxylates, and their non-halogenated precursors in water and sludge by liquid chromatography electrospray tandem mass spectrometry. J. Am. Soc. Mass Spectrom. 14, 516–527. https://doi.org/10.1016/S1044- 0305(03)00139-9 Pietrzak, D., Kania, J., Malina, G., Kmiecik, E., Wątor, K., 2019. Pesticides from the EU First and Second Watch Lists in the Water Environment. Clean - Soil, Air, Water 47. https://doi.org/10.1002/clen.201800376 Posthuma, L., Altenburger, R., Backhaus, T., Kortenkamp, A., Müller, C., Focks, A., de Zwart, D., Brack, W., 2019. Improved component-based methods for mixture risk assessment are key to characterize complex chemical pollution in surface waters. Environ. Sci. Eur. 31. https://doi.org/10.1186/s12302-019-0246-5 Quednow, K., Püttmann, W., 2009. Temporal concentration changes of DEET, TCEP, terbutryn, and nonylphenols in freshwater streams of Hesse, Germany: Possible influence of mandatory regulations and voluntary environmental agreements. Environ. Sci. Pollut. Res. 16, 630–640. https://doi.org/10.1007/s11356-009-0169-6 Rainieri, S., Conlledo, N., Langerholc, T., Madorran, E., Sala, M., Barranco, A., 2017. Toxic effects of perfluorinated compounds at human cellular level and on a model vertebrate. Food Chem. Toxicol. 104, 14–25. https://doi.org/10.1016/j.fct.2017.02.041 Reichenberg, F., Mayer, P., 2006. Two complementary sides of bioavailability: accessibility and chemical activity of organic contaminants in sediments and soils. Environ. Toxicol. Chem. 25, 1239–45. Ren, Y., Geng, J., Li, F., Ren, H., Ding, L., Xu, K., 2016. The oxidative stress in the liver of Carassius auratus exposed to acesulfame and its UV irradiance products. Sci. Total Environ. 571, 755–762. https://doi.org/10.1016/j.scitotenv.2016.07.047 Richardson, S.D., Ternes, T.A., 2018. Water Analysis: Emerging Contaminants and Current Issues. Anal. Chem. 90, 398–428. https://doi.org/10.1021/acs.analchem.7b04577 Riva, F., Zuccato, E., Davoli, E., Fattore, E., Castiglioni, S., 2019. Risk assessment of a mixture of emerging contaminants in surface water in a highly urbanized area in Italy. J. Hazard. Mater. 361, 103–110. https://doi.org/10.1016/j.jhazmat.2018.07.099 Rojickova, R., Dvorakova, D., 1998. The Use of Miniaturized Algal Bioassays in Comparison to the Standard Flask Assay. Env. Toxicol Water Qual 13, 235–241. Rusina, T.P., Smedes, F., Koblizkova, M., Klanova, J., 2010. Calibration of silicone rubber passive samplers: Experimental and modeled relations between sampling rate and compound properties. Environ. Sci. Technol. 44, 362–367. Rykowska, I., Wasiak, W., 2015. Research trends on emerging environment pollutants - A review. Open Chem. 13, 1353–1370. https://doi.org/10.1515/chem-2015-0151 Salimi, M., Esrafili, A., Gholami, M., Jonidi Jafari, A., Rezaei Kalantary, R., Farzadkia, M., Kermani, M., Sobhi, H.R., 2017. Contaminants of emerging concern: a review of new approach in AOP technologies. Environ. Monit. Assess. 189. https://doi.org/10.1007/s10661-017-6097-x Saucedo-Vence, K., Elizalde-Velázquez, A., Dublán-García, O., Galar-Martínez, M., Islas-Flores, H., SanJuan-Reyes, N., García-Medina, S., Hernández-Navarro, M.D., Gómez-Oliván, L.M., 2017. Toxicological hazard induced by sucralose to environmentally relevant concentrations in common carp (Cyprinus carpio). Sci. Total Environ. 575, 347–357. REFERENCES 57 https://doi.org/10.1016/j.scitotenv.2016.09.230 Scheurer, M., Brauch, H.J., Lange, F.T., 2009. Analysis and occurrence of seven artificial sweeteners in German waste water and surface water and in soil aquifer treatment (SAT). Anal. Bioanal. Chem. 394, 1585–1594. https://doi.org/10.1007/s00216-009-2881-y Schirmer, K., Chan, A.G.J., Greenberg, B.M., Dixon, D.G., Bols, N.C., 1998. Ability of 16 priority PAHs to be photocytotoxic to a cell line from the rainbow trout gill. Toxicology 127, 143–155. https://doi.org/10.1016/S0300-483X(98)00031-6 Schriks, M., Van Leerdam, J.A., Van Der Linden, S.C., Van Der Burg, B., Van Wezel, A.P., De Voogt, P., 2010. High-resolution mass spectrometric identification and quantification of glucocorticoid compounds in various wastewaters in the Netherlands. Environ. Sci. Technol. 44, 4766–4774. https://doi.org/10.1021/es100013x Schulze, Tobias, Ahel, M., Ahlheim, J., Aït-Aïssa, S., Brion, F., Di Paolo, C., Froment, J., Hidasi, A.O., Hollender, J., Hollert, H., Hu, M., Kloß, A., Koprivica, S., Krauss, M., Muz, M., Oswald, P., Petre, M., Schollée, J.E., Seiler, T.-B., Shao, Y., Slobodnik, J., Sonavane, M., Suter, M.J.-F., Tollefsen, K.E., Tousova, Z., Walz, K.-H., Brack, W., 2017. Assessment of a novel device for onsite integrative large-volume solid phase extraction of water samples to enable a comprehensive chemical and effect-based analysis. Sci. Total Environ. 581–582, 350–358. https://doi.org/10.1016/j.scitotenv.2016.12.140 Schulze, T., Ahel, M., Ahlheim, J., Aït-Aïssa, S., Brion, F., Di Paolo, C., Froment, J., Hidasi, A.O., Hollender, J., Hollert, H., Hu, M., Kloß, A., Koprivica, S., Krauss, M., Muz, M., Oswald, P., Petre, M., Schollée, J.E., Seiler, T.-B., Shao, Y., Slobodnik, J., Sonavane, M., Suter, M.J.-F., Tollefsen, K.E., Tousova, Z., Walz, K.-H., Brack, W., 2017. Assessment of a novel device for onsite integrative large-volume solid phase extraction of water samples to enable a comprehensive chemical and effect-based analysis. Sci. Total Environ. 581–582. https://doi.org/10.1016/j.scitotenv.2016.12.140 Schymanski, E.L., Singer, H.P., Longrée, P., Loos, M., Ruff, M., Stravs, M. a., Ripollés Vidal, C., Hollender, J., 2014. Strategies to characterize polar organic contamination in wastewater: Exploring the capability of high resolution mass spectrometry. Environ. Sci. Technol. 48, 1811–1818. https://doi.org/10.1021/es4044374 Schymanski, E.L., Singer, H.P., Slobodnik, J., Ipolyi, I.M., Oswald, P., Krauss, M., Schulze, T., Haglund, P., Letzel, T., Grosse, S., Thomaidis, N.S., Bletsou, A., Zwiener, C., Ibáñez, M., Portolés, T., de Boer, R., Reid, M.J., Onghena, M., Kunkel, U., Schulz, W., Guillon, A., Noyon, N., Leroy, G., Bados, P., Bogialli, S., Stipaničev, D., Rostkowski, P., Hollender, J., 2015. Non-target screening with high-resolution mass spectrometry: critical review using a collaborative trial on water analysis. Anal. Bioanal. Chem. 6237–6255. https://doi.org/10.1007/s00216- 015-8681-7 Shi, Z.Q., Liu, Y.S., Xiong, Q., Cai, W.W., Ying, G.G., 2019. Occurrence, toxicity and transformation of six typical benzotriazoles in the environment: A review. Sci. Total Environ. 661, 407–421. https://doi.org/10.1016/j.scitotenv.2019.01.138 Sima, L., Amador, J., Silva, A.D.K., Miller, S.M., Morse, A.N., Pellegrin, M.-L., Rock, C., Wells, M.J.M., 2014. Emerging Pollutants – Part I: Occurrence, Fate and Transport. Water Environ. Res. 86, 1994–2035. https://doi.org/10.2175/106143014X14031280668731 Slobodnik, J., Mrafkova, L., Carere, M., Ferrara, F., Pennelli, B., Schüürmann, G., von der Ohe, P.C., 2012. Identification of river basin specific pollutants and derivation of environmental quality standards: A case study in the Slovak Republic. TrAC Trends Anal. Chem. 41, 133– 145. https://doi.org/10.1016/j.trac.2012.08.008 Smital, T., Terzić, S., Lončar, J., Senta, I., Žaja, R., Popović, M., Mikac, I., Tollefsen, K.-E., Thomas, K. V., Ahel, M., 2013. Prioritisation of organic contaminants in a river basin using chemical analyses and bioassays. Environ. Sci. Pollut. Res. 20, 1384–1395. https://doi.org/10.1007/s11356-012-1059-x Sonneveld, E., Jansen, H.J., Riteco, J.A.C., Brouwer, A., van der Burg, B., 2005. Development of androgen- and estrogen-responsive bioassays members of a panel of human cell line-based highly selective steroid-responsive bioassays. Toxicol. Sci. 83, 136–148. https://doi.org/10.1093/toxsci/kfi005 Tang, J.Y.M., Escher, B.I., 2014. Realistic environmental mixtures of micropollutants in surface, 58 drinking, and recycled water: Herbicides dominate the mixture toxicity toward algae. Environ. Toxicol. Chem. 33, 1427–1436. https://doi.org/10.1002/etc.2580 Tangtian, H., Bo, L., Wenhua, L., Shin, P.K.S., Wu, R.S.S., 2012. Estrogenic potential of benzotriazole on marine medaka (Oryzias melastigma). Ecotoxicol. Environ. Saf. 80, 327–332. https://doi.org/10.1016/j.ecoenv.2012.03.020 Taylor, A.C., Fones, G.R., Vrana, B., Mills, G.A., 2019. Applications for Passive Sampling of Hydrophobic Organic Contaminants in Water—A Review. Crit. Rev. Anal. Chem. 0, 1–35. https://doi.org/10.1080/10408347.2019.1675043 Thomaidis, N.S., Asimakopoulos, A.G., Bletsou, A.A., 2012. Emerging contaminants: A tutorial mini-review. Glob. Nest J. 14, 72–79. https://doi.org/10.30955/gnj.000823 Tijani, J.O., Fatoba, O.O., Babajide, O.O., Petrik, L.F., 2016. Pharmaceuticals, endocrine disruptors, personal care products, nanomaterials and perfluorinated pollutants: a review. Environ. Chem. Lett. 14, 27–49. https://doi.org/10.1007/s10311-015-0537-z Tousova, Z., Oswald, P., Slobodnik, J., Blaha, L., Muz, M., Hu, M., Brack, W., Krauss, M., Di Paolo, C., Tarcai, Z., Seiler, T.-B., Hollert, H., Koprivica, S., Ahel, M., Schollée, J.E., Hollender, J., Suter, M.J.-F., Hidasi, A.O., Schirmer, K., Sonavane, M., Ait-Aissa, S., Creusot, N., Brion, F., Froment, J., Almeida, A.C., Thomas, K., Tollefsen, K.E., Tufi, S., Ouyang, X., Leonards, P., Lamoree, M., Torrens, V.O., Kolkman, A., Schriks, M., Spirhanzlova, P., Tindall, A., Schulze, T., 2017. European demonstration program on the effect-based and chemical identification and monitoring of organic pollutants in European surface waters. Sci. Total Environ. 601–602, 1849–1868. https://doi.org/10.1016/j.scitotenv.2017.06.032 US-EPA, 2015. ECOTOX Database [WWW Document]. URL http://cfpub.epa.gov/ecotox/ (accessed 4.20.15). US-EPA, 2012. ECOSAR 1.11 [WWW Document]. URL http://www.epa.gov/oppt/newchems/tools/21ecosar.htm Venzmer, H.V., 2008. Forum Altbausanierung 2. Biofilme und funktionale Baustoffoberflächen.: 8. Dahlberg-Kolloquium vom 25. bis 26. September 2008 im Zeughaus Wismar. Fraunhofer IRB Verlag. Vera-Candioti, L., Gil García, M.D., Martínez Galera, M., Goicoechea, H.C., 2008. Chemometric assisted solid-phase microextraction for the determination of anti-inflammatory and antiepileptic drugs in river water by liquid chromatography-diode array detection. J. Chromatogr. A 1211, 22–32. https://doi.org/10.1016/j.chroma.2008.09.093 Vermeirssen, E.L.M., Hollender, J., Bramaz, N., Van Der Voet, J., Escher, B.I., 2010. Linking toxicity in algal and bacterial assays with chemical analysis in passive samplers deployed in 21 treated sewage effluents. Environ. Toxicol. Chem. 29, 2575–2582. https://doi.org/10.1002/etc.311 Von der Ohe, P.C., Dulio, V., Slobodnik, J., De Deckere, E., Kühne, R., Ebert, R.U., Ginebreda, A., De Cooman, W., Schüürmann, G., Brack, W., 2011. A new risk assessment approach for the prioritization of 500 classical and emerging organic microcontaminants as potential river basin specific pollutants under the European Water Framework Directive. Sci. Total Environ. 409, 2064–2077. https://doi.org/10.1016/j.scitotenv.2011.01.054 Vonberg, D., Vanderborght, J., Cremer, N., Pütz, T., Herbst, M., Vereecken, H., 2014. 20 years of long-term atrazine monitoring in a shallow aquifer in western Germany. Water Res. 50. https://doi.org/10.1016/j.watres.2013.10.032 Vrana, B., Allan, I.J., Greenwood, R., Mills, G.A., Dominiak, E., Svensson, K., Knutsson, J., Morrison, G., 2005. Passive sampling techniques for monitoring pollutants in water. TrAC - Trends Anal. Chem. 24, 845–868. Vrana, Branislav, Klučárová, V., Benická, E., Abou-Mrad, N., Amdany, R., Horáková, S., Draxler, A., Humer, F., Gans, O., 2014. Passive sampling: An effective method for monitoring seasonal and spatial variability of dissolved hydrophobic organic contaminants and metals in the Danube river. Environ. Pollut. 184, 101–112. https://doi.org/10.1016/j.envpol.2013.08.018 Vrana, B, Klučárová, V., Benická, E., Abou-Mrad, N., Amdany, R., Horáková, S., Draxler, A., Humer, F., Gans, O., 2014. Passive sampling: An effective method for monitoring seasonal and spatial variability of dissolved hydrophobic organic contaminants and metals in the REFERENCES 59 Danube river. Environ. Pollut. 184, 101–112. Waldmeier, P.C., Baumann, P.A., Wicki, P., Feldtrauer, J.-J., Stierlin, C., Schmutz, M., 1995. Similar potency of carbamazepine, oxcarbazepine, and lamotrigine in inhibiting the release of glutamate and other neurotransmitters. Neurology 45, 1907 LP – 1913. https://doi.org/10.1212/WNL.45.10.1907 Weissinger, R.H., Blackwell, B.R., Keteles, K., Battaglin, W.A., Bradley, P.M., 2018. Bioactive contaminants of emerging concern in National Park waters of the northern Colorado Plateau, USA. Sci. Total Environ. 636, 910–918. https://doi.org/10.1016/j.scitotenv.2018.04.332 Wilson, V.S., Bobseine, K., Lambright, C.R., Gray, L.E., 2002. A novel cell line, MDA-kb2, that stably expresses an androgen- and glucocorticoid-responsive reporter for the detection of hormone receptor agonists and antagonists. Toxicol. Sci. 66, 69–81. Wode, F., van Baar, P., Dünnbier, U., Hecht, F., Taute, T., Jekel, M., Reemtsma, T., 2015. Search for over 2000 current and legacy micropollutants on a wastewater infiltration site with a UPLC-high resolution MS target screening method. Water Res. 69, 274–283. https://doi.org/10.1016/J.WATRES.2014.11.034 Zedda, M., Zwiener, C., 2012. Is nontarget screening of emerging contaminants by LC-HRMS successful? A plea for compound libraries and computer tools. Anal. Bioanal. Chem. 403, 2493–2502. https://doi.org/10.1007/s00216-012-5893-y Zhao, J.-L., Ying, G.-G., Yang, B., Liu, S., Zhou, L.-J., Chen, Z.-F., Lai, H.-J., 2011. Screening of multiple hormonal activities in surface water and sediment from the Pearl River system, South China, using effect-directed in vitro bioassays. Environ. Toxicol. Chem. 30, 2208–15. https://doi.org/10.1002/etc.625 60 7 ANNEXES List of annexes: ANNEX I Tousova, Z., Oswald, P., Slobodnik, J., Blaha, L., Muz, M., Hu, M., Brack, W., Krauss, M., Di Paolo, C., Tarcai, Z., Seiler, T.-B., Hollert, H., Koprivica, S., Ahel, M., Schollée, J.E., Hollender, J., Suter, M.J.-F., Hidasi, A.O., Schirmer, K., Sonavane, M., Ait-Aissa, S., Creusot, N., Brion, F., Froment, J., Almeida, A.C., Thomas, K., Tollefsen, K.E., Tufi, S., Ouyang, X., Leonards, P., Lamoree, M., Torrens, V.O., Kolkman, A., Schriks, M., Spirhanzlova, P., Tindall, A., Schulze, T., 2017. European demonstration program on the effect-based and chemical identification and monitoring of organic pollutants in European surface waters. Sci. Total Environ. 601–602. https://doi.org/10.1016/j.scitotenv.2017.06.032 ANNEX II Tousova, Z., Froment, J., Oswald, P., Slobodník, J., Hilscherova, K., Thomas, K.V., Tollefsen, K.E., Reid, M., Langford, K., Blaha, L., 2018. Identification of algal growth inhibitors in treated waste water using effect-directed analysis based on non-target screening techniques. J. Hazard. Mater. 358. https://doi.org/10.1016/j.jhazmat.2018.05.031 ANNEX III Toušová, Z., Vrana, B., Smutná, M., Novák, J., Klučárová, V., Grabic, R., Slobodník, J., Giesy, J.P., Hilscherová, K., 2019. Analytical and bioanalytical assessments of organic micropollutants in the Bosna River using a combination of passive sampling, bioassays and multi-residue analysis. Sci. Total Environ. 650. https://doi.org/10.1016/j.scitotenv.2018.08.336 ANNEX IV Schulze, T., Ahel, M., Ahlheim, J., Aït-Aïssa, S., Brion, F., Di Paolo, C., Froment, J., Hidasi, A.O., Hollender, J., Hollert, H., Hu, M., Kloß, A., Koprivica, S., Krauss, M., Muz, M., Oswald, P., Petre, M., Schollée, J.E., Seiler, T.-B., Shao, Y., Slobodnik, J., Sonavane, M., Suter, M.J.-F., Tollefsen, K.E., Tousova, Z., Walz, K.-H., Brack, W., 2017. Assessment of a novel device for onsite integrative large-volume solid phase extraction of water samples to enable a comprehensive chemical and effect-based analysis. Sci. Total Environ. 581–582. https://doi.org/10.1016/j.scitotenv.2016.12.140 LIST OF ANNEXES 61 ANNEX V Ouyang, X., Leonards, P.E.G., Tousova, Z., Slobodnik, J., De Boer, J., Lamoree, M.H., 2016. Rapid Screening of Acetylcholinesterase Inhibitors by Effect- Directed Analysis Using LC × LC Fractionation, a High Throughput in Vitro Assay, and Parallel Identification by Time of Flight Mass Spectrometry. Anal. Chem. 88. https://doi.org/10.1021/acs.analchem.5b04311 ANNEX I Science of the Total Environment 601–602 (2017) 1849–1868 Contents lists available at ScienceDirect Science of the Total Environment journal homepage: www.elsevier.com/locate/scitotenv European demonstration program on the effect-based and chemical identification and monitoring of organic pollutants in European surface waters Zuzana Tousova a,b,PeterOswalda, Jaroslav Slobodnik a, Ludek Blaha b, Melis Muz c,d,MengHuc,d, Werner Brack c,d, Martin Krauss c, Carolina Di Paolo d, Zsolt Tarcai d, Thomas-Benjamin Seiler d, Henner Hollert d, Sanja Koprivica e,MarijanAhele,JenniferE.Scholléef,g, Juliane Hollender f,g, Marc J.-F. Suter f,g, Anita O. Hidasi f,h, Kristin Schirmer f,g,h, Manoj Sonavane i, Selim Ait-Aissa i, Nicolas Creusot i, Francois Brion i, Jean Froment c,j, Ana Catarina Almeida j, Kevin Thomas j,k, Knut Erik Tollefsen j,l,SaraTufi m,XiyuOuyangm, Pim Leonards m,MarjaLamoreem, Victoria Osorio Torrens n, Annemieke Kolkman n, Merijn Schriks n,o, Petra Spirhanzlova p,AndrewTindallp, Tobias Schulze c,⁎ a Environmental Institute (EI), Okruzna 784/42, 972 41 Kos, Slovak Republic b Masaryk University, Faculty of Science, RECETOX, Kamenice 753/5, 625 00 Brno, Czech Republic c UFZ Helmholtz Centre for Environmental Research GmbH, Permoserstrasse 15, 04318 Leipzig, Germany d RWTH Aachen University, Institute for Environmental Research (Biology V), Department of Ecosystem Analysis, Worringerweg 1, 52074 Aachen, Germany e Rudjer Boskovic Institute, Bijenicka cesta 54, 10000 Zagreb, Croatia f Eawag, Überlandstrasse 133, 8600 Dübendorf, Switzerland g Institute of Biogeochemistry and Pollutant Dynamics, ETH Zürich, 8092 Zürich, Switzerland h EPF Lausanne, School of Architecture, Civil and Environmental Engineering, 1015 Lausanne, Switzerland i Institut National de l'Environnement Industriel et des Risques (INERIS), Unité ECOT, Parc ALATA - BP2, 60550 Verneuil-en-Halatte, France j Norwegian Institute for Water Research (NIVA), Ecotoxicology and Risk Assessment, Gaustadallèen 21, NO-0349 Oslo, Norway k Queensland Alliance for Environmental Health Sciences (QAEHS), The University of Queensland, 39 Keesels Road, Coopers Plains 4108, Australia l Norwegian University of Life Sciences (NMBU), Faculty of Environmental Science & Technology, Dept. for Environmental Sciences, Post Box 5003, N-1432 Ås, Norway m Vrije Universiteit Amsterdam, Department Environment & Health, De Boelelaan 1087, 1081 HV Amsterdam, The Netherlands n KWR, Watercycle Research Institute, Department of Chemical Water, Quality and Health, P.O. Box 1072, 3430 BB Nieuwegein, The Netherlands o Vitens drinking water company, P.O Box 1205, 8001 BE Zwolle, The Netherlands p WatchFrog S. A., 1 rue Pierre Fontaine, 91000 Evry, France HIGHLIGHTS GRAPHICAL ABSTRACT • Simplified protocol for effect-based monitoring of micropollutants was ap- plied. • Large volume solid phase extraction de- vice was applied to European rivers. • Target compounds did not explain ma- jor portions of observed biological ef- fects • 21 micropollutants were prioritized based on risk assessment. http://dx.doi.org/10.1016/j.scitotenv.2017.06.032 0048-9697/© 2017 Elsevier B.V. All rights reserved. 1850 Z. Tousova et al. / Science of the Total Environment 601–602 (2017) 1849–1868 article info abstract Article history: Growing concern about the adverse environmental and human health effects of a wide range of micropollutants Received 12 April 2017 requires the development of novel tools and approaches to enable holistic monitoring of their occurrence, fate Received in revised form 4 June 2017 and effects in the aquatic environment. A European-wide demonstration program (EDP) for effect-based moni- Accepted 4 June 2017 toring of micropollutants in surface waters was carried out within the Marie Curie Initial Training Network EDA- Available online xxxx EMERGE. The main objectives of the EDP were to apply a simplified protocol for effect-directed analysis, to Editor: D. Barcelo link biological effects to target compounds and to estimate their risk to aquatic biota. Onsite large volume solid phase extraction of 50 L of surface water was performed at 18 sampling sites in four European river Keywords: basins. Extracts were subjected to effect-based analysis (toxicity to algae, fish embryo toxicity, neurotoxicity, Adverse effects (anti-)estrogenicity, (anti-)androgenicity, glucocorticoid activity and thyroid activity), to target analysis (151 or- Large volume solid phase extraction ganic micropollutants) and to nontarget screening. The most pronounced effects were estrogenicity, toxicity to EDA-EMERGE algae and fish embryo toxicity. In most bioassays, major portions of the observed effects could not be explained fi Simpli ed effect-directed analysis protocol by target compounds, especially in case of androgenicity, glucocorticoid activity and fish embryo toxicity. Estrone Environmental health and nonylphenoxyacetic acid were identified as the strongest contributors to estrogenicity, while herbicides, Human health with a minor contribution from other micropollutants, were linked to the observed toxicity to algae. Fipronil and nonylphenol were partially responsible for the fish embryo toxicity. Within the EDP, 21 target compounds were prioritized on the basis of their frequency and extent of exceedance of predicted no effect concentrations. The EDP priority list included 6 compounds, which are already addressed by European legislation, and 15 micropollutants that may be important for future monitoring of surface waters. The study presents a novel sim- plified protocol for effect-based monitoring and draws a comprehensive picture of the surface water status across Europe. © 2017 Elsevier B.V. All rights reserved. 1. Introduction sampling techniques like passive sampling or solid phase extraction have been developed (Dimpe and Nomngongo, 2016; Jones et al., Environmental quality monitoring of surface waters is fundamental 2015). Onsite large volume solid phase extraction (LVSPE), introduced to the sustainable management of water resources and to reducing risks by Schulze et al. (2017), enables collection and pre-concentration of posed by multiple anthropogenic stressors (Geissen et al., 2015). Cur- large samples for effect assessment and chemical analyses and allows rent EU legislation requires that member states regularly monitor a direct quantification of water concentrations of the identified com- wide range of biological and chemical parameters with the objective pounds. Chemical analyses of CECs are based on highly sensitive analyt- of achieving good chemical and ecological status of all European ical methods, which have become available with the advance of gas waterbodies (Directive 2000/60/EC). Forty-five priority substances chromatography (GC) and (ultra)-high performance liquid chromatog- have been established by the European Commission in the Water raphy [(U)HPLC] coupled either to tandem mass spectrometry (MS/ Framework Directive (WFD) for chemical monitoring and results of tar- MS), or a wide range of high resolution mass spectrometers (HR-MS), get analyses are compared against the environmental quality standards the latter allowing also for identification of unknowns in nontarget (EQSs), while biological monitoring evaluates flora, benthic inverte- screening workflows (Schymanski et al., 2015). brates and fish communities (Directive 2000/60/EC; Directive The biological effects of CECs have been assessed by means of in vitro 2008/1005/EC; Directive 2013/39/EC). However, adverse effects on and in vivo methods at various levels of biological organization, among aquatic biota cannot be easily associated with the analysis of priority which the use of cell-based methods, bacterial algal and fish models substances and the assessment of chemical and ecological status often such as zebrafish (Danio rerio), medaka (Oryzias latipes) or fathead min- results in contradictory outcomes (Altenburger et al., 2015). Therefore, now (Pimephales promelas) have been well established. For assessing ef- novel holistic monitoring approaches combining sensitive effect-based fects of CECs in surface waters and in addition to the traditional tools with sophisticated chemical analysis are needed to address a ecotoxicological endpoints like survival, growth and reproduction, spe- wide range of contaminants of emerging concerns (CEC) occurring in cific assays on neurotoxicity, immunotoxicity, oxidative stress or complex environmental mixtures. The growing list of CECs currently in- genotoxicity and hormonal activity such as (anti-)estrogenicity, (anti- cludes personal care products, human and veterinary pharmaceuticals, )androgenicity, glucocorticoid activity and thyroid activity are available surfactants and surfactant-derived compounds, X-ray contrast media, (Connon et al., 2012). These specific biological effects were repeatedly pesticides, disinfection by-products, algal toxins, flame retardants, plas- observed for surface waters, drinking waters, raw and treated wastewa- ticizers, UV-filters, industrial compounds and transformation products ters and technical systems (Dizer et al., 2002; Loos et al., 2013; Jálová et (Sima et al., 2014). Due to the broad range of physicochemical proper- al., 2013; Osman et al., 2015; Zhao et al., 2011) and some can be associ- ties, toxic modes of action and usage patterns, CECs constitute a major ated with particular compounds as the main toxicity drivers, e.g., challenge for analytical chemists, ecotoxicologists and environmental estrogenicity due to E2 and EE2 (Aerni et al., 2004; Koenig et al., 2017; regulators (Ginebreda et al., 2014). Miège et al., 2009) or phytotoxicity to herbicides (Tang and Escher, CECs in surface waters tend to occur at rather low (down to sub- 2014). However, the cause of many effects observed in bioassays re- ng L−1)andalsofluctuating concentrations, and hence integrative mains unexplained by the chemicals detected in the same samples Abbreviations: AR, androgen receptor; Dex-EQ, dexamethasone equivalent; CEC, contaminants of emerging concern; DHT, dihydrotestosterone; DHT-EQ, dihydrotestosterone equivalent; E1, estrone; E2, 17β-estradiol; E2-EQ, 17-β-estradiol equivalent; E3, estriol; EC50, concentration at which the effect reaches 50% of the effect in untreated control; EDA, effect directed analysis; EDP, European demonstration program; EE2, 17α-ethinylestradiol; EEQ-SSE, concentration of E2-EQ which is safe regarding major steroid estrogens E1, E2, E3 or EE2; EQS, environmental quality standard; ER, estrogen receptor; Flu-EQ, flutamide equivalent; GC, gas chromatography; GR, glucocorticoid receptor; HPLC, high performance liquid chromatography; (HR)MS,(highresolution)massspectrometry;LOEC, lowest observed effect concentration; LOD, limit of detection; LOQ, limit of quantification; LVSPE, large volume solid phase extraction; MEC95, 95th percentile of maximum environmental concentration; MeOH, methanol; OH-Tam-EQ, hydroxytamoxifen equivalent; PAH, polycyclic aromatic hydrocarbons,; (P-)PNEC, (provisional) predicted no effect concentration; RB, river basin; REF, relative enrichment factor; SPE, solid phase extraction; T, testosterone; TOF, time of flight; TU, toxic unit. ⁎ Corresponding author at: UFZ Helmholtz Centre for Environmental Research GmbH, Department of Effect-Directed Analysis, Permoserstrasse 15, 04318 Leipzig, Germany. E-mail address: [email protected] (T. Schulze). Z. Tousova et al. / Science of the Total Environment 601–602 (2017) 1849–1868 1851 (Burgess et al., 2013; Neale et al., 2015). Effect-directed analysis (EDA), anthropogenic pressures is given in Table 1. The main parameters of se- a method of combining biotesting, fractionation and chemical analysis lected RBs differed; however, in general all were affected by wastewater of bioactive fractions, is a powerful tool to facilitate identification of un- discharges and to a lesser extent by industrial and/or agricultural pollu- known toxicants (Brack, 2003; Brack et al., 2016). EDA studies are very tion. The impact of WWTP effluents were assessed for the Saale and costly, tedious and time consuming. However, novel sampling, Sava basin as well as for Swiss sites in SD (Table S12). The overall bioanalytical, chemical and software tools make EDA studies increasing- scheme of the EDP is given in Fig. 2. ly feasible through miniaturized and automated high-throughput for- mats, lowered detection limits, and optimized multi-target and non target screening (Guijarro et al., 2015; Brack et al., 2013). 2.2. Sampling In the presented study, a newly developed simplified EDA protocol was applied within a European demonstration program (EDP), which Single sampling was carried out in the different RBs between July included effect-based monitoring of CECs at 18 sampling sites in 4 Euro- 2013 and August 2014 using both grab sampling and LVSPE sampling. pean river basins, as a part of the EU-funded EDA-EMERGE project. The The 2 L-grab samples were kept at 4 °C and shipped immediately to simplified EDA protocol was based on existing EDA methods and includ- WatchFrog laboratory (Evry, France) for an in vivo thyroid activity ed a broad extraction procedure, simple clean-up and fractionation assay. The LVSPE sampling was performed according to a method de- techniques, followed by a set of sensitive, rapid and low volume bioas- scribed (Schulze et al., 2017). In brief, 50 L of river water were extracted says and extensive chemical screening analyses designed to add signif- with an onsite LVSPE device (UFZ, Leipzig, Germany; Maxx Mess- und icant value to classical chemical target monitoring. The specific goals of Probenahmetechnik GmbH, Rangendingen, Germany) containing this study were to i) combine the LVSPE sampling with a set of bioas- three different sorbents in sequence, designed to capture neutral says, chemical target analysis and nontarget screening in a novel simpli- (Chromabond® HR-X, Macherey-Nagel, Düren, Germany), weakly acid- fied EDA protocol; ii) identify target compounds responsible for ic (Chromabond® HR-XAW) and weakly basic (Chromabond® HR- biological effects and quantify their contribution to the observed effects; XCW) organic compounds. Each sorbent was eluted separately and iii) prioritize target compounds according to their estimated ecological the extracts were then divided into 13 aliquots according to the volume risk. This study is a result of the joint research activities of the EDA- requirements of each bioassay and chemical analysis. The aliquots were EMERGE fellows as a training event. dried and shipped to the respective laboratories for bioassays or chem- ical analyses (supplementary data (SD) - Table S3). Two blank samples 2. Materials and methods were collected with the LVSPE sampler. Extraction cartridges filled with clean sorbents were conditioned and extracted accordingly to obtain 2.1. River basins and sampling sites fabrication blanks. Laboratory blanks were prepared by percolation of 2 L of mineralized LC-grade distilled water through the LVSPE device Four case studies in four river basins (RBs) in six European countries for 100 cycles, which is equivalent to sampling of 50 L. Both blanks (i.e., Germany, Czech Republic, Slovakia, Hungary, Croatia and Switzer- were otherwise treated identically to the field samples. The detailed de- land) were selected for the EDP (Fig. 1) based on results of earlier mon- scription of the LVSPE sampling method and individual samples is pro- itoring campaigns. Overview of the sampling sites and major vided in SD (Table S1 and S2). Fig. 1. Overview map of the EDP sampling locations (from top right to bottom left) in the Saale RB (Germany), Danube RB (Czech Rep., Slovakia and Hungary), Sava RB (Croatia) and Emme RB (Switzerland). 1852 Z. Tousova et al. / Science of the Total Environment 601–602 (2017) 1849–1868 Table 1 Description of the EDP sampling sites. For more details see supplementary data (Table S1). Code Site Coordinates Site characterization and major pressures Danube RB 1-1 WWTP effluent - Svratka, 49.12447N 16.62697E Site at the effluent of large scale municipal WWTP (400,000 citizens); estrogenicity, Brno Modrice, Czech Republic dioxin-like activity, cytotoxicity, pharmaceuticals, polar pesticides, bisphenol A, PFOA and steroid hormones were previously observed (Vrana et al., 2011, Jedličková et al., 2011) 1-2 Danube - Morava Tributary, 48.17864N 16.97597E Site 300 m upstream of Morava confluence to Danube; site affected by automobile Devinska Nová Ves, Slovakia industy and agriculture; polluted sediments and SPM - DDT,DDE, DDD, PAHs, petroleum hydrocarbons and DEHP were previously observed (“Danube River Basin Water Quality Database, https://www.icpdr.org/wq-db” 2012) 1-3 Danube - Bratislava, Slovakia 48.11834N 17.14435E Site in the city centre of Bratislava, site affected by urban pressures and navigation; nonylphenols, triazine pesticides, 2,4-D, pharmaceuticals, polluted sediments and SPM with PAHs, nonylphenols, DDT, dioxins, PBDEs, PCBs and DEHP were previously observed (“Danube River Basin Water Quality Database,” 2012) 1-4 Danube - Szob, Hungary 47.59890N 19.08361E Site downstream of towns - Estergom and Sturovo and tributary of Hron; river bank filtration facility for drinking water production in the vicinity; PFCs, triazine pesticides, DEHP and pharmaceuticals were observed (“Danube River Basin Water Quality Database,” 2012); high dilution by the Danube River expected, reference to the upstream sites Sava RB 2-1 Sava - Otok Samoborski, 45.843083N 15.729167E 10 km upstream of the city of Zagreb, upstream reference location to the downstream Samobor, Croatia sites affected by various pollution sources in the area of the city of Zagreb; a checkpoint for transboundary pollution from Slovenia (Smital et al., 2013) 2-2 Sava - Podsused, Zagreb, Croatia 45.793583N 15.852783E Situated in the western part of the city of Zagreb; about 5 km downstream from the discharge point of the wastewaters from the WWTP of the city of Zapresic, including significant contribution of effluents from pharmaceutical industry; marked input of macrolide antimicrobials (Terzic and Ahel, 2011) 2-3 Sava - Oborovo, 45.686450N 16.246900E Several km downstream of the main wastewater outlets of WWTPs of the cities of Velika Gorica, Croatia Zagreb (800,000 inh.) and Velika Gorica (60,000 inh.); the wastewaters of both cities are of the mixed type, including significant contribution of industrial WWs; both WWs receive full mechanical and biological treatment before the discharge 2-4 Sava - Crnac, Sisak, Croatia 45.445267N 16.419267E 2 km downstream from the city of Sisak (50,000 inhabitants), affected by urban WW, iron works and oil refinery; notable petroleum hydrocarbon pollution in water and sediments (Smital et al., 2013) Emme RB 3-1 Limpach – Messen, Switzerland 47.10458N 7.44811E Upstream of WWTP Messen effluent to the Limpach Creek, site unaffected by WW 3-2 Limpach - Messen, Switzerland 47.10989N 7.46659E Downstream of WWTP Messen effluent to the Limpach Creek, no other pressures 3-3 Urtenen - Kernenried, Switzerland 47.05544N 7.53750E Upstream of WWTP effluent Kernenried to the Urtenen Creek, site unaffected by WW 3-4 Urtenen - Kernenried, Switzerland 47.05864N 7.54021E Downstream of WWTP effluent Kernenried to the Urtenen Creek, no other pressures 3-5 Emme - Wiler bei Utzenstorf, 47.16023N 7.54703E Confluence of Limpach Creek to Emme, downstream of both Messen and Kernenried WWTPs Switzerland Saale RB 4-1 Holtemme - Steinerne Renne, Germany 51.81794 N 10.72889E Pristine water in Harz mountains; reference site to 4-2, 4-5 4-2 Holtemme - Derenburg, Germany 51.86645N 10.879617E Site downstream of WWTP in Silstedt, impact of hospitals and wellness clinics 4-3 Saale - Rudolstadt, Germany 50.71865N 11.3984E Site with industrial impact from paper and chemical industry 4-4 Saale - Klein Rosenburg, Germany 51.931944N 11.888889E Integrative site on the Saale before entering into the Elbe river in Groß Rosenburg, sampling station of the Watch List Study 2012 4-5 Holtemme - Nienhagen, Germany 51.94153N 11.15856E Integrative site of Holtemme near the confluence with Bode; downstream of WWTP Halberstadt; urban and agricultural impact Abbreviations: DDT - 1,1,1-trichlor-2,2-bis(4-chlorfenyl)ethan; DDD - dichlorodiphenyldichloroethane; DDE - dichlorodiphenyldichloroethylene; DEHP - Bis(2-ethylhexyl) phthalate; PAHs - polycyclic aromatic hydrocarbons; PBDEs - polybrominated diphenyl ethers or PBDEs- polybrominated diphenyl ethers; PCBs - polychlorinate bifenyls; PFCs - perfluorocarbons; SPM - suspended particulate matter; WW - waste water; WWTP - water water treatment plant. 2.3. Effect assessment available RegTox Microsoft Excel™ Macro (http://www.normalesup. org/~vindimian/fr_index.html). A basic description of the assays is pro- A set of seven bioassays performed in six different laboratories was vided below, while more details on the bioassays methods and test con- applied to screen for both nonspecificandspecifictoxicityofthe ditions are provided in the SD (Table S4). water samples, including both in vitro and in vivo methods (SD - Table S3). The effective concentrations of surface water samples were 2.3.1. ER-mediated activity (MELN cells) assay expressed in relative enrichment factors (REFs) as proposed by Escher The estrogen receptor mediated bioassay based on human et al. (2006). For active samples in the receptor-mediated assays, effect breast cancer cell line (MELN cells transfected with a promoter contain- equivalents of standard agonists or antagonists were calculated. Prior to ing estrogen responsive elements driving expression of luciferase testing, the extracts of all sample fractions were reconstituted in a sol- (Balaguer et al., 1999)) was used for assessment of estrogenicity and vent and mixed with the test media to reach safe solvent concentration antiestrogenicity of the sample. The receptor mediated activity of a sam- related to the bioassay used and final REFs from 1 to 100, meaning that ple was assessed in 96-well plates by measurement of luminescence the tested range covered the original river concentrations (REF = 1) as after 16 h of exposure at 37 °C. Effective concentrations and equivalents well as concentrations up to 100-times higher (REF = 100). In excep- of estrogenic potency (using the reference estrogen estradiol, E2) and tional cases, REF 500 were used for the upper limit (e.g., in the AChE in- anti-estrogenic potency (with OH-tamoxifen as the reference) were de- hibition assay, where no effects were observed up to REF 100). One of termined as previously detailed (Creusot et al., 2016). the bioassays (in vivo thyroid activity assay with Xenopus)wasper- formed with whole raw water samples. Laboratory blank, fabrication 2.3.2. AR-mediated activity (MDA-kb2 cells) assay blank, solvent control and positive control with standard reference The androgen receptor mediated bioassay based on human breast compounds were tested in parallel with the samples in each bioassay. cancer cell line (MDA-kb2 cells transfected with a promoter containing The results of bioassays were evaluated using nonlinear regression androgen responsive elements driving expression of luciferase models to derive the effective concentrations (EC50 or EC20)in (Wilson et al., 2002)) was used for assessment of androgenicity and GraphPad Prism 6 (GraphPad Software, Inc., La Jolla, USA) or freely antiandrogenicity of the sample. The receptor mediated activity of the Z. Tousova et al. / Science of the Total Environment 601–602 (2017) 1849–1868 1853 Fig. 2. Scheme of the European demonstration programme. For more details see supplementary data (Table S1–S5). sample was assessed in 96-well plates by measurement of lumines- Potencies were expressed as equivalents (EQ) of standard reference com- cence after 16 h of exposure at 37 °C. Potencies were expressed as pound dexamethasone Dex-EQ as previously described in Macikova et al. equivalents (EQ) of standard reference compounds, i.e., DHT-EQ for (2014). androgenicity and flutamide-EQ for anti-androgenicity as described previously (Creusot et al., 2014). The 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide 2.3.4. Zebrafish embryo acute toxicity assay (FET) (MTT) test, based on the color change of tetrazolium dye induced by ox- Fish embryo acute toxicity was assessed in zebrafish (Danio rerio) idoreductase enzymes of viable cells as previously described by Creusot according to the OECD guideline 236 (OECD, 2013). Freshly fertilized et al. (2015) was used to assess cytotoxicity of the samples in the cell zebrafish eggs were exposed to six sample concentrations in 24-well based bioassays (MELN and MDA-kb2 cells). plates for 96 h, at 26 °C and 14 h light : 10 h dark cycle. Four apical end- points were observed daily as indicators of lethality: embryo coagula- tion, lack of somite formation, nondetachment of the tail and absence 2.3.3. GR-CALUX® assay of heartbeat. Additionally, the occurrence of sublethal morphological ef- GR-CALUX is a patented glucocorticoid receptor mediated bioassay fects was also recorded. Effects are described either as the occurrence of based on human osteoblastic cells (U2-OS) transfected with a promoter lethality, or as the cumulative occurrence of any endpoint for lethality or containing glucocorticoid responsive elements driving expression of lucif- sublethal toxicity (i.e., any lethal or sublethal effect). In each test, ten erase. Glucocorticoid activity of the samples was assessed in 96-well embryos were exposed per exposure concentration. 3,4-Dichloraniline plates by measurement of luminescence after 24 h of exposure at 37 °C. was used as the standard reference compound. 1854 Z. Tousova et al. / Science of the Total Environment 601–602 (2017) 1849–1868 Table 2 Summary results of target analysis. For more details see supplementary data (Table S5 and S7). Compound name CAS Compound LOQ PNEC PNEC MEC95 Median Frequency of PNEC Frequency of group/usage pattern [ng/L] [ng/L] type [ng/L] [ng/L] occurrence [%] exceedance PNEC exceedance N=18 [%] N = 18 17-alpha-Ethinylestradiol 57-63-6 Pharmaceutical, 0.40 0.007 PNEC ND ND 0 no NA estrogen, WFD watch list 17-beta-Estradiol 50-28-2 Hormone, estrogen, 0.30 0.08 PNEC ND ND 0 no NA WFD watch list 1H-Benzotriazole 95-14-7 Corrosion inhibitor, 1.0 4840 P-PNEC 1026 135 94 no NA industrial compound 2,4-Dichlorophenoxyacetic acid 94-75-7 Herbicide 1.0 200 PNEC 4.2 bLOQ 44 no NA 2,4-Dinitrophenol 51-28-5 Herbicide, 1.0 1 PNEC 7.7 1.9 72 7.70 56 transformation product 2,6-Dichlorobenzamide 2008-58-4 Herbicide, 1.0 NA NA 20.1 bLOQ 50 NA NA transformation product 4-Androstene-3,17-dione 63-05-8 Hormone, androgen 1.50 NA NA 2.99 bLOQ 11 NA NA 4-Hydroxytamoxifen 68392-35-8 Pharmaceutical, 0.60 NA NA ND ND 0 NA NA anti-estrogen 4-Nitrophenol 100-02-7 Industrial compound 5.0 NA NA bLOQ bLOQ 22 NA NA 4-Toluenesulfonamide 70-55-3 Industrial compound 10 NA NA 308 56 94 NA NA 5-Methyl-1H-benzotriazole 136-85-6 Corrosion inhibitor 1.0 5285 P-PNEC 955 53.7 94 no NA 6-alpha-Methylprednisolone 83-43-2 Pharmaceutical, 0.63 NA NA 1.44 bLOQ 17 NA NA glucocorticoid Acesulfame 55589-62-3 Artificial sweetener, 1.0 1731 P-PNEC 19.8 2.3 72 no NA Marker compound Acetyl-Sulfamethoxazole 21312-10-7 Pharmaceutical, 1.0 178177 PNEC 5.9 1.7 61 no NA transformation product Aclonifen 74070-46-5 Herbicide, WFD 0.10 12.0 PNEC ND ND 0 no NA priority Alachlor 15972-60-8 Herbicide, WFD 1.0 300 PNEC ND ND 0 no NA priority Aldosterone 52-39-1 Hormone, 14 NA NA ND ND 0 NA NA glucocorticoid Amcinonide 51022-69-6 Pharmaceutical, 0.88 NA NA ND ND 0 NA NA glucocorticoid Anastrozole 120511-73-1 Pharmaceutical, 0.15 NA NA 0.24 bLOQ 6 NA NA anti-estrogen Androsterone 53-41-8 Hormone, androgen 2.0 NA NA ND ND 0 NA NA Anthracene 120-12-7 PAH, WFD priority 1.0 100 PNEC 5.9 bLOQ 67 no NA Atrazine 1912-24-9 Herbicide, WFD 1.0 600 PNEC 18.0 5.0 94 no NA priority Atrazine-desethyl 6190-85-4 Herbicide, 1.0 300 PNEC 23.3 9.0 72 no NA transformation product Azithromycin 83905-01-5 Pharmaceutical, 1.0 90 PNEC 1022 bLOQ 17 11.36 11 antibiotic Azoxystrobin 131860-33-8 Fungicide 1.0 56 PNEC 5.0 b LOQ 67 no NA Bentazone 25057-89-0 Herbicide 1.0 60 PNEC 52.2 4.5 89 no 6 Benzo(a)pyrene 50-32-8 PAH, WFD priority 1.0 0.17 PNEC ND ND 0 no NA Benzo(b)fluoranthene 205-99-2 PAH, WFD priority 1.0 0.17 PNEC ND ND 0 no NA Benzo(g,h,i)perylene 191-24-2 PAH, WFD priority 1.0 0.17 PNEC ND ND 0 no NA Benzo(k)fluoranthene 207-08-9 PAH, WFD priority 1.0 0.17 PNEC ND ND 0 no NA Benzophenone-3 131-57-7 Sunscreen agent 1.0 345 P-PNEC 26.8 bLOQ 56 no NA Benzophenone-4 4065-45-6 UV filter 1.0 NA NA 205 40.5 89 NA NA Betamethasone 378-44-9 Pharmaceutical, 0.52 658425 P-PNEC bLOQ bLOQ 6 no NA glucocorticoid Bezafibrate 41859-67-0 pharmaceutical, 1.0 460 PNEC 5.1 bLOQ 28 no NA lipid regulator Bifenox 42576-02-3 Herbicide, WFD 20 1.2 PNEC ND ND 0 no NA priority Bisphenol A 80-05-7 Plasticizer, 4.0 100 PNEC 43.3 5.3 89 no NA xenoestrogen Caffeine 58-08-2 Stimulant in 1.0 100 PNEC 237 12.8 72 2.37 17 beverages, marker compound Canrenone 976-71-6 Pharmaceutical, 0.70 NA NA 2.64 bLOQ 11 NA NA diuretic, anti-androgen Carbamazepine 298-46-4 Pharmaceutical, 1.0 500 PNEC 393 25.2 94 no 6 antiepileptic, marker compound Carbendazim 10605-21-7 Fungicide 1.0 150 PNEC 2.8 bLOQ 44 no NA Ciprofloxacin 85721-33-1 Pharmaceutical, 10 89 PNEC bLOQ bLOQ 6 no NA Z. Tousova et al. / Science of the Total Environment 601–602 (2017) 1849–1868 1855 Table 2 (continued) Compound name CAS Compound LOQ PNEC PNEC MEC95 Median Frequency of PNEC Frequency of group/usage pattern [ng/L] [ng/L] type [ng/L] [ng/L] occurrence [%] exceedance PNEC exceedance N=18 [%] N = 18 antibiotic Chlorfenvinphos 470-90-6 Insecticide, WFD 1.0 100 PNEC ND ND 0 no NA priority Chlorotoluron 15545-48-9 Herbicide 1.0 100 PNEC 2.4 bLOQ 56 no NA Chlorpyrifos 2921-88-2 Insecticide, WFD 1.0 30 PNEC 1.7 bLOQ 11 no NA priority Clarithromycin 81103-11-9 Pharmaceutical, 1.0 60 PNEC 57.6 3.0 94 no 6 antibiotic Clobetasol propionate 25122-46-7 Pharmaceutical, 0.93 NA NA ND ND 0 NA NA glucocorticoid Clobetasone butyrate 25122-57-0 Pharmaceutical, 0.39 NA NA ND ND 0 NA NA glucocorticoid Clopidogrel 113665-84-2 Pharmaceutical, 0.20 NA NA ND ND 0 NA NA anti-coagulant Clozapine 5786-21-0 Pharmaceutical, 0.25 NA NA 11.98 0.35 61 NA NA anti-psychotic Cortisone 53-06-5 Hormone, 1.0 NA NA bLOQ bLOQ 28 NA NA glucocorticoid Cyproterone 2098-66-0 Pharmaceutical, 0.70 NA NA ND ND 0 NA NA Anti-androgen Desoximethasone 382-67-2 Pharmaceutical, 0.42 NA NA ND ND 0 NA NA glucocorticoid Dexamethasone 50-02-2 Pharmaceutical, 0.54 663141 P-PNEC bLOQ bLOQ 6 no NA glucocorticoid Diazinon 333-41-5 Insecticide 1.0 1.0 PNEC 10.8 bLOQ 50 10.75 50 Diclofenac 15307-86-5 Pharmaceutical, 1.0 10 PNEC 443 2.9 83 44.34 39 analgesic, WFD watch list Dicofol 115-32-2 Insecticide, WFD 10 0.03 PNEC ND ND 0 no NA priority Didecyldimethylammonium 2390-68-3 Surfactant 1.0 NA NA bLOQ bLOQ 6 NA NA Diethyltoluamid 134-62-3 Repellent 1.0 360 P-PNEC 191 16.7 94 no 6 Difluprednate 23674-86-4 Pharmaceutical, 0.73 NA NA bLOQ bLOQ 6 NA NA glucocorticoid Diglyme 111-96-6 Industrial solvent 1.0 NA NA 5.3 bLOQ 22 NA NA Dihydrotestosterone 521-18-6 Hormone, androgen 1.5 NA NA ND ND 0 NA NA Dichlorvos 62-73-7 Insecticide, WFD 1.0 0.6 PNEC 3.8 bLOQ 11 6.40 11 priority Dimethenamid-p 87674-68-8 Herbicide 1.0 200 PNEC 1.9 bLOQ 22 no NA Dimethoate 60-51-5 insecticide 1.0 70 PNEC bLOQ bLOQ 6 no NA Diuron 330-54-1 Herbicide, WFD 1.0 200 PNEC 29.6 1.8 72 no NA priority Drospirenone 67392-87-4 Pharmaceutical, 1.0 NA NA ND ND 0 NA NA progesterone Endosulfan aplha and beta 115-29-7 Insecticide, WFD 1.0 5 PNEC bLOQ bLOQ 6 no NA priority Epi-Androsterone 481-29-8 Hormone, androgen 2.0 NA NA 4.6 bLOQ 33 NA NA Erythromycin 114-07-8 Pharmaceutical, 10 40 PNEC 54 bLOQ 39 1.34 11 antibiotic Estriol 50-27-1 Hormone, estrogen 1.0 232530 P-PNEC ND ND 0 no NA Estrone 53-16-7 Hormone, estrogen 0.10 3.6 PNEC 1.04 0.27 78 no NA Fipronil 120068-37-3 Insecticide 0.10 12.0 PNEC 1.51 0.22 72 no NA Flunisolide 3385-03-3 Pharmaceutical, 0.65 NA NA ND ND 0 NA NA glucocorticoid Fluoranthene 206-44-0 PAH, WFD priority 1.0 6.3 PNEC 13.8 1.6 61 2.19 17 Fluorometholone 426-13-1 Pharmaceutical, 0.41 NA NA ND ND 0 NA NA glucocorticoid Gestoden 60282-87-3 Pharmaceutical, 1.0 NA NA ND ND 0 NA NA progesterone Heptachlor 76-44-8 Insecticide, WFD 1.0 0.00001 PNEC ND ND 0 no NA priority Heptachlor epoxide 1024-57-3 Insecticide, 1.0 0.00001 PNEC ND ND 0 no NA transformation product Hexachlorobenzene 118-74-1 Fungicide, WFD 1.0 10 PNEC bLOQ bLOQ 11 no NA priority Hexachlorocyclohexane 608-73-1 Insecticide, WFD 1.0 NA NA bLOQ bLOQ 6 NA NA priority Hexamethoxymethylmelamine 3089-11-0 Industrial 5.0 1350 PNEC 152 27.8 78 no NA compound, marker compound Hydrocortisone 50-23-7 Hormone, 0.86 NA NA 2.47 bLOQ 11 NA NA glucocorticoid (continued on next page) 1856 Z. Tousova et al. / Science of the Total Environment 601–602 (2017) 1849–1868 Table 2 (continued) Compound name CAS Compound LOQ PNEC PNEC MEC95 Median Frequency of PNEC Frequency of group/usage pattern [ng/L] [ng/L] type [ng/L] [ng/L] occurrence [%] exceedance PNEC exceedance N=18 [%] N = 18 Ibuprofen 15687-27-1 Pharmaceutical, 10 10 PNEC 20 bLOQ 28 2.01 22 analgesic Indeno(1,2,3-cd)pyrene 193-39-5 PAH, WFD priority 1.0 0.17 PNEC ND ND 0 no NA Irgarol (cybutryn) 28159-98-0 Biocide, anti-fouling 1.0 2.5 PNEC bLOQ bLOQ 11 no NA agent Isoproturon 34123-59-6 Herbicide, WFD 1.0 300 PNEC 21.8 bLOQ 89 no NA priority Ketoprofen 22071-15-4 Pharmaceutical, 1.0 3452 P-PNEC 5.6 bLOQ 28 no NA analgesic Lauryl diethanolamide 120-40-1 Surfactant 9.4 NA NA bLOQ bLOQ 6 NA NA Levo-norgestrel 797-63-7 Pharmaceutical, 1.0 NA NA ND ND 0 NA NA progesterone Mecoprop 93-65-2 Herbicide 1.0 18000 PNEC 13.9 1.5 94 no NA Medroxyprogesterone 520-85-4 Pharmaceutical, 0.40 NA NA bLOQ bLOQ 6 NA NA progesterone Metoprolol 37350-58-6 Pharmaceutical, beta 1.0 64000 PNEC 507 4.4 83 no NA blocker N,N-Dimethyldodecylamine-N-oxide 1643-20-5 Surfactant 1.0 NA NA 10.3 bLOQ 11 NA NA Naphthalene 91-20-3 PAH, WFD priority 1.0 2000 PNEC 2.8 bLOQ 50 no NA Naproxen 22204-53-1 Pharmaceutical, 1.0 1700 PNEC 38.6 1.4 50 no NA analgesic Nonylphenol 25154-52-3 Surfactant, WFD 10 21 PNEC 148 17 78 7.17 39 priority, xenoestrogen Nonylphenoxyacetic acid 3115-49-9 Surfactant, 1.0 176 P-PNEC 276 12.0 94 1.57 6 xenoestrogen Nonylphenoxydiethoxyacetic acid 106807–78–7 Surfactant, 1.0 NA NA 69.6 2.3 72 NA NA xenoestrogen Norethindrone 68-22-4 Pharmaceutical, 1.0 NA NA ND ND 0 NA NA progesterone Norfloxacin 70458-96-7 Pharmaceutical, 10 16600 PNEC ND ND 0 no NA antibiotic Norgestimate 35189-28-7 Pharmaceutical, 0.30 NA NA ND ND 0 NA NA progesterone Octylphenol 140-66-9 Surfactant, 3.0 100 PNEC 16.1 bLOQ 22 no NA xenoestrogen Octylphenoxyacetic acid 15234-85-2 Surfactant, 1.0 319 P-PNEC ND ND 0 no NA xenoestrogen Pentabromodiphenylether 32534-81-9 Flame retadant, WFD 1.0 NA NA ND ND 0 NA NA (congener numbers 28, 47, 99, priority 100, 153 and 154) Pentachlorobenzene 608-93-5 Industrial 1.0 7 PNEC 1.8 bLOQ 28 no NA compound, WFD priority Pentachlorophenol 87-86-5 Fungicide, WFD 100 400 PNEC ND ND 0 no NA priority Perfluorooctanesulfonic acid 1763-23-1 Surfactant, WFD 10 0.13 PNEC 31 bLOQ 67 239.23 67 priority Perfluorooctanoic acid 335-67-1 Surfactant 1.0 10000 PNEC 71.5 17.5 83 no NA Phenazone 60-80-0 Pharmaceutical, 1.0 2969 P-PNEC 20.8 bLOQ 67 no NA analgesic Phenylbenzimidazolesulfonic acid 27503-81-7 Sunscreen agent 1.0 NA NA 1026 58.9 94 NA NA Pirimicarb 23103-98-2 Insecticide 1.0 NA NA bLOQ bLOQ 17 NA NA Prednicarbate 73771-04-7 Pharmaceutical, 0.41 NA NA ND ND 0 NA NA glucocorticoid Prednisolone 50-24-8 Pharmaceutical, 0.88 2450000 PNEC ND ND 0 no NA glucocorticoid Prednisone 53-03-2 pharmaceutical, 0.76 NA NA bLOQ bLOQ 6 NA NA glucocorticoid Progesterone 57-83-0 Hormone, 1.0 NA NA bLOQ bLOQ 6 NA NA progesterone Propiconazole 60207-90-1 Fungicide 1.0 230 PNEC 4.3 bLOQ 72 no NA Prothioconazole-desthio 120983-64-4 Fungicide, 1.0 NA NA 1.2 bLOQ 39 NA NA transformation product Quinoxyfen 124495-18-7 Fungicide, WFD 10 15 PNEC ND ND 0 no NA priority Raloxifene 84449-90-1 Pharmaceutical, 0.40 NA NA bLOQ bLOQ 6 NA NA anti-/estrogen Rimexolone 49697-38-3 Pharmaceutical, 0.39 NA NA ND ND 0 NA NA glucocorticoid Roxithromycin 80214-83-1 Pharmaceutical, 1.0 1000 PNEC 31.9 bLOQ 22 no NA antibiotic Z. Tousova et al. / Science of the Total Environment 601–602 (2017) 1849–1868 1857 Table 2 (continued) Compound name CAS Compound LOQ PNEC PNEC MEC95 Median Frequency of PNEC Frequency of group/usage pattern [ng/L] [ng/L] type [ng/L] [ng/L] occurrence [%] exceedance PNEC exceedance N=18 [%] N = 18 Simazine 122-34-9 Herbicide, WFD 1.0 1000 PNEC 6.2 1.7 72 no NA priority Sotalol 3930-20-9 Pharmaceutical, beta 1.0 1958 P-PNEC 82.4 2.9 78 no NA blocker Spiroxamine 118134-30-8 Fungicide 1.0 NA NA bLOQ bLOQ 11 NA NA Sucralose 56038-13-2 Artificial sweetener, 10 346561 P-PNEC 1579 100 89 no NA marker compound Sulfadimidin/sulfamethazine 57-68-1 Pharmaceutical, 1.0 1735 P-PNEC 4.6 bLOQ 39 no NA sulfonamide antimicrobial Sulfamethoxazole 723-46-6 Pharmaceutical, 1.0 600 PNEC 92.2 8.0 83 no NA sulfonamide antimicrobial Sulfapyridine 144-83-2 Pharmaceutical, 1.0 10500 P-PNEC 88.5 1.8 61 no NA sulfonamide antimicrobial Sulfathiazole 72-14-0 Pharmaceutical, 1.0 NA NA 1.5 bLOQ 6 NA NA sulfonamide antimicrobial Tamoxifen 10540-29-1 Pharmaceutical, 0.70 NA NA ND ND 0 NA NA anti-estrogen Tebuconazole 107534-96-3 Fungicide 1.0 100 PNEC 7.4 bLOQ 56 no NA Terbuthylazine 5915-41-3 Herbicide 1.0 190 PNEC 17.7 4.7 83 no NA Terbuthylazine-desethyl 30125-63-4 Herbicide, 1.0 2.4 PNEC 8.1 2.8 94 3.38 61 transformation product Terbutryn 886-50-0 Herbicide, WFD 1.0 6.5 PNEC 4.9 1.5 94 no 6 priority Testosterone 58-22-0 Hormone, androgen 0.40 NA NA ND ND 0 NA NA Tetraglyme 143-24-8 Industrial solvent 0.10 NA NA 66.20 3.90 67 NA NA Thiacloprid 111988-49-9 Insecticide 1.0 NA NA bLOQ bLOQ 17 NA NA Trenbolone 10161-33-8 Pharmaceutical, 2.0 NA NA bLOQ bLOQ 6 NA NA growth promoter Triamcihexacetonide 5611-51-8 Pharmaceutical, 18 NA NA ND ND 0 NA NA glucocorticoid Triamcinolone 83474-03-7 Pharmaceutical, 1.2 NA NA bLOQ bLOQ 6 NA NA glucocorticoid Triamcinolone acetonide 8054-16-8 Pharmaceutical, 0.59 NA NA 1.73 bLOQ 33 NA NA glucocorticoid Triclosan 3380-34-5 Biocide 1.0 0.7 PNEC 2.4 bLOQ 39 3.43 28 Triethylcitrate 77-93-0 Plasticizer 1.0 240700 P-PNEC 126 bLOQ 39 no NA Trifluralin 1582-09-8 Herbicide, WFD 1.0 30 PNEC ND ND 0 no NA priority Triglyme 112-49-2 Industrial solvent 0.10 88221 P-PNEC 4.6 0.2 61 no NA Trimethoprim 738-70-5 Pharmaceutical, 1.0 60000 PNEC 9.4 1.2 83 no NA antibiotic Trimethyloctylammonium 2083-68-3 Surfactant 1.0 NA NA bLOQ bLOQ 6 NA NA Triphenylphosphate 115-86-6 Plasticizer, flame 5.0 30 PNEC 51.7 bLOQ 11 1.72 6 retadant Triphenylphosphine oxide 791-28-6 Industrial 5.0 17092 P-PNEC 294 bLOQ 72 no NA compound, marker compound Tris(2-butoxyethyl)phosphate 78-51-3 Flame retardant 5.0 6800 PNEC 140 bLOQ 50 no NA Verapamil 152-11-4 Pharmaceutical, 1.0 30 P-PNEC 14.5 bLOQ 56 no 6 antihypertensive MEC95 - maximum environmental concentration 95 or 95th percentile of measured concentration; ND - non-detected; NA - not available; LOQ - limit of quantification; PNEC - predicted no effect concentration based on experimental value or existing environmental quality standard (EQS); P-PNEC - predicted-PNEC based on in silico prediction. 2.3.5. Algal growth inhibition assay coding sequence under the control of thyroid hormone responsive re- Growth of a population of unicellular green alga (Raphidocelis gion of the TH/bZIP-eGFP promoter) were exposed to the sample for subcapitata) was assessed after 72 h of exposure to the sample by mea- 72 h to assess thyroid axis activity by fluorescence microscopy (Fini et suring absorbance at 680 nm. The assay was performed in transparent al., 2007). The assay was performed with unfiltered whole water 96-well plates, at 24 °C and under permanent illumination (1600 lx). samples. Potassium dichromate was used as the standard reference compound. The method was based on OECD guideline 201 (OECD, 2011)andmod- ified according to Rojickova and Dvorakova (1998). 2.3.7. Acetylcholine esterase (AChE) inhibition assay Changes in the activity of the AChE were determined after 1 h of exposure at room temperature using the color change of 2.3.6. In vivo thyroid activity assay dithiobisnitrobenzoate to 5-thio-2-nitrobenzoic acid, which is a direct Tadpoles of a stable line of transgenic Xenopus (Xenopus laevis), har- measure of hydrolysis catalyzed by AChE (Ellman et al., 1961; Galgani boring the TH/bZIP-eGFP genetic construct (green fluorescent protein and Bocquene, 1991). The assay was performed in 96-well plates, 1858 Z. Tousova et al. / Science of the Total Environment 601–602 (2017) 1849–1868 Table 3 Frequency of effects in bioassays (expressed as percentage) for all sampling sites and different LVSPE sorbents. For more details see supplementary data (Table S6). Bioassay Endpoint Sampling sites Neutral sorbent (HR-X) Weak anionic exchanger (HR-XAW) Weak cationic exchanger (HR-XCW) N=18 N=18 N=18 N=17 ER-mediated activity Estrogenicity 83 83 28 12 MELN cells Anti-estrogenicity 46 46 0 (N = 13) 0 (N = 12) AR-mediated activity Androgenicity 22 22 0 0 MDA-kb2 cells Anti-androgenicity 17 17 0 0 GR CALUX Glucocorticoid activity 28 28 0 (N = 13) 0 (N = 12) Zebrafish embryo acute toxicity Survival 83 83 28 13 (N = 16) Sublethal endpoints 73 73 28 13 (N = 16) Algal growth inhibition Growth rate inhibition 83 78 6 18 In vivo thyroid activity Thyroid activity 6 –– – AChE inhibition AChE inhibition 61 17 22 41 using an automated liquid handling system. Dichlorvos was used as a detected compounds in a sample was used to quantify the overall con- standard reference compound. tribution of target compounds to the observed toxicity.