Seas and oceans are a global reservoir of antibiotic resistance genes (ARGs). Only a few studies investigated the dynamics of ARGs along the water column of the Black Sea, a unique environment, with a peculiar geology, biology and history of anthropogenic pollution. In this study, we analyzed metagenomic data from two sampling campaigns (2013 and 2019) collected across three different sites in the Western Black Sea at depths ranging from 5 to 2000 m. The data were processed to annotate ARGs, metal resistance genes (MRGs) and integron integrase genes. The ARG abundance was significantly higher in the deep water layers and depth was the main driver of beta-diversity both for ARGs and MRGs. Moreover, ARG and MRG abundances strongly correlated (r = 0.95). The integron integrase gene abundances and composition were not influenced by the water depth and did not correlate with ARGs. The analysis of the obtained MAGs showed that some of them harbored intI gene together with several ARGs and MRGs, suggesting the presence of multidrug resistant bacteria and that MRGs and integrons could be involved in the selection of ARGs. These results demonstrate that the Black Sea is not only an important reservoir of ARGs, but also that they accumulate in the deep water layers where co-selection with MRGs could be assumed as a relevant mechanism of their persistence.

Understanding the processes that underlay an ecological disaster represents a major scientific challenge. Here, we investigated phytoplankton and zooplankton community changes before and during a fauna mass kill in a European protected wetland. Evidence on gradual development and collapse of harmful phytoplankton blooms, allowed us to delineate the biotic and abiotic interactions that led to this ecological disaster. Before the mass fauna kill, mixed blooms of known harmful cyanobacteria and the killer alga Prymnesium parvum altered biomass flow and minimized zooplankton resource use efficiency. These blooms collapsed under high nutrient concentrations and inhibitory ammonia levels, with low phytoplankton biomass leading to a dramatic drop in photosynthetic oxygenation and a shift to a heterotrophic ecosystem phase. Along with the phytoplankton collapse, extremely high numbers of red planktonic crustaceans-Daphnia magna, visible through satellite images, indicated low oxygen conditions as well as a decrease or absence of fish predation pressure. Our findings provide clear evidence that the mass episode of fish and birds kill resulted through severe changes in phytoplankton and zooplankton dynamics, and the alternation on key abiotic conditions. Our study highlights that plankton-related ecosystem functions mirror the accumulated heavy anthropogenic impacts on freshwaters and could reflect a failure in conservation and restoration measures.

Environmental exposure to 4-nonylphenol (4-NP) is extensive, and studies related to human risk assessment must continue. Especially, prediction of toxicodynamics (TDs) related to reproductive toxicity in males is very important in risk-level assessment and management of 4-NP. This study aimed to develop a physiologically-based-toxicokinetic-toxicodynamic (PBTK-TD) model that added a TD prostate model to the previously reported 4-n-nonylphenol (4-n-NP) physiologically-based-pharmacokinetic (PBPK) model. Modeling was performed under the assumption of similar TKs between 4-n-NP and 4-NP because TK experiments on 4-NP, a random-mixture, are practically difficult. This study was very important to quantitatively predict the TKs and TDs of 4-NP by age at exposure using an advanced PBTK-TD model that reflected physiological-changes according to age. TD-modeling was performed based on the reported toxic effects of 4-NP on RWPE-1 cells, a human-prostate-epithelial-cell-line. Through a meta-analysis of reported human physiological data, body weight, tissue volume, and blood flow rate patterns according to age were mathematically modeled. These relationships were reflected in the PBTK-TD model for 4-NP so that the 4-NP TK and TD changes according to age and their differences could be confirmed. Differences in TK and TD parameters of 4-NP at various ages were not large, within 3.61-fold. Point-of-departure (POD) and reference-doses for each age estimated using the model varied as 426.37–795.24 and 42.64–79.52 μg/kg/day, but the differences (in POD or reference doses between ages) were not large, at less than 1.87-times. The PBTK-TD model simulation predicted that even a 100-fold 4-NP PODₘₐₙ dose would not have large toxicity to the prostate. With a focus on TDs, the predicted maximum possible exposure of 4-NP was as high as 6.06–23.60 mg/kg/day. Several toxicity-related values estimated by the dose-response curve were higher than those calculated, depending upon the PK or TK, which would be useful as a new exposure limit for prostate toxicity of 4-NP.

Reservoirs located in middle and high latitudes freeze for months in winter, where the accumulation characteristics of pollutants are changed by superimposed influence of salt exclusion from ice on the surface and pollution release from sediments at the bottom. Taking total nitrogen (TN) of Biliuhe reservoir in Northeast China as an example, we developed a model to simulate TN accumulation characteristics influenced by ice and sediments during the freezing period (NACISF), and quantified contributions of TN from ice and sediments. Model parameters of ice and sediments were determined by laboratory freeze-up simulation experiment and sediment release flux simulation experiment, and water quality data were obtained from field investigations. Results showed that the annual average amount of TN input during the ice-covered period from 2015 to 2020 was 220.77 t, the output was 400.11 t, and the accumulated amount was 589.52 t. TN excluded from ice and released from sediments contributed 8.12% and 7.17% of the total TN inputs in winter, respectively. Analysis showed that the TN excluded from ice was positively correlated with ice thickness and initial TN concentration. The maximum ice thickness of Biliuhe reservoir had a 13 year cyclic feature, and the proportion of TN excluded from ice to the total TN inputs in different periods ranged from 10.68% to 17.30% (mean 13.18%). Meanwhile, TN accumulated seasonally as summer > autumn > winter > spring. The TN exclusion effect in 2050 would be weakened when considering the combined effects of climate change and human activities, with a reduction of about 40.85% compared to the current. It is concluded that the NACISF model took into account the influences of both ice and sediments, which provided a detailed understanding of the accumulation characteristics of TN during freezing period, and had important reference significance for water quality management in winter.

Under the climate change context, warming Southern Ocean waters may allow mercury (Hg) to become more bioavailable to the Antarctic marine food web (i.e., ice-stored Hg release and higher methylation rates by microorganisms), whose biomagnification processes are poorly documented. Biomagnification of Hg in the food web of the Antarctic Peninsula, one of the world's fastest-warming regions, was examined using carbon (δ¹³C) and nitrogen (δ¹⁵N) stable isotope ratios for estimating feeding habitat and trophic levels, respectively. The stable isotope signatures and total Hg (T-Hg) concentrations were measured in Antarctic krill Euphausia superba and several Antarctic predator species, including seabirds (gentoo penguins Pygoscelis papua, chinstrap penguins Pygoscelis antarcticus, brown skuas Stercorarius antarcticus, kelp gulls Larus dominicanus, southern giant petrels Macronectes giganteus) and marine mammals (southern elephant seals Mirounga leonina). Significant differences in δ¹³C values among species were noted with a great overlap between seabird species and M. leonina. As expected, significant differences in δ¹⁵N values among species were found due to interspecific variations in diet-related to their trophic position within the marine food web. The lowest Hg concentrations were registered in E. superba (0.007 ± 0.008 μg g⁻¹) and the highest values in M. giganteus (12.090 ± 14.177 μg g⁻¹). Additionally, a significant positive relationship was found between Hg concentrations and trophic levels (reflected by δ¹⁵N values), biomagnifying nearly 2 times its concentrations at each level. Our results support that trophic interaction is the major pathway for Hg biomagnification in Southern Ocean ecosystems and warn about an increase in the effects of Hg on long–lived (and high trophic level) Antarctic predators under climate change in the future.

Bacteria play crucial roles in the biogeochemical cycle of arsenic (As) and selenium (Se) as these elements are metabolized via detoxification, energy generation (anaerobic respiration) and biosynthesis (e.g. selenocysteine) strategies. To date, arsenic and selenium biomineralization in bacteria were studied separately. In this study, the anaerobic metabolism of As and Se in Shewanella sp. O23S was investigated separately and mixed, with an emphasis put on the biomineralization products of this process. Multiple analytical techniques including ICP-MS, TEM-EDS, XRD, Micro-Raman, spectrophotometry and surface charge (zeta potential) were employed. Shewanella sp. O23S is capable of reducing selenate (SeO₄²⁻) and selenite (SeO₃²⁻) to red Se(-S)⁰, and arsenate (AsO₄³⁻) to arsenite (AsO₃³⁻). The release of H₂S from cysteine led to the precipitation of AsS minerals: nanorod AsS and granular As₂S₃. When As and Se oxyanions were mixed, both As–S and Se(-S)⁰ biominerals were synthesized. All biominerals were extracellular, amorphous and presented a negative surface charge (−24 to −38 mV). Kinetic analysis indicated the following reduction yields: SeO₃²⁻ (90%), AsO₄³⁻ (60%), and SeO₄²⁻ (<10%). The mix of SeO₃²⁻ with AsO₄³⁻ led to a decrease in As removal to 30%, while Se reduction yield was unaffected (88%). Interestingly, SeO₄²⁻ incubated with AsO₄³⁻ boosted the Se removal (71%). The exclusive extracellular formation of As and Se biominerals might indicate an extracellular respiratory process characteristic of various Shewanella species and strains. This is the first study documenting a complex interplay between As and Se oxyanions: selenite decreased arsenate reduction, whereas arsenate stimulated selenate reduction. Further investigation needs to clarify whether Shewanella sp. O23S employs multi-substrate respiratory enzymes or separate, high affinity enzymes for As and Se oxyanion respiration.

Periphyton is considered important for removal of organic pollutants from water bodies, but knowledge of the impacts of antibiotics on the community structure and ecological function of waterbodies remains limited. In this study, the effects of oxytetracycline hydrochloride (OTC) on the communities of photoautotrophic epilithon and epipelon and its effect on nitrogen and phosphorus concentrations in the water column were studied in a 12-day mesocosm experiment. The dynamics of nitrogen and phosphorus concentrations in the epipelon and epilithon experiment showed similar patterns. The concentrations of total nitrogen, dissolved total nitrogen, ammonium nitrogen, total phosphorus and dissolved total phosphorus in the water column increased rapidly during the initial days of exposure, after which a downward trend occurred. In the epilithon experiment, we found that the photosynthesis (Fv/Fm) and biomass of epilithon were significantly (P < 0.05) stimulated in the low concentration group. Contrarily, growth and photosynthesis (Fv/Fm) were significantly (P < 0.05) reduced in the medium and high concentration group. We further found that the photosynthetic efficiency of photoautotrophic epilithon was negatively correlated with the concentrations of nitrogen and phosphorus in the water column (P < 0.05). Principal coordinate analysis (PCoA) showed that the communities of epilithic algae in the control group and in the low concentration group were significantly (P < 0.05) different from that of the high concentration group during the initial 4 days. After 8 days’ exposure, all groups tended to be similar, indicating that epilithon showed rapid adaptability and/or resilience. Similar results were found for the relative abundance of some epilithic algae. Our findings indicate that the biofilm system has strong tolerance and adaptability to OTC as it recovered fast after an initial suppression, thus showing the important role of periphyton in maintaining the dynamic balance of nutrients with other processes in aquatic ecosystems.

Per- and polyfluoroalkyl substances (PFAS) have become a target of rigorous scientific research due to their ubiquitous nature and adverse health effects. However, there are still gaps in knowledge about their environmental fate and health implications. More attention is needed for remote locations with source exposures. This study focuses on assessing PFAS exposure in Gustavus, a small Alaska community, located near a significant PFAS source from airport operations and fire training sites. Residential water (n = 25) and serum (n = 40) samples were collected from Gustavus residents and analyzed for 39 PFAS compounds. In addition, two water samples were collected from the previously identified PFAS source near the community. Fourteen distinct PFAS were detected in Gustavus water samples, including 6 perfluorinated carboxylic acids (PFCAs), 7 perfluorosulfonic acids (PFSAs), and 1 fluorotelomer sulfonate (FTS). ΣPFAS concentrations in residential drinking water ranged from not detected to 120 ng/L. High ΣPFAS levels were detected in two source samples collected from the Gustavus Department of Transportation (14,600 ng/L) and the Gustavus Airport (228 ng/L), confirming these two locations as a nearby major source of PFAS contamination. Seventeen PFAS were detected in serum and ΣPFAS concentrations ranged from 0.0170 to 13.1 ng/mL (median 0.0823 ng/mL). Perfluorooctanesulfonic acid (PFOS) and perfluorohexanesulfonic acid (PFHxS) were the most abundant PFAS in both water and serum samples and comprised up to 70% of ΣPFAS concentrations in these samples. Spearman's correlation analysis revealed PFAS concentrations in water and sera were significantly and positively correlated (r = 0.495; p = 0.0192). Our results confirm a presence of a significant PFAS source near Gustavus, Alaska and suggest that contaminated drinking water from private wells contributes to the overall PFAS body burden in Gustavus residents.

We examined ingestion and retention rates of microplastics (MPs) by two freshwater (Japanese medaka and zebrafish) and two marine fish species (Indian medaka and clown anemonefish) to determine their color preferences and gastrointestinal-tract retention times. In our ingestion experiments, clown anemonefish ingested the most MP particles, followed by zebrafish, and then Japanese and Indian medaka. Next, we investigated color preferences among five MP colors. Red, yellow, and green MP were ingested at higher rates than gray and blue MPs for all tested fish species. To test whether these differences truly reflect a recognition of and preference for certain colors based on color vision, we investigated the preferences of clown anemonefish for MP colors under light and dark conditions. Under dark conditions, ingestion of MP particles was reduced, and color preferences were not observed. Finally, we assessed gastrointestinal-tract retention times for all four fish species. Some individuals retained MP particles in their gastrointestinal tracts for over 24 h after ingestion. Our results show that fish rely on color vision to recognize and express preferences for certain MP colors. In addition, MP excretion times varied widely among individuals. Our results provide new insights into accidental MP ingestion by fishes.

Nitrated and oxygenated polycyclic aromatic hydrocarbons (NPAHs and OPAHs) are receiving attention because of their high toxicity compared with parent PAHs. However, the experimental data of their physicochemical properties has been limited. This study proposed the gas chromatographic retention time (GC-RT) technique as an effective alternative one to determine octanol-air partition coefficients (KOA) and sub-cooled liquid vapor pressures (PL) for 11 NPAHs, 10 OPAHs, and 19 parent PAHs. The slopes and intercepts of the linear regressions between temperature versus KOA and PL were provided and can be used to estimate KOA and PL for the 40 targeted compounds at any temperature. The internal energies of phase transfer (ΔUOA) and enthalpies of vaporization (ΔHL) for all targeted compounds were also calculated using the GC-RT technique. High-molecular-weight compounds may release or absorb higher heat energy to transform between different phases. NPAHs and OPAHs had a non-ideal solution behavior with activity in octanol (γₒcₜ) in the range of 19–53 and 18–1,078, respectively, which is larger than the unity threshold. A comparison among four groups of PAH derivatives showed that a functional group (nitro-, oxygen-, chloro-, and bromo-) in PAH derivatives increased γₒcₜ for corresponding parent PAHs by tens (mono-group) to hundreds of times (di-group). This study suggests that the GC-RT method is applicable for indirectly measuring the physicochemical properties of various groups of organic compounds.