Department Environmental Chemistry

Organic contaminants enter aquatic ecosystems from various sources, including wastewater treatment plant effluent. Freshwater biofilms play a major role in the removal of organic contaminants from receiving water bodies, but knowledge of the molecular mechanisms driving contaminant biotransformations in complex stream biofilm (periphyton) communities remains limited. Previously, we demonstrated that biofilms in experimental flume systems grown at higher ratios of treated wastewater (WW) to stream water displayed an increased biotransformation potential for a number of organic contaminants. We identified a positive correlation between WW percentage and biofilm biotransformation rates for the widely-used insect repellent, N,N-diethyl-meta-toluamide (DEET) and a number of other wastewater-borne contaminants with hydrolyzable moieties. Here, we conducted deep shotgun sequencing of flume biofilms and identified a positive correlation between WW percentage and metagenomic read abundances of DEET hydrolase (DH) homologs. To test the causality of this association, we constructed a targeted metagenomic library of DH homologs from flume biofilms. We screened our complete metagenomic library for activity with four different substrates, including DEET, and a subset thereof with 183 WW-related organic compounds. The majority of active hydrolases in the metagenomic library preferred aliphatic and aromatic ester substrates while, remarkably, only a single reference enzyme was capable of DEET hydrolysis. Of the 626 total enzyme-substrate combinations tested, approximately 5% were active enzyme-substrate pairs. Metagenomic DH family homologs revealed a broad substrate promiscuity spanning 22 different compounds when summed across all enzymes tested. We biochemically characterized the most promiscuous and active enzymes identified based on metagenomic analysis from uncultivated Rhodospirillaceae and Planctomycetaceae. In addition to characterizing new DH family enzymes, we exemplified a framework for linking metagenome-guided hypothesis generation with experimental validation. Overall, this study expands the scope of known enzymatic contaminant biotransformations for metagenomic hydrolases from WW-receiving stream biofilm communities.

Biodegradation plays a key role in the fate of chemicals in the environment. The variability of biodegradation in time can cause uncertainty in evaluating the environmental persistence and risk of chemicals. However, the seasonality of biodegradation in rivers has not yet been the subject of environmentally relevant testing and systematic investigation for large numbers of chemicals. In this work, we studied the biodegradation of 96 compounds during four seasons at four locations (up- and downstream of WWTPs located on two Swedish rivers). Significant seasonality (ANOVA, p < 0.05) of the first-order rate constant for primary biodegradation was observed for most compounds. Variations in pH and total bacterial cell count were not the major factors explaining the seasonality of biodegradation. Deviation from the classical Arrhenius-type behavior was observed for most of the studied compounds, which calls into question the application of this relationship to correct biodegradation rate constants for differences in environmental temperature. Similarities in magnitude and seasonality of biodegradation rate constants were observed for some groups of chemicals possessing the same functional groups. Moreover, reduced seasonality of biodegradation was observed downstream of WWTPs, while biodegradation rates of most compounds were not significantly different between up- and downstream.

The use of agrochemical and pharmaceutical active ingredients is essential in our modern society. Given the increased concern and awareness of the potential risks of some chemicals, there is a growing need to align with 'green chemistry' and 'safe and sustainable by design' principles and thus to evaluate the hazards of agrochemical and pharmaceutical active ingredients in early stages of R&D. We give an overview of the current challenges and opportunities to assess the principle of biodegradability in the environment. Development of new medium/high-throughput methodologies, combining predictive tools and wet-lab experimentation are essential to design biodegradable chemicals early in the active ingredient discovery and selection process.

The assessment of persistence (P), bioaccumulation (B), and toxicity (T) of a chemical is a crucial first step at ensuring chemical safety and is a cornerstone of the European Union's chemicals regulation REACH (Registration, Evaluation, Authorization, and Restriction of Chemicals). Existing methods for PBT assessment are overly complex and cumbersome, have produced incorrect conclusions, and rely heavily on animal-intensive testing. We explore how new-approach methodologies (NAMs) can overcome the limitations of current PBT assessment. We propose two innovative hazard indicators, termed cumulative toxicity equivalents (CTE) and persistent toxicity equivalents (PTE). Together they are intended to replace existing PBT indicators and can also accommodate the emerging concept of PMT (where M stands for mobility). The proposed "toxicity equivalents" can be measured with high throughput in vitro bioassays. CTE refers to the toxic effects measured directly in any given sample, including single chemicals, substitution products, or mixtures. PTE is the equivalent measure of cumulative toxicity equivalents measured after simulated environmental degradation of the sample. With an appropriate panel of animal-free or alternative in vitro bioassays, CTE and PTE comprise key environmental and human health hazard indicators. CTE and PTE do not require analytical identification of transformation products and mixture components but instead prompt two key questions: is the chemical or mixture toxic, and is this toxicity persistent or can it be attenuated by environmental degradation? Taken together, the proposed hazard indicators CTE and PTE have the potential to integrate P, B/M and T assessment into one high-throughput experimental workflow that sidesteps the need for analytical measurements and will support the Chemicals Strategy for Sustainability of the European Union.

The assessment of environmental hazard indicators such as persistence, mobility, toxicity, or bioaccumulation of chemicals often results in highly variable experimental outcomes. Persistence is particularly affected due to a multitude of influencing environmental factors, with biodegradation experiments resulting in half-lives spanning several orders of magnitude. Also, half-lives may lie beyond the limits of reliable half-life quantification, and the number of available data points per substance may vary considerably, requiring a statistically robust approach for the characterization of data. Here, we apply Bayesian inference to address these challenges and characterize the distributions of reported soil half-lives. Our model estimates the mean, standard deviation, and corresponding uncertainties from a set of reported half-lives experimentally obtained for a single substance. We apply our inference model to 893 pesticides and pesticide transformation products with experimental soil half-lives of varying data quantity and quality, and we infer the half-life distribution for each compound. By estimating average half-lives, their experimental variability, and the uncertainty of the estimations, we provide a reliable data source for building predictive models, which are urgently needed by regulatory authorities to manage existing chemicals and by industry to design benign, nonpersistent chemicals. Our approach can be readily adapted for other environmental hazard indicators.

Identifying a chemical's potential for biotransformation in the aquatic environment is crucial to predict its fate and manage its potential hazards. Due to the complexity of natural water bodies, especially river networks, biotransformation is often studied in laboratory experiments, assuming that study outcomes can be extrapolated to compound behavior in the field. Here, we investigated to what extent outcomes of laboratory simulation studies indeed reflect biotransformation kinetics observed in riverine systems. To determine in-field biotransformation, we measured loads of 27 wastewater treatment plant effluent-borne compounds along the Rhine and its major tributaries during two seasons. Up to 21 compounds were detected at each sampling location. Measured compound loads were used in an inverse model framework of the Rhine river basin to derive k'_bio,field values – a compound-specific parameter describing the compounds' average biotransformation potential during the field studies. To support model calibration, we performed phototransformation and sorption experiments with all the study compounds, identifying 5 compounds that are susceptible towards direct phototransformation and determining K_oc values covering four orders of magnitude. On the laboratory side, we used a similar inverse model framework to derive k'_bio,lab values from water-sediment experiments run according to a modified OECD 308-type protocol. The comparison of k'_bio,lab and k'_bio,field revealed that their absolute values differed, pointing towards faster transformation in the Rhine river basin. Yet, we could demonstrate that relative rankings of biotransformation potential and groups of compounds with low, moderate and high persistence agree reasonably well between laboratory and field outcomes. Overall, our results provide evidence that laboratory-based biotransformation studies using the modified OECD 308 protocol and k'_bio values derived thereof bear considerable potential to reflect biotransformation of micropollutants in one of the largest European river basins.

Motivation: Transformation products (TPs) of man-made chemicals, formed through microbially mediated transformation in the environment, can have serious adverse environmental effects, yet the analytical identification of TPs is challenging. Rule-based prediction tools are successful in predicting TPs, especially in environmental chemistry applications that typically have to rely on small datasets, by imparting the existing knowledge on enzyme-mediated biotransformation reactions. However, the rules extracted from biotransformation reaction databases usually face the issue of being over/under-generalized and are not flexible to be updated with new reactions.
Results: We developed an automatic rule extraction tool called enviRule. It clusters biotransformation reactions into different groups based on the similarities of reaction fingerprints, and then automatically extracts and generalizes rules for each reaction group in SMARTS format. It optimizes the genericity of automatic rules against the downstream TP prediction task. Models trained with automatic rules outperformed the models trained with manually curated rules by 30% in the area under curve (AUC) scores. Moreover, automatic rules can be easily updated with new reactions, highlighting enviRule's strengths for both automatic extraction of optimized reactions rules and automated updating thereof.
Availability and implementation: enviRule code is freely available at github.com/zhangky12/enviRule.

For many polar organic micropollutants, biotransformation by activated sludge microorganisms is a major removal process during wastewater treatment. However, our current understanding of how wastewater treatment operations influence microbial communities and their micropollutant biotransformation potential is limited, leaving major parts of observed variability in biotransformation rates across treatment facilities unexplained. Here, we present biotransformation rate constants for 42 micropollutants belonging to different chemical classes along a gradient of solids retention time (SRT). The geometric mean of biomass-normalized first-order rate constants shows a clear increase between 3 d and 15 d SRT by 160% and 87%, respectively, in two experiments. However, individual micropollutants show a variety of trends. Rate constants of oxidative biotransformation reactions mostly increased with SRT. Yet, nitrifying activity could be excluded as primary driver. For substances undergoing other than oxidative reactions, i.e. mostly substitution-type reactions, more diverse dependencies on SRT were observed. Most remarkably, characteristic trends were observed for groups of substances undergoing similar types of initial transformation reaction, suggesting that shared enzymes or enzyme systems that are conjointly regulated catalyze biotransformation reactions within such groups. These findings open up opportunities for correlating rate constants with measures of enzyme abundance, which in turn should help to identify genes or gene products associated with the respective biotransformation reactions.

Aerobic biodegradation half-lives (half-lives) are key parameters used to evaluate pesticide persistence in soil. However, half-life estimates for individual pesticides often span several orders of magnitude, reflecting the impact that various environmental or experimental parameters have on half-lives in soil. In this work, we collected literature-reported half-lives for eleven pesticides along with associated metadata describing the environmental or experimental conditions under which they were derived. We then developed a multivariable framework to discover relationships between the half-lives and associated metadata. We first compared data for the herbicide atrazine collected from 95 laboratory and 65 field studies. We discovered that atrazine application history and soil texture were the parameters that have the largest influence on the observed half-lives in both types of studies. We then extended the analysis to include ten additional pesticides with data collected exclusively from laboratory studies. We found that, when data were available, pesticide application history and biomass concentrations were always positively associated with half-lives. The relevance of other parameters varied among the pesticides, but in some cases the variability could be explained by the physicochemical properties of the pesticides. For example, we found that the relative significance of the organic carbon content of soil for determining half-lives depends on the relative solubility of the pesticide. Altogether, our analyses highlight the reciprocal influence of both environmental parameters and intrinsic physicochemical properties for determining half-lives in soil.

The presence of antibiotics in treated wastewater and consequently in surface and groundwater resources raises concerns about the formation and spread of antibiotic resistance. Improving the removal of antibiotics during wastewater treatment therefore is a prime objective of environmental engineering. Here we obtained a detailed picture of the fate of sulfonamide antibiotics during activated sludge treatment using a combination of analytical methods. We show that pterin-sulfonamide conjugates, which are formed when sulfonamides interact with their target enzyme to inhibit folic acid synthesis, represent a major biotransformation route for sulfonamides in laboratory batch experiments with activated sludge. The same major conjugates were also present in the effluents of nine Swiss wastewater treatment plants. The demonstration of this biotransformation route, which is related to bacterial growth, helps explain seemingly contradictory views on optimal conditions for sulfonamide removal. More importantly, since pterinsulfonamide conjugates show retained antibiotic activity, our findings suggest that risk from exposure to sulfonamide antibiotics may be less reduced during wastewater treatment than previously assumed. Our results thus further emphasize the inadequacy of focusing on parent compound removal and the importance of investigating biotransformation pathways and removal of bioactivity to properly assess contaminant removal in both engineered and natural systems.

To optimize removal of organic micropollutants from the water cycle, understanding the processes during activated sludge treatment is essential. In this study, we hypothesize that aliphatic amines, which are highly abundant among organic micropollutants, are partly removed from the water phase in activated sludge through ion trapping in protozoa. In ion trapping, which has been extensively investigated in medical research, the neutral species of amine-containing compounds diffuse through the cell membrane and further into acidic vesicles present in eukaryotic cells such as protozoa. There they become trapped because diffusion of the positively charged species formed in the acidic vesicles is strongly hindered. We tested our hypothesis with two experiments. First, we studied the distribution of the fluorescent amine acridine orange in activated sludge by confocal fluorescence imaging. We observed intense fluorescence in distinct compartments of the protozoa, but not in the bacterial biomass. Second, we investigated the distribution of 12 amine-containing and eight control micropollutants in both regular activated sludge and sludge where the protozoa had been inactivated. In contrast to most control compounds, the amine-containing micropollutants displayed a distinctly different behavior in the noninhibited sludge compared to the inhibited one: (i) more removal from the liquid phase; (ii) deviation from first-order kinetics for the removal from the liquid phase; and (iii) higher amounts in the solid phase. These results provide strong evidence that ion trapping in protozoa occurs and that it is an important removal mechanism for amine-containing micropollutants in batch experiments with activated sludge that has so far gone unnoticed. We expect that our findings will trigger further investigations on the importance of this process in full-scale wastewater treatment systems, including its relevance for accumulation of ammonium.

Developing models for the prediction of microbial biotransformation pathways and half-lives of trace organic contaminants in different environments requires as training data easily accessible and sufficiently large collections of respective biotransformation data that are annotated with metadata on study conditions. Here, we present the Eawag-Soil package, a public database that has been developed to contain all freely accessible regulatory data on pesticide degradation in laboratory soil simulation studies for pesticides registered in the EU (282 degradation pathways, 1535 reactions, 1619 compounds and 4716 biotransformation half-life values with corresponding metadata on study conditions). We provide a thorough description of this novel data resource, and discuss important features of the pesticide soil degradation data that are relevant for model development. Most notably, the variability of half-life values for individual compounds is large and only about one order of magnitude lower than the entire range of median half-life values spanned by all compounds, demonstrating the need to consider study conditions in the development of more accurate models for biotransformation prediction. We further show how the data can be used to find missing rules relevant for predicting soil biotransformation pathways. From this analysis, eight examples of reaction types were presented that should trigger the formulation of new biotransformation rules, e.g., Ar-OH methylation, or the extension of existing rules, e.g., hydroxylation in aliphatic rings. The data were also used to exemplarily explore the dependence of half-lives of different amide pesticides on chemical class and experimental parameters. This analysis highlighted the value of considering initial transformation reactions for the development of meaningful quantitative-structure biotransformation relationships (QSBR), which is a novel opportunity offered by the simultaneous encoding of transformation reactions and corresponding half-lives in Eawag-Soil. Overall, Eawag-Soil provides an unprecedentedly rich collection of manually extracted and curated biotransformation data, which should be useful in a great variety of applications.

The University of Minnesota Biocatalysis/Biodegradation Database and Pathway Prediction System (UM-BBD/PPS) has been a unique resource covering microbial biotransformation pathways of primarily xenobiotic chemicals for over 15 years. This paper introduces the successor system, enviPath (The Environmental Contaminant Biotransformation Pathway Resource), which is a complete redesign and reimplementation of UM-BBD/PPS. enviPath uses the database from the UM-BBD/PPS as a basis, extends the use of this database, and allows users to include their own data to support multiple use cases. Relative reasoning is supported for the refinement of predictions and to allow its extensions in terms of previously published, but not implemented machine learning models. User access is simplified by providing a REST API that simplifies the inclusion of enviPath into existing workflows. An RDF database is used to enable simple integration with other databases. enviPath is publicly available at envipath.org with free and open access to its core data.

The main removal process for polar organic micropollutants during activated sludge treatment is biotransformation, which often leads to the formation of stable transformation products (TPs). Because the analysis of TPs is challenging, the use of pathway prediction systems can help by generating a list of suspected TPs. To complete and refine pathway prediction, comprehensive biotransformation studies for compounds exhibiting pertinent functional groups under environmentally relevant conditions are needed. Because many polar organic micropollutants present in wastewater contain one or several amine functional groups, we systematically explored amine biotransformation by conducting experiments with 19 compounds that contained 25 structurally diverse primary, secondary, and tertiary amine moieties. The identification of 144 TP candidates and the structure elucidation of 101 of these resulted in a comprehensive view on initial amine biotransformation reactions. The reactions with the highest relevance were N-oxidation, N-dealkylation, N-acetylation, and N-succinylation. Whereas many of the observed reactions were similar to those known for the mammalian metabolism of amine-containing xenobiotics, some N-acylation reactions were not previously described. In general, different reactions at the amine functional group occurred in parallel. Finally, recommendations on how these findings can be implemented to improve microbial pathway prediction of amine-containing micropollutants are given.