Local and Regional Diversity Reveals Dispersal Limitation and Drift as Drivers for Groundwater Bacterial Communities from a Fractured Granite Formation

Spatial and temporal groundwater biogeochemical variability help inform subsurface connectivity within a high-altitude Alpine catchment (Riale di Ronco, Switzerland)

Accessing the deep terrestrial subsurface (greater than 1 km below the surface) presents significant practical challenges, leaving these ecosystems largely uncharacterized despite their extensive presence beneath Earth's landmasses. In this study, we introduce the BedrettoLab Deep Life Observatory (DELOS), a new underground laboratory to study the biogeochemical diversity of groundwater in a high-altitude Alpine catchment tens of meters to 1.6 km underground. Biogeochemical monitoring of DELOS over spatial and temporal scales highlight three dominant ecotypes throughout DELOS: (1) Shallow groundwater with low electrical conductivity enriched in Leptospirillia; (2) High-inflow fault zones enriched in ultra-small bacteria and archaea; (3) Bicarbonate-enriched waters that are enriched in Candidatus Kryptonia and Spirochaetota. Despite a consistent lithology throughout DELOS, groundwater from fractures that are spatially near each other are not always represented by the same ecotype and can be more similar to groundwater emitted from fractures thousands of meters away. Despite this heterogeneity, the biological and hydrochemical compositions of the groundwater of individual fractures remained relatively stable throughout the course of a 1-year monitoring period. An exception to this trend occurred after a series of seismic events near one groundwater-bearing fracture. Here, the microbial community and hydrochemical composition of the groundwater changed after the seismic events but returned to the site's "baseline" composition within 3 weeks. Taken together, these findings provide new insights into the spatial and temporal heterogeneity of deep subsurface ecosystems and the subsurface connectivity of an Alpine subsurface environment.

Eukaryotic plankton community and assembly processes in a large-scale water diversion project in China

The Middle Route of the South to North Water Diversion Project (MRP) and its water source, the Danjiangkou Reservoir (DJK), play a pivotal role in mitigating the chronic water scarcity challenges faced by northern China. Eukaryotic plankton are widespread in aquatic ecosystems, which are crucial for the water quality stability of DJK and MRP, yet comparative studies on their contemporaneous dynamics and assembly processes are scarce. In this study, amplicon sequencing was used to investigate the eukaryotic plankton communities. The results revealed that the similarity in community composition of DJK is significantly higher than that of MRP, exhibiting distance-decay patterns. Environmental heterogeneity exhibits significant differences between DJK and MRP, and it significantly influences community composition and alpha diversity. Additionally, the assembly processes of eukaryotic plankton in both DJK and MRP are predominantly influenced by stochastic processes. However, in comparison to DJK, deterministic processes have a more pronounced impact on MRP, accounting for 39.29% and 1.82%, respectively. The variations in total nitrogen (TN), chlorophy IIa (Chl.a), and conductivity (Spc) have led to a transition in the assembly of eukaryotic phytoplankton communities in MRP from a stochastic process to a deterministic process. This study extends insights into the dynamics and assembly processes of eukaryotic plankton communities in the large, engineered drinking water diversion project and its water source, which is also useful for the management and regulation of the DJK and MRP.

Metabolic adaptations underpin high productivity rates in relict subsurface water

Groundwater aquifers are ecological hotspots with diverse microbes essential for biogeochemical cycles. Their ecophysiology has seldom been studied on a basin scale. In particular, our knowledge of chemosynthesis in the deep aquifers where temperatures reach 60 °C, is limited. Here, we investigated the diversity, activity, and metabolic potential of microbial communities from nine wells reaching ancient groundwater beneath Israel's Negev Desert, spanning two significant, deep (up to 1.5 km) aquifers, the Judea Group carbonate and Kurnub Group Nubian sandstone that contain fresh to brackish, hypoxic to anoxic water. We estimated chemosynthetic productivity rates ranging from 0.55 ± 0.06 to 0.82 ± 0.07 µg C L-1 d-1 (mean ± SD), suggesting that aquifer productivity may be underestimated. We showed that 60% of MAGs harbored genes for autotrophic pathways, mainly the Calvin-Benson-Bassham cycle and the Wood-Ljungdahl pathway, indicating a substantial chemosynthetic capacity within these microbial communities. We emphasize the potential metabolic versatility in the deep subsurface, enabling efficient carbon and energy use. This study set a precedent for global aquifer exploration, like the Nubian Sandstone Aquifer System in the Arabian and Western Deserts, and reconsiders their role as carbon sinks.

Environmental stress mediates groundwater microbial community assembly

Community assembly describes how different ecological processes shape microbial community composition and structure. How environmental factors impact community assembly remains elusive. Here we sampled microbial communities and >200 biogeochemical variables in groundwater at the Oak Ridge Field Research Center, a former nuclear waste disposal site, and developed a theoretical framework to conceptualize the relationships between community assembly processes and environmental stresses. We found that stochastic assembly processes were critical (>60% on average) in shaping community structure, but their relative importance decreased as stress increased. Dispersal limitation and 'drift' related to random birth and death had negative correlations with stresses, whereas the selection processes leading to dissimilar communities increased with stresses, primarily related to pH, cobalt and molybdenum. Assembly mechanisms also varied greatly among different phylogenetic groups. Our findings highlight the importance of microbial dispersal limitation and environmental heterogeneity in ecosystem restoration and management.

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.

Insights into microbial diversity on plastisphere by multi-omics

Plastic pollution is a major concern in marine environment as it takes many years to degrade and is one of the greatest threats to marine life. Plastic surface, referred to as plastisphere, provides habitat for growth and proliferation of various microorganisms. The discovery of these microbes is necessary to identify significant genes, enzymes and bioactive compounds that could help in bioremediation and other commercial applications. Conventional culture techniques have been successful in identifying few microbes from these habitats, leaving majority of them yet to be explored. As such, to recognize the vivid genetic diversity of microbes residing in plastisphere, their structure and corresponding ecological roles within the ecosystem, an emerging technique, called metagenomics has been explored. The technique is expected to provide hitherto unknown information on microbes from the plastisphere. Metagenomics along with next generation sequencing provides comprehensive knowledge on microbes residing in plastisphere that identifies novel microbes for plastic bioremediation, bioactive compounds and other potential benefits. The following review summarizes the efficiency of metagenomics and next generation sequencing technology over conventionally used methods for culturing microbes. It attempts to illustrate the workflow mechanism of metagenomics to elucidate diverse microbial profiles. Further, importance of integrated multi-omics techniques has been highlighted in discovering microbial ecology residing on plastisphere for wider applications.

Seasonal variation significantly affected bacterioplankton and eukaryoplankton community composition in Xijiang River, China

Both bacterioplankton and eukaryoplankton communities play important roles in the geochemical cycles and energy flows of river ecosystems. However, whether a seasonal change in bacterioplankton and eukaryoplankton communities is synchronous remains unclear. To test the synchronicity and analyze how physical and chemical environmental factors affect these communities, we compared bacterioplankton and eukaryoplankton communities in surface water samples between March (dry season) and June (rainfall season) considering water environmental factors. Our results showed that there was no significant difference in operational taxonomic unit number, Shannon index, and Chao1 index in bacterioplankton and eukaryoplankton communities between March and June. However, principal component analysis showed that the communities were significantly different between the sampling times and sampling sites. Water temperature (WT), oxidation-reduction potential (ORP), water transparency (SD), NO₃-N, and NH₃ significantly influenced bacterioplankton communities, and WT, SD, ORP, and NH₄-N significantly influenced eukaryoplankton communities in the river. These results implied that compared with the sampling sites, sampling times more significantly affected the bacterioplankton and eukaryoplankton river communities by influencing WT, ORP, SD, and nitrogen forms.

Microbial Communities in a Serpentinizing Aquifer Are Assembled through Strong Concurrent Dispersal Limitation and Selection

In recent years, our appreciation of the extent of habitable environments in Earth’s subsurface has greatly expanded, as has our understanding of the biodiversity contained within. Most studies have relied on single sampling points, rather than considering the long-term dynamics of subsurface environments and their microbial populations. One such habitat are aquifers associated with the aqueous alteration of ultramafic rocks through a process known as serpentinization. Ecological modeling performed on a multiyear time series of microbiology, hydrology, and geochemistry in an ultrabasic aquifer within the Coast Range Ophiolite reveals that community assembly is governed by undominated assembly (i.e., neither stochastic [random] nor deterministic [selective] processes alone govern assembly). Controls on community assembly were further assessed by characterizing aquifer hydrogeology and microbial community adaptations to the environment. These analyses show that low permeability rocks in the aquifer restrict the transmission of microbial populations between closely situated wells. Alpha and beta diversity measures and metagenomic and metatranscriptomic data from microbial communities indicate that high pH and low dissolved inorganic carbon levels impose strong environmental selection on microbial communities within individual wells. Here, we find that the interaction between strong selection imposed by extreme pH and enhanced ecological drift due to dispersal limitation imposed by slow fluid flow results in the undominated assembly signal observed throughout the site. Strong environmental selection paired with extremely low dispersal in the subsurface results in low diversity microbial communities that are well adapted to extreme pH conditions and subject to enhanced stochasticity introduced by ecological drift over time.

IMPORTANCE Microbial communities existing under extreme or stressful conditions have long been thought to be structured primarily by deterministic processes. The application of macroecology theory and modeling to microbial communities in recent years has spurred assessment of assembly processes in microbial communities, revealing that both stochastic and deterministic processes are at play to different extents within natural environments. We show that low diversity microbial communities in a hard-rock serpentinizing aquifer are assembled under the influence of strong selective processes imposed by high pH and enhanced ecological drift that occurs as the result of dispersal limitation due to the slow movement of water in the low permeability aquifer. This study demonstrates the important roles that both selection and dispersal limitation play in terrestrial serpentinites, where extreme pH assembles a microbial metacommunity well adapted to alkaline conditions and dispersal limitation drives compositional differences in microbial community composition between local communities in the subsurface.

Direct quantification of ecological drift at the population level in synthetic bacterial communities

In community ecology, drift refers to random births and deaths in a population. In microbial ecology, drift is estimated indirectly via community snapshots but in this way, it is almost impossible to distinguish the effect of drift from the effect of other ecological processes. Controlled experiments where drift is quantified in isolation from other processes are still missing. Here we isolate and quantify drift in a series of controlled experiments on simplified and tractable bacterial communities. We detect drift arising randomly in the populations within the communities and resulting in a 1.4–2% increase in their growth rate variability on average. We further use our experimental findings to simulate complex microbial communities under various conditions of selection and dispersal. We find that the importance of drift increases under high selection and low dispersal, where it can lead to ~5% of species loss and to ~15% increase in β-diversity. The species extinct by drift are mainly rare, but they become increasingly less rare when selection increases, and dispersal decreases. Our results provide quantitative insights regarding the properties of drift in bacterial communities and suggest that it accounts for a consistent fraction of the observed stochasticity in natural surveys.

Selection imposed by local environmental conditions drives differences in microbial community composition across geographically distinct groundwater aquifers

Several studies have analyzed biogeographic distribution patterns of microbial communities across broad spatial scales. However, it is often unclear to what extent differences in community composition across different regions are caused by dispersal limitation or selection, and if selection is caused by local environmental conditions alone or additional broad-scale region-specific factors. This is especially true for groundwater environments, which have been understudied in this context relative to other non-subsurface habitats. Here, we analyzed microbial community composition based on exact 16S rRNA amplicon sequence variants (ASVs) from four geographically separated aquifers located in different regions along a latitudinal transect of ∼700 km across Germany. Using a combination of variation partitioning and ecological null models revealed that differences in microbial community composition were mainly the product of selection imposed by local environmental conditions and to a smaller but still significant extent dispersal limitation and drift across regions. Only ∼23% of the total variation in microbial community composition remained unexplained, possibly due to underestimated effects of dispersal limitation among local communities within regions and temporal drift. No evidence was found for selection due to region-specific factors independent of local environmental conditions.

Comamonadaceae OTU as a Remnant of an Ancient Microbial Community in Sulfidic Waters

Intraterrestrial waters harbor microbial communities being extensively studied to understand microbial processes underlying subsurface ecosystem functioning. This paper provides the results of an investigation on the microbiomes of unique, subsurface sulfidic waters associated with Upper Jurassic, Cretaceous, and Miocene sediments. We used high-throughput 16S rDNA amplicon sequencing to reveal the structure of bacterial and archaeal communities in water samples differing in sulfide content (20-960 mg/dm3), salinity (1.3-3.2%), and depth of extraction (60-660 m below ground level). Composition of the bacterial communities strongly varied across the samples; however, the bacteria participating in the sulfur cycle were common in all sulfidic waters. The shallowest borehole water (60 m bgl) was dominated by sulfur-oxidizing Epsilonproteobacteria (Sulfurimonas, Sulfurovum). In the waters collected from greater depths (148-300 m bgl), the prevalence of Betaproteobacteria (Comamonadaceae) and sulfate/sulfur-reducing Deltaproteobacteria (Desulfopila, Desulfomicrobium, MSBL7) was observed. Sulfate reducers (members of Clostridia: Candidatus Desulforudis) were the most abundant bacteria in the deepest borehole water (660 m bgl). Out of 850 bacterial OTUs, only one, affiliated with the Comamonadaceae family, was found abundant (> 1% of total bacterial sequences) in all samples. Contribution of Archaea to the whole microbial communities was lower than 0.5%. Archaeal communities did not differ across the samples and they consisted of Halobacteriaceae. Out of 372 archaeal OTUs, five, belonging to the four genera Natronomonas, Halorubrum, Halobellus, and Halorhabdus, were the most numerous.

Changes in bacterial and archaeal communities during the concentration of brine at the graduation towers in Ciechocinek spa (Poland)

This study evaluates the changes in bacterial and archaeal community structure during the gradual evaporation of water from the brine (extracted from subsurface Jurassic deposits) in the system of graduation towers located in Ciechocinek spa, Poland. The communities were assessed with 16S rRNA gene sequencing (MiSeq, Illumina) and microscopic methods. The microbial cell density determined by direct cell count was at the order of magnitude of 10⁷ cells/mL. It was found that increasing salt concentration was positively correlated with both the cell counts, and species-level diversity of bacterial and archaeal communities. The archaeal community was mostly constituted by members of the phylum Euryarchaeota, class Halobacteria and was dominated by Halorubrum-related sequences. The bacterial community was more diverse, with representatives of the phyla Proteobacteria and Bacteroidetes as the most abundant. The proportion of Proteobacteria decreased with increasing salt concentration, while the proportion of Bacteroidetes increased significantly in the more concentrated samples. Representatives of the genera Idiomarina, Psychroflexus, Roseovarius, and Marinobacter appeared to be tolerant to changes of salinity. During the brine concentration, the relative abundances of Sphingobium and Sphingomonas were significantly decreased and the raised contributions of genera Fabibacter and Fodinibius were observed. The high proportion of novel (not identified at 97% similarity level) bacterial reads (up to 42%) in the 16S rRNA gene sequences indicated that potentially new bacterial taxa inhabit this unique environment.

Spatial Autocorrelation, Source Water and the Distribution of Total and Viable Microbial Abundances within a Crystalline Formation to a Depth of 800 m

Proposed radioactive waste repositories require long residence times within deep geological settings for which we have little knowledge of local or regional subsurface dynamics that could affect the transport of hazardous species over the period of radioactive decay. Given the role of microbial processes on element speciation and transport, knowledge and understanding of local microbial ecology within geological formations being considered as host formations can aid predictions for long term safety. In this relatively unexplored environment, sampling opportunities are few and opportunistic. We combined the data collected for geochemistry and microbial abundances from multiple sampling opportunities from within a proposed host formation and performed multivariate mixing and mass balance (M3) modeling, spatial analysis and generalized linear modeling to address whether recharge can explain how subsurface communities assemble within fracture water obtained from multiple saturated fractures accessed by boreholes drilled into the crystalline formation underlying the Chalk River Laboratories site (Deep River, ON, Canada). We found that three possible source waters, each of meteoric origin, explained 97% of the samples, these are: modern recharge, recharge from the period of the Laurentide ice sheet retreat (ca. ∼12000 years before present) and a putative saline source assigned as Champlain Sea (also ca. 12000 years before present). The distributed microbial abundances and geochemistry provide a conceptual model of two distinct regions within the subsurface associated with bicarbonate – used as a proxy for modern recharge – and manganese; these regions occur at depths relevant to a proposed repository within the formation. At the scale of sampling, the associated spatial autocorrelation means that abundances linked with geochemistry were not unambiguously discerned, although fine scale Moran’s eigenvector map (MEM) coefficients were correlated with the abundance data and suggest the action of localized processes possibly associated with the manganese and sulfate content of the fracture water.