High biodiversity and distinct assembly patterns of microbial communities in groundwater compared with surface water

Exploring microbial diversity in relation to hydrological signatures and anthropic pressures in the Beauce Aquifer, France

Groundwater is one of the largest reservoirs of liquid fresh water on earth. Aquifers are complex environments where water quality can be affected by geological origins, geochemistry and biological activity. These environments are under pressure due to climate change and pollution. Strongly involved in the biogeochemistry of groundwater ecosystems, microorganisms can also contribute to reducing pollution and potentially be used as indicators of water quality. However, the diversity and role of microorganisms in groundwater remains largely unknown. The present work aimed to characterize bacterial diversity in 60 wells over 8000 km2 in the Beauce Aquifer, France, determine the baseline diversity expected in relation to the hydro geochemical signatures and identify variations that could be due to anthropogenic pressures. We analysed amplicon sequence variants (ASVs) of the 16S rRNA gene as well as the abundance of functional genes associated to arsenic and nitrate (aioA, arsB, narG), organic carbon contents and the geochemical composition. Bacterial diversity and bacterial co-occurrences were found to change according to the hydrogeochemical signatures identified in the aquifer. A lower diversity was measured in confined aquifers, where low oxygen and reducing conditions occurred with high levels of dissolved iron and manganese. This could indicate a higher vulnerability in the face of pressure. Finally, a baseline bacterial diversity composition was found to be linked to each hydrogeochemical signature and several wells that differed from the expected variations of this diversity were identified.

Running against the clock: exploring microbial diversity in an extremely endangered microbial oasis in the Chihuahuan Desert

The Cuatro Ciénegas Basin is a biodiversity hotspot known for its unique biodiversity. However, this ecosystem is facing severe anthropogenic threats that are drying its aquatic systems. We investigated microbial communities at three sites with different physicochemical and environmental characteristics (Pozas Rojas, Archean Domes, and the Churince system) within the basin to explore potential connections to deep aquifers and determine if the sites shared microorganisms. Utilizing 16S rRNA gene data, we identified a core microbiota between Pozas Rojas (PR) and Archean Domes (AD). Sulfur reduction appears to shape the microbial connectivity among sites, since sulfur-reducing bacteria has the highest prevalence between samples from PR and AD: Halanaerobium sp. (88.46%) and Desulfovermiculus halophilus (65%); and between the Churince system and AD: Halanaerobium sp. (63%) and D. halophilus (60%). Furthermore, metagenome-assembled genomes from Ectothiorhodospira genus were found in both AD and Churince, suggesting microbial dispersal. An important finding is that microbial diversity in the AD system declined, from 2016 to 2023 the ecosystem lost 29 microbial phyla. If this trend continues, the basin will lose most of its water, resulting in the loss of various prokaryotic lineages and potential biotechnological solutions, such as enzymes or novel antibiotics. Our findings highlighting the need for water extraction regulations to preserve the basin's biodiversity.

Spatio-temporal dynamics of bacterial community composition in a Western European watershed, the Meuse River watershed

This study aimed to identify factors influencing bacterial diversity in the Meuse River watershed by analyzing 42 locations sampled in spring and summer 2019, combined with biweekly sampling of one mid-stream location for a year. Bacterial community composition (BCC) was assessed in the small (SF; <5 µm) and large fractions (LF; ≥5 µm,), alongside physico-chemical parameters. LF consistently exhibited greater alpha diversity than SF. During the spatial campaigns, alpha diversity increased downstream in spring with high discharge, and BCC differed significantly between headwaters and the main river. Along this axis, several genera, Flavobacterium, Limnohabitans, and Aquirufa stood out as indicators of good water quality. Rhodoferax, another taxon indicative of good water quality, prevailed in the headwaters and during winter. In contrast, two cyanobacteria genera indicators of poor river quality, Microcystis PCC 7914 and Cyanobium PCC 6307, peaked in summer. BCC in spring and summer temporal samples aligned with spatial ones, while winter and autumn samples had distinct BCC. Finally, season, temperature, and distance from river mouth were the main driving parameters of beta diversity, outweighing the effect of fraction size on the BCC. These findings reinforce the notion that local conditions exert significant influence on bacterial communities in rivers.

Contribution of the microbial community to operational stability in an anammox reactor: Neutral theory and functional redundancy perspectives

A comprehensive understanding of microbial assembly is essential for achieving stable performance in biological wastewater treatment. Nevertheless, few studies have quantified these phenomena in detail, particularly in anammox-based processes. This study integrated mathematical and microbial approaches to analyze a 330-day anammox reactor with stable nitrogen removal efficiency (97 - 99%) despite changes in the high nitrogen loading rate, nitrogen concentration, and hydraulic retention time. A high value of functional redundancy (0.82) was obtained, with 84.6% of the microbial species following the neutral community model in stochastic processes, thus maintaining the stability of the dominant species and function in the microbial community. This study represents an initial attempt to quantify and evaluate the importance of functional redundancy in an anammox reactor. Based on these findings, engineering strategies have also been proposed to preserve high functional redundancy in stabilizing system performance under varying operating conditions.

Physicochemical and isotopic similarity between well water and intruding surface water is not synonymous with similarity in prokaryotic diversity and community composition

Intruding surface water can impact the physicochemical and microbiological quality of groundwater. Understanding these impacts is important because groundwater provides much of the world's potable water, and reduced quality is a potential public health risk. In this study, we monitored six shallow groundwater wells and three surface water bodies in the North Ostrobothnia region of Finland twice monthly for 12 months (October 2021-October 2022) via (i) on-site and off-site measurements of physicochemical water quality parameters, (ii) determination of stable water isotope compositions, and (iii) analysis of microbial communities (via amplicon sequencing of the V3-V4 16S rRNA gene sub-regions). Water from one well showed clear overall physicochemical and isotopic similarity with a nearby pond, as well as temporal fluctuations in water temperature and isotopes that mirrored those of the pond. Isotope mixing analyses suggested that about 80-95 % of the well water comes from the pond. Such large-scale intrusion might be expected to reduce prokaryotic diversity and composition in the aquifer, either by strong influx of surface water taxa or changes to aquifer physicochemistry. Compared to the pond, however, prokaryotic communities from the well showed significantly higher alpha diversity and a composition more similar to a nearby well unaffected by intrusion. The finding that physicochemical and isotopic similarity between well water and intruding surface water is not synonymous with similarity in prokaryotic diversity and community composition makes clear the need for a multi-method approach when studying the impact of surface water intrusion on shallow wells.

Metabarcoding reveals ecologically distinct fungal assemblages in river and groundwater along an Austrian alpine to lowland gradient

Biodiversity, the source of origin, and ecological roles of fungi in groundwater are to this day a largely neglected field in fungal and freshwater ecology. We used DNA-based Illumina high-throughput sequence analysis of both fungal gene markers 5.8S and internal transcribed spacers region 2 (ITS2), improving taxonomic classification. This study focused on the groundwater and river mycobiome along an altitudinal and longitudinal transect of a pre-alpine valley in Austria in two seasons. Using Bayesian network modeling approaches, we identified patterns in fungal community assemblages that were mostly shaped by differences in landscape (climatic, topological, and geological) and environmental conditions. While river fungi were comparatively more diverse, unique fungal assemblages could be recovered from groundwater, including typical aquatic lineages such as Rozellomycota and Olpidiomycota. The most specious assemblages in groundwater were not linked to the input of organic material from the surface, and as such, seem to be sustained by characteristic groundwater conditions. Based on what is known from closely related fungi, our results suggest that the present fungal communities potentially contribute to mineral weathering, carbon cycling, and denitrification in groundwater. Furthermore, we were able to observe the effects of varying land cover due to agricultural practices on fungal biodiversity in groundwater ecosystems. This study contributes to improving our understanding of fungi in the subsurface aquatic biogeosphere.

Dispersal limitation determines the ecological processes that regulate the seasonal assembly of bacterial communities in a subtropical river

Bacteria play a crucial role in pollutant degradation, biogeochemical cycling, and energy flow within river ecosystems. However, the underlying mechanisms governing bacterial community assembly and their response to environmental factors at seasonal scales in subtropical rivers remain poorly understood. In this study, we conducted 16S rRNA gene amplicon sequencing on water samples from the Liuxi River to investigate the composition, assembly processes, and co-occurrence relationships of bacterial communities during the wet season and dry season. The results demonstrated that seasonal differences in hydrochemistry significantly influenced the composition of bacterial communities. A more heterogeneous community structure and increased alpha diversity were observed during the dry season. Water temperature emerged as the primary driver for seasonal changes in bacterial communities. Dispersal limitation predominantly governed community assembly, however, during the dry season, its contribution increased due to decreased immigration rates. Co-occurrence network analysis reveals that mutualism played a prevailing role in shaping bacterial community structure. Compared to the wet season, the network of bacterial communities exhibited higher modularity, competition, and keystone species during the dry season, resulting in a more stable community structure. Although keystone species displayed distinct seasonal variations, Proteobacteria and Actinobacteria were consistently abundant keystone species maintaining network structure in both seasons. Our findings provide insights into how bacterial communities respond to seasonal environmental changes, uncovering underlying mechanisms governing community assembly in subtropical rivers, which are crucial for the effective management and conservation of riverine ecosystems.

Response of microbial communities in aquifers with multiple organic solvent contamination: Implications for MNA remedy

The application of Monitored Natural Attenuation (MNA) technology has been widespread, while there is a paucity of data on groundwater with multiple co-contaminants. This study focused on high permeability, low hydraulic gradient groundwater with co-contamination of benzene, toluene, ethylbenzene, and xylenes (BTEX), chlorinated aliphatic hydrocarbons (CAHs), and chlorinated aromatic hydrocarbons (CPs). The objective was to investigate the responses of microbial communities during natural attenuation processes. Results revealed greater horizontal variation in groundwater microbial community composition compared to vertical variation. The variation was strongly correlated with the total contaminant quantity (r = 0.722, p < 0.001) rather than individual contaminants. BTEX exerted a more significant influence on community diversity than other contaminants. The assembly of groundwater microbial communities was primarily governed by deterministic processes (βNTI < -2) in high contaminant concentration zones, while stochastic processes (|βNTI| < 2) dominated in low-concentration zones. Moreover, the microbial interactions shifted at different depths indicating the degradation rate variation in the vertical. This study makes fundamental contribution to the understanding for the effects of groundwater flow and material fields on indigenous microbial communities, which will provide a scientific basis for more precise adoption of microbial stimulation/augmentation to accelerate the rate of contaminant removal.

Using the Potential Transformation of Dissolved Organic Matter to Understand Carbon Emissions from Inland Rivers

Understanding the transformation of river dissolved organic matter (DOM) is important for assessing the emissions of greenhouse gases (GHGs) in inland waters. However, the relationships between the variations in DOM components and GHGs remain largely unknown. Here, parallel factor analysis (PARAFAC) was applied to investigate the DOM components in 46 inland rivers in China. We found that the GHG emissions in peri-urban rivers were 1.10-2.15 times greater than those in urban rivers. Microbial and environmental factors (e.g., living cell numbers, microbial activity and pH) explained more than 70% of the total variance in GHG emissions in rivers. DOM variations relationships between different components ware revealed based on compositional data principal component analysis (CoDA-PCA). Microbial-mediated DOM production and degradation were quantified, and the degradation levels in peri-urban rivers were 11.8-25.2% greater than those in urban rivers. Differences in carbon emission potential between urban and peri-urban rivers were related to DOM variances and transformations and were affected by water chemistry (e.g., NH4-N and As). This study clarifies the regulatory effects of DOM composition variations and transformations on GHG emissions, and enhances the understanding of the DOM biogeochemical cycle.

Insight into spatial variations of DOM fractions and its interactions with microbial communities of shallow groundwater in a mesoscale lowland river watershed

Dissolved organic matter (DOM) plays a crucial role in driving biogeochemical processes and determining water quality in shallow groundwater systems, where DOM could be susceptible to dynamic influences of surface water influx. This study employed fluorescence excitation-emission matrix (EEM) spectroscopy combined with principal component coefficients, parallel factor analysis (PARAFAC), co-occurrence network analysis and structural equation modeling (SEM) to examine changes of DOM fractions from surface water to shallow groundwater in a mesoscale lowland river basin. Combining stable isotope and hydrochemical parameters, except for surface water (SW), two groups of groundwater samples were defined, namely, deeply influenced by surface water (IGW) and groundwater nearly non-influenced by surface water (UGW), which were 50.34 % and 19.39 % recharged by surface water, respectively. According to principal component coefficients, reassembled EEM data of these categories highlighted variations of the tyrosine-like peak in DOM. EEMs coupled with PARAFAC extracted five components (C1-C5), i.e. C1, protein-like substances, C2 and C4, humic-like substances, and C3 and C5, microbial-related substances. The abundance of the protein-like was SW > IGW > UGW, while the order of the humic-like was opposite. The bacterial communities exhibited an obvious cluster across three regions, which hinted their sensitivity to variations in environmental conditions. Based on co-occurrence, SW represented the highest connectivity between bacterial OTUs and DOM fractions, followed by IGW and UGW. SEM revealed that microbial activities increased bioavailability of the humic-like in the SW and IGW, whereas microbial compositions promoted the evolution of humic-like substances in the UGW. Generally, these results could be conducive to discern dissimilarity in DOM fractions across surface water and shallow groundwater, and further trace their interactions in the river watershed.

Biostimulation with oxygen and electron donors supports micropollutant biodegradation in an experimentally simulated nitrate-reducing aquifer

The availability of suitable electron donors and acceptors limits micropollutant natural attenuation in oligotrophic groundwater. This study investigated how electron donors with different biodegradability (humics, dextran, acetate, and ammonium), and different oxygen concentrations affect the biodegradation of 15 micropollutants (initial concentration of each micropollutant = 50 μg/L) in simulated nitrate reducing aquifers. Tests mimicking nitrate reducing field conditions showed no micropollutant biodegradation, even with electron donor amendment. However, 2,4-dichlorophenoxyacetic acid and mecoprop were biodegraded under (micro)aerobic conditions with and without electron donor addition. The highest 2,4-dichlorophenoxyacetic acid and mecoprop biodegradation rates and removal efficiencies were obtained under fully aerobic conditions with amendment of an easily biodegradable electron donor. Under microaerobic conditions, however, amendment with easily biodegradable dissolved organic carbon (DOC) inhibited micropollutant biodegradation due to competition between micropollutants and DOC for the limited oxygen available. Microbial community composition was dictated by electron acceptor availability and electron donor amendment, not by micropollutant biodegradation. Low microbial community richness and diversity led to the absence of biodegradation of the other 13 micropollutants (such as bentazon, chloridazon, and carbamazepine). Finally, adaptation and potential growth of biofilms interactively determined the location of the micropollutant removal zone relative to the point of amendment. This study provides new insight on how to stimulate in situ micropollutant biodegradation to remediate oligotrophic groundwaters as well as possible limitations of this process.

Spatial variations of fungal community assembly and soil enzyme activity in rhizosphere of zonal Stipa species in inner Mongolia grassland

Rhizosphere soil fungal and enzyme activities affect the nutrient cycling of terrestrial ecosystems, and rhizosphere fungi are also important participants in the ecological process of vegetation succession, responding to changes in plant communities. Stipa is an excellent forage grass with important ecological and economic value, and has the spatial distribution pattern of floristic geographical substitution. In order to systematically investigate the synergistic response strategies of fungal communities and enzyme activities in the rhizosphere under the vegetation succession. Here we explored the turnover and assembly mechanisms of Stipa rhizosphere fungal communities and the spatial variation of metabolic activity under the succession of seven Stipa communities in northern China grassland under large scale gradients. The results indicated that the composition, abundance and diversity of fungal communities and microbial enzyme activities in rhizosphere soil differed among different Stipa species and were strikingly varied along the Stipa community changes over the geographic gradient. As the geographical distribution of Stipa community changed from east to west in grassland transect, Mortierellomycetes tended to be gradually replaced by Dothideomycetes. The null models showed that the rhizosphere fungal communities were governed primarily by the dispersal limitation of stochastic assembly processes, which showed decreased relative importance from S. grandis to S. gobica. Moreover, the MAT and MAP were the most important factors influencing the changes in the fungal community (richness, β-diversity and composition) and fungal community assembly, while SC and NP also mediated fungal community assembly processes. These findings deepen our understanding of the responses of the microbial functions and fungal community assembly processes in the rhizosphere to vegetation succession.

Bacterial diversity in a continuum from supraglacial habitats to a proglacial lake on the Tibetan Plateau

Mountain glaciers are frequently assessed for their hydrological connectivity from glaciers to proglacial lakes. Ecological process on glacier surfaces and downstream ecosystems have often been investigated separately, but few studies have focused on the connectivity between the different glacial habitats. Therefore, it remains a limited understanding of bacterial community assembly across different habitats along the glacier hydrological continuum. In this study, we sampled along a glacial catchment from supraglacial snow, cryoconite holes, supraglacial runoff, ice-marginal moraine and proglacial lake on the Tibetan Plateau. The bacterial communities in these habitats were analyzed using high-throughput DNA sequencing of the 16S rRNA gene to determine the bacterial composition and assembly. Our results showed that each habitat hosted unique bacterial communities, with higher bacterial α-diversity in transitional habitats (e.g. runoff and ice-marginal moraine). Null model analysis indicated that deterministic processes predominantly shaped bacterial assembly in snow, cryoconite holes and lake, while stochastic process dominantly governed bacterial community in transitional habitats. Collectively, our findings suggest that local environment play a critical role in filtering bacterial community composition within glacier habitats. This study enhances our understanding of microbial assembly process in glacier environments and provides valuable insights into the factors governing bacterial community compositions across different habitats along the glacial hydrological continuum.

Vertical and seasonal dynamics of bacterial pathogenic communities at an aged organic contaminated site: Insights into microbial diversity, composition, interactions, and assembly processes

Under the background of the Coronavirus Disease 2019 (COVID-19) pandemic, research on pathogens deserves greater attention in the natural environment, especially in the widely distributed contaminated sites with complicated and severe organic pollution. In this study, the community composition and assembly of soil pathogens identified by the newly-developed 16S-based pipeline of multiple bacterial pathogen detection (MBPD) have been investigated on spatiotemporal scales in the selected organic polluted site. We demonstrated that the richness and diversity of the pathogenic communities were primarily controlled by soil depth, while the structure and composition of pathogenic communities varied pronouncedly with seasonal changes, which were driven by the alterations in both physiochemical parameters and organic contaminants over time. Network analysis revealed that the overwhelmingly positive interactions, identified multiple keystone species, and a well-organized modular structure maintained the stability and functionality of the pathogenic communities under environmental pressures. Additionally, the null-model analysis showed that deterministic processes dominated the pathogenic community assembly across soil profiles. In three seasons, stochasticity-dominated processes in spring and summer changed into determinism-dominated processes in winter. These findings extend our knowledge of the response of the bacterial pathogenic community to environmental disruptions brought on by organic contaminated sites.

A global perspective on compositions, risks, and ecological genesis of antibiotic resistance genes in biofilters of drinking water treatment plants

Antibiotic resistance genes (ARGs) in biofilters of drinking water treatment plants (DWTPs) are regarded to be a remarkable potential health risk to human. A global survey on ARGs in biofilters may help evaluate their risk features as a whole. This study aims to explore the compositions, risks, and ecological genesis of ARGs in the biofilters of DWTPs. In total, 98 metagenomes of DWTP biofilters were collected from Sequence Read Archive  (SRA)  of National center for Biotechnology Information (NCBI), and the main ARG types were recognized, with multidrug, bacitracin, and beta-lactam as the first three types. Source water types (surface water vs. groundwater) were found to significantly influence antibiotic resistome, overpassing biofilter media and locations. Although ARG abundances of surface water biofilters were approximately five times higher than that of groundwater biofilters, the risk pattern of ARGs was highly similar between surface water biofilters and groundwater biofilters, and up to 99.61% of the ARGs on average belong to the least risk and unassessed ranks, and only 0.23% the highest risk rank. Monobactam biosynthesis pathway and prodigiosin biosynthesis pathway, two antibiotics biosynthesis pathways, were observed to be positively correlated with several ARG types and total ARG abundance in samples of surface water and groundwater biofilters, respectively, suggesting their potential roles in ecological genesis of ARGs. Overall, the results of this study would deepen our understanding of the ARG risks in biofilters of DWTPs and shed light on their ecological genesis inside.

From the Mountain to the Valley: Drivers of Groundwater Prokaryotic Communities along an Alpine River Corridor

Rivers are the “tip of the iceberg”, with the underlying groundwater being the unseen freshwater majority. Microbial community composition and the dynamics of shallow groundwater ecosystems are thus crucial, due to their potential impact on ecosystem processes and functioning. In early summer and late autumn, samples of river water from 14 stations and groundwater from 45 wells were analyzed along a 300 km transect of the Mur River valley, from the Austrian alps to the flats at the Slovenian border. The active and total prokaryotic communities were characterized using high-throughput gene amplicon sequencing. Key physico-chemical parameters and stress indicators were recorded. The dataset was used to challenge ecological concepts and assembly processes in shallow aquifers. The groundwater microbiome is analyzed regarding its composition, change with land use, and difference to the river. Community composition and species turnover differed significantly. At high altitudes, dispersal limitation was the main driver of groundwater community assembly, whereas in the lowland, homogeneous selection explained the larger share. Land use was a key determinant of the groundwater microbiome composition. The alpine region was more diverse and richer in prokaryotic taxa, with some early diverging archaeal lineages being highly abundant. This dataset shows a longitudinal change in prokaryotic communities that is dependent on regional differences affected by geomorphology and land use.

Spatial distributions of microbial diversity in the contaminated deep groundwater: A case study of the Huaibei coalfield

The impact of coal mining activities on the structure of groundwater microbial communities in coal mining areas has gradually received academic attention. In this study, hydrochemical analysis and sequencing of the V4 region of the 16S rRNA gene were carried out using groundwater samples from the fourth aquifer in the loose layer of Cenozoic, the sandstone fissure aquifer in the coal measure strata of Permian, the Carboniferous Taiyuan Formation limestone aquifer, and the Ordovician limestone aquifer, at depths of 250 m, 600 m, 750 m, and 1000 m in monitoring wells. Results showed that the operational taxonomy units (OTUs) in the deep groundwater ecosystem were clustered distinguishably between the contaminated and the uncontaminated aquifers. The microbial community alpha-diversity of groundwater was significantly correlated with depth, and the richness of microbial community composition decreased with increasing depth. Proteobacteria (34.41%-97.41%), was found to be the dominant phylum, Gammaproteobacteria (10.05%-92.06%) was the dominant class and "Unassigned" (4.12%-64.72%) was dominant at the genus level. The number of endemic bacteria in the four aquifers was 1, 33, 99 and 11, respectively. It was also found that F^-, oxidation-reduction potential (ORP), and TOC were the main environmental variables affecting the groundwater all OTUs, abundant OTUs, and rare OTUs, respectively. These results indicate that the activity of rare OTU subcommunities increases gradually with increasing aquifer depth and that mining significantly enriched Thiovirga in deep groundwater. In addition, it was found that with the increase of depth, the effect of ORP on microbial community abundance decreased. This study deepens our understanding of the evolution characteristics of microbial communities in deep groundwater in coal mining areas. The unique characteristics of microbial communities characteristics of four aquifers with different depths provide a microbial perspective for understanding the characteristics of deep aquifers.

Environmental and Anthropogenic Factors Shape the Skin Bacterial Communities of a Semi-Arid Amphibian Species

The amphibian skin microbiome is important in maintaining host health, but is vulnerable to perturbation from changes in biotic and abiotic conditions. Anthropogenic habitat disturbance and emerging infectious diseases are both potential disrupters of the skin microbiome, in addition to being major drivers of amphibian decline globally. We investigated how host environment (hydrology, habitat disturbance), pathogen presence, and host biology (life stage) impact the skin microbiome of wild Dhofar toads (Duttaphrynus dhufarensis) in Oman. We detected ranavirus (but not Batrachochytrium dendrobatidis) across all sampling sites, constituting the first report of this pathogen in Oman, with reduced prevalence in disturbed sites. We show that skin microbiome beta diversity is driven by host life stage, water source, and habitat disturbance, but not ranavirus infection. Finally, although trends in bacterial diversity and differential abundance were evident in disturbed versus undisturbed sites, bacterial co-occurrence patterns determined through network analyses revealed high site specificity. Our results therefore provide support for amphibian skin microbiome diversity and taxa abundance being associated with habitat disturbance, with bacterial co-occurrence (and likely broader aspects of microbial community ecology) being largely site specific.