Groundwater bacterial communities evolve over time in response to recharge

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.

Assessing long-term ecological impacts of PCE contamination in groundwater using a flow cytometric fingerprint approach

This study aims to develop and validate a comprehensive method for assessing ecological disturbances in groundwater ecosystems caused by tetrachloroethylene (PCE) contamination, utilizing flow cytometry (FCM) fingerprint approach. We hypothesized that the ecological disturbance resulting from PCE contamination would exhibit 'press disturbance', persisting over extended periods, and inducing notable phenotypic differences in the microbial community compared to undisturbed groundwater. We collected 40 groundwater samples from industrial district with a history of over twenty years of PCE contamination, along with 56 control groundwater from the national surveillance groundwater system. FCM revealed significant alterations in the phenotypic diversity of microbial communities in PCE-contaminated groundwater, particularly during the dry season. The presence of specific dechlorinating bacteria (Dehalococcoides, Dehalogenimonas, and Geobacter) and their syntrophic partners was identified as an indicator of contamination. Phenotypic diversity measures provided clearer and more direct reflections of contamination impact compared to taxonomic diversity measures. This study establishes FCM fingerprinting as a simple, robust, and accurate method for evaluating ecological disturbances, with potential applications in early warning systems and continuous monitoring of groundwater contamination. The findings not only underscore the sensitivity of FCM in detecting phenotypic variations induced by environmental stressors but also highlight its utility in understanding the complex dynamics of microbial communities in contaminated groundwater ecosystems.

Revealing the hierarchical structure of microbial communities

Measuring the dynamics of microbial communities results in high-dimensional measurements of taxa abundances over time and space, which is difficult to analyze due to complex changes in taxonomic compositions. This paper presents a new method to investigate and visualize the intrinsic hierarchical community structure implied by the measurements. The basic idea is to identify significant intersection sets, which can be seen as sub-communities making up the measured communities. Using the subset relationship, the intersection sets together with the measurements form a hierarchical structure visualized as a Hasse diagram. Chemical organization theory (COT) is used to relate the hierarchy of the sets of taxa to potential taxa interactions and to their potential dynamical persistence. The approach is demonstrated on a data set of community data obtained from bacterial 16S rRNA gene sequencing for samples collected monthly from four groundwater wells over a nearly 3-year period (n = 114) along a hillslope area. The significance of the hierarchies derived from the data is evaluated by showing that they significantly deviate from a random model. Furthermore, it is demonstrated how the hierarchy is related to temporal and spatial factors; and how the idea of a core microbiome can be extended to a set of interrelated core microbiomes. Together the results suggest that the approach can support developing models of taxa interactions in the future.

Iron coatings on carbonate rocks shape the attached bacterial aquifer community

Most studies of groundwater ecosystems target planktonic microbes, which are easily obtained via water samples. In contrast, little is known about the diversity and function of microbes adhering to rock surfaces, particularly to consolidated rocks. To investigate microbial attachment to rock surfaces, we incubated rock chips from fractured aquifers in limestone-mudstone alternations in bioreactors fed with groundwater from two wells representing oxic and anoxic conditions. Half of the chips were coated with iron oxides, representing common secondary mineralization in fractured rock. Our time-series analysis showed bacteria colonizing the chips within two days, reaching cell numbers up to 4.16 × 105 cells/mm2 after 44 days. Scanning electron microscopy analyses revealed extensive colonization but no multi-layered biofilms, with chips from oxic bioreactors more densely colonized than from anoxic ones. Estimated attached-to-planktonic cell ratios yielded values of up to 106: 1 and 103: 1, for oxic and anoxic aquifers, respectively. We identified distinct attached and planktonic communities with an overlap between 17 % and 42 %. Oxic bioreactors were dominated by proteobacterial genera Aquabacterium and Rhodoferax, while Rheinheimera and Simplicispira were the key players of anoxic bioreactors. Motility, attachment, and biofilm formation traits were predicted in major genera based on groundwater metagenome-assembled genomes and reference genomes. Early rock colonizers appeared to be facultative autotrophs, capable of fixing CO2 to synthesize biomass and a biofilm matrix. Late colonizers were predicted to possess biofilm degrading enzymes such as beta-glucosidase, beta-galactosidase, amylases. Fe-coated chips of both bioreactors featured more potential iron reducers and oxidizers than bare rock chips. As secondary minerals can also serve as energy source, they might favor primary production and thus contribute to subsurface ecosystem services like carbon fixation. Since most subsurface microbes seem to be attached, their contribution to ecosystem services should be considered in future studies.

Temporal profiling of sediment microbial communities in the Three Gorges Reservoir Area discovered time-dissimilarity patterns and multiple stable states

Microbial communities play vital roles in cycling nutrients and maintaining water quality in aquatic ecosystems. To better understand the dynamics of microbial communities and to pave way to effective ecological remediation, it's essential to reveal the temporal patterns of the communities and to identify their states. However, research exploring the dynamic changes of microbial communities needs a large amount of time-series data, which could be an extravagant requirement for a single study. In this research, we overcame this challenge by conducting a meta-analysis of years of accumulations of 16S rRNA high-throughput sequencing data from the Three Gorges Reservoir Area (TGRA), an ecological and environmental hotspot. For better understanding the microbial communities time-dissimilarity dynamics, three microbial communities time-dissimilarity patterns were hypothesized, and the linear pattern in the TGRA was validated. In addition, to explore the stability of microbial communities in the TGRA, two alternative stable states were revealed, and their differences in community richness, alpha diversity indices, community composition, ecological network topological properties, and metabolic functions were demonstrated. In short, two states of microbial communities showed distinct richness and alpha diversity indices, and the communities in one state were more dominated by Halomonas and Nitrosopumilaceae genera, facilitating nitrogen cycling metabolic processes; whilst the main genera of the other state were Bathyarchaeia and Methanosaeta, which favored methane-related metabolism. Moreover, different studies and environmental differences between mainstream and tributaries were attributed as the potential inducing factors of the state division. Our study provides a comprehensive insight into the dynamics and stability of microbial communities in the TGRA, and a reference for future studies on microbial community dynamics.

Microbial-induced MnO2 precipitation in a carbonate coastal aquifer

A study was carried out to identify biogeochemical reactions along a transect of a coastal dolomitic aquifer. In this transect, the physicochemical parameters of the groundwater as well as the microbial composition of samples taken at different depths and salinities were measured. Many of the dissolved ions measured in the groundwater follow a pattern that reflects the distribution of the water masses (fresh, interface and salt) in the aquifer, while others such as Ca and Mg ions deviate from this trend by identifying the zones of maximum dissolution of the carbonate matrix. The concentrations of minor ions, such as Fe and Mn, also follow a singular pattern, with maximum concentrations in the reducing zones of the aquifer and lower values in the oxidizing zones. Precipitates of Mn oxides along with other metals, such as Fe, Ba, Zn and Ni, were observed in the saline zone displaying oxidizing conditions close to the coastline, where a continuous core was recovered. This zone, which is located below the freshwater-seawater mixing zone and features percentages of seawater higher than 80 %, is characterized by the presence of Marinobacter as the predominant genus. These bacteria are also related to the formation of Mn-rich polymetallic oxides in other contexts such as the ocean floor (Wang et al., 2012; Cao et al., 2021). All in all, a biogeochemical reaction model is proposed that describes the formation of these oxides in areas close to the discharge zone of coastal aquifers. To do this, it has been necessary to integrate the results obtained from geochemical, hydrogeological and microbiological information.

Rapid Consumption of Dihydrogen Injected into a Shallow Aquifer by Ecophysiologically Different Microbes

The envisaged future dihydrogen (H2) economy requires a H2 gas grid as well as large deep underground stores. However, the consequences of an unintended spread of H2 through leaky pipes, wells, or subterranean gas migrations on groundwater resources and their ecosystems are poorly understood. Therefore, we emulated a short-term leakage incident by injecting gaseous H2 into a shallow aquifer at the TestUM test site and monitored the subsequent biogeochemical processes in the groundwater system. At elevated H2 concentrations, an increase in acetate concentrations and a decrease in microbial α-diversity with a concomitant change in microbial β-diversity were observed. Additionally, microbial H2 oxidation was indicated by temporally higher abundances of taxa known for aerobic or anaerobic H2 oxidation. After H2 concentrations diminished below the detection limit, α- and β-diversity approached baseline values. In summary, the emulated H2 leakage resulted in a temporally limited change of the groundwater microbiome and associated geochemical conditions due to the intermediate growth of H2 consumers. The results confirm the general assumption that H2, being an excellent energy and electron source for many microorganisms, is quickly microbiologically consumed in the environment after a leakage.

Redox gradients drive microbial community assembly patterns and molecular ecological networks in the hyporheic zone of effluent-dominated rivers

The impacts of effluent discharge on receiving waterbodies have been a research hotspot. Nonetheless, limited information is available on the microbial community assembly patterns in the hyporheic zones (HZ) responding to the changes in the microenvironments, e.g., solute transport and redox gradient variations. Using two representative effluent-dominated rivers as model systems, the spatio-temporal bacterial community dynamics and assembly patterns in oxic and suboxic zones in the shallow riverbed sediments were disentangled via null model- and neutral model-based approaches. Bacterial dynamics in community composition were observed driven by environmental filtering, i.e., impacts of environmental variables, more than geographic distances, i.e., the depths of sediments. The communities in samples collected in summer were largely shaped by stochasticity, in which homogeneous selection occupied a higher proportion in oxic (∼39%) than in suboxic zone (∼23%). Deterministic processes contributed to a more complex community structure for samples from oxic zones, whereas weakened the interspecies interactions in suboxic zones. The richness and abundances of non-neutral community were confirmed governing the deterministic assembly in oxic zones. Key species ascribed to 'connectors' and 'network hubs' dominated the community assembly variations in samples collected in winter, and in oxic zones, respectively. Significant positive relationships between β-nearest taxon index and dissolved organic nitrogen (DON) and nitrate highlighted their vital roles in community assembly via deterministic selective pressures in oxic zones. The significance thresholds of nitrogen species for community transition in winter (ΔDON: 2.81 mg-N/L, ΔNO3-: 1.09 mg-N/L) were lower than in summer, probably implying that stricter effluent quality standards should be established in colder seasons. Combined, our work poses first insights on the roles of redox zonation in driving microbial community assembly in HZ, which is of significance in guiding ecological remediation processes in effluent-dominated rivers.

Microbial dark matter sequences verification in amplicon sequencing and environmental metagenomics data

Although microorganisms constitute the most diverse and abundant life form on Earth, in many environments, the vast majority of them remain uncultured. As it is based on information gleaned mainly from cultivated microorganisms, our current body of knowledge regarding microbial life is partial and does not reflect actual microbial diversity. That diversity is hidden in the uncultured microbial majority, termed by microbiologists as "microbial dark matter" (MDM), a term borrowed from astrophysics. Metagenomic sequencing analysis techniques (both 16S rRNA gene and shotgun sequencing) compare gene sequences to reference databases, each of which represents only a small fraction of the existing microorganisms. Unaligned sequences lead to groups of "unknown microorganisms" that are usually ignored and rarefied from diversity analysis. To address this knowledge gap, we analyzed the 16S rRNA gene sequences of microbial communities from four different environments-a living organism, a desert environment, a natural aquatic environment, and a membrane bioreactor for wastewater treatment. From those datasets, we chose representative sequences of potentially unknown bacteria for additional examination as "microbial dark matter sequences" (MDMS). Sequence existence was validated by specific amplification and re-sequencing. These sequences were screened against databases and aligned to the Genome Taxonomy Database to build a comprehensive phylogenetic tree for additional sequence classification, revealing potentially new candidate phyla and other lineages. These putative MDMS were also screened against metagenome-assembled genomes from the explored environments for additional validation and for taxonomic and metabolic characterizations. This study shows the immense importance of MDMS in environmental metataxonomic analyses of 16S rRNA gene sequences and provides a simple and readily available methodology for the examination of MDM hidden behind amplicon sequencing results.

Karst spring microbiome: Diversity, core taxa, and community response to pathogens and antibiotic resistance gene contamination

Karst aquifers are important water resources for drinking water supplies worldwide. Although they are susceptible to anthropogenic contamination due to their high permeability, there is a lack of detailed knowledge on the stable core microbiome and how contamination may affect these communities. In this study, eight karst springs (distributed across three different regions in Romania) were sampled seasonally for one year. The core microbiota was analysed by 16S rRNA gene amplicon sequencing. To identify bacteria carrying antibiotic resistance genes and mobile genetic elements, an innovative method was applied, consisting of high-throughput antibiotic resistance gene quantification performed on potential pathogen colonies cultivated on Compact Dry™ plates. A taxonomically stable bacterial community consisting of Pseudomonadota, Bacteroidota, and Actinomycetota was revealed. Core analysis reaffirmed these results and revealed primarily freshwater-dwelling, psychrophilic/psychrotolerant species affiliated to Rhodoferax, Flavobacterium, and Pseudomonas genera. Both sequencing and cultivation methods indicated that more than half of the springs were contaminated with faecal bacteria and pathogens. These samples contained high levels of sulfonamide, macrolide, lincosamide and streptogramins B, and trimethoprim resistance genes spread primarily by transposase and insertion sequences. Differential abundance analysis found Synergistota, Mycoplasmatota, and Chlamydiota as suitable candidates for pollution monitoring in karst springs. This is the first study highlighting the applicability of a combined approach based on high-throughput SmartChip™ antibiotic resistance gene quantification and Compact Dry™ pathogen cultivation for estimating microbial contaminants in karst springs and other challenging low biomass environments.

Metabolic versatility enables sulfur-oxidizers to dominate primary production in groundwater

High rates of CO2 fixation and the genetic potential of various groundwater microbes for autotrophic activity have shown that primary production is an important source of organic C in groundwater ecosystems. However, the contribution of specific chemolithoautotrophic groups such as S-oxidizing bacteria (SOB) to groundwater primary production and their adaptation strategies remain largely unknown. Here, we stimulated anoxic groundwater microcosms with reduced S and sampled the microbial community after 1, 3 and 6 weeks. Genome-resolved metaproteomics was combined with 50at-% 13CO2 stable isotope probing to follow the C flux through the microbial food web and infer traits expressed by active SOB in the groundwater microcosms. Already after 7 days, 90% of the total microbial biomass C in the microcosms was replaced by CO2-derived C, increasing to 97% at the end of incubation. Stable Isotope Cluster Analysis revealed active autotrophs, characterized by a uniform 13C-incorporation of 45% in their peptides, to dominate the microbial community throughout incubation. Mixo- and heterotrophs, characterized by 10 to 40% 13C-incorporation, utilized the primarily produced organic C. Interestingly, obligate autotrophs affiliated with Sulfuricella and Sulfuritalea contained traits enabling the storage of elemental S in globules to maintain primary production under energy limitation. Others related to Sulfurimonas seemed to rapidly utilize substrates for fast proliferation, and most autotrophs further maximized their energy yield via efficient denitrification and the potential for H2 oxidation. Mixotrophic SOB, belonging to Curvibacter or Polaromonas, enhanced metabolic flexibility by using organic compounds to satisfy their C requirements. Time series data spanning eight years further revealed that key taxa of our microcosms composed up to 15% of the microbial groundwater community, demonstrating their in-situ importance. This showed that SOB, by using different metabolic strategies, are able to account for high rates of primary production in groundwater, especially at sites limited to geogenic nutrient sources. The widespread presence of SOB with traits such as S storage, H2 oxidation, and organic C utilization in many aquatic habitats further suggested that metabolic versatility governs S-fueled primary production in the environment.

Ecological differentiation and assembly processes of abundant and rare bacterial subcommunities in karst groundwater

The ecological health of karst groundwater has been of global concern due to increasing anthropogenic activities. Bacteria comprising a few abundant taxa (AT) and plentiful rare taxa (RT) play essential roles in maintaining ecosystem stability, yet limited information is known about their ecological differentiation and assembly processes in karst groundwater. Based on a metabarcoding analysis of 64 groundwater samples from typical karst regions in southwest China, we revealed the environmental drivers, ecological roles, and assembly mechanisms of abundant and rare bacterial communities. We found a relatively high abundance of potential functional groups associated with parasites and pathogens in karst groundwater, which might be linked to the frequent regional anthropogenic activities. Our study confirmed that AT was dominated by Proteobacteria and Campilobacterota, while Patescibacteria and Chloroflexi flourished more in the RT subcommunity. The node-level topological features of the co-occurrence network indicated that AT might share similar niches and play more important roles in maintaining bacterial community stability. RT in karst groundwater was less environmentally constrained and showed a wider environmental threshold response to various environmental factors than AT. Deterministic processes, especially homogeneous selection, tended to be more important in the community assembly of AT, whereas the community assembly of RT was mainly controlled by stochastic processes. This study expanded our knowledge of the karst groundwater microbiome and was of great significance to the assessment of ecological stability and drinking water safety in karst regions.

Effect of Snowmelt on Groundwater Bacterial Community Composition and Potential Role of Surface Environments as Microbial Seed Bank in Two Distinct Aquifers from the Region of Quebec, Canada

Events of groundwater recharge are associated with changes in the composition of aquifer microbial communities but also abiotic conditions. Modification in the structure of the community can be the result of different environmental condition favoring or hindering certain taxa, or due to the introduction of surface-derived taxa. Yet, in both cases, the local hydrogeochemical settings of the aquifer is likely to affect the amount of variation observed. Therefore, in our study, we used 16S rRNA gene sequencing to assess how microbial communities change in response to snowmelt and the potential connectivity between subsurface and surface microbiomes in two distinct aquifers located in the region of Vaudreuil-Soulanges (Québec, Canada). At both sites, we observed an increase in groundwater level and decrease in temperature following the onset of snow melt in March 2019. Bacterial community composition of each aquifer was significantly different (p < 0.05) between samples collected prior and after groundwater recharge. Furthermore, microbial source tracking results suggested a low contribution of surface environments to the groundwater microbiome except for in the months associated with recharge (March 2019 and April 2019). Overall, despite differences in soil permeability between both sites, the period of snow melt was followed by important changes in the composition of microbial communities from aquifers.

Phylogenetic group-based assembly and co-occurrence pattern of the microbial community in full-scale wastewater treatment plants during the Chinese spring festival

The quality and quantity of domestic sewage discharge vary significantly during the Chinese Spring Festival due to the huge population shift. The dynamics of microbial community traits during the Spring Festival, particularly the phylogenetic group-based assembly and co-occurrence patterns, are however little understood. Here, influent and activated sludge samples from 2 full-scale wastewater treatment plants were collected bi-daily throughout a 20-day Spring Festival period and subjected to high-throughput Illumina-MiSeq sequencing. The findings revealed that the microbial communities in the activated sludge displayed a comparatively stable pattern, and that the influent communities experienced significant temporal fluctuations in terms of diversity and composition. The characterization by "Infer Community Assembly Mechanisms by Phylogenetic-bin based null model" demonstrated that for Competibacter glycogen-accumulating organisms, the assembly mechanism shifted from deterministic process (HoS = 69.5%) before the Spring Festival to stochastic process (DR = 65.9%) after the Spring Festival. The network analysis revealed that the network structure of sludge communities was more stable before the Spring Festival than that after the Spring Festival. Additionally, sludge communities had no keystone species in common with the influent before the Spring Festival, while the sludge and influent communities shared two keystone taxa after the Spring Festival (Sebaldella and Candidatus Competibacter). This study would deepen our understanding of the microbial ecology in biological wastewater treatment systems, which also aids in managing wastewater treatment plants.

Meta-analysis addressing the characterization and risk identification of antibiotics and antibiotic resistance genes in global groundwater

Antimicrobial resistance (AMR) is one of the significant global issues to public health. Compared to other aquatic environments, research on AMR in groundwater is scarce. In the study, a meta-analysis was conducted to explore the characteristics and risks of antibiotics and antibiotic resistance genes (ARGs) in global groundwater, using a data set of antibiotic concentrations collected from publications during 2000-2021 and a large-scale metagenomes of groundwater samples (n = 330). The ecotoxicological risks of antibiotics in the global groundwater were evaluated using mixture risk quotient with concentration addition model to consider the synergistic effects of multiple antibiotics. Bioinformatic annotations identified 1413 ARGs belonging to 37 ARG types in the global groundwater, dominated by rifamycin, polyketide, and quinolone resistance genes and including some emerging ARGs such as mcr-family and carbapenem genes. Relatively, the level of ARGs in the groundwater from spring was significantly higher (ANOVA, p < 0.01) than those from the riparian zone, sand and deep aquifer. Similarly, metal resistance genes (MRGs) were prevalent in the global groundwater, and network analysis suggested the MRGs presented non-random co-occurrence with the ARGs in such environments. Taxonomic annotations showed Proteobacteria, Actinobacteria, Eukaryota, Acidobacteria and Thaumarchaeota were the dominant phylum in the groundwater, and the microbial community largely shaped profile of ARGs in the environment. Notably, the ARGs presented co-occurrence with mobile genetic elements, virulence factors and human bacterial pathogens, indicating potential dissemination risk of ARGs in the groundwater. Furthermore, an omics-based approach was used for health risk assessment of antibiotic resistome and screened out 152 risk ARGs in the global groundwater. Comparatively, spring and cold creek presented higher risk index, which deserves more attention to ensure the safety of water supply.

Subsurface microbial communities as a tool for characterizing regional-scale groundwater flow

Subsurface microbial community distribution patterns are influenced by biogeochemical and groundwater fluxes and may inform hydraulic connections along groundwater-flow paths. This study examined the regional-scale microbial community of the Death Valley Regional Flow System and evaluated whether subsurface communities can be used to identify groundwater-flow paths between recharge and discharge areas. Samples were collected from 36 sites in three groundwater basins: Pahute Mesa-Oasis Valley (PMOV), Ash Meadows (AM), and Alkali Flat-Furnace Creek Ranch (AFFCR). Microbial diversity within and between communities varied by location, and communities were separated into two overall groups that affiliated with the AM and PMOV/AFFCR basins. Network analysis revealed patterns between clusters of common microbes that represented groundwaters with similar geochemical conditions and largely corroborated hydraulic connections between recharge and discharge areas. Null model analyses identified deterministic and stochastic ecological processes contributing to microbial community assemblages. Most communities were more different than expected and governed by dispersal limitation, geochemical differences, or undominating processes. However, certain communities from sites located within or near the Nevada National Security Site were more similar than expected and dominated by homogeneous dispersal or selection. Overall, the (dis)similarities between the microbial communities of DVRFS recharge and discharge areas supported previously documented hydraulic connections between: (1) Spring Mountains and Ash Meadows; (2) Frenchman and Yucca Flat and Amargosa Desert; and (3) Amargosa Desert and Death Valley. However, only a portion of the flow path between Pahute Mesa and Oasis Valley could be supported by microbial community analyses, likely due to well-associated artifacts in samples from the two Oasis Valley sites. This study demonstrates the utility of combining microbial data with hydrologic, geologic, and water-chemistry information to comprehensively characterize groundwater systems, highlighting both strengths and limitations of this approach.

Microbiome characterization of alpine water springs for human consumption reveals site- and usage-specific microbial signatures

The microbiome of water springs is gaining increasing interest, especially in water intended for human consumption. However, the knowledge about large-scale patterns in water springs microbiome is still incomplete. The presence of bacteria in water sources used for human consumption is a major concern for health authorities; nonetheless, the standard microbiological quality checks are focused only on pathogenic species and total microbial load. Using 16S rRNA high throughput sequencing, we characterized the microbiome from 38 water springs in Trentino (Northern Italy) for 2 consecutive years in order to gain precious insights on the microbiome composition of these unexplored yet hardly exploited environments. The microbiological studies were integrated with standard measurements of physico-chemical parameters performed by the Provincial Office for Environmental Monitoring in order to highlight some of the dynamics influencing the microbial communities of these waters. We found that alpha diversity showed consistent patterns of variation overtime, and showed a strong positive correlation with the water nitrate concentration and negatively with fixed residue, electrical conductivity, and calcium concentration. Surprisingly, alpha diversity did not show any significant correlation with neither pH nor temperature. We found that despite their remarkable stability, different water springs display different coefficients of variation in alpha diversity, and that springs used for similar purposes showed similar microbiomes. Furthermore, the springs could be grouped according to the number of shared species into three major groups: low, mid, and high number of shared taxa, and those three groups of springs were consistent with the spring usage. Species belonging to the phyla Planctomycetes and Verrucomicrobia were prevalent and at relatively high abundance in springs classified as low number of shared species, whereas the phylum Lentisphaerae and the Candidate Phyla radiation were prevalent at higher abundance in the mineral and potable springs. The present study constitutes an example for standard water spring monitoring integrated with microbial community composition on a regional scale, and provides information which could be useful in the design and application of future water management policies in Trentino.

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

The differences in bacterial community assembly mechanism between surface water and groundwater, as well as the driving factors of environmental factors, are still unknown. Here we aimed to answer these questions by analyzing microbial community samples from surface water and groundwater. We observed a strong connection between microbial communities in surface water and groundwater and several human pathogens are shared between surface water and groundwater; however, the richness and diversity of groundwater microbial communities were greater than those of surface water, regardless of the season. Additionally, bacterial community compositions of surface water and groundwater differed significantly between seasons. Most importantly, the groundwater community exhibited a highly deterministic assembly process (56% contributed by deterministic process, with neutral community model R² = 0.277) compared with surface water (51% contributed by deterministic process, with R² = 0.526). This study provides a deep understanding of the effects of environmental factors on surface water and groundwater microbial communities, to better protect water resources.

Biocarriers facilitated gravity-driven membrane filtration of domestic wastewater in cold climate: Combined effect of temperature and periodic cleaning

In this study, two lava stone biocarrier facilitated gravity-driven membrane (GDM) reactors were operated at ~8 °C and ~22 °C in parallel for treating primary wastewater effluent. Although the biocarrier reactor at 8 °C displayed less efficient removals of biodegradable organics than that at 22 °C, both GDM systems (without cleaning) showed comparable fouling resistance distribution patterns, accompanying with similar cake filtration constants and pore constriction constants by modelling simulation. Compared to the GDM at 8 °C, more foulants were accumulated on the GDM at 22 °C, but they presented similar soluble organics/inorganics contents and specific cake resistances. This indicated the cake layers at 22 °C may contain greater-sized foulants due to proliferation of both prokaryotes and eukaryotes, leading to a relatively less-porous nature. In the presence of periodic cleaning (at 50 °C), the cleaning effectiveness followed a sequence as ultrasonication-enhanced physical cleaning > two-phase flow cleaning > chemical-enhanced physical cleaning > physical cleaning, regardless of GDM operation temperature. However, significantly higher cake resistances were observed in the GDM system at 22 °C than those at 8 °C, because shear force tended to remove loosely-attached foulant layers and may compress the residual dense cake layer. The presence of periodic cleaning led to dissimilar dominant prokaryotic and eukaryotic communities in the cake layers as those without cleaning and in the lava stone biocarriers. Nevertheless, operation temperature did not influence GDM permeate quality, which met EU discharge standards.

Bacterial Necromass Is Rapidly Metabolized by Heterotrophic Bacteria and Supports Multiple Trophic Levels of the Groundwater Microbiome

Pristine groundwater is a highly stable environment with microbes adapted to dark, oligotrophic conditions. Input events like heavy rainfalls can introduce the excess particulate organic matter, including surface-derived microorganisms, thereby disturbing the groundwater microbiome. Some surface-derived bacteria will not survive this translocation, leading to an input of necromass to the groundwater. Here, we investigated the effects of necromass addition to the microbial community in fractured bedrock groundwater, using groundwater mesocosms as model systems. We followed the uptake of ¹³C-labeled necromass by the bacterial and eukaryotic groundwater community quantitatively and over time using a complementary protein-stable and DNA-stable isotope probing approach. Necromass was rapidly depleted in the mesocosms within 4 days, accompanied by a strong decrease in Shannon diversity and a 10-fold increase in bacterial 16S rRNA gene copy numbers. Species of Flavobacterium, Massilia, Rheinheimera, Rhodoferax, and Undibacterium dominated the microbial community within 2 days and were identified as key players in necromass degradation, based on a ¹³C incorporation of >90% in their peptides. Their proteomes comprised various proteins for uptake and transport functions and amino acid metabolization. After 4 and 8 days, the autotrophic and mixotrophic taxa Nitrosomonas, Limnohabitans, Paucibacter, and Acidovorax increased in abundance with a ¹³C incorporation between 0.5% and 23%. Likewise, eukaryotes assimilated necromass-derived carbon either directly or indirectly. Our data point toward a fast and exclusive uptake of labeled necromass by a few specialists followed by a concerted action of groundwater microorganisms, including autotrophs presumably fueled by released, reduced nitrogen and sulfur compounds generated during necromass degradation.

IMPORTANCE Subsurface microbiomes provide essential ecosystem services, like the generation of drinking water. These ecosystems are devoid of light-driven primary production, and microbial life is adapted to the resulting oligotrophic conditions. Modern groundwater is most vulnerable to anthropogenic and climatic impacts. Heavy rainfalls, which will increase with climate change, can result in high surface inputs into shallow aquifers by percolation or lateral flow. These inputs include terrestrial organic matter and surface-derived microbes that are not all capable to flourish in aquatic subsurface habitats. Here, we investigated the response of groundwater mesocosms to the addition of bacterial necromass, simulating event-driven surface input. We found that the groundwater microbiome responds with a rapid bloom of only a few primary degraders, followed by the activation of typical groundwater autotrophs and mixotrophs, as well as eukaryotes. Our results suggest that this multiphase strategy is essential to maintain the balance of the groundwater microbiome to provide ecosystem services.

Linkage and driving mechanisms of antibiotic resistome in surface and ground water: Their responses to land use and seasonal variation

Antibiotic resistance in the environment, mostly mediated by antibiotic resistance genes (ARGs), has posed a threat to ecological and human health. Contamination of surface water and groundwater with ARGs has become a serious environmental concern. However, the distinctions and similarities across ARG profiles, the various ecological processes associated with ARGs, the driving mechanisms for ARG profiles in surface water and groundwater, and how they respond to land use and seasonal variation remain unknown. To tackle these issues, the contamination of ARGs in surface water and groundwater in central China was investigated using metagenomic technology. The results indicated that seasonal changes in ARG abundance and diversity were inconsistent across surface water and groundwater, and that the relationship between ARGs in surface water and groundwater was stronger during the rainy season. Land use had a greater effect on ARGs in surface water than in groundwater and was stronger during the dry season than during the rainy season. More interestingly, the ideal buffer zones with the greatest impact of land use on the ARGs of surface water and groundwater had distinct radii: 1500 m for both dry and rainy seasons in surface water, and 1000 m for dry season and 500 m for rainy season in groundwater. Furthermore, stochastic mechanisms mediated by mobile gene elements (MGEs) contribute significantly more to ARG assemblages than deterministic processes, particularly in groundwater. Furthermore, our results also showed that ARG enrichment in microbial communities was host- dependent, and the risk of ARGs in groundwater was greater both during the rainy season and dry season. In conclusion, the findings have improved our understanding of the relationship and driving mechanisms of ARGs in surface and ground water, as well as their responses to land use and seasonal variation, which may be beneficial for limiting ARG pollution in a watershed with high levels of anthropogenic activity.

Decreased levels and ecological risks of disinfection by-product chloroform in a field-scale artificial groundwater recharge project by colloid supplement

To bolster freshwater supply, artificial groundwater recharge with recycled water has increasingly attracted research attentions and interests. However, artificial groundwater recharge has potential risks to groundwater quality, as recharge water disinfection is frequently used for pathogen inactivation and causes the concerns of disinfection by-products (DBPs). Colloid supplement is a good approach solving this problem, but its roles in mitigating DBPs remains unclear. In this study, we collected 20 groundwater and soil samples from a field-scale groundwater recharge project, and explored the impacts of silica colloids on chloroform migration and groundwater bacterial communities during the recharge process. Water physicochemical variables changed along the recharge time, and colloid supplement significantly reduced chloroform formation and slowed its migration in groundwater. Bacterial communities in groundwater, river water and recharge water were significantly different. Gammaproteobacteria in recharge water (71.7%) was more abundant than in river water (30.5%) and groundwater (33.5%), while Actinobacteria dominated groundwater (40.6%). After recharge, Gammaproteobacteria increased more with colloid supplement (75.7%) than without (52.6%), attributing to its dominance in soils (74.6%). Our results suggested more bacterial lineages released from soils into aquifer by silica colloid supplement, owing to the competitive adsorption encouraging microbial transfer, especially Gram-negative bacteria. Our findings unraveled the effects of colloid supplement on chloroform formation and migration during artificial groundwater recharge, which consequently altered groundwater bacterial communities, and offered valuable suggestions for the safety management of DBPs in aquifer recharge.

The economical lifestyle of CPR bacteria in groundwater allows little preference for environmental drivers

BACKGROUND

The highly diverse Cand. Patescibacteria are predicted to have minimal biosynthetic and metabolic pathways, which hinders understanding of how their populations differentiate in response to environmental drivers or host organisms. Their mechanisms employed to cope with oxidative stress are largely unknown. Here, we utilized genome-resolved metagenomics to investigate the adaptive genome repertoire of Patescibacteria in oxic and anoxic groundwaters, and to infer putative host ranges.

RESULTS

Within six groundwater wells, Cand. Patescibacteria was the most dominant (up to 79%) super-phylum across 32 metagenomes sequenced from DNA retained on 0.2 and 0.1 µm filters after sequential filtration. Of the reconstructed 1275 metagenome-assembled genomes (MAGs), 291 high-quality MAGs were classified as Cand. Patescibacteria. Cand. Paceibacteria and Cand. Microgenomates were enriched exclusively in the 0.1 µm fractions, whereas candidate division ABY1 and Cand. Gracilibacteria were enriched in the 0.2 µm fractions. On average, Patescibacteria enriched in the smaller 0.1 µm filter fractions had 22% smaller genomes, 13.4% lower replication measures, higher proportion of rod-shape determining proteins, and of genomic features suggesting type IV pili mediated cell-cell attachments. Near-surface wells harbored Patescibacteria with higher replication rates than anoxic downstream wells characterized by longer water residence time. Except prevalence of superoxide dismutase genes in Patescibacteria MAGs enriched in oxic groundwaters (83%), no major metabolic or phylogenetic differences were observed. The most abundant Patescibacteria MAG in oxic groundwater encoded a nitrate transporter, nitrite reductase, and F-type ATPase, suggesting an alternative energy conservation mechanism. Patescibacteria consistently co-occurred with one another or with members of phyla Nanoarchaeota, Bacteroidota, Nitrospirota, and Omnitrophota. Among the MAGs enriched in 0.2 µm fractions,, only 8% Patescibacteria showed highly significant one-to-one correlation, mostly with Omnitrophota. Motility and transport related genes in certain Patescibacteria were highly similar to genes from other phyla (Omnitrophota, Proteobacteria and Nanoarchaeota).

CONCLUSION

Other than genes to cope with oxidative stress, we found little genomic evidence for niche adaptation of Patescibacteria to oxic or anoxic groundwaters. Given that we could detect specific host preference only for a few MAGs, we speculate that the majority of Patescibacteria is able to attach multiple hosts just long enough to loot or exchange supplies.