Strong impact of micropollutants on prokaryotic communities at the horizontal but not vertical scales in a subtropical reservoir, China

Urbanization enhances consumer protist-driven ARGs dissemination in riverine ecosystems

Despite the emergence of antibiotic resistance genes (ARGs) and antibiotic resistant bacteria (ARBs), how biological inter-trophic interactions, modulated by watershed urbanization, shape the resistome remains unexplored. We collected water samples from the highly urbanized (western: 65 % built land, sewage-affected) and lesser-urbanized (northern: 25 % built land, drinking water source) downstream tributaries of the Jiulong River in southeast China over dry and wet seasons. We utilized metagenomic and amplicon (16S and 18S rDNA) sequencing to investigate the relationships among microeukaryotic algae, consumer protists, bacterial communities, and the resistome. Metagenomic results showed that ARG-MGE-carrying contigs (mobile ARGs), rather than ARG-carrying contigs (non-mobile ARGs), exhibited more pronounced discrepancies between tributaries. A higher total abundance of ARGs and a greater number of co-shared ARGs between pathogen and non-pathogen bacteria were observed in the more urbanized western tributary. Structural equation modeling revealed that consumer protist-bacteria and algae-bacteria cohesions predominantly influenced the resistome in the western and northern tributaries, respectively. Additionally, consumer protists had more significant associations (511 out of 634) with bacteria carrying mobile ARGs in western tributary, while algae had more significant associations (73 out of 105) in northern tributary. These results highlight the distinct inter-trophic driving factors of the resistome modulated by watershed urbanization.

Combined effects of micropollutants and their degradation on prokaryotic communities at the sediment-water interface

Pesticides and pharmaceuticals enter aquatic ecosystems as complex mixtures. Various processes govern their dissipation and effect on the sediment and surface waters. These micropollutants often show persistence and can adversely affect microorganisms even at low concentrations. We investigated the dissipation and effects on procaryotic communities of metformin (antidiabetic drug), metolachlor (agricultural herbicide), and terbutryn (herbicide in building materials). These contaminants were introduced individually or as a mixture (17.6 µM per micropollutant) into laboratory microcosms mimicking the sediment-water interface. Metformin and metolachlor completely dissipated within 70 days, whereas terbutryn persisted. Dissipation did not differ whether the micropollutants were introduced individually or as part of a mixture. Sequence analysis of 16S rRNA gene amplicons evidenced distinct responses of prokaryotic communities in both sediment and water. Prokaryotic community variations were mainly driven by matrix composition and incubation time. Micropollutant exposure played a secondary but influential role, with pronounced effects of recalcitrant metolachlor and terbutryn within the micropollutant mixture. Antagonistic and synergistic non-additive effects were identified for specific taxa across taxonomic levels in response to the micropollutant mixture. This study underscores the importance of considering the diversity of interactions between micropollutants, prokaryotic communities, and their respective environments when examining sediment-water interfaces affected by multiple contaminants.

Biogeographic Patterns and Community Assembly Processes of Bacterioplankton and Potential Pathogens in Subtropical Estuaries in China

Microbial communities in coastal waters are diverse and dynamic and play important roles in ecosystem functions and services. Despite the ecological impact of bacterioplankton or pathogens, little is known about whether bacterioplankton and pathogen communities exhibit similar patterns. Here, using 16S RNA gene amplicon sequencing, the geographic patterns and assembly processes of bacterioplankton and pathogen communities in 30 subtropical estuaries were studied. Results showed that the estuarine bacterioplankton communities mainly consisted of Proteobacteria (49.06%), Actinobacteria (17.62%), and Bacteroidetes (16.33%), among which 31 pathogen genera (186 amplicon sequence variants [ASVs]) were identified. Under the influence of salinity, bacterioplankton and pathogens showed similar biogeographic patterns. Redundancy and correlation analyses indicated that the bacterioplankton communities were strongly correlated with estuarine environmental factors, but potential pathogens were less influenced. Co-occurrence network analysis revealed a close relationship between bacterioplankton and potential pathogens, with two pathogens identified as connectors (i.e., ASV340 [Clostridium perfringens] and ASV1624 [Brevundimonas diminuta]), implying potential impacts of pathogens on structure, function, and stability of estuarine bacterioplankton communities. Null-model analysis revealed that deterministic processes (heterogeneous selection) dominated bacterioplankton community assembly, while stochastic processes (undominated effect) shaped the potential pathogen community. Our findings illustrate the biogeographic patterns and community assembly mechanisms of bacterioplankton and pathogens in estuaries, which should provide guidance and a reference for the control of potential pathogenic bacteria. IMPORTANCE Bacterioplankton play an important role in estuarine ecosystem functions and services; however, potentially pathogenic bacteria may exhibit infectivity and pose a serious threat to environmental and human health. In this study, geographic patterns and assembly processes of bacterioplankton communities in 30 subtropical estuaries were explored, and potential pathogenic bacteria in the estuaries were detected and profiled. Our results demonstrate here that bacterioplankton and pathogens show similar biogeographic patterns under the influence of salinity. Interestingly, heterogeneous selection dominated bacterioplankton assembly, while stochasticity dominated pathogen assembly. This study provides important information for future risk assessment of potential pathogenic bacteria as well as management in estuarine ecosystems.

Network analysis reveals significant joint effects of microplastics and tetracycline on the gut than the gill microbiome of marine medaka

Microplastics could accumulate and enrich antibiotics in the aquatic environment. Despite this, the joint effects of microplastics and antibiotics on aquatic organisms are not clear. Here, we investigated the changes of microbial interactions in both gill and gut of marine medaka exposed to polystyrene microbeads (PS) and/or tetracycline for 30 days by using co-occurrence network analysis based on 16S rRNA gene amplicon sequences. We found that the single and combined effects of PS and tetracycline were more profound on the gut than on the gill microbiome. SourceTracker analysis showed that the relative contributions from the gill microbiome to the gut microbiome increased under combined exposure. Moreover, the combined exposure reduced the complexity and stability of the gut microbial network more than those induced by any single exposure, suggesting the synergistic effects of PS and tetracycline on the gut microbiome. The PS and tetracycline combined exposure also caused a shift in the keystone taxa of the gut microbial network. However, no similar pattern was found for gill microbial networks. Furthermore, single and combined exposure to PS and/or tetracycline altered the associations between the gut network taxa and indicator liver metabolites. Altogether, these findings enhanced our understanding of the hazards of the co-occurring environmental microplastics and antibiotics to the fish commensal microbiome.

Seasonal hydrological dynamics govern lifestyle preference of aquatic antibiotic resistome

Antibiotic resistance genes (ARGs) are a well-known environmental concern. Yet, limited knowledge exists on the fate and transport of ARGs in deep freshwater reservoirs experiencing seasonal hydrological changes, especially in the context of particle-attached (PA) and free-living (FL) lifestyles. Here, the ARG profiles were examined using high-throughput quantitative PCR in PA and FL lifestyles during four seasons representing two hydrological phenomena (vertical mixing and thermal stratification) in the Shuikou Reservoir (SR), Southern China. The results indicated that seasonal hydrological dynamics were critical for influencing the ARGs in PA and FL and the transition of ARGs between the two lifestyles. ARG profiles both in PA and FL were likely to be shaped by horizontal gene transfer. However, they exhibited distinct responses to the physicochemical (e.g., nutrients and dissolved oxygen) changes under seasonal hydrological dynamics. The particle-association niche (PAN) index revealed 94 non-conservative ARGs (i.e., no preferences for PA and FL) and 23 and 16 conservative ARGs preferring PA and FL lifestyles, respectively. A sharp decline in conservative ARGs under stratified hydrologic suggested seasonal influence on the ARGs transition between PA and FL lifestyles. Remarkably, the conservative ARGs (in PA or FL lifestyle) were more closely related to bacterial OTUs in their preferred lifestyle than their counterparts, indicating lifestyle-dependent ARG enrichment. Altogether, these findings enhanced our understanding of the ARG lifestyles and the role of seasonal hydrological changes in governing the ARG transition between the lifestyles in a typical deep freshwater ecosystem.

Continuous antibiotic attenuation in algal membrane photobioreactor: Performance and kinetics

The attenuation of 10 mixed antibiotics along with nutrients in a continuous flow mode by four freshwater microalgae (Haematococcus pluvialis, Selenastrum capricornutum, Scenedesmus quadricauda, and Chlorella vulgaris) was examined in membrane photobioreactors (MPBRs). At lab-scale, consistent removal of both antibiotic and nutrient was shown by H. pluvialis and S. quadricauda, respectively. The system exhibited better performance with enhanced removal at HRT 24 h compared to 12 h and 48 h. The highest removal efficiency of antibiotics was observed in H. pluvialis MPBR, with the mean antibiotic removal values of 53.57%- 96.33%. Biodegradation was the major removal pathway of the antibiotics in the algal-MPBR (AMPBR), while removal by bioadsorption, bioaccumulation, membrane rejection, and abiotic was minor. Then, the bacterial feature was studied and showed significant influence from system hydrodynamics. The kinetics of continuous flow antibiotic removal followed Stover-Kincannon and Grau second-order models, which revealed great potential of AMPBR to withstand antibiotic load. The latter coupled with the computational fluid dynamic simulation was successfully applied for the residual antibiotic prediction and potential system optimization. Overall, these results provide an important reference for continuous flow antibiotic removal using AMPBR.

Repeated introduction of micropollutants enhances microbial succession despite stable degradation patterns

The increasing-volume release of micropollutants into natural surface waters has raised great concern due to their environmental accumulation. Persisting micropollutants can impact multiple generations of organisms, but their microbially-mediated degradation and their influence on community assembly remain understudied. Here, freshwater microbes were treated with several common micropollutants, alone or in combination, and then transferred every 5 days to fresh medium containing the same micropollutants to mimic the repeated exposure of microbes. Metabarcoding of 16S rRNA gene makers was chosen to study the succession of bacterial assemblages following micropollutant exposure. The removal rates of micropollutants were then measured to assess degradation capacity of the associated communities. The degradation of micropollutants did not accelerate over time but altered the microbial community composition. Community assembly was dominated by stochastic processes during early exposure, via random community changes and emergence of seedbanks, and deterministic processes later in the exposure, via advanced community succession. Early exposure stages were characterized by the presence of sensitive microorganisms such as Actinobacteria and Planctomycetes, which were then replaced by more tolerant bacteria such as Bacteroidetes and Gammaproteobacteria. Our findings have important implication for ecological feedback between microbe-micropollutants under anthropogenic climate change scenarios.

Domestic wastewater causes nitrate pollution in an agricultural watershed, China

Excessive quantities of nitrates in the aquatic environment can cause eutrophication and raise water safety concerns. Therefore, identification of the sources of nitrate is crucial to mitigate nitrate pollution and for better management of the water resources. Here, the spatiotemporal variations and sources of nitrate were investigated by stable isotopes (δ¹⁵N and δ¹⁸O), hydrogeochemical variables (e.g., NO₃^- and Cl^-), and exogenous microbial signals (i.e., sediments, soils, domestic and swine sewage) in an agricultural watershed (Changle River watershed) in China. The concentration ranges of δ¹⁵N- and δ¹⁸O-NO₃^- between 3.03‰-18.97‰ and -1.55‰-16.47‰, respectively, suggested that soil nitrogen, chemical fertilizers, and manure and sewage (M&S) were the primary nitrate sources. Bayesian isotopic mixing model suggested that the major proportion of nitrate within the watershed (53.12 ± 10.40% and 63.81 ± 15.08%) and tributaries (64.43 ± 5.03% and 76.20 ± 4.34%) were contributed by M&S in dry and wet seasons, respectively. Community-based microbial source tracking (MST) showed that untreated and treated domestic wastewater was the major source (>70%) of river microbiota. Redundancy analysis with the incorporation of land use, hydrogeochemical variables, dual stable isotope, and exogenous microbial signals revealed domestic wastewater as the dominant cause of nitrate pollution. Altogether, this study not only identifies and quantifies the spatiotemporal variations in nitrate sources in the study area but also provides a new analytical framework by combining nitrate isotopic signatures and community-based MST approaches for source appointment of nitrate in other polluted watersheds.

Spatial autocorrelation and temporal variation of contaminants of emerging concern in a typical urbanizing river

The distribution and fate of contaminants of emerging concern (CECs) was studied in relation to hydrological conditions, land use characteristics, and spatial contiguity in Houxi River. Thirty-four CECs were detected in the surface water during a three-year sampling campaign. Caffeine was most prevalent (99% frequency), while bisphenol A had the highest median concentration (78.2 ng/L) among the detected CECs. Caffeine and the other prevalent CECs lincomycin and bisphenol A, with median concentrations of 3.89 ng/L, 0.26 ng/L, and 78.2 ng/L, respectively, were positively correlated with land use types related to anthropogenic activities (grass, barren, built up, and cropland areas and landscape indexes for human activities). The analysis of similarities revealed significant annual variations, with increasing trends in both the concentrations and detection frequencies of CECs. Spatial variations were demonstrated by higher concentrations and detection frequencies downstream compared to upstream. The singular value decomposition analysis revealed that the downstream sites were the major contributors (55.6%-100%) to the spatial variability of most CECs. Moran's I analysis based on downstream contiguity indicated strong spatial autocorrelation among the connected sites for most CECs. This was further supported by longer correlation lengths for 18 CECs than the average distance between the sampling sites. The spatial autocorrelation can be attributed to the physicochemical properties of CECs and local hydrological dynamics, including temperature, wind speed, and sunshine hours. For most CECs, local contribution predominated over neighbor influence with an average value of 75.5%. The results of this study provide new insight to evaluate CEC distributions, which will be beneficial to policymakers for the management and prioritization of CEC contaminants in the Houxi watershed.

A Comprehensive Profile of Antibiotic Resistance Genes in the Water Column of a Shallow-Sea Hydrothermal Vent Ecosystem

Antibiotic resistance genes (ARGs) became an emerging contaminant, and were found to accumulate in natural and man-made environments. A comprehensive understanding of the diversity and abundance of ARGs in pristine environments is critical for defining the baseline levels of environmental ARGs. However, there is limited information available on the ARG profiles of pristine environments, especially for shallow-sea hydrothermal vents ecosystems. Here, we combined 16S rRNA gene full-length amplicon sequencing and high-throughput quantitative PCR (HT-qPCR) to study the bacterial communities, and ARG abundance and diversity in the shallow-sea hydrothermal vent ecosystem of the Kueishantao Islet. The results of the 16S rRNA gene amplicon sequencing showed that several sulfur-cycling related bacterial genera, including Thiomicrorhabdus, Thioreductor, Sulfurovum, Sulfurimonas and Lebetimonas, dominated in the water column of the shallow-sea system. Temperature was the significant factor shaping the bacterial communities. The results of HT-qPCR analysis showed that the Kueishantao shallow-sea system harbored the lowest diversity (average 10 ARG subtypes) and abundance (average 1.0 × 10−3 copy per bacterial cell) of ARGs compared with other pristine (i.e., Tibet lake sediments, marine water and sediments) and anthropogenic-disturbed (i.e., drinking water reservoirs, urban ponds and wastewater treatment plants) environments. Procrustes analysis demonstrated a concordant pattern between the compositions of bacterial communities and ARGs in the shallow-sea system, while variation partition analysis revealed that the shared effects of physicochemical and bacterial communities explained >80% of the variation in the composition of ARGs. These results suggest that the vent bacterial communities and local environmental factors played an important role in shaping the distribution of the ARG profiles. Our study provides the first comprehensive overview of the background level of ARGs in a shallow-sea hydrothermal vent ecosystem.

Effects of Recreational Boating on Microbial and Meiofauna Diversity in Coastal Shallow Ecosystems of the Baltic Sea

Recreational boating can impact benthic ecosystems in coastal waters. Reduced height and cover of aquatic vegetation in shallow Baltic Sea inlets with high boat traffic have raised concerns about cascading effects on benthic communities in these ecosystems. Here, we characterized the diversity and composition of sediment-associated microbial and meiofaunal communities across five bays subjected to low and high degrees of boating activity and examined the community-environment relationships and association with bay morphometry. We found that recreational boating activity altered meiofauna alpha diversity and the composition of both micro- and meiobenthic communities, and there were strong correlations between community structure and morphometric variables like topographic openness, wave exposure, water surface area, and total phosphorous concentrations. Inlets with high boat traffic showed an increase of bacterial taxa like Hydrogenophilaceae and Burkholderiaceae. Several meiofauna taxa previously reported to respond positively to high levels of suspended organic matter were found in higher relative abundances in the bays with high boat traffic. Overall, our results show that morphometric characteristics of inlets are the strongest drivers of benthic diversity in shallow coastal environments. However, while the effects were small, we found significant effects of recreational boating on benthic community structure that should be considered when evaluating the new mooring projects.

IMPORTANCE With the increase of recreational boating activity and development of boating infrastructure in shallow, wave-protected areas, there is growing concern for their impact on coastal ecosystems. In order to properly assess the effects and consider the potential for recovery, it is important to investigate microbial and meiofaunal communities that underpin the functioning of these ecosystems. Here, we present the first study that uses DNA metabarcoding to assess how benthic biodiversity in shallow coastal areas is impacted by recreational boating. Our study shows a relatively small, but significant, effect of recreational boating both on meiofauna alpha diversity and meiofauna and bacterial community composition. However, both meiofauna and bacterial community composition in shallow benthic habitats is mediated to a higher degree by abiotic variables, such as topographic openness, area or size of the inlets, and wave exposure. Despite the fact that the effects were small, such impacts on benthic biodiversity should be considered in the management of coastal shallow habitats.

Abundance-Occupancy Relationships Along Taxonomic Ranks Reveal a Consistency of Niche Differentiation in Marine Bacterioplankton With Distinct Lifestyles

Abundance-occupancy relationships (AORs) are an important determinant of biotic community dynamics and habitat suitability. However, little is known about their role in complex bacterial communities, either within a phylogenetic framework or as a function of niche breadth. Based on data obtained in a field study in the St. Lawrence Estuary, we used 16S rRNA gene sequencing to examine the vertical patterns, strength, and character of AORs for particle-attached and free-living bacterial assemblages. Free-living communities were phylogenetically more diverse than particle-attached communities. The dominant taxa were consistent in terms of their presence/absence but population abundances differed in surface water vs. the cold intermediate layer. Significant, positive AORs characterized all of the surveyed communities across all taxonomic ranks of bacteria, thus demonstrating an ecologically conserved trend for both free-living and particle-attached bacteria. The strength of the AORs was low at the species level but higher at and above the genus level. These results demonstrate that an assessment of the distributions and population densities of finely resolved taxa does not necessarily improve determinations of apparent niche differences in marine bacterioplankton communities at regional scales compared with the information inferred from a broad taxonomic classification.

Urban ponds as hotspots of antibiotic resistome in the urban environment

The occurrence, dissemination and assembly processes of antibiotic resistance genes (ARGs) in urban water ecosystems are far from being understood. Here, we examined the diversity and abundance of ARGs in urban water ecosystems including landscape ponds, drinking water reservoirs, influents (IFs) and effluents (EFs) of wastewater treatment plants of a coastal city, China through high-throughput quantitative PCR. A total of 237 ARGs were identified, where multidrug, aminoglycoside and beta-lactamase resistance genes were the most abundant. Urban ponds had a comparatively high diversity and large numbers of shared ARGs with IFs and EFs. The average absolute abundance of ARGs (1.38 × 10⁷ copies/mL) and mobile genetic elements (MGEs) (4.19 × 10⁶ copies/mL) in ponds were only one order of magnitude lower than those of IFs, but higher than those of EFs and reservoirs. Stochastic processes dominated the ARG community assembly in IFs and ponds due to the random horizontal gene transfer caused by MGEs. These results imply that urban ponds are hotspots of ARGs. We further identified 25, 3, and 11 indicator ARGs for tracing the ARG contamination from IFs, EFs and ponds, respectively. Our study represents the first to highlight the role of urban ponds in the dissemination of ARGs.

Distinct mechanisms underlying the assembly of microeukaryotic generalists and specialists in an anthropogenically impacted river

Microeukaryotic communities are sensitive to environmental changes; and are considered essential for microbial food webs and biogeochemical cycles. Therefore, understanding the community responses of microeukaryotes to environmental changes is of great ecological significance. Very little is known about the assembly mechanisms underlying the microeukaryotic communities, especially for the key ecological groups (e.g., habitat generalists/specialists) in the riverine ecosystems. Here, we employed 18S rDNA amplicon sequencing to study the assembly processes governing the microeukaryotic communities and their habitat generalists and specialists across three hydrological seasons in a subtropical river in China. The results showed that deterministic and stochastic processes jointly shaped the microeukaryotic communities, where the relative importance of stochastic processes decreased in the following order: wet > normal > dry seasons. However, deterministic processes played more important role in shaping the microeukaryotic communities than those of prokaryotes. Meanwhile, stochastic and deterministic processes were responsible for structuring the microeukaryotic habitat generalists and specialists, respectively. Generally, the pure effects of physicochemical factors on the microeukaryotic communities and their key ecological groups were ~ 1.7 folds than those of the micropollutants. However, several micropollutants (e.g., acetaminophen, benzophenone-3, bisphenol A, sulfadiazine, triclocarban and triclosan) were found to have a strong effect on the habitat specialists. Altogether, we suggested that the ecological responses of the riverine microeukaryotes to micropollutants may vary among species due to the intrinsic differences in their environmental plasticity.