Archaeal Diversity and CO2 Fixers in Carbonate-/Siliciclastic-Rock Groundwater Ecosystems

Microbial degradation and assimilation of veratric acid in oxic and anoxic groundwaters

Microbial communities are key players in groundwater ecosystems. In this dark environment, heterotrophic microbes rely on biomass produced by the activity of lithoautotrophs or on the degradation of organic matter seeping from the surface. Most studies on bacterial diversity in groundwater habitats are based on 16S gene sequencing and full genome reconstructions showing potential metabolic pathways used in these habitats. However, molecular-based studies do not allow for the assessment of population dynamics over time or the assimilation of specific compounds and their biochemical transformation by microbial communities. Therefore, in this study, we combined DNA-, phospholipid fatty acid-, and metabolomic-stable isotope probing to target and identify heterotrophic bacteria in the groundwater setting of the Hainich Critical Zone Exploratory (CZE), focusing on 2 aquifers with different physico-chemical conditions (oxic and anoxic). We incubated groundwater from 4 different wells using either 13C-labeled veratric acid (a lignin-derived compound) (single labeling) or a combination of 13CO2 and D-labeled veratric acid (dual labeling). Our results show that heterotrophic activities dominate all groundwater sites. We identified bacteria with the potential to break down veratric acid (Sphingobium or Microbacterium). We observed differences in heterotrophic activities between the oxic and anoxic aquifers, indicating local adaptations of bacterial populations. The dual labeling experiments suggested that the serine pathway is an important carbon assimilation pathway and that organic matter was an important source of hydrogen in the newly produced lipids. These experiments also yielded different labeled taxa compared to the single labeling experiments, showing that there exists a complex interaction network in the groundwater habitats.

Deep Isolated Aquifer Brines Harbor Atypical Halophilic Microbial Communities in Quebec, Canada

The deep terrestrial subsurface, hundreds of meters to kilometers below the surface, is characterized by oligotrophic conditions, dark and often anoxic settings, with fluctuating pH, salinity, and water availability. Despite this, microbial populations are detected and active, contributing to biogeochemical cycles over geological time. Because it is extremely difficult to access the deep biosphere, little is known about the identity and metabolisms of these communities, although they likely possess unknown pathways and might interfere with deep waste deposits. Therefore, we analyzed rock and groundwater microbial communities from deep, isolated brine aquifers in two regions dating back to the Ordovician and Devonian, using amplicon and whole genome sequencing. We observed significant differences in diversity and community structure between both regions, suggesting an impact of site age and composition. The deep hypersaline groundwater did not contain typical halophilic bacteria, and genomes suggested pathways involved in protein and hydrocarbon degradation, and carbon fixation. We identified mainly one strategy to cope with osmotic stress: compatible solute uptake and biosynthesis. Finally, we detected many bacteriophage families, potentially indicating that bacteria are infected. However, we also found auxiliary metabolic genes in the viral genomes, probably conferring an advantage to the infected hosts.

Response of Prokaryotic Communities to Freshwater Salinization

Each year, millions of tons of sodium chloride are dumped on roads, contributing to the salinization of freshwater environments. Thus, we sought to understand the effect of sodium chloride (NaCl) on freshwater lake prokaryotic communities, an important and understudied component of food webs. Using mesocosms with 0.01–2.74 ppt NaCl (0.27–1110.86 mg/L Cl−), we evaluated the effect generated on the diversity and absolute abundance of prokaryotic populations after three and six weeks. A positive relationship between Cl− values and absolute bacterial abundance was found after three weeks. The influence of eukaryotic diversity variation was observed as well. Significant differentiation of bacterial communities starting at 420 mg/L Cl− was observed after three weeks, levels lower than the Canadian and US recommendations for acute chloride exposure. The partial recovery of a “pre-disturbance” community was observed following a drop in salinity at the threshold level of 420 mg/L Cl−. A gradual transition of dominance from Betaproteobacteria and Actinobacteria to Bacteroidia and Alphaproteobacteria was observed and is overall similar to the natural transition observed in estuaries.

Diverse ecophysiological adaptations of subsurface Thaumarchaeota in floodplain sediments revealed through genome-resolved metagenomics

The terrestrial subsurface microbiome contains vastly underexplored phylogenetic diversity and metabolic novelty, with critical implications for global biogeochemical cycling. Among the key microbial inhabitants of subsurface soils and sediments are Thaumarchaeota, an archaeal phylum that encompasses ammonia-oxidizing archaea (AOA) as well as non-ammonia-oxidizing basal lineages. Thaumarchaeal ecology in terrestrial systems has been extensively characterized, particularly in the case of AOA. However, there is little knowledge on the diversity and ecophysiology of Thaumarchaeota in deeper soils, as most lineages, particularly basal groups, remain uncultivated and underexplored. Here we use genome-resolved metagenomics to examine the phylogenetic and metabolic diversity of Thaumarchaeota along a 234 cm depth profile of hydrologically variable riparian floodplain sediments in the Wind River Basin near Riverton, Wyoming. Phylogenomic analysis of the metagenome-assembled genomes (MAGs) indicates a shift in AOA population structure from the dominance of the terrestrial Nitrososphaerales lineage in the well-drained top ~100 cm of the profile to the typically marine Nitrosopumilales in deeper, moister, more energy-limited sediment layers. We also describe two deeply rooting non-AOA MAGs with numerous unexpected metabolic features, including the reductive acetyl-CoA (Wood-Ljungdahl) pathway, tetrathionate respiration, a form III RuBisCO, and the potential for extracellular electron transfer. These MAGs also harbor tungsten-containing aldehyde:ferredoxin oxidoreductase, group 4f [NiFe]-hydrogenases and a canonical heme catalase, typically not found in Thaumarchaeota. Our results suggest that hydrological variables, particularly proximity to the water table, impart a strong control on the ecophysiology of Thaumarchaeota in alluvial sediments.

From Surface to Subsurface: Diversity, Composition, and Abundance of Sessile and Endolithic Bacterial, Archaeal, and Eukaryotic Communities in Sand, Clay and Rock Substrates in the Laurentians (Quebec, Canada)

Microbial communities play an important role in shallow terrestrial subsurface ecosystems. Most studies of this habitat have focused on planktonic communities that are found in the groundwater of aquifer systems and only target specific microbial groups. Therefore, a systematic understanding of the processes that govern the assembly of endolithic and sessile communities is still missing. This study aims to understand the effect of depth and biotic factors on these communities, to better unravel their origins and to compare their composition with the communities detected in groundwater. To do so, we collected samples from two profiles (~0–50 m) in aquifer sites in the Laurentians (Quebec, Canada), performed DNA extractions and Illumina sequencing. The results suggest that changes in geological material characteristics with depth represent a strong ecological and phylogenetical filter for most archaeal and bacterial communities. Additionally, the vertical movement of water from the surface plays a major role in shallow subsurface microbial assembly processes. Furthermore, biotic interactions between bacteria and eukaryotes were mostly positive which may indicate cooperative or mutualistic potential associations, such as cross-feeding and/or syntrophic relationships in the terrestrial subsurface. Our results also point toward the importance of sampling both the geological formation and groundwater when it comes to studying its overall microbiology.

Influence of Environmental Factors on the Variability of Archaeal Communities in a Karst Wetland

Archaea are ubiquitous and play an important role in elemental cycles in Earth’s biosphere; but little is known about their diversity, distribution, abundance, and impact in karst environments. The present study investigated the effect of environmental factors on the variability of archaeal communities in the sediment of the Huixian karst wetland, the largest karst wetland in South China. Sediment cores were obtained from four sampling sites with different water depths and macrophyte inhabitants in both the winter of 2016 and the summer of 2018. The community analysis was based on PacBio sequencing and quantitative PCR of the archaeal 16S rRNA gene. The results showed that Euryarchaeota (57.4%) and Bathyarchaeota (38.7%) were dominant in all the samples. Methanogenic Methanosarcinales (25.1%) and Methanomicrobiales (13.7%), and methanotrophic archaea ANME-2d (9.0%) were the dominant Euryarchaeota; MCG-11 (16.5%), MCG-6 (9.1%), and MCG-5b (5.5%) were the dominant Bathyarchaeota. The community composition remained stable between summer and winter, and the vertical distributions of the archaeal phyla conformed to two patterns among the four sampling sites. In the winter samples, the archaeal 16S rRNA gene abundance was approximately 1.0E+10 copies/g of wet sediment and the Shannon index was 7.3±5, which were significantly higher than in the summer samples and in other karst environments. A correlation analysis showed that the moisture content and pH were the factors that mostly affected the archaeal communities. The prevalence of nitrate in the summer may be a key factor causing a significant decrease in archaeal abundance and diversity. Two features specific to karst environments, calcium-richness and weak alkalescence of the water supplies, may benefit the prevalence of bathyarchaeotal subgroups MCG-11, MCG-5b, and MCG-6. These results suggest that in karst wetlands, most of the archaea belong to clades that have significant roles in carbon turnover; their composition remains stable, but their abundance and diversity vary significantly from season to season.

Side Effects of Pesticides and Metabolites in Groundwater: Impact on Denitrification

The impact of two pesticides (S-metolachlor and propiconazole) and their respective main metabolites (ESA-metolachlor and 1,2,4-triazole) on bacterial denitrification in groundwater was studied. For this, the denitrification activity and the bacterial diversity of a microbial community sampled from a nitrate-contaminated groundwater were monitored during 20 days in lab experiments in the presence or absence of pesticides or metabolites at 2 or 10 μg/L. The kinetics of nitrate reduction along with nitrite and N₂O production all suggested that S-metolachlor had no or only little impact, whereas its metabolite ESA-metolachlor inhibited denitrification by 65% at 10 μg/L. Propiconazole and 1,2,4-triazole also inhibited denitrification at both concentrations, but to a lesser extent (29–38%) than ESA-metolachlor. When inhibition occurred, pesticides affected the reduction of nitrate into nitrite step. However, no significant differences were detected on the abundance of nitrate reductase narG and napA genes, suggesting an impact of pesticides/metabolites at the protein level rather than on denitrifying bacteria abundance. 16S rRNA gene Illumina sequencing indicated no major modification of bacterial diversity in the presence or absence of pesticides/metabolites, except for ESA-metolachlor and propiconazole at 10 μg/L that tended to increase or decrease Shannon and InvSimpson indices, respectively. General growth parameters suggested no impact of pesticides, except for propiconazole at 10 μg/L that partially inhibited acetate uptake and induced a decrease in microbial biomass. In conclusion, pesticides and metabolites can have side effects at environmental concentrations on microbial denitrification in groundwater and may thus affect ecosystem services based on microbial activities.

The microbial dimension of submarine groundwater discharge: current challenges and future directions

Despite the relevance of submarine groundwater discharge (SGD) for ocean biogeochemistry, the microbial dimension of SGD remains poorly understood. SGD can influence marine microbial communities through supplying chemical compounds and microorganisms, and in turn, microbes at the land–ocean transition zone determine the chemistry of the groundwater reaching the ocean. However, compared with inland groundwater, little is known about microbial communities in coastal aquifers. Here, we review the state of the art of the microbial dimension of SGD, with emphasis on prokaryotes, and identify current challenges and future directions. Main challenges include improving the diversity description of groundwater microbiota, characterized by ultrasmall, inactive and novel taxa, and by high ratios of sediment-attached versus free-living cells. Studies should explore microbial dynamics and their role in chemical cycles in coastal aquifers, the bidirectional dispersal of groundwater and seawater microorganisms, and marine bacterioplankton responses to SGD. This will require not only combining sequencing methods, visualization and linking taxonomy to activity but also considering the entire groundwater–marine continuum. Interactions between traditionally independent disciplines (e.g. hydrogeology, microbial ecology) are needed to frame the study of terrestrial and aquatic microorganisms beyond the limits of their presumed habitats, and to foster our understanding of SGD processes and their influence in coastal biogeochemical cycles.

Metabolic versatility of freshwater sedimentary archaea feeding on different organic carbon sources

Members of the phylum Bathyarchaeota and the class Thermoplasmata are widespread in marine and freshwater sediments where they have been recognized as key players in the carbon cycle. Here, we tested the responsiveness of archaeal communities on settled plant debris and sediment from a karstic lake to different organic carbon amendments (amino acids, plant-derived carbohydrates, and aromatics) using a lab-scale microcosm. Changes in the composition and abundance of sediment and biofilm archaeal communities in both DNA and RNA fractions were assessed by 16S rRNA gene amplicon sequencing and qPCR, respectively, after 7 and 30 days of incubation. Archaeal communities showed compositional changes in terms of alpha and beta diversity in relation to the type of carbon source (amino acids vs. plant-derived compounds), the nucleic acid fraction (DNA vs. RNA), and the incubation time (7 vs. 30 days). Distinct groups within the Bathyarchaeota (Bathy-15 and Bathy-6) and the Thermoplasmata (MBG-D) differently reacted to carbon supplements as deduced from the analysis of RNA libraries. Whereas Bathyarchaeota in biofilms showed a long-term positive response to humic acids, their counterparts in the sediment were mainly stimulated by the addition of tryptophan, suggesting the presence of different subpopulations in both habitats. Overall, our work presents an in vitro assessment of the versatility of archaea inhabiting freshwater sediments towards organic carbon and introduces settled leaf litter as a new habitat for the Bathyarchaeota and the Thermoplasmata.

Complex food webs coincide with high genetic potential for chemolithoautotrophy in fractured bedrock groundwater

Groundwater ecosystems face the challenge of energy limitation due to the absence of light-driven primary production. Lack of space and low oxygen availability might further contribute to generally assumed low food web complexity. Chemolithoautotrophy provides additional input of carbon within the subsurface, however, we still do not understand how abundances of chemolithoautotrophs, differences in surface carbon input, and oxygen availability control subsurface food web complexity. Using a molecular approach, we aimed to disentangle the different levels of potential trophic interactions in oligotrophic groundwater along a hillslope setting of alternating mixed carbonate-/siliciclastic bedrock with contrasting hydrochemical conditions and hotspots of chemolithoautotrophy. Across all sites, groundwater harbored diverse protist communities including Ciliophora, Cercozoa, Centroheliozoa, and Amoebozoa but correlations with hydrochemical parameters were less pronounced for eukaryotes compared to bacteria. Ciliophora-affiliated reads dominated the eukaryotic data sets across all sites. DNA-based evidence for the presence of metazoan top predators such as Cyclopoida (Arthropoda) and Stenostomidae (Platyhelminthes) was only found at wells where abundances of functional genes associated with chemolithoautotrophy were 10-100 times higher compared to wells without indications of these top predators. At wells closer to recharge areas with presumably increased inputs of soil-derived substances and biota, fungi accounted for up to 85% of the metazoan-curated eukaryotic sequence data, together with a low potential for chemolithoautotrophy. Although we did not directly observe higher organisms, our results point to the existence of complex food webs with several trophic levels in oligotrophic groundwater. Chemolithoautotrophy appears to provide strong support to more complex trophic interactions, feeding in additional biomass produced by light-independent CO₂-fixation.

Optimal Growth Temperature and Intergenic Distances in Bacteria, Archaea, and Plastids of Rhodophytic Branch

The lengths of intergenic regions between neighboring genes that are convergent, divergent, or unidirectional were calculated for plastids of the rhodophytic branch and complete archaeal and bacterial genomes. Statistically significant linear relationships between any pair of the medians of these three length types have been revealed in each genomic group. Exponential relationships between the optimal growth temperature and each of the three medians have been revealed as well. The leading coefficients of the regression equations relating all pairs of the medians as well as temperature and any of the medians have the same sign and order of magnitude. The results obtained for plastids, archaea, and bacteria are also similar at the qualitative level. For instance, the medians are always low at high temperatures. At low temperatures, the medians tend to statistically significant greater values and scattering. The original model was used to test our hypothesis that the intergenic distances are optimized in particular to decrease the competition of RNA polymerases within the locus that results in transcribing shortened RNAs. Overall, this points to an effect of temperature for both remote and close genomes.

The endolithic bacterial diversity of shallow bedrock ecosystems

Terrestrial subsurface microbial communities are not restricted to the fluid-filled void system commonly targeted during groundwater sampling but are able to inhabit and dwell in rocks. However, compared to the exploration of the deep biosphere, endolithic niches in shallow sedimentary bedrock have received little interest so far. Despite the potential contribution of rock matrix dwellers to matter cycling and groundwater resource quality, their identity and diversity patterns are largely unknown. Here, we investigated the bacterial diversity in twenty-two rock cores in common limestone-mudstone alternations that differed in rock permeabilities and other geostructural and petrological factors. 16S rRNA gene analysis showed the existence of a unique rock matrix microbiome compared to surrounding groundwater. Typically, shallow weathered limestones contained bacterial groups most likely originating from soil habitats. In low-permeable mudstones, we found similar communities of oligotrophic heterotrophs, and thiosulfate-oxidizing autotrophs, without relation to depth, rock type and bulk rock permeability. In fractured limestone, the bacterial communities of fracture surfaces were distinct from their matrix counterparts and ranged from organic matter decomposers in outcrop areas to autotrophs in downdip positions that receive limited surface input. Contrastingly, rock matrices from lithologically corresponding, but highly isolated environments, were dominated by spore-forming bacteria, oligotrophic heterotrophs and hydrogen-oxidizing autotrophs. Neither depth, matrix permeability nor major mineralogy dominantly controlled the endolithic bacterial diversity. Instead, a combination of subsurface factors drives the supply of niches by fluids, matter and energy as well as the (re)dispersal conditions that likely shape bacterial diversity.

Archaeal Communities in Deep Terrestrial Subsurface Underneath the Deccan Traps, India

Archaeal community structure and potential functions within the deep, aphotic, oligotrophic, hot, igneous provinces of ∼65 Myr old basalt and its Archean granitic basement was explored through archaeal 16S rRNA gene amplicon sequencing from extracted environmental DNA of rocks. Rock core samples from three distinct horizons, basaltic (BS), transition (weathered granites) (TZ) and granitic (GR) showed limited organic carbon (4–48 mg/kg) and varied concentrations (<1.0–5000 mg/kg) of sulfate, nitrate, nitrite, iron and metal oxides. Quantitative PCR estimated the presence of nearly 10³–10⁴ archaeal cells per gram of rock. Archaeal communities within BS and GR horizons were distinct. The absence of any common OTU across the samples indicated restricted dispersal of archaeal cells. Younger, relatively organic carbon- and Fe₂O₃-rich BS rocks harbor Euryarchaeota, along with varied proportions of Thaumarchaeota and Crenarchaeota. Extreme acid loving, thermotolerant sulfur respiring Thermoplasmataceae, heterotrophic, ferrous-/H-sulfide oxidizing Ferroplasmaceae and Halobacteriaceae were more abundant and closely interrelated within BS rocks. Samples from the GR horizon represent a unique composition with higher proportions of Thaumarchaeota and uneven distribution of Euryarchaeota and Bathyarchaeota affiliated to Methanomicrobia, SAGMCG-1, FHMa11 terrestrial group, AK59 and unclassified taxa. Acetoclastic methanogenic Methanomicrobia, autotrophic SAGMCG-1 and MCG of Thaumarcheaota could be identified as the signature groups within the organic carbon lean GR horizon. Sulfur-oxidizing Sulfolobaceae was relatively more abundant in sulfate-rich amygdaloidal basalt and migmatitic gneiss samples. Methane-oxidizing ANME-3 populations were found to be ubiquitous, but their abundance varied greatly between the analyzed samples. Changes in diversity pattern among the BS and GR horizons highlighted the significance of local rock geochemistry, particularly the availability of organic carbon, Fe₂O₃ and other nutrients as well as physical constraints (temperature and pressure) in a niche-specific colonization of extremophilic archaeal communities. The study provided the first deep sequencing-based illustration of an intricate association between diverse extremophilic groups (acidophile-halophile-methanogenic), capable of sulfur/iron/methane metabolism and thus shed new light on their potential role in biogeochemical cycles and energy flow in deep biosphere hosted by hot, oligotrophic igneous crust.

Divergent microbial communities in groundwater and overlying soils exhibit functional redundancy for plant-polysaccharide degradation

Light driven primary production by plants is the main source of biomass in terrestrial ecosystems. But also in subsurface habitats like aquifers, life is fueled largely by this plant-derived biomass. Here, we investigate the degradation of plant-derived polysaccharides in a groundwater microbiome to identify the microbial key players involved, and compare them to those from soil of the groundwater recharge area. We quantified the activities of enzymes degrading the abundant plant polymers starch, cellulose and hemicellulose in oligotrophic groundwater samples, despite the low cell numbers present. Normalized to 16S rRNA gene copy numbers, these activities were only one order of magnitude lower than in soil. Stimulation of the groundwater microbiome with either starch or cellulose and hemicellulose led to changes of the enzymatic activity ratios, indicating autochthonous production of enzymes in response to the plant polymers. Furthermore, DNA stable isotope probing with 13C labelled plant polymers allowed us to identify microbes involved in the degradation of these compounds. In (hemi)cellulose microcosms, Bacteroidia and Candidatus Parcubacteria were active, while the active community in starch microcosms mostly comprised Candidatus Saccharibacteria, Cytophagia, and Actinobacteria. Not a single one of the active OTUs was also found to be labelled in soil microcosms. This indicates that the degradation of plant-derived polysaccharides in groundwater is driven by organisms completely distinct from those active in soil. The involvement of members of the candidate phyla Cand. Parcubacteria and Cand. Saccharibacteria, organisms known to be abundant in groundwater, in plant-derived organic matter degradation might strongly impact subsurface carbon cycling.

14C-Free Carbon Is a Major Contributor to Cellular Biomass in Geochemically Distinct Groundwater of Shallow Sedimentary Bedrock Aquifers

Despite the global significance of the subsurface biosphere, the degree to which it depends on surface organic carbon (OC) is still poorly understood. Here, we compare stable and radiogenic carbon isotope compositions of microbial phospholipid fatty acids (PLFAs) with those of in situ potential microbial C sources to assess the major C sources for subsurface microorganisms in biogeochemical distinct shallow aquifers (Critical Zone Exploratory, Thuringia Germany). Despite the presence of younger OC, the microbes assimilated 14C-free OC to varying degrees; ~31% in groundwater within the oxic zone, ~47% in an iron reduction zone, and ~70% in a sulfate reduction/anammox zone. The persistence of trace amounts of mature and partially biodegraded hydrocarbons suggested that autochthonous petroleum-derived hydrocarbons were a potential 14C-free C source for heterotrophs in the oxic zone. In this zone, Δ14C values of dissolved inorganic carbon (-366 ± 18‰) and 11MeC16:0 (-283 ± 32‰), an important component in autotrophic nitrite oxidizers, were similar enough to indicate that autotrophy is an important additional C fixation pathway. In anoxic zones, methane as an important C source was unlikely since the 13C-fractionations between the PLFAs and CH4 were inconsistent with kinetic isotope effects associated with methanotrophy. In the sulfate reduction/anammox zone, the strong 14C-depletion of 10MeC16:0 (-942 ± 22‰), a PLFA common in sulfate reducers, indicated that those bacteria were likely to play a critical part in 14C-free sedimentary OC cycling. Results indicated that the 14C-content of microbial biomass in shallow sedimentary aquifers results from complex interactions between abundance and bioavailability of naturally occurring OC, hydrogeology, and specific microbial metabolisms.

Impact of hydrologic boundaries on microbial planktonic and biofilm communities in shallow terrestrial subsurface environments

Subsurface environments contain a large proportion of planetary microbial biomass and harbor diverse communities responsible for mediating biogeochemical cycles important to groundwater used by human society for consumption, irrigation, agriculture and industry. Within the saturated zone, capillary fringe and vadose zones, microorganisms can reside in two distinct phases (planktonic or biofilm), and significant differences in community composition, structure and activity between free-living and attached communities are commonly accepted. However, largely due to sampling constraints and the challenges of working with solid substrata, the contribution of each phase to subsurface processes is largely unresolved. Here, we synthesize current information on the diversity and activity of shallow freshwater subsurface habitats, discuss the challenges associated with sampling planktonic and biofilm communities across spatial, temporal and geological gradients, and discuss how biofilms may be constrained within shallow terrestrial subsurface aquifers. We suggest that merging traditional activity measurements and sequencing/-omics technologies with hydrological parameters important to sediment biofilm assembly and stability will help delineate key system parameters. Ultimately, integration will enhance our understanding of shallow subsurface ecophysiology in terms of bulk-flow through porous media and distinguish the respective activities of sessile microbial communities from more transient planktonic communities to ecosystem service and maintenance.

First insights into the living groundwater mycobiome of the terrestrial biogeosphere

Although fungi play important roles in biogeochemical cycling in aquatic ecosystems and have received a great deal of attention, much remains unknown about the living fractions of fungal communities in aquifers of the terrestrial subsurface in terms of diversity, community dynamics, functional roles, the impact of environmental factors and presence of fungal pathogens. Here we address this gap in knowledge by using RNA-based high throughput pair-end illumina sequencing analysis of fungal internal transcribed spacer (ITS) gene markers, to target the living fractions of groundwater fungal communities from fractured alternating carbonate-/siliciclastic-rock aquifers of the Hainich Critical Zone Exploratory. The probed levels of the hillslope multi-storey aquifer system differ primarily in their oxygen and nitrogen content due to their different connections to the surface. We discovered highly diverse living fungal communities (384 Operational Taxonomic Units, OTUs) with different taxonomic affiliations and ecological functions. The observed fungal communities primarily belonged to three phyla: Ascomycota, Basidiomycota and Chytridiomycota. Perceived dynamics in the composition of living fungal communities were significantly shaped by the concentration of ammonium in the moderately agriculturally impacted aquifer system. Apart from fungal saprotrophs, we also detected living plant and animal pathogens for the first time in this aquifer system. This work also demonstrates that the RNA-based high throughput pair-end illumina sequencing method can be used in future for water quality monitoring in terms of living fungal load and subsequent risk assessments. In general, this study contributes towards the growing knowledge of aquatic fungi in terrestrial subsurface biogeosphere.