Microbial diversity in the deep-subsurface hydrothermal aquifer feeding the giant gypsum crystal-bearing Naica Mine, Mexico

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.

Ribosomal proteins can hold a more accurate record of bacterial thermal adaptation compared to rRNA

Ribosomal genes are widely used as 'molecular clocks' to infer evolutionary relationships between species. However, their utility as 'molecular thermometers' for estimating optimal growth temperature of microorganisms remains uncertain. Previously, some estimations were made using the nucleotide composition of ribosomal RNA (rRNA), but the universal application of this approach was hindered by numerous outliers. In this study, we aimed to address this problem by identifying additional indicators of thermal adaptation within the sequences of ribosomal proteins. By comparing sequences from 2021 bacteria with known optimal growth temperature, we identified novel indicators among the metal-binding residues of ribosomal proteins. We found that these residues serve as conserved adaptive features for bacteria thriving above 40°C, but not at lower temperatures. Furthermore, the presence of these metal-binding residues exhibited a stronger correlation with the optimal growth temperature of bacteria compared to the commonly used correlation with the 16S rRNA GC content. And an even more accurate correlation was observed between the optimal growth temperature and the YVIWREL amino acid content within ribosomal proteins. Overall, our work suggests that ribosomal proteins contain a more accurate record of bacterial thermal adaptation compared to rRNA. This finding may simplify the analysis of unculturable and extinct species.

Depth wide distribution and metabolic potential of chemolithoautotrophic microorganisms reactivated from deep continental granitic crust underneath the Deccan Traps at Koyna, India

Characterization of inorganic carbon (C) utilizing microorganisms from deep crystalline rocks is of major scientific interest owing to their crucial role in global carbon and other elemental cycles. In this study we investigate the microbial populations from the deep [up to 2,908 meters below surface (mbs)] granitic rocks within the Koyna seismogenic zone, reactivated (enriched) under anaerobic, high temperature (50°C), chemolithoautotrophic conditions. Subsurface rock samples from six different depths (1,679-2,908 mbs) are incubated (180 days) with CO₂ (+H₂) or HCO₃ ^- as the sole C source. Estimation of total protein, ATP, utilization of NO₃ ^- and SO₄ ^2- and 16S rRNA gene qPCR suggests considerable microbial growth within the chemolithotrophic conditions. We note a better response of rock hosted community towards CO₂ (+H₂) over HCO₃ ^-. 16S rRNA gene amplicon sequencing shows a depth-wide distribution of diverse chemolithotrophic (and a few fermentative) Bacteria and Archaea. Comamonas, Burkholderia-Caballeronia-Paraburkholderia, Ralstonia, Klebsiella, unclassified Burkholderiaceae and Enterobacteriaceae are reactivated as dominant organisms from the enrichments of the deeper rocks (2335-2,908 mbs) with both CO₂ and HCO₃ ^-. For the rock samples from shallower depths, organisms of varied taxa are enriched under CO₂ (+H₂) and HCO₃ ^-. Pseudomonas, Rhodanobacter, Methyloversatilis, and Thaumarchaeota are major CO₂ (+H₂) utilizers, while Nocardioides, Sphingomonas, Aeromonas, respond towards HCO₃ ^-. H₂ oxidizing Cupriavidus, Hydrogenophilus, Hydrogenophaga, CO₂ fixing Cyanobacteria Rhodobacter, Clostridium, Desulfovibrio and methanogenic archaea are also enriched. Enriched chemolithoautotrophic members show good correlation with CO₂, CH₄ and H₂ concentrations of the native rock environments, while the organisms from upper horizons correlate more to NO₃ ^-, SO₄ ^2- _, Fe and TIC levels of the rocks. Co-occurrence networks suggest close interaction between chemolithoautotrophic and chemoorganotrophic/fermentative organisms. Carbon fixing 3-HP and DC/HB cycles, hydrogen, sulfur oxidation, CH₄ and acetate metabolisms are predicted in the enriched communities. Our study elucidates the presence of live, C and H₂ utilizing Bacteria and Archaea in deep subsurface granitic rocks, which are enriched successfully. Significant impact of depth and geochemical controls on relative distribution of various chemolithotrophic species enriched and their C and H₂ metabolism are highlighted. These endolithic microorganisms show great potential for answering the fundamental questions of deep life and their exploitation in CO₂ capture and conversion to useful products.

Seeds of imperfection rule the mesocrystalline disorder in natural anhydrite single crystals

In recent years, we have come to appreciate the astounding intricacies associated with the formation of minerals from ions in aqueous solutions. In this context, a number of studies have revealed that the nucleation of calcium sulfate systems occurs nonclassically, involving the aggregation and reorganization of nanosized prenucleation species. In recent work, we have shown that this particle-mediated nucleation pathway is actually imprinted in the resultant micrometer-sized CaSO₄ crystals. This property of CaSO₄ minerals provides us with the unique opportunity to search for evidence of nonclassical nucleation pathways in geological environments. In particular, we focused on large anhydrite crystals extracted from the Naica Mine in Mexico. We were able to shed light on this mineral's growth history by mapping defects at different length scales. Based on this, we argue that the nanoscale misalignment of the structural subunits, observed in the initial calcium sulfate crystal seeds, propagates through different length scales both in morphological, as well as in strictly crystallographic aspects, eventually causing the formation of large mesostructured single crystals of anhydrite. Hence, the nonclassical nucleation mechanism introduces a "seed of imperfection," which leads to a macroscopic "single" crystal whose fragments do not fit together at different length scales in a self-similar manner. Consequently, anisotropic voids of various sizes are formed with very well-defined walls/edges. However, at the same time, the material retains in part its single crystal nature.

Evolution of the Astonishing Naica Giant Crystals in Chihuahua, Mexico

Calcium sulfate (CaSO4) is one of the most common evaporites found in the earth’s crust. It can be found as four main variations: gypsum (CaSO4∙2H2O), bassanite (CaSO4∙0.5H2O), soluble anhydrite, and insoluble anhydrite (CaSO4), being the key difference the hydration state of the sulfate mineral. Naica giant crystals’ growth starts from a supersaturated solution in a delicate thermodynamic balance close to equilibrium, where gypsum can form nanocrystals able to grow up to 11–12 m long. The growth rates are reported to be as slow as (1.4 ± 0.2) × 10−5 nm/s, taking thousands of years to form crystals with a unique smoothness and diaphaneity, which may or may not include solid or liquid inclusions. Conservation efforts can be traced back to other gypsum structures found prior to Naica’s. Furthermore, in the last two decades, several authors have explored the unique requirements in which these crystals grow, the characterization of their environment and microclimatic conditions, and the prediction of deterioration scenarios. We present a state-of-the-art review on the mentioned topics. Beyond the findings on the origin, in this work we present the current state and the foreseeable future of these astounding crystals.

Revisiting the Relationships Between Genomic G + C Content, RNA Secondary Structures, and Optimal Growth Temperature

Over twenty years ago Galtier and Lobry published a manuscript entitled "Relationships between Genomic G + C Content, RNA Secondary Structure, and Optimal Growth Temperature" in the Journal of Molecular Evolution that showcased the lack of a relationship between genomic G + C content and optimal growth temperature (OGT) in a set of about 200 prokaryotes. Galtier and Lobry also assessed the relationship between RNA secondary structures (rRNA stems, tRNAs) and OGT, and in this case a clear relationship emerged. Increasing structured RNA G + C content (particularly in regions that are double-stranded) correlates with increased OGT. Both of these fundamental relationships have withstood test of many additional sequences and spawned a variety of different applications that include prediction of OGT from rRNA sequence and computational ncRNA identification approaches. In this work, I present the motivation behind Galtier and Lobry's original paper and the larger questions addressed by the work, how these questions have evolved over the last two decades, and the impact of Galtier and Lobry's manuscript in fields beyond these questions.

Microorganisms and Their Metabolic Capabilities in the Context of the Biogeochemical Nitrogen Cycle at Extreme Environments

Extreme microorganisms (extremophile) are organisms that inhabit environments characterized by inhospitable parameters for most live beings (extreme temperatures and pH values, high or low ionic strength, pressure, or scarcity of nutrients). To grow optimally under these conditions, extremophiles have evolved molecular adaptations affecting their physiology, metabolism, cell signaling, etc. Due to their peculiarities in terms of physiology and metabolism, they have become good models for (i) understanding the limits of life on Earth, (ii) exploring the possible existence of extraterrestrial life (Astrobiology), or (iii) to look for potential applications in biotechnology. Recent research has revealed that extremophilic microbes play key roles in all biogeochemical cycles on Earth. Nitrogen cycle (N-cycle) is one of the most important biogeochemical cycles in nature; thanks to it, nitrogen is converted into multiple chemical forms, which circulate among atmospheric, terrestrial and aquatic ecosystems. This review summarizes recent knowledge on the role of extreme microorganisms in the N-cycle in extremophilic ecosystems, with special emphasis on members of the Archaea domain. Potential implications of these microbes in global warming and nitrogen balance, as well as their biotechnological applications are also discussed.

Microbial Composition and Diversity Patterns in Deep Hyperthermal Aquifers from the Western Plain of Romania

A limited number of studies have investigated the biodiversity in deep continental hyperthermal aquifers and its influencing factors. Here, we present the first description of microbial communities inhabiting the Pannonian and Triassic hyperthermal aquifers from the Western Plain of Romania, the first one being considered a deposit of "fossilized waters," while the latter is embedded in the hydrological cycle due to natural refilling. The 11 investigated drillings have an open interval between 952 and 3432 m below the surface, with collected water temperatures ranging between 47 and 104 °C, these being the first microbial communities characterized in deep continental water deposits with outflow temperatures exceeding 80 °C. The abundances of bacterial 16S rRNA genes varied from approximately 10⁵-10⁶ mL^-1 in the Pannonian to about 10²-10⁴ mL^-1 in the Triassic aquifer. A 16S rRNA gene metabarcoding analysis revealed distinct microbial communities in the two water deposits, especially in the rare taxa composition. The Pannonian aquifer was dominated by the bacterial genera Hydrogenophilus and Thermodesulfobacterium, together with archaeal methanogens from the Methanosaeta and Methanothermobacter groups. Firmicutes was prevalent in the Triassic deposit with a large number of OTUs affiliated to Thermoanaerobacteriaceae, Thermacetogenium, and Desulfotomaculum. Species richness, evenness, and phylogenetic diversity increased alongside with the abundance of mesophiles, their presence in the Triassic aquifer being most probably caused by the refilling with large quantities of meteoric water in the Carpathian Mountains. Altogether, our results show that the particular physico-cheminal characteristics of each aquifer, together with the water refilling possibilities, seem to determine the microbial community structure.

On the role, ecology, phylogeny, and structure of dual-family immunophilins

The novel class of dual-family immunophilins (henceforth abbreviated as DFI) represents naturally occurring chimera of classical FK506-binding protein (FKBP) and cyclophilin (CYN), connected by a flexible linker that may include a three-unit tetratricopeptide (TPR) repeat. Here, I report a comprehensive analysis of all current DFI sequences and their host organisms. DFIs are of two kinds: CFBP (cyclosporin- and FK506-binding protein) and FCBP (FK506- and cyclosporin-binding protein), found in eukaryotes. The CFBP type occurs in select bacteria that are mostly extremophiles, such as psychrophilic, thermophilic, halophilic, and sulfur-reducing. Essentially all DFI organisms are unicellular. I suggest that DFIs are specialized bifunctional chaperones that use their flexible interdomain linker to associate with large polypeptides or multisubunit megacomplexes to promote simultaneous folding or renaturation of two clients in proximity, essential in stressful and denaturing environments. Analysis of sequence homology and predicted 3D structures of the FKBP and CYN domains as well as the TPR linkers upheld the modular nature of the DFIs and revealed the uniqueness of their TPR domain. The CFBP and FCBP genes appear to have evolved in parallel pathways with no obvious single common ancestor. The occurrence of both types of DFI in multiple unrelated phylogenetic clades supported their selection in metabolic and environmental niche roles rather than a traditional taxonomic relationship. Nonetheless, organisms with these rare immunophilins may define an operational taxonomic unit (OTU) bound by the commonality of chaperone function.

Metagenomic potential for and diversity of N-cycle driving microorganisms in the Bothnian Sea sediment

The biological nitrogen cycle is driven by a plethora of reactions transforming nitrogen compounds between various redox states. Here, we investigated the metagenomic potential for nitrogen cycle of the in situ microbial community in an oligotrophic, brackish environment of the Bothnian Sea sediment. Total DNA from three sediment depths was isolated and sequenced. The characterization of the total community was performed based on 16S rRNA gene inventory using SILVA database as reference. The diversity of diagnostic functional genes coding for nitrate reductases (napA;narG), nitrite:nitrate oxidoreductase (nxrA), nitrite reductases (nirK;nirS;nrfA), nitric oxide reductase (nor), nitrous oxide reductase (nosZ), hydrazine synthase (hzsA), ammonia monooxygenase (amoA), hydroxylamine oxidoreductase (hao), and nitrogenase (nifH) was analyzed by blastx against curated reference databases. In addition, Polymerase chain reaction (PCR)‐based amplification was performed on the hzsA gene of anammox bacteria. Our results reveal high genomic potential for full denitrification to N₂, but minor importance of anaerobic ammonium oxidation and dissimilatory nitrite reduction to ammonium. Genomic potential for aerobic ammonia oxidation was dominated by Thaumarchaeota. A higher diversity of anammox bacteria was detected in metagenomes than with PCR‐based technique. The results reveal the importance of various N‐cycle driving processes and highlight the advantage of metagenomics in detection of novel microbial key players.

Stable and Variable Parts of Microbial Community in Siberian Deep Subsurface Thermal Aquifer System Revealed in a Long-Term Monitoring Study

The goal of this work was to study the diversity of microorganisms inhabiting a deep subsurface aquifer system in order to understand their functional roles and interspecies relations formed in the course of buried organic matter degradation. A microbial community of a deep subsurface thermal aquifer in the Tomsk Region, Western Siberia was monitored over the course of 5 years via a 2.7 km deep borehole 3P, drilled down to a Palaeozoic basement. The borehole water discharges with a temperature of ca. 50°C. Its chemical composition varies, but it steadily contains acetate, propionate, and traces of hydrocarbons and gives rise to microbial mats along the surface flow. Community analysis by PCR-DGGE 16S rRNA genes profiling, repeatedly performed within 5 years, revealed several dominating phylotypes consistently found in the borehole water, and highly variable diversity of prokaryotes, brought to the surface with the borehole outflow. The major planktonic components of the microbial community were Desulfovirgula thermocuniculi and Methanothermobacter spp. The composition of the minor part of the community was unstable, and molecular analysis did not reveal any regularity in its variations, except some predominance of uncultured Firmicutes. Batch cultures with complex organic substrates inoculated with water samples were set in order to enrich prokaryotes from the variable part of the community. PCR-DGGE analysis of these enrichments yielded uncultured Firmicutes, Chloroflexi, and Ignavibacteriae. A continuous-flow microaerophilic enrichment culture with a water sample amended with acetate contained Hydrogenophilus thermoluteolus, which was previously detected in the microbial mat developing at the outflow of the borehole. Cultivation results allowed us to assume that variable components of the 3P well community are hydrolytic organotrophs, degrading buried biopolymers, while the constant planktonic components of the community degrade dissolved fermentation products to methane and CO2, possibly via interspecies hydrogen transfer. Occasional washout of minor community components capable of oxygen respiration leads to the development of microbial mats at the outflow of the borehole where residual dissolved fermentation products are aerobically oxidized. Long-term community analysis with the combination of molecular and cultivation techniques allowed us to characterize stable and variable parts of the community and propose their environmental roles.

Microbiome composition and geochemical characteristics of deep subsurface high-pressure environment, Pyhäsalmi mine Finland

Pyhäsalmi mine in central Finland provides an excellent opportunity to study microbial and geochemical processes in a deep subsurface crystalline rock environment through near-vertical drill holes that reach to a depth of more than two kilometers below the surface. However, microbial sampling was challenging in this high-pressure environment. Nucleic acid yields obtained were extremely low when compared to the cell counts detected (1.4 × 10(4) cells mL(-1)) in water. The water for nucleic acid analysis went through high decompression (60-130 bar) during sampling, whereas water samples for detection of cell counts by microscopy could be collected with slow decompression. No clear cells could be identified in water that went through high decompression. The high-pressure decompression may have damaged part of the cells and the nucleic acids escaped through the filter. The microbial diversity was analyzed from two drill holes by pyrosequencing amplicons of the bacterial and archaeal 16S rRNA genes and from the fungal ITS regions from both DNA and RNA fractions. The identified prokaryotic diversity was low, dominated by Firmicute, Beta- and Gammaproteobacteria species that are common in deep subsurface environments. The archaeal diversity consisted mainly of Methanobacteriales. Ascomycota dominated the fungal diversity and fungi were discovered to be active and to produce ribosomes in the deep oligotrophic biosphere. The deep fluids from the Pyhäsalmi mine shared several features with other deep Precambrian continental subsurface environments including saline, Ca-dominated water and stable isotope compositions positioning left from the meteoric water line. The dissolved gas phase was dominated by nitrogen but the gas composition clearly differed from that of atmospheric air. Despite carbon-poor conditions indicated by the lack of carbon-rich fracture fillings and only minor amounts of dissolved carbon detected in formation waters, some methane was found in the drill holes. No dramatic differences in gas compositions were observed between different gas sampling methods tested. For simple characterization of gas composition the most convenient way to collect samples is from free flowing fluid. However, compared to a pressurized method a relative decrease in the least soluble gases may appear.

Bacterial gene import and mesophilic adaptation in archaea

It is widely believed that the archaeal ancestor was hyperthermophilic, but during archaeal evolution, several lineages - including haloarchaea and their sister methanogens, the Thaumarchaeota, and the uncultured Marine Group II and Marine Group III Euryarchaeota (MGII/III) - independently adapted to lower temperatures. Recent phylogenomic studies suggest that the ancestors of these lineages were recipients of massive horizontal gene transfer from bacteria. Many of the acquired genes, which are often involved in metabolism and cell envelope biogenesis, were convergently acquired by distant mesophilic archaea. In this Opinion article, we explore the intriguing hypothesis that the import of these bacterial genes was crucial for the adaptation of archaea to mesophilic lifestyles.

New mode of energy metabolism in the seventh order of methanogens as revealed by comparative genome analysis of “Candidatus methanoplasma termitum”

The recently discovered seventh order of methanogens, the Methanomassiliicoccales (previously referred to as “Methanoplasmatales”), so far consists exclusively of obligately hydrogen-dependent methylotrophs. We sequenced the complete genome of “Candidatus Methanoplasma termitum” from a highly enriched culture obtained from the intestinal tract of termites and compared it with the previously published genomes of three other strains from the human gut, including the first isolate of the order. Like all other strains, “Ca. Methanoplasma termitum” lacks the entire pathway for CO2 reduction to methyl coenzyme Mand produces methane by hydrogen-dependent reduction of methanol or methylamines, which is consistent with additional physiological data. However, the shared absence of cytochromes and an energy-converting hydrogenase for the reoxidation of the ferredoxin produced by the soluble heterodisulfide reductase indicates that Methanomassiliicoccales employ a new mode of energy metabolism, which differs from that proposed for the obligately methylotrophic Methanosphaera stadtmanae. Instead, all strains possess a novel complex that is related to the F420:methanophenazine oxidoreductase (Fpo) of Methanosarcinales butlacks an F420-oxidizing module, resembling the apparently ferredoxin-dependent Fpo-like homolog in Methanosaeta thermophila. Since all Methanomassiliicoccales also lack the subunit E of the membrane-bound heterodisulfide reductase (HdrDE), wepropose that the Fpo-like complex interacts directly with subunit D, forming an energy-converting ferredoxin: heterodisulfideoxidoreductase. The dual function of heterodisulfide in Methanomassiliicoccales, which serves both in electron bifurcation and as terminal acceptor in a membrane-associated redox process, may be a unique characteristic of the novel order.

Oxygen availability and distance to surface environments determine community composition and abundance of ammonia-oxidizing prokaroytes in two superimposed pristine limestone aquifers in the Hainich region, Germany

We followed the abundance and compared the diversity of ammonia-oxidizing archaea (AOA) and bacteria (AOB) in the groundwater of two superimposed pristine limestone aquifers located in the Hainich region (Thuringia, Germany) over 22 months. Groundwater obtained from the upper aquifer (12 m depth) was characterized by low oxygen saturation (0-20%) and low nitrate concentrations (0-20 μM), contrasting with 50-80% oxygen saturation and 40-200 μM nitrate in the lower aquifer (48 m and 88 m depth). Quantitative PCR targeting bacterial and archaeal amoA and 16S rRNA genes suggested a much higher ammonia oxidizer fraction in the lower aquifer (0.4-7.8%) compared with the upper aquifer (0.01-0.29%). In both aquifers, AOB communities were dominated by one phylotype related to Nitrosomonas ureae, while AOA communities were more diverse. Multivariate analysis of amoA DGGE profiles revealed a stronger temporal variation of AOA and AOB community composition in the upper aquifer, pointing to a stronger influence of surface environments. Parallel fluctuations of AOA, AOB, and total microbial abundance suggested that hydrological factors (heavy rain falls, snow melt) rather than specific physicochemical parameters were responsible for the observed community dynamics.

Pangenome evidence for extensive interdomain horizontal transfer affecting lineage core and shell genes in uncultured planktonic thaumarchaeota and euryarchaeota

Horizontal gene transfer (HGT) is an important force in evolution, which may lead, among other things, to the adaptation to new environments by the import of new metabolic functions. Recent studies based on phylogenetic analyses of a few genome fragments containing archaeal 16S rRNA genes and fosmid-end sequences from deep-sea metagenomic libraries have suggested that marine planktonic archaea could be affected by high HGT frequency. Likewise, a composite genome of an uncultured marine euryarchaeote showed high levels of gene sequence similarity to bacterial genes. In this work, we ask whether HGT is frequent and widespread in genomes of these marine archaea, and whether HGT is an ancient and/or recurrent phenomenon. To answer these questions, we sequenced 997 fosmid archaeal clones from metagenomic libraries of deep-Mediterranean waters (1,000 and 3,000 m depth) and built comprehensive pangenomes for planktonic Thaumarchaeota (Group I archaea) and Euryarchaeota belonging to the uncultured Groups II and III Euryarchaeota (GII/III-Euryarchaeota). Comparison with available reference genomes of Thaumarchaeota and a composite marine surface euryarchaeote genome allowed us to define sets of core, lineage-specific core, and shell gene ortholog clusters for the two archaeal lineages. Molecular phylogenetic analyses of all gene clusters showed that 23.9% of marine Thaumarchaeota genes and 29.7% of GII/III-Euryarchaeota genes had been horizontally acquired from bacteria. HGT is not only extensive and directional but also ongoing, with high HGT levels in lineage-specific core (ancient transfers) and shell (recent transfers) genes. Many of the acquired genes are related to metabolism and membrane biogenesis, suggesting an adaptive value for life in cold, oligotrophic oceans. We hypothesize that the acquisition of an important amount of foreign genes by the ancestors of these archaeal groups significantly contributed to their divergence and ecological success.

Thaumarchaeal ammonium oxidation and evidence for a nitrogen cycle in a subsurface radioactive thermal spring in the Austrian Central Alps

Previous studies had suggested the presence of ammonium oxidizing Thaumarchaeota as well as nitrite oxidizing Bacteria in the subsurface spring called Franz Josef Quelle (FJQ), a slightly radioactive thermal mineral spring with a temperature of 43.6-47°C near the alpine village of Bad Gastein, Austria. The microbiological consortium of the FJQ was investigated for its utilization of nitrogen compounds and the putative presence of a subsurface nitrogen cycle. Microcosm experiments made with samples from the spring water, containing planktonic microorganisms, or from biofilms, were used in this study. Three slightly different media, enriched with vitamins and trace elements, and two incubation temperatures (30 and 40°C, respectively) were employed. Under aerobic conditions, high rates of conversion of ammonium to nitrite, as well as nitrite to nitrate were measured. Under oxygen-limited conditions nitrate was converted to gaseous compounds. Stable isotope probing with (15)NH4Cl or ((15)NH4)2SO4as sole energy sources revealed incorporation of (15)N into community DNA. Genomic DNA as well as RNA were extracted from all microcosms. The following genes or fragments of genes were successfully amplified, cloned and sequenced by standard PCR from DNA extracts: Ammonia monooxygenase subunit A (amoA), nitrite oxidoreductase subunits A and B (nxrA and nxrB), nitrate reductase (narG), nitrite reductase (nirS), nitric oxide reductases (cnorB and qnorB), nitrous oxide reductase (nosZ). Reverse transcription of extracted total RNA and real-time PCR suggested the expression of each of those genes. Nitrogen fixation (as probed with nifH and nifD) was not detected. However, a geological origin of NH(+) 4 in the water of the FJQ cannot be excluded, considering the silicate, granite and gneiss containing environment. The data suggested the operation of a nitrogen cycle in the subsurface environment of the FJQ.