Background  Antimicrobial resistance is a major global health concern, with the environment playing a key role in its emergence and spread. Understanding the relationships between environmental factors, microbial communities, and resistance mechanisms is vital for elucidating environmental resistome dynamics. In this study, we characterized the environmental resistome of the Baltic Sea and evaluated how environmental gradients and spatial variability, alongside its microbial communities and associated functional genes, influence resistome diversity and composition across geographic regions.

Results  We analyzed the metagenomes of benthic sediments from 59 monitoring stations across a 1,150 km distance of the Baltic Sea, revealing an environmental resistome comprised of predicted antimicrobial resistance genes (ARGs) associated with resistance against 26 antibiotic classes. We observed spatial variation in its resistance profile, with higher resistome diversity in the northern regions and a decline in the dead zones and the southern areas. The combined effects of salinity and temperature gradients, alongside nutrient availability, created a complex environmental landscape that shaped the diversity and distribution of the predicted ARGs. Salinity predominantly influenced microbial communities and predicted ARG composition, leading to clear distinctions between high-saline regions and those with lower to mid-level salinity. Furthermore, our analysis suggests that microbial community composition and mobile genetic elements might be crucial in shaping ARG diversity and composition.

Conclusions  We presented that salinity and temperature were identified as the primary environmental factors influencing resistome diversity and distribution across geographic regions, with nutrient availability further shaping these patterns in the Baltic Sea. Our study also highlighted the interplay between microbial communities, resistance, and associated functional genes in the benthic ecosystem, underscoring the potential role of microbial and mobile genetic element composition in ARG distribution. Understanding how environmental factors and microbial communities modulate environmental resistomes will help predict the impact of future environmental changes on resistance mechanisms in complex aquatic ecosystems.

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.

Biodegradation is one of the most important processes influencing the fate of organic contaminants in the environment. Quantitative understanding of the spatial variability in environmental biodegradation is still largely uncharted territory. Here, we conducted modified OECD 309 tests to determine first-order biodegradation rate constants for 97 compounds in 18 freshwater river segments in five European countries: Sweden, Germany, Switzerland, Spain, and Greece. All but two of the compounds showed significant spatial variability in rate constants across European rivers (ANOVA, P < 0.05). The median standard deviation of the biodegradation rate constant between rivers was a factor of 3. The spatial variability was similar between pristine and contaminated river segments. The longitude, total organic carbon, and clay content of sediment were the three most significant explanatory variables for the spatial variability (redundancy analysis, P < 0.05). Similarities in the spatial pattern of biodegradation rates were observed for some groups of compounds sharing a given functional group. The pronounced spatial variability presents challenges for the use of biodegradation simulation tests to assess chemical persistence. To reflect the variability in the biodegradation rate, the modified OECD 309 test would have to be repeated with water and sediment from multiple sites.

There is a design-to-function knowledge gap regarding how engineered stream restoration structures can maximize hyporheic contaminant attenuation. Surface and subsurface structures have each been studied in isolation as techniques to restore hyporheic exchange, but surface-subsurface structures have not been investigated or optimized in an integrated manner. Here, we used a numerical model to systematically evaluate key design variables for combined surface (i.e., weir height and length) and subsurface (i.e., upstream and downstream baffle plate spacing) structures. We also compared performance metrics that place differing emphasis on hyporheic flux versus transit times. We found that surface structures tended to create higher flux, shorter transit time flowpaths, whereas subsurface structures promoted moderate-flux, longer transit time flowpaths. Optimal combined surface-subsurface structures could increase fluxes and transit times simultaneously, thus providing conditions for contaminant attenuation that were many times more effective than surface or subsurface structures alone. All performance metrics were improved by the presence of an upstream plate and the absence of a downstream plate. Increasing the weir length tended to improve all metrics, whereas the optimal weir height varied based on metrics. These findings may improve stream restoration by better aligning specific restoration goals with appropriate performance metrics and hyporheic structure designs. 

The environmental relevance of standard biodegradation tests such as OECD 309 has been questioned. Challenges include the interpretation of changing degradation kinetics over the 60–90 incubation days and the effects of chemical spiking on the microbial community. To ameliorate these weaknesses, we evaluated a modified OECD 309 test using water and sediment from three Swedish rivers. For each river, we had three treatments (no spiking, 0.5 μg L^–1 spiking, and 5 μg L^–1 spiking). The dissipation of a mixture of 56–80 spiked chemicals was followed over 14 days. Changes in dissipation kinetics during the incubation were interpreted as a departure of the microbial community from its initial (natural) state. The biodegradation kinetics were first-order throughout the incubation in the no spiking and 0.5 μg L^–1 spiking treatments for almost all chemicals, but for the 5 μg L^–1 treatment, more chemicals showed changes in kinetics. The rate constants in the no spiking and 0.5 μg L^–1 treatments agreed within a factor of 2 for 35 of 37 cases. We conclude that the environmental relevance of OECD 309 is improved by considering only the initial biodegradation phase and that it is not compromised by spiking multiple chemicals at 0.5 μg L^–1. 

Artificial tracers are often used for quantitative estimates of solute transport properties in rivers. However, single-injection tracer tests give insights in transport characteristics limited to the ecohydrological conditions at the testing time. Series of time-consuming and laborious tracer tests would be required to properly capture seasonal changes. The present study uses intrinsic diurnal fluctuations of electrical conductivity (EC) caused by discharge of treated wastewater as a tracer to evaluate solute transport processes along a 4.7-km reach of the River Erpe, Germany. By reproducing the fluctuations recorded along the river using the solute transport model one-dimensional transport with inflow and storage, this study investigated the long-term dynamics in solute transport properties. Individual 48-hr curves of EC were used in the steady state configuration of the model to gain 48-hr-integrated estimates of selected transport parameters. Using a sliding window approach in 1-hr steps along the 2,270-hr time series of EC the temporal variability of solute transport between April and June 2016 was assessed. To test the identifiability of parameters using the proposed method, sensitivity analyses and a breakthrough curve analysis of selected 48-hr windows were implemented. With time advancing into the summer, a significant rising trend (Mann-Kendall test p-value < 0.05) of the cross sectional area of the channel was observed and attributed to the growth of macrophytes and a significant slightly decreasing trend for the storage rate was found. The presented method is of high value for river management, as promoting transient storage enhances biogeochemical cycling and benefits water quality.

The hyporheic zone (HZ) is considered a hydrodynamically-driven bioreactor with significant pollutant removal capacities and can therefore not only improve wholestream water quality but also preserve human and ecosystem health. Microbial metabolism is hypothesized to play a key role in pollutant transformation in hyporheic sediments of natural streams. However, previous work investigating the influence of microbial metabolism on pollutant transformation has been predominantly laboratory studies. The key challenge for field studies is the appropriate determination of net microbial metabolism, i.e. information on the actual exposure times to specific microbial processes in the investigated system. The present study uses reactive fluorescent tracers to determine microbial metabolism and ultimately its influence on pollutant transformation, e.g. for trace organic compounds, in hyporheic sediments under natural conditions. In particular, the reactive fluorescent tracers resazurin and its main transformation product resorufin were used to determine the microbial metabolism of facultative or obligate aerobes. The influence of the derived microbial metabolism on the transformation of 20 trace organic compounds, such as pharmaceuticals, including 3 parent–daughter pairs, was examined. The present findings validate laboratory results on the microbially-mediated transformation of the anticonvulsant gabapentin to its main transformation product gabapentin lactam under natural conditions. All other TrOCs investigated did not show a clear link between TrOC reactivity to the microbial metabolism informed by the resazurin–resorufin-system. Overall, the present study not only demonstrates the use of the fluorescent tracer-system resazurin and resorufin for determining microbial metabolism of facultative or obligate aerobes but also generally highlights the potential of reactive fluorescent tracers to disentangle specific reactive properties and ultimately their influence on the fate of pollutants in natural HZs.

Toxicometabolomics and biotransformation product (bioTP) elucidation were carried out in single zebrafish (ZF) embryos exposed to carbamazepine (CBZ). Exposures were conducted in 96-well plates containing six CBZ concentrations ranging from 0.5 μg/L to 50 mg/L (n = 12 embryos per dose). In the 50 mg/L dose group, 33% of embryos developed edema during the exposure (120 hpf), while hatching was significantly delayed in three of the lower-dose groups (0.46, 3.85, and 445 μg/L) compared to the control at 48 hpf. Toxicometabolomic analysis together with random forest modeling revealed a total of 80 significantly affected metabolites (22 identified via targeted lipidomics and 58 via nontarget analysis). The wide range of doses enabled the observation of both monotonic and nonmonotonic dose responses in the metabolome, which ultimately produced a unique and comprehensive biochemical picture that aligns with existing knowledge on the mode of action of CBZ. The combination of high dose exposures and apical endpoint assessment in single embryos also enabled hypothesis generation regarding the target organ for the morphologically altering insult. In addition, two CBZ bioTPs were identified without additional exposure experiments. Overall, this work showcases the potential of toxicometabolomics and bioTP determination in single ZF embryos for rapid and comprehensive chemical hazard assessment.

Metoprolol is widely used as a beta-blocker and considered an emerging contaminant of environmental concern due to pseudo persistence in wastewater effluents that poses a potential ecotoxicological threat to aquatic ecosystems. Microbial removal of metoprolol in the redox-delineated hyporheic zone (HZ) was investigated using streambed sediments supplemented with 15 or 150 mu M metoprolol in a laboratory microcosm incubation under oxic and anoxic conditions. Metoprolol disappeared from the aqueous phase under oxic and anoxic conditions within 65 and 72 days, respectively. Metoprolol was refed twice after initial depletion resulting in accelerated disappearance under both conditions. Metoprolol disappearance was marginal in sterile control microcosms with autoclaved sediment. Metoprolol was transformed mainly to metoprolol acid in oxic microcosms, while metoprolol acid and alpha-hydroxymetoprolol were formed in anoxic microcosms. Transformation products were transient and disappeared within 30 days under both conditions. Effects of metoprolol on the HZ bacterial community were evaluated using DNA- and RNA-based time-resolved amplicon Illumina MiSeq sequencing targeting the 16S rRNA gene and 16S rRNA, respectively, and were prominent on 16S rRNA rather than 16S rRNA gene level suggesting moderate metoprolol-induced activity-level changes. A positive impact of metoprolol on Sphingomonadaceae and Enterobacteriaceae under oxic and anoxic conditions, respectively, was observed. Nitrifiers were impaired by metoprolol under oxic and anoxic conditions. Collectively, our findings revealed high metoprolol biodegradation potentials in the hyporheic zone under contrasting redox conditions associated with changes in the active microbial communities, thus contributing to the attenuation of micropollutants.