Potential for hydrogen-oxidizing chemolithoautotrophic and diazotrophic populations to initiate biofilm formation in oligotrophic, deep terrestrial subsurface waters
© The Author(s). 2017
Received: 5 July 2016
Accepted: 12 March 2017
Published: 23 March 2017
Deep terrestrial biosphere waters are separated from the light-driven surface by the time required to percolate to the subsurface. Despite biofilms being the dominant form of microbial life in many natural environments, they have received little attention in the oligotrophic and anaerobic waters found in deep bedrock fractures. This study is the first to use community DNA sequencing to describe biofilm formation under in situ conditions in the deep terrestrial biosphere.
In this study, flow cells were attached to boreholes containing either “modern marine” or “old saline” waters of different origin and degree of isolation from the light-driven surface of the earth. Using 16S rRNA gene sequencing, we showed that planktonic and attached populations were dissimilar while gene frequencies in the metagenomes suggested that hydrogen-fed, carbon dioxide- and nitrogen-fixing populations were responsible for biofilm formation across the two aquifers. Metagenome analyses further suggested that only a subset of the populations were able to attach and produce an extracellular polysaccharide matrix. Initial biofilm formation is thus likely to be mediated by a few bacterial populations which were similar to Epsilonproteobacteria, Deltaproteobacteria, Betaproteobacteria, Verrucomicrobia, and unclassified bacteria.
Populations potentially capable of attaching to a surface and to produce extracellular polysaccharide matrix for attachment were identified in the terrestrial deep biosphere. Our results suggest that the biofilm populations were taxonomically distinct from the planktonic community and were enriched in populations with a chemolithoautotrophic and diazotrophic metabolism coupling hydrogen oxidation to energy conservation under oligotrophic conditions.
Keywords16S rRNA gene Metagenome Binning Terrestrial Deep biosphere Biofilm formation Metabolism
Microbial life in natural environments typically occurs in biofilms adhered to surfaces via a matrix of extracellular polymeric substances (EPS). This is despite the higher energetic costs compared to the free-living state, e.g., up-regulation of genes involved in motility during attachment . Biofilm formation involves a number of sequential steps. These include an establishment phase characterized by motility, cell-to-cell communication, initial adhesion, and EPS production [2, 3]. Living in a biofilm can confer many ecological advantages including more efficient nutrient recycling , genetic exchange , resistance to grazing , stress tolerance , and facilitation of syntrophy plus metabolite exchange .
The deep biosphere is separated from the light-driven surface by the time needed for waters to penetrate to these depths and is also set apart by other environmental conditions such as high pressure, stable reducing conditions, extensive exposure to mineral surfaces, and increasing temperatures with depth. The deep biome is the largest microbial ecosystem on earth , and the deep continental biosphere is estimated to host 2 to 19% of earth’s total biomass . Although a large proportion of the earth’s microbial cells reside in the deep biosphere, our knowledge of the biology and identity of these organisms is scarce . One reason for this is the difficulty to obtain uncontaminated samples. The present study was carried out at the Äspö Hard Rock Laboratory (Äspö HRL). This 460-m-deep underground laboratory circumvents many problems associated with contamination of the low biomass samples of the deep biosphere . These advantages include that the groundwaters are physically separated from the oxidizing environment in the tunnel and that the flow of water into the boreholes is by gravity rather than pumping.
Intrusive igneous rocks make up the vast majority of the earth’s crust. These rocks are typically considerably fractured, providing space for water transport and habitats for microbial life. Despite being highly oligotrophic, at least a portion of the microbes in fracture groundwaters are proposed to be active. This is consistent with the presence of up to ~10 μM concentrations of hydrogen  suggested to support lithotrophic growth [13, 14] and the occurrence of bacteriophages that depend on active microorganisms to survive [15, 16]. Most knowledge of the terrestrial deep biosphere is from studies of free-living cells, revealing communities comprised of anaerobes capable of reducing nitrate, ferric iron, and sulfur or sulfate alongside methanogens and acetogens. All of these energy conservation strategies are supported by carbon dioxide and hydrogen as carbon and energy sources [12, 17–21]. Recent metagenomic studies also point to extensive metabolic versatility with heterotrophic, mixotrophic, and autotrophic metabolic strategies [22, 23]. However, as a substantial fraction of the microbiome in the deep biosphere is likely to live in biofilms [24, 25], the earlier work may provide an incomplete and possibly misleading description of this biome.
Despite biofilms in rock fractures being present in the deep biosphere [24–26], the factors controlling biofilm formation are poorly understood, especially in deep subsurface environments. This study is the first study to use metagenomic sequencing to address this knowledge gap.
Results and discussion
Geochemical characteristics of the borehole waters
Chemical and isotopic composition of the groundwaters. Concentrations of chemical elements and ions are given in mg L−1 except for Cs which is in μg L−1. The 18O/16O ratio (δ18O) is reported via the δ notation “per mil” as a deviation from the Standard Mean Ocean Water (SMOW)
−7.8 to −7.4
−12.4 to −12.4
KA2198A had high magnesium and potassium concentrations, which are tracers of marine waters , and also had slightly lower values for chloride and δ18O compared to modern brackish Baltic Sea water  (Table 1). This implies that this groundwater mainly consisted of infiltrated brackish Baltic Sea water, to some extent diluted with meteoric water, and was classified as “modern marine” (sample name defined as “MM”). The precise infiltration age of the groundwater was <20 years and probably even more recent. The groundwater cesium concentrations (4.5–4.8 μg L−1) were much higher than in surface waters in the region , indicating that the trace-metal hydrochemistry had been altered during and after the marine water infiltration . The latest chemical measurements of this groundwater (from 2009) predate the microbiological experiments (Additional file 2: Figure S2). However, the chemical characteristics were unlikely to have changed as, e.g., the chloride and sulfate concentrations in this type of groundwater within the Äspö HRL have not altered over time scales of decades . KF0069A01 had typically high chloride concentrations and relatively low δ18O values for saline groundwater with a hydrological residence time in the order of millions of years  (Table 1). This groundwater was thus termed “old saline” (defined as “OS”). Its high concentrations of calcium, cesium, and sulfate result from mineral weathering and dissolution over millions of years , and the water also contains relatively low amounts of dissolved organic carbon, bicarbonate, and ammonium. This type of groundwater within the Äspö HRL has had overall stable chloride concentrations since the early 21st century , and therefore, the geochemical measurements of this groundwater (2003–2008) can be considered representative of the prevailing geochemistry in the borehole during the microbiological experiments.
Sampling and sequencing data
Flow cells were directly connected to boreholes and the groundwater was allowed to pass under in situ temperature and pressure for 33 days. This time was sufficient to allow investigation of initial biofilm formation. The flow cells were loaded with garnet grains and glass beads as solid support for biofilm growth. These surfaces were chosen as (i) they were sterile, DNA-free, and RNAse/DNAse-free and could be sterilized by heating to 450 °C, respectively; (ii) in testing biofilm formation at the Äspö HRL, it was found that bedrock from the same environment was porous and that unattached minerals and particles disrupted DNA extractions from the formed biofilms [13, 14, 32]; and (iii) the fracture surfaces in the bedrock are mineralogically very heterogeneous consisting of various proportions of primary minerals (e.g., quartz, feldspar, plagioclase, and mica) and secondary precipitates (e.g., calcite, pyrite, epidote, Fe-oxides, and clay minerals) [33, 34]. Therefore, the flow cell system was simplified by using a single silicate mineral (garnet) that in terms of biofilm formation can be considered as representative for the silicate rocks (granites) in the Äspö HRL. However, the drawback of these solid supports was that they have different characteristics to the rock surfaces that potentially affected the initial biofilm forming populations.
The volume of groundwater that passed through the flow cells and details of the 16S rRNA gene sequencing data are given in Additional file 3: Table S1. The Illumina sequencing yielded 1.41 × 108 to 1.70 × 108 raw reads and 1.38 × 108 to 1.65 × 108 trimmed reads per sample that assembled into 5671 to 48857 contigs ≥1000 bp in length (Additional file 4: Table S2). The contigs were binned into 33 garnet and 35 glass near-complete metagenome-assembled genomes (MAGs) from the modern marine water and 11 garnet and 9 glass MAGs from the old saline, representing 44.4 ± 5.3% of the reads. The MAGs contained ≥31 of the 36 CONCOCT single-copy genes, with an estimated bin completeness of ≥86% (Additional file 5: Table S3).
Estimation of biofilm and planktonic cell numbers
Based upon normalized ATP measurements for the different conditions in the two contrasting groundwaters, the modern marine water flow cells had consistently higher estimated cell abundances with 4.6 × 106 and 2.0 × 108 cells cm2 on the garnet (sample name defined as “MMR”) and glass surfaces (defined as “MMG”), respectively. This compared to 1.9 × 104 and 4.3 × 103 cells cm2 on the corresponding solid surfaces from the old saline water (defined as “OSR” and “OSG”; Additional file 6: Table S4). The higher cell abundance in the biofilms of the modern marine water was consistent with the greater DOC in this water type (Table 1). For comparison, calculated planktonic cell concentrations for the modern marine borehole was 2.8 × 104 cells mL−1 . The planktonic cell concentration was not measured for the old saline water, but water from an adjacent borehole with similar chemical characteristics was estimated to hold approximately 100 cells mL−1 . Hence, the biofilms were quantitatively significant components of these deep aquifer systems and should be considered if we are to understand the interplay between the biosphere and geosphere.
Even with the extremely oligotrophic conditions in the Äspö HRL fracture waters, earlier data from cell growth in response to addition of putative substrates such as hydrogen and carbon dioxide, suggests that deep biosphere microorganisms are viable [13, 14, 36]. In the present study, direct measurements of ATP showed that the investigated microbial communities were metabolically active (Additional file 6: Table S4). This is further corroborated by the cells attaching to and colonizing solid surfaces within a 33-day period without any experimental manipulation of resources and with respective numbers of cells per square centimeter in the modern marine garnet and glass biofilms being ~164 and ~7000 fold greater than the number of planktonic cells in a milliliter of water.
Microbial diversity from 16S rRNA gene sequencing and metagenome binning
Relative abundance (%) of OTUs in bacterial 16S rRNA gene v4v6 sequence libraries. Sequences with ≥1% abundance frequency are shown in detail. The data are from three PCR-generated amplicon samples that were pooled before sequencing
Relative abundance (%)
Candidate division OP3f
The biofilm MAGs were compared with planktonic cell MAGs from borehole groundwaters in the Äspö HRL . This suggested MMG_Bin_14 (Candidatus OD1) was similar to a population present in borehole SA1229A also containing modern marine water. In addition, OSR_Bin_45, OSG_Bin_23, and OSG_Bin_24 aligning with Thiobacillus denitrificans were similar to a planktonic population from old saline water borehole KA3385A:1. The suggested overlap of just two MAG populations between planktonic and biofilm populations from separate studies supported the observation that these communities were dissimilar. In addition, the identified community was different from an established biofilm previously recovered from an in situ rock surface , suggesting that the biofilm community would further develop as it matures, such as by the inclusion of species incapable of initiating the biofilm formation.
Partition of genes between biofilm and planktonic populations
Potential for biofilm formation
The flow cells were run for a total of 33 days and generation times for attached microbial populations in the studied deep hard rock fractures have been extrapolated from cell counts to be between 16 and 90 days [14, 42]. Therefore, the cells attached to the solid surfaces had likely not proliferated more than at most one generation and largely represent cells imported from the aqueous phase. This experimental design allowed the investigation of populations with the genetic potential to initiate biofilm formation prior to the attachment of later colonizing populations that lack such traits.
Summary of the metabolic pathways for the reconstructed genomes. The table provides gene homologs for chemotaxis, motility, and EPS production and secretion for biofilm formation
Metabolic pathways and energy conservation
% mapped readsa
Modern marine water
EPS production, secretion and motility
H2 oxidation, sulfate and nitrate reduction and Rnf
MMR_Bins_36, 41 and 98/MMG_Bins_22, 48 and 93
EPS production and secretion
(Fermentation)b, (formate and H2 oxidation), (Sulfur reduction), (nitrate reduction), (denitrification) and (Rnf)
Formate and H2 oxidation, sulfate/sulfur reduction and Rnf
CO2 and N2
Old saline water
EPS production, secretion and motility
Formate oxidation, Sulfate sulfur reduction
CO2 and N2
EPS production and secretion
Fermentation, H2 oxidation, Rnf
OSG_Bins_16, 23 and 24
(CO2) and (N2)
In the modern marine water, only one MAG (MMR_Bin_58; assigned to Group A) featured all traits for biofilm formation capacity (i.e., chemotaxis, flagellum-assisted motility, and EPS production and export). This MAG represented 0.13% of the mapped reads from modern marine water garnet surface biofilm. In contrast, none of the MAGs from the glass biofilm contained all these traits for surface attachment and biofilm formation. This may have been due to syntrophic interactions between populations able to e.g., attach to the surface allowing subsequent populations that lack this ability to attach and build a viable biofilm.
The second group of potential biofilm-producing populations (able to produce EPS and transport it to the cell surface) included MMR_Bin_58 was most closely related to Sulfuricella denitrificans, a facultative anaerobe that couples sulfur oxidation with nitrate reduction . However, the metagenome assembly also suggests that MMR_Bin_58 grows via anaerobic hydrogen oxidation coupled to sulfate reduction, potentially producing an ion motive force via the ferredoxin:NAD+ oxidoreductase Rnf complex and fixing carbon dioxide via the Calvin-Benson-Bassham (CBB) cycle. Three MAGs from the modern marine water garnet surface (totaling 1.2% of the mapped reads) and three MAGs from the modern marine water glass surface (13.4% of the mapped reads) contained genes implicated in EPS production and export (assigned to Group B). MMR_Bin_36, MMR_Bin_41, and MMG_Bin_48 were affiliated with the Verrucomicrobia and MMG_Bin_22 was from candidate division OP3. All of these MAGs were suggested to generate an ion motive force via the Rnf complex. In addition, MMR_Bin_41 and MMG_Bin_48 also had the potential to utilize hydrogen via type III anaerobic hydrogen oxidation. MMR_Bin_98 (assigned to a poorly defined clade within the Deltaproteobacteria) and MMG_Bin_93 (most related to Sulfurimonas denitrificans) were suggested to have the potential to oxidize formate (fdh gene), convert nitrate to nitrite or nitrogen gas (nar/nap and nar/nir/nor/nos genes, respectively), and fix nitrogen for cellular growth. In modern marine water, only one MAG (MMG_Bin_17; 3.10% of the mapped reads) was suggested to be able to sense chemical stimuli and move via a flagellum and was defined as Group C. This MAG contained genes assigned to hydrogen and formate oxidation, sulfate and sulfur reduction, carbon dioxide and nitrogen fixation, and the Rnf complex. Although the modern marine water was ostensibly anaerobic, both MMG_Bin_17 and OSG_Bin_16 were most closely related to Sideroxydans lithotrophicus, an iron-oxidizing bacterium that grows at oxic-anoxic interfaces .
The old saline water contained one MAG (OSR_Bin_1; 0.22% of the reads) that was suggested to code for all genes required for biofilm formation (group A). This population affiliated most closely to D. aespoeensis  and was suggested to reduce sulfate to sulfide during anaerobic respiration. However, OSR_Bin_1 also contained gene homologs suggested to encode enzymes involved in formate oxidation, carbon dioxide fixation via the CBB cycle, and nitrogen assimilation. The old saline water biofilms also contained OSR_Bin_39 and OSG_Bin_9 from Group B (0.65% and 0.34% of mapped reads, respectively) that contain genes to produce and export EPS. These bins were potentially able to metabolize pyruvate to either acetate or ethanol (all genes in one of the eight fermentation pathways to acetate or two out of three genes for the three pyruvate fermentation pathways to ethanol), oxidize hydrogen, and generate an ion motive force via the Rnf complex. Group C, suggested to sense chemical stimuli and move via a flagellum, was represented by one MAG on the garnet surface of the old saline water (OSR_Bin_45; 1.7% of reads) and three MAGs on the glass surfaces (OSG_Bin_16, OSG_Bin_23 and OSG_Bin_24; totaling 2.8% of reads). OSR_Bin_45 and OSG_Bin_24 from the old saline water were both phylogenetically similar to T. denitrificans . These MAGs contained genes assigned to carbon dioxide fixation by the CBB cycle, nitrogen assimilation, formate oxidation to carbon dioxide, sulfate and nitrate reduction, and generation of an ion motive force via the Rnf complex. OSG_Bin_16 showed the potential ability to reduce sulfate to sulfide, fix carbon dioxide and nitrogen, and generate a proton motive force via the Rnf complex.
Potential metabolism in deep biosphere biofilms
All MAGs were scrutinized for genes encoding electron donor and acceptor use as well as carbon dioxide and nitrogen fixation (Table 3 and Additional file13: Table S7). There were no major differences in the predicted metabolic pathways for MAGs implicated in initial biofilm formation between either the water types or surfaces.
The metabolic potential of the surface attached communities’ implied organic carbon (e.g., MMR_Bin_67 and OSG_Bin_4) and/or hydrogen oxidation (e.g., MMG_Bin_17 and OSR_Bin_28) coupled to reduction of sulfur or sulfate for energy acquisition. Additional populations were suggested to link organic carbon (e.g., MMG_Bin_16 and OSR_Bin_36) and/or hydrogen oxidation (MMR_Bin_58 and OSR_Bin_6) to nitrate reduction. Both of these observations are in agreement with previous studies from the Äspö HRL [22, 51]. The dominant populations in both the modern marine and old saline water biofilm formers were putative carbon- and nitrogen-fixers (MMG_Bin_93 representing 11.73% of the reads and OSG_Bin_24 representing 2.15% of the reads). The importance of previously unknown microorganisms with a simple, fermentative growth strategy has recently been recognized at the Äspö HRL  as well as at another site . These earlier findings are consistent with the metabolic features of several biofilm populations identified here (e.g., MMR_Bin_23 and OSR_Bin_0). Hence, the biofilm community seems able to utilize a range of carbon and energy sources.
Populations predicted to carry out all stages of biofilm formation in the deep terrestrial biosphere were taxonomically distinct from the planktonic community. The biofilm contained a mixed community dominated by hydrogen-fed autotrophs able to fix nitrogen that reflected the oligotrophic conditions in the waters of deep terrestrial aquifers.
Description of the Äspö HRL site
The Swedish Nuclear Fuel and Waste Management Company (SKB) operated Äspö HRL is situated on the southeastern coast of Sweden (Lat N 57° 26′ 4′′ Lon E 16° 39′ 36′′) in a bedrock dominated by 1800 Ma granite and quartz monzodiorite [12, 52]. The structure of the Äspö HRL tunnel has been shown in a previous study , and the geology, chemistry, and hydrology of the boreholes extending from the tunnel have been described [28, 54, 55]. The boreholes were sampled for planktonic cells, and flow cells attached to investigate for biofilm development (described below).
Chemical and isotopic composition of the water
Descriptions of the analytical techniques and precision of the variables are given elsewhere: Cl− and δ18O ; Na, K, Ca, Mg, Cs, and NH4 + ; Fe2+, Fe-total, dissolved organic carbon (DOC), and HS− ; and Mn, SO4 2−, HCO3 −, and PO4 3− . Multiple measurements were taken over a maximum of 7 years with median plus minimum/maximum values presented.
Biofilm growth, sampling, DNA preparation, and sequencing
The flow cells were directly attached to the boreholes in the Äspö HRL tunnel. They had a stainless steel shell (length 300 mm, diameter 65 mm), were lined with polyvinyldifloride plastic, and were equipped with manometers and a pressure relief valve to enable biofilms to form under the high in situ pressure and low redox conditions prevailing in the fissures intersected by the boreholes (Additional file 1: Figure S1). Each flow cell had a 120-mm-long polyvinyldifloride insert with a 22 × 32 mm opening that supported ~100 g of sterile, DNA-free and RNAse/DNAse-free garnet grains (0.7 mm in diameter; MOBIO Laboratories (USA)) and glass beads (1 mm in diameter; VWR International). Glass beads were sterilized by heating to 450 °C for 5 h in a muffle furnace. The flow cells were connected to boreholes KA2198A and KF0069A01 from 23 May 2013 until 24 June 2013 (33 days). Immediately after disconnection from the borehole, the capped flow cells were transported at 4 °C to the laboratory (transport time <8 h), and DNA extraction was carried out on the same day. Approximately, 6 g of each support material was collected from the flow cells and used for DNA preparation.
Planktonic cells from the boreholes were collected on 2 February 2012 for KA2198A and 20 July 2012 for KF0069A01 by filtration under pressure through 47-mm diameter and 0.22-μm pore size membrane filters (supplied with the PowerWater® DNA Isolation Kit, (MO BIO Laboratories, Immuno diagnostics, Hämeenlinna, Finland) contained within a stainless steel filter holder (Millipore) at a flow rate of 0.2 L min−1 for 16 and 16.5 h, respectively. The filters were aseptically removed, placed in sample tubes provided with the PowerWater Kit, and stored frozen at −20 °C until DNA extraction. Genomic DNA from groundwater filters plus garnet grains and glass beads was extracted using the PowerWater Kit according to the protocol provided by the manufacturer.
A Bacterial 16S rDNA v4v6 amplicon library for sequencing was generated by using the degenerative forward (518F, CCAGCAGCYGCGGTAAN)  and reverse primer (1064R, CGACRRCCATGCANCACCT) . Conditions for the PCR reaction were 1× Platinum HiFi Taq polymerase buffer, 1.6 units Platinum HiFi polymerase, 3.7 mM MgSO4, 200 μM dNTPs (PurePeak polymerization mix, ThermoFisher), and 400 nM primers. Between 5 and 25 ng of sample DNA was added to a master mix to a final volume of 100 μL, and this was divided into three replicate 33-μL reactions. Cycling conditions included an initial denaturation at 94 °C for 3 min; 30 cycles of 94 °C for 30 s, 57–60 °C for 45 s, and 72 °C for 1 min; and a final extension at 72 °C for 2 min using a Bio-Rad mycycler. The quality and concentration of the amplicon library was evaluated by using the Agilent Tapestation 2000 instrument according to the manufacturer’s protocol. The reactions were cleaned and products under 300 bp were removed using AMpure beads at 0.75 × volume (Beckman Coulter, Brea CA). The final products were re-suspended in 100 μL of 10 mM Tris-EDTA + 0.05% Tween-20, quantified using PicoGreen Quant-IT assay (Life Technologies), and assayed once again on the Tapestation 2000 instrument. Amplicons were further titrated in equimolar concentration before emulsion-PCR based on their dsDNA concentrations. A GS-FLX Sequencer was used to generate pyrotag sequence reads with the Roche Titanium reagents.
Metagenome libraries of extracted DNA from garnet grains and glass beads were prepared using the ThruPlex DNA-seq Kit with 96 dual indexes (Rubicon Genomics, MI, USA) using an Agilent NGS workstation (Agilent, CA, USA) and purified [59, 60]. The libraries were sequenced on an Illumina HiSeq (2 × 150 bp) in rapid mode at the Science for Life Laboratory in Stockholm, Sweden.
The 16S rRNA gene amplicon sequencing data was first trimmed to remove primer bases, barcodes, and low quality sequences . The trimmed sequences were screened for chimeras by using the UCHIME algorithm . Clustering into operational taxonomic units (OTUs) was by an open reference OTU-picking methodology with the USEARCH algorithm which uses both de novo and reference-based approaches . Representative sequences were chosen from each OTU that were phylogenetically classified against the SILVA database (SILVA123_QIIME-release). Rarefaction to the lowest number of sequences (n = 10 000) was used to normalize sample count before analysis of the datasets.
Metagenome analysis was carried out as previously described . In brief, adapters were removed with Seqprep before the sequences were trimmed and assembled using Sickle , Ray (version 2.3.1) , and Newbler (version 2.6). The assembled contigs were then binned to individual near-complete genomes using CONCOCT (version 0.3.0) . CONCOCT uses 36 single-copy genes to evaluate the coverage of the assembled genome bins. The bins which had ≥31 single-copy genes and ≤2 duplicated single-copy genes were chosen for individual taxonomic and functional annotation using Phylosift v1.0.  and Prokka v1.10 . Phylosift uses a suite of 37 marker genes for phylogenetic classification that are available at https://phylosift.wordpress.com/tutorials/scripts-markers/ . Gene frequency comparisons between planktonic and biofilm samples were calculated by an in-house pipeline. In brief, all trimmed reads were first co-assembled using MEGAHIT v1.0.3  with a minimum kmer size of 31 and a maximum kmer size of 81 at a step of 10. The assembled contigs (of ≥1000 bp in length) were then annotated using Prokka v1.10 . Coverage (% of a locus represented in the assembly) of each gene predicted and annotated by Prokka was calculated using bedtools v2.17.0  and the scripts prokkagff2bed.sh and get_coverage_for_genes.py, after the raw reads were mapped back to the assembly using the script map-bowtie2-markduplicates.sh. These scripts, developed by the Environmental Genomics group at SciLifeLab Stockholm, are available at http://metagenomics-workshop.readthedocs.io/. Gene coverage was calculated and normalized by dividing the coverage values by the total coverage of the sample. Frequencies of the genes of interest were computed and normalized by dividing the coverage values by the total coverage of the sample. The one-way ANOVA values were then calculated on the genes used to define the presence or absence of a pathway in the MAGs by averaging the normalized values and dividing the coverage values by the total coverage of the sample. The one-way ANOVAs were based on the biofilm values for garnet and glass (number of replicates =1 each) compared to large and small planktonic cells (number of replicates =2 each).
ATP measurement and calculation of cell numbers
ATP was measured from the cells using the ATP Biomass Kit HS (BioThema, Sweden). The garnet grains and glass beads were added to 1 mL of reagent BS, vortexed for 30 s, and placed in the dark at room temperature for 30 min. Thereafter, the ATP was analyzed as previously described . After analysis, the garnets and glass beads were washed, dried, and weighed for calculation of their total surface area. The total sampled surface area ranged from 8.5 to 15.0 cm2. The number of cells was then divided by the total area of the respective surface to get the approximate number of cells per cm2. ATP measurements were carried out in triplicate and averages ± standard deviation are presented.
- Äspö HRL:
Äspö Hard Rock Laboratory
Dissolved organic carbon
Extracellular polymeric substances
All authors thank the Swedish Nuclear Fuel and Waste Management Company for access to the Äspö HRL tunnel, laboratories, and Sicada database. Pyrotag sequencing was performed at the Marine Biological Laboratory (Woods Hole, MA, USA), where we received excellent assistance from Sharon Grim, Hilary Morrison, Susan Huse, Mitch Sogin, Joseph Vineis, and Andrew Voorhis. Metagenome sequencing was performed at the National Genomics Infrastructure within SciLifeLab. Bioinformatics utilized the Uppsala Multidisciplinary Center for Advanced Computational Science resource at Uppsala University (project b2013127). We acknowledge Johannes Alneberg, John Sundh, Luisa Hugerth, and Ino de Bruijn at SciLifeLab Stockholm for developing the scripts for quantitative metagenome analysis.
This work was supported by the grant from the Swedish Research Council, Vetenskapsrådet (contracts 2014-4398 and 2012-3892) to MD, MÅ, AA, and SB; The Crafoord Foundation (contract 20130557) to MD; Nova Center for University Studies Research and Development, and Familjen Hellmans stiftelse to MD; and the Swedish Research Council Formas to SB. Pyrotag sequencing was funded by the Deep Carbon Observatory’s Census of Deep Life program supported by the Alfred P. Sloan Foundation.
Availability of data and materials
The amplicon sequencing data supporting the results of this article are available in the Short Read Archive (SRA) with accession numbers SRP041926 for biofilms from KF0069A01 and KA2198A and SRP041904 for the corresponding groundwaters. The metagenome sequencing data are deposited in the Sequence Read Archive (NCBI) with the accession numbers: SRR2540949 (MMR), SRR2544004 (MMG), SRR2544011 (OSR), and SRR2544012 (OSG).
KP, AA, SB, and MD designed the study; XW, JE, LE and AA performed the experimental work and carried out the in silico analysis; XW, KP, SB, and MD wrote the paper. All authors read and approved the final manuscript.
The authors KP, JE, and LE are employed by the company, Microbial Analytics Sweden AB. However, this study was carried out in the absence of a commercial contract. Therefore, the authors declare that they have no competing financial interests.
Consent for publication
Ethics approval and consent to participate
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Kalmokoff M, Lanthier P, Tremblay TL, Foss M, Lau PC, Sanders G, Austin J, Kelly J, Szymanski CM. Proteomic analysis of Campylobacter jejuni 11168 biofilms reveals a role for the motility complex in biofilm formation. J Bacteriol. 2006;188:4312–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Hori K, Matsumoto S. Bacterial adhesion: from mechanism to control. Biochem Engin J. 2010;48:424–34.View ArticleGoogle Scholar
- Petrova OE, Sauer K. Sticky situations: key components that control bacterial surface attachment. J Bacteriol. 2012;194:2413–25.View ArticlePubMedPubMed CentralGoogle Scholar
- van Gestel J, Vlamakis H, Kolter R. Division of labor in biofilms: the ecology of cell differentiation. Microbiol Spectrum. 2015;3:MB-0002-2014.
- Merod RT, Wuertz S. Extracellular polymeric substance architecture influences natural genetic transformation of Acinetobacter baylyi in biofilms. Appl Environ Microbiol. 2014;80:7752–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Sun SY, Tay QXM, Kjelleberg S, Rice SA, McDougald D. Quorum sensing-regulated chitin metabolism provides grazing resistance to Vibrio cholerae biofilms. ISME J. 2015;9:1812–20.View ArticlePubMedPubMed CentralGoogle Scholar
- Periasamy S, Nair HAS, Lee KWK, Ong J, Goh JQJ, Kjelleberg S, Rice SA. Pseudomonas aeruginosa PA01 exopolysaccharides are important for mixed species biofilm community development and stress tolerance. Front Microbiol. 2015;6: doi:10.3389/fmicb.2015.00851.
- Morris BE, Henneberger R, Huber H, Moissl-Eichinger C. Microbial syntrophy: interaction for the common good. FEMS Microbiol Rev. 2013;37:384–406.View ArticlePubMedGoogle Scholar
- Edwards KJ, Becker K, Colwell F. The deep, dark energy biosphere: Intraterrestrial life on earth. In Annu Rev Earth Planet Sci. Volume 40. Edited by Jeanloz R. Palo Alto; 2012: 551-568
- McMahon S, Parnell J. Weighing the deep continental biosphere. FEMS Microbiol Ecol. 2014;87:113–20.View ArticlePubMedGoogle Scholar
- Wilkins MJ, Daly RA, Mouser PJ, Trexler R, Sharma S, Cole DR, Wrighton KC, Biddle JF, Denis EH, Fredrickson JK, et al. Trends and future challenges in sampling the deep terrestrial biosphere. Front Microbiol. 2014;5:481.PubMedPubMed CentralGoogle Scholar
- Hallbeck L, Pedersen K. Characterization of microbial processes in deep aquifers of the Fennoscandian Shield. Appl Geochem. 2008;23:1796–819.View ArticleGoogle Scholar
- Pedersen K. Influence of H2 and O2 on sulphate-reducing activity of a subterranean community and the coupled response in redox potential. FEMS Microbiol Ecol. 2012;82:653–65.View ArticlePubMedGoogle Scholar
- Pedersen K. Subterranean microbial populations metabolize hydrogen and acetate under in situ conditions in granitic groundwater at 450 m depth in the Äspö Hard Rock Laboratory, Sweden. FEMS Microbiol Ecol. 2012;81:217–29.View ArticlePubMedGoogle Scholar
- Eydal HS, Jagevall S, Hermansson M, Pedersen K. Bacteriophage lytic to Desulfovibrio aespoeensis isolated from deep groundwater. ISME J. 2009;3:1139–47.View ArticlePubMedGoogle Scholar
- Kyle JE, Eydal HS, Ferris FG, Pedersen K. Viruses in granitic groundwater from 69 to 450 m depth of the Äspö hard rock laboratory, Sweden. ISME J. 2008;2:571–4.View ArticlePubMedGoogle Scholar
- Pedersen K. Metabolic activity of subterranean microbial communities in deep granitic groundwater supplemented with methane and H2. ISME J. 2013;7:839–49.View ArticlePubMedGoogle Scholar
- Itavaara M, Nyyssonen M, Kapanen A, Nousiainen A, Ahonen L, Kukkonen I. Characterization of bacterial diversity to a depth of 1500 m in the Outokumpu deep borehole, Fennoscandian Shield. FEMS Microbiol Ecol. 2011;77:295–309.View ArticlePubMedGoogle Scholar
- Nyyssonen M, Hultman J, Ahonen L, Kukkonen I, Paulin L, Laine P, Itavaara M, Auvinen P. Taxonomically and functionally diverse microbial communities in deep crystalline rocks of the Fennoscandian Shield. ISME J. 2014;8:126–38.View ArticlePubMedGoogle Scholar
- Purkamo L, Bomberg M, Nyyssonen M, Kukkonen I, Ahonen L, Itavaara M. Heterotrophic communities supplied by ancient organic carbon predominate in deep Fennoscandian bedrock fluids. Microb Ecol. 2015;69:319–32.View ArticlePubMedGoogle Scholar
- Dzaugis ME, Spivack AJ, Dunlea AG, Murray RW, D'Hondt S. Radiolytic hydrogen production in the subseafloor basaltic aquifer. Front Microbiol. 2016;7. 10.3389/fmicb.2016.00076.
- Wu X, Holmfeldt K, Hubalek V, Lundin D, Åström M, Bertilsson S, Dopson M: Microbial metagenomes from three aquifers in the Fennoscandian Shield terrestrial deep biosphere reveal metabolic partitioning among populations. ISME J. 2015;10:1192–203.
- Bagnoud A, Chourey K, Hettich RL, de Bruijn I, Andersson AF, Leupin OX, Schwyn B, Bernier-Latmani R. Reconstructing a hydrogen-driven microbial metabolic network in Opalinus Clay rock. Nat Commun. 2016;7:12770.View ArticlePubMedPubMed CentralGoogle Scholar
- Wanger G, Southam G, Onstott TC. Structural and chemical characterization of a natural fracture surface from 2.8 kilometers below land surface: biofilms in the deep subsurface. Geomicrobiol J. 2006;23:443–52.View ArticleGoogle Scholar
- Pfiffner SM, Cantu JM, Smithgall A, Peacock AD, White DC, Moser DP, Onstott TC, van Heerden E. Deep subsurface microbial biomass and community structure in Witwatersrand Basin mines. Geomicrobiol J. 2006;23:431–42.View ArticleGoogle Scholar
- Jägevall S, Rabe L, Pedersen K. Abundance and diversity of biofilms in natural and artificial aquifers of the Äspö Hard Rock Laboratory, Sweden. Microb Ecol. 2011;61:410–22.View ArticlePubMedGoogle Scholar
- Gimeno MJ, Auqué LF, Acero P, Gómez JB. Hydrogeochemical characterisation and modelling of groundwaters in a potential geological repository for spent nuclear fuel in crystalline rocks (Laxemar, Sweden). Appl Geochem. 2014;45:50–71.View ArticleGoogle Scholar
- Mathurin FA, Astrom ME, Laaksoharju M, Kalinowski BE, Tullborg EL. Effect of tunnel excavation on source and mixing of groundwater in a coastal granitoidic fracture network. Environ Sci Technol. 2012;46:12779–86.View ArticlePubMedGoogle Scholar
- Mathurin FA, Drake H, Tullborg E-L, Berger T, Peltola P, Kalinowski BE, Åström ME. High cesium concentrations in groundwater in the upper 1.2 km of fractured crystalline rock—Influence of groundwater origin and secondary minerals. Geochim Cosmochim Acta. 2014;132:187–213.View ArticleGoogle Scholar
- Mathurin FA, Åström ME, Drake H, Maskenskaya OM, Kalinowski BE. REE and Y in groundwater in the upper 1.2 km of Proterozoic granitoids (Eastern Sweden)—Assessing the role of composition and origin of groundwaters, geochemistry of fractures, and organic/inorganic aqueous complexation. Geochim Cosmochim Acta. 2014;144:342–78.View ArticleGoogle Scholar
- Louvat D, Michelot JL, Aranyossy JF. Origin and residence time of salinity in the Äspö groundwater system. Appl Geochem. 1999;14:917–25.View ArticleGoogle Scholar
- Pedersen K, Bengtsson AF, Edlund JS, Eriksson LC. Sulphate-controlled diversity of subterranean microbial communities over depth in deep groundwater with opposing gradients of sulphate and methane. Geomicrobiol J. 2014;31:617–31.View ArticleGoogle Scholar
- Drake H, Tullborg E-L, Page L. Distinguishing multiple events of fracture mineralisations related to far-field orogenic effects in Paleoproterozoic crystalline rocks, Simpevarp area, SE Sweden. Lithos. 2009;110:37–49.View ArticleGoogle Scholar
- Drake H, Tullborg EL. Detecting the near-surface redox front in crystalline bedrock using fracture mineral distribution, geochemistry and U-series disequilibrium. Appl Geochem. 2009;24:1023–39.View ArticleGoogle Scholar
- Andersen CR, James RE, Fru EC, Kennedy MA, Pedersen K. In situ ecological development of a bacteriogenic iron oxide-producing microbial community from a subsurface granitic rock environment. Geobiology. 2006;4:29–42.View ArticleGoogle Scholar
- Pedersen K, Ekendahl S. Assimilation of CO2 and introduced organic compounds by bacterial communities in groundwater from southeastern Sweden deep crystalline bedrock. Microb Ecol. 1992;23:1–14.View ArticlePubMedGoogle Scholar
- Hall-Stoodley L, Costerton JW, Stoodley P. Bacterial biofilms: from the natural environment to infectious diseases. Nat Rev Microbiol. 2004;2:95–108.View ArticlePubMedGoogle Scholar
- Edlund J, Rabe L, Bengtsson A, Hallbeck B, Eriksson L, Johansson J, Johansson L, Pedersen K. Understanding microbial reduction of sulphate to sulphide in deep Olkiluoto groundwater. Compilation and interpretation of three consecutive Sulphate Reduction Experiments (SURE) performed during 2010 - 2014. In Posiva Working Report 2016-48. pp. 1-164; 2016.
- Pedersen K. The MICROBE project. Achievements of a 10-year research program. SKB report R-13-49. Stockholm Sweden. pp. 1-40; 2013.
- Hubalek V, Wu X, Eiler A, Buck M, Heim C, Dopson M, Bertilsson S, Ionescu D. Connectivity to the surface determines diversity patterns in subsurface aquifers of the Fennoscandian Shield. ISME J. 2016;10:2447–58.View ArticlePubMedPubMed CentralGoogle Scholar
- Rocca JD, Hall EK, Lennon JT, Evans SE, Waldrop MP, Cotner JB, Nemergut DR, Graham EB, Wallenstein MD. Relationships between protein-encoding gene abundance and corresponding process are commonly assumed yet rarely observed. ISME J. 2015;9:1693–9.View ArticlePubMedGoogle Scholar
- Ekendahl S, Pedersen K. Carbon transformations by attached bacterial populations in granitic groundwater from deep crystalline bed-rock of the Stripa research mine. Microbiology. 1994;140:1565–73.View ArticlePubMedGoogle Scholar
- Holscher T, Bartels B, Lin YC, Gallegos-Monterrosa R, Price-Whelan A, Kolter R, Dietrich LE, Kovacs AT. Motility, chemotaxis and aerotaxis contribute to competitiveness during bacterial pellicle biofilm development. J Mol Biol. 2015;427:3695–708.View ArticlePubMedPubMed CentralGoogle Scholar
- Bogino PC, de las Mercedes Oliva M, Sorroche FG, Giordano W. The role of bacterial biofilms and surface components in plant-cacterial associations. Int J Mo Sci. 2013;14:15838–59.View ArticleGoogle Scholar
- Flemming H-C, Wingender J. The biofilm matrix. Nat Rev Microbiol. 2010;8:623–33.PubMedGoogle Scholar
- Perez-Rodriguez I, Bolognini M, Ricci J, Bini E, Vetriani C. From deep-sea volcanoes to human pathogens: a conserved quorum-sensing signal in Epsilonproteobacteria. ISME J. 2015;9:1222–34.View ArticlePubMedGoogle Scholar
- Kojima H, Fukui M. Sulfuricella denitrificans gen. nov., sp. nov., a sulfur-oxidizing autotroph isolated from a freshwater lake. Int J Syst Evol Bacteriol. 2010;60:2862–6.View ArticleGoogle Scholar
- Liu J, Wang Z, Belchik SM, Edwards MJ, Liu C, Kennedy DW, Merkley ED, Lipton MS, Butt JN, Richardson DJ, et al: Identification and characterization of MtoA: A decaheme c-type cytochrome of the neutrophilic Fe(II)-oxidizing bacterium Sideroxydans lithotrophicus ES-1. Front Microbiol. 2012;3. doi:10.3389/fmicb.2012.00037.
- Pedersen K, Bengtsson A, Edlund J, Rabe L, Hazen T, Chakraborty R, Goodwin L, Shapiro N. Complete genome sequence of the subsurface, mesophilic sulfate-reducing bacterium Desulfovibrio aespoeensis Aspo-2. Genome Announc. 2014;2. doi:10.1128/genomeA.00509-00514.
- Beller HR, Chain PS, Letain TE, Chakicherla A, Larimer FW, Richardson PM, Coleman MA, Wood AP, Kelly DP. The genome sequence of the obligately chemolithoautotrophic, facultatively anaerobic bacterium Thiobacillus denitrificans. J Bacteriol. 2006;188:1473–88.View ArticlePubMedPubMed CentralGoogle Scholar
- Nielsen ME, Fisk MR, Istok JD, Pedersen K. Microbial nitrate respiration of lactate at in situ conditions in ground water from a granitic aquifer situated 450 m underground. Geobiology. 2006;4:43–52.View ArticleGoogle Scholar
- Ström A, Andersson J, Skagius K, Winberg A. Site descriptive modelling during characterization for a geological repository for nuclear waste in Sweden. Appl Geochem. 2008;23:1747–60.View ArticleGoogle Scholar
- Pedersen K. Microbial life in deep granitic rock. FEMS Microbiol Rev. 1997;20:399–414.View ArticleGoogle Scholar
- Laaksoharju M, Gascoyne M, Gurban I. Understanding groundwater chemistry using mixing models. Appl Geochem. 2008;23:1921–40.View ArticleGoogle Scholar
- Smellie JAT, Laaksoharju M, Wikberg P. Äspö, SE Sweden—a natural groundwater-flow model derived from hydrogeological observations. J Hydrol. 1995;172:147–69.View ArticleGoogle Scholar
- Alakangas LJ, Mathurin FA, Faarinen M, Wallin B, Åström ME. Sampling and characterizing rare earth elements in groundwater in deep-lying fractures in granitoids under in situ high-pressure and low-redox conditions. Aquat Geochem. 2014;20:405–18.View ArticleGoogle Scholar
- Marteinsson VT, Runarsson A, Stefansson A, Thorsteinsson T, Johannesson T, Magnusson SH, Reynisson E, Einarsson B, Wade N, Morrison HG, Gaidos E. Microbial communities in the subglacial waters of the Vatnajokull ice cap, Iceland. ISME J. 2013;7:427–37.View ArticlePubMedGoogle Scholar
- Huber JA, Mark Welch DB, Morrison HG, Huse SM, Neal PR, Butterfield DA, Sogin ML. Microbial population structures in the deep marine biosphere. Science. 2007;318:97–100.View ArticlePubMedGoogle Scholar
- Lundin S, Stranneheim H, Pettersson E, Klevebring D, Lundeberg J. Increased throughput by parallelization of library preparation for massive sequencing. PLoS ONE. 2010;5:e10029.View ArticlePubMedPubMed CentralGoogle Scholar
- Borgström E, Lundin S, Lundeberg J. Large scale library generation for high throughput sequencing. PLoS ONE. 2011;6:e19119.View ArticlePubMedPubMed CentralGoogle Scholar
- Huse SM, Huber JA, Morrison HG, Sogin ML, Welch DM. Accuracy and quality of massively parallel DNA pyrosequencing. Genome Biol. 2007;8:R143.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC, Haas BJ, Clemente JC, Quince C, Knight R. UCHIME improves sensitivity and speed of chimera detection. Bioinformatics. 2011;27:2194–200.View ArticlePubMedPubMed CentralGoogle Scholar
- Edgar RC. Search and clustering orders of magnitude faster than BLAST. Bioinformatics. 2010;26:2460–1.View ArticlePubMedGoogle Scholar
- Joshi NA, Fass JN: Sickle. A sliding-window, adaptive, quality-based trimming tool for FastQ files. Version 1.33 edition. 2011. https://github.com/najoshi/sickle.
- Boisvert S, Raymond F, Godzaridis E, Laviolette F, Corbeil J. Ray Meta: scalable de novo metagenome assembly and profiling. Genome Biol. 2012;13:R122.View ArticlePubMedPubMed CentralGoogle Scholar
- Alneberg J, Bjarnason BS, de Bruijn I, Schirmer M, Quick J, Ijaz UZ, Lahti L, Loman NJ, Andersson AF, Quince C. Binning metagenomic contigs by coverage and composition. Nat Methods. 2014;11:1144–6.View ArticlePubMedGoogle Scholar
- Darling AE, Jospin G, Lowe E, Matsen FA, Bik HM, Eisen JA. PhyloSift: phylogenetic analysis of genomes and metagenomes. PeerJ. 2014;2:e243.View ArticlePubMedPubMed CentralGoogle Scholar
- Seemann T. Prokka: rapid prokaryotic genome annotation. Bioinformatics. 2014;30:2068–9.View ArticlePubMedGoogle Scholar
- Li D, Liu CM, Luo R, Sadakane K, Lam TW. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph. Bioinformatics. 2015;31:1674–6.View ArticlePubMedGoogle Scholar
- Quinlan AR, Hall IM. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics. 2010;26:841–2.View ArticlePubMedPubMed CentralGoogle Scholar
- Eydal HS, Pedersen K. Use of an ATP assay to determine viable microbial biomass in Fennoscandian Shield groundwater from depths of 3-1000 m. J Microbiol Methods. 2007;70:363–73.View ArticlePubMedGoogle Scholar