Main

High-precision sampling was conducted during cruises of the research vessel Urania equipped with the Modus–Scipack system (http://www.geo.unimib.it/BioDeep/Project.html; Fig. 1a). The vehicle Modus, connected by cable to the research vessel, held a second instrument, the Scipack, with a 10-m data transmission cable. The Scipack, consisting of a Rosette sampler equipped with a CTD (conductivity–temperature–depth probe) and a series of Niskin bottles, was connected to the Modus through the Sciskid, a module equipped with a pressure sensor for recording the pressure at which the Niskin bottles were closed (Fig. 1c). A camera on the Modus could provide an image of the Scipack entering the brine lake (Fig. 1b, and Supplementary Fig. S1). Immediately after sampling, the Modus–Scipack was raised, the Niskin bottles were retrieved and their contents were carefully fractionated on board ship by slowly recovering 0.5-litre, 1-litre or 2-litre fractions from the bottom tap. These were then immediately analysed for salinity (Fig. 1d). The reconstructed interface salinity profile was strongly positively correlated (r = 0.98, P < 0.001) with the CTD conductivity profile recorded in independent non-sampling casts (Fig. 2d), indicating that little or no mixing had occurred.

Figure 1: Procedures for sample recovery from the seawater–Bannock brine interface.
figure 1

a, Schematic drawing of the Modus–Scipack system illustrating the Modus vehicle holding the Scipack–Sciskid unit consisting of a Rosette sampler equipped with a CTD and a series of Niskin bottles (Scipack), and a pressure sensor (Sciskid) for recording the pressure at which the Niskin bottles were closed. b, Photographs of the Scipack Rosette sampler at 0.5-s intervals while going into and out of the seawater–Bannock brine interface obtained during a video survey from the Modus vehicle (from left to right). The disappearance of the Scipack is caused by the transmittance change resulting from the high salinity of the brine. c, This sampling system yielded the exact time and depth at which each Niskin bottle (represented by cylinders over the curve) was closed. In this way it was possible to determine at which point of the conductivity curve during the cast the bottles were closed. The circled numbers below selected bottles identify those Niskin bottles (shaded) that were subsampled on board ship for the characterization of the gradient. d, Example of conductivity measurements of the different water fractions sub-sampled from Niskin bottle no. 4 reported in c, using the conductivity probe from the Scipack.

Figure 2: Microbiological and geochemical profiling of the seawater–Bannock brine interface.
figure 2

a, Filled circles, DAPI microbial counts; open circles, Shannon–Weaver indices calculated from ARISA data. b, Filled circles, prokaryote 16S rRNA gene abundance; open squares, Archaea 16S rRNA gene abundance; open circles, ATP. c, Filled circles, SRR; open circles, methane production rate (MPR). d, Filled circles, conductivity; open circles, salinity. e, Filled circles, sodium; open circles, redox potential (Eh). f, Filled circles, dissolved manganese; open circles, nitrate. g, Filled circles, ammonium; open circles, sulphate. Shading indicates zones with relatively high biological activity. Conductivity data in c were obtained during a non-sampling cast by continuous measurements with a conductivity probe mounted on the Scipack system. All the data were from station AB27SCI except values for deep sea water and Bannock brine that were from station AB29SCI. Error bars indicate standard deviations. For sodium, dissolved manganese, nitrate, ammonium and sulphate, standard deviations were 2–5% of each value.

The interface halocline was about 2.5 m deep, in agreement with previous estimates that employed alternative sampling strategies1. Although biomass values fluctuated along the halocline, there were significantly greater numbers of microbial cells in the interface (about 106 cells ml-1) than in either the deep sea water or the underlying hypersaline brine, both of which had about 104 cells ml-1 (Fig. 2a). Comparison of archaeal and total prokaryotic 16S ribosomal RNA gene abundance revealed that Bacteria and Archaea were present at equivalent levels in the oxic seawater above the hypersaline brine, as reported previously for the deep sea8. In contrast, the interface samples contained 20.5 and 1.8 pg ml-1 bacterial and archaeal 16S rRNA gene, respectively (Fig. 2b), contrary to the notion that extreme environments are dominated by Archaea but consistent with the low level of Archaea-mediated methanogenesis measured in the upper 2 m of the halocline (Fig. 2c). The high concentration of bacteria in the interface indicates that it acts as a barrier to particulates that sink through the marine water column, leading to an accumulation of organic carbon, nutrients and reactive surfaces that support microbial growth9. In addition, methane originating in the hypersaline brine will rise and enter the interface, where it can be oxidized under more energetically favourable conditions.

Levels of ATP, an important indicator of metabolic activity10, were significantly higher in the interface than in the overlying sea water and the brine lake (Fig. 2b), as were extracellular aminopeptidase and alkaline phosphatase activities (Supplementary Table S1). At certain depths of the interface, sulphate reduction rates (SRRs) showed appreciably higher values than those of the deep oxic sea water11, but the highest SRRs were measured in the underlying anoxic brine (Fig. 2c).

Down the halocline, nitrate concentration decreased from 6.6 to 0.3 µM, ammonium increased from 5 to 3,450 µM (refs 2, 12), sulphate, the most abundant electron acceptor along the halocline, increased from 31 to 84 mM, manganese increased from 0.4 to 8.3 µM and redox potential decreased from 210 mV to less than zero, indicating a total depletion of oxygen in the lower part of the chemocline (Fig. 2e–g). Electron acceptors and donors potentially available for metabolic processes in the interface are thus variously distributed along the gradient. Whereas Na+ showed chemically conservative behaviour along the halocline, dissolved manganese, nitrate, ammonium and sulphate concentrations exhibited nonlinear slopes, demonstrating non-conservative behaviour presumably reflecting biologically mediated redox cycling13 (Supplementary Fig. S2).

Detailed vertical profiles of biological parameters showed two maxima in diversity at about 8% and 15–22% salinity, which corresponded to maxima in biomass, determined by total cell counts and prokaryote 16S rRNA gene abundance analysis (r = 0.93, P < 0.001; Fig. 2a), and in ATP concentration (r = 0.74, P < 0.01 and r = 0.64, P < 0.05; Fig. 2b). The decreases in these biological parameters in the 10–15% salinity range, and the low SRR in this region, can be explained by a depletion of nutrients and organic substrates from above and below, and/or a localized shortage of energetically favourable redox couplings (for example nitrate and oxygen are depleted at this point; Fig. 2e, f).

NaCl is well known for generating high osmotic potentials, reducing the mole fraction of water, and its kosmotropic (order-generating) activity stabilizes cellular macromolecules and stiffens membranes (Supplementary Fig. S3): such parameters influence the environmental windows of growth for any microbial species. In contrast with other hypersaline environments, such as coastal solar salt pans4, bacterial diversity in the Bannock interface was higher than that found in the adjacent low-salt habitat (the overlying deep seawater), as shown by Shannon–Weaver indices derived from amplified ribosomal intergenic spacer analysis (ARISA) of Bacteria at 18 different salinities (Fig. 2a). To identify the Bacteria present at different locations in the halocline, clone libraries of RT–PCR amplicons of 16S rRNA (complementary DNA) and of PCR amplicons of 16S rRNA genes were generated; 539 clones, including 179 cDNA clones representing the more metabolically active microbial fraction, were sequenced. This analysis confirmed the high bacterial phylogenetic diversity suggested by the Shannon–Weaver indices and identified four new candidate divisions designated Mediterranean Sea Brine Lake groups 3–6 (MSBL3–6) and a dominant candidate division, newly named MSBL2 (Supplementary Fig. S4). Candidate division MSBL2 is phylogenetically related to the candidate division SB1, previously identified in the interface of the Shaban Deep in the Red Sea14. MSBL2-related sequences were also detected in the interfaces of other hypersaline basins, but were absent from the underlying brines or the overlying sea water3, which indicates that bacteria belonging to the MSBL2 division are specifically adapted to the seawater–hypersaline brine interfaces.

The deeper part of the halocline (22.2–25.0% salinity) was inhabited by MSBL2, by Sphingobacteriaceae and by candidate divisions KB1 (refs 14, 15) and MSBL3–6 (Fig. 3a). The candidate division KB1 includes species related to Thermotogales that have been found in brine lakes and their sediments in the Red Sea14,15, and in other anaerobic hypersaline sediments16. Despite being less abundant than those of MSBL2, bacteria of MSBL3–6 candidate divisions were found only in the deeper part of the halocline, indicating that they are adapted to high salinity. Archaeal clones related to members of candidate divisions MSBL1, previously found to dominate in the brine3, and ANME-1, associated with anaerobic oxidation of methane17, were found in cDNA libraries generated from samples from the deeper part of the halocline, namely the anoxic zone where methane and sulphate coexist. Sulphate-reducing δ-Proteobacteria were abundant, comprising 6–15% of the clones in the libraries from all anoxic levels of the interface. Desulfobacteriaceae and Desulfobulbaceae were found mainly at lower and higher salinities, respectively (Fig. 3b), indicating specific adaptation of these two groups to the different prevailing environmental factors of salinity, redox potential, oxygen18,19, sulphate concentration and availability of organic substrates.

Figure 3: Distribution of prokaryotic diversity at different salinities along the seawater–Bannock brine interface.
figure 3

a, Percentage distribution of 16S rRNA gene clones belonging to different bacterial phylogenetic groups. New divisions MSBL3–6 are included under ‘new taxa’. b, Percentage distribution of different families of the δ-Proteobacteria.

Eighty-four bacterial isolates were cultured from the Bannock interface; on the basis of partial 16S rRNA sequences they were assigned to the Firmicutes, Bacteroidetes, α-Proteobacteria, γ-Proteobacteria and ɛ-Proteobacteria (Supplementary Table S2). The isolates belonged to divisions widely represented in the rRNA libraries; however, they did not correspond to any specific sequence found in the libraries (Supplementary Fig. S4). Most isolates were moderately halophilic facultative anaerobes, able to grow aerobically or by fermentation or denitrification, and therefore well adapted to the seawater–brine interface. In particular, a Halothiobacillus species was able to aerobically oxidize thiosulphate, with CO2 as sole carbon source, over a NaCl range of 0.5–23%. In addition, two obligately anaerobic, fermentative, halophilic isolates had 99.8% 16S rRNA sequence similarity to Halanaerobium sp. KT-2/3-3, which was isolated from a similar seawater–brine interface in the Kebrit deep20.

Two obligately anaerobic strains, one from the Bacteroidetes and one from the ɛ-Proteobacteria from the same enrichment, had less than 92% 16S rRNA sequence similarity to known organisms. Both were moderate halophiles, growing between 2% and 12% NaCl and between 0.1% and 9.0% NaCl, respectively. The Bacteroidetes isolate fermented a variety of sugars and biopolymers as sole sources of carbon and energy, whereas the ɛ-Proteobacteria member fermented a variety of organic acids, and used nitrate, nitrite, thiosulphate, trimethylamine N-oxide, dimethylsulphoxide and elemental sulphur as terminal electron acceptors with formate and acetate as electron donors. It can therefore be envisaged that these organisms, which are highly suited to life in Bannock interface, would benefit mutually from the degradation of organic polymers, which are likely to accumulate at the interface9, coupled to the respiration of sulphur compounds or nitrate.

Sulphate reducers, despite being abundantly represented in clone libraries, escaped cultivation. Nevertheless, analysis of our isolates confirmed that a range of biological processes—such as fermentation and the oxidation of reduced sulphur species and organic matter coupled to a variety of terminal electron acceptors—are likely to occur in the interface.

In conclusion, the Bannock brine lake–seawater interface is revealed here to be a deep-sea microbial oasis, consisting of a biomass-dense, metabolically active microbial ecosystem that receives methane from the hypersaline brine below and organic particulates sedimenting through the water column from above, recycles this organic matter and interconverts diverse ionic species, which are used and reused as electron acceptors and donors. Although functionally analogous to the deep-sea floor, which also collects and recycles sinking particulates, the interfaces of deep hypersaline basin lakes differ significantly in the stress imposed by their steep halocline, and consequently have selected for stratified microbial communities containing new organisms, including divisions of bacteria not previously described. Recent work detailing the retrieval of novel enzymes from the analogous Urania basin interface21 would suggest that this new phylogenetic diversity reflects new functional activity. Such biomass-rich environments are stable over hundreds of years, constituting localized sources of organic matter that may contribute to the geological record, for example in the form of sapropels22,23.

Methods

Sampling of seawater–brine interface in Bannock basin

Sampling of the Bannock basin was conducted from the research vessel Urania at location 34° 21.640′ N, 20° 02.260′ E in 2001, location 34° 17.949′ N, 20° 00.985′ E (station BD29CT) in 2002, and at locations 34° 17.488′ N, 20° 00.692′ E (AB27SCI) and 34° 17.397′ N, 20° 00.709′ E (AB29SCI) in 2003. The halocline water fractions, sampled with the Modus–Scipack system, were analysed as reported below for microbial activity, microbial abundance and diversity, physicochemical properties, and isolation of microbial strains.

Activity measurements

Aminopeptidase and phosphatase activities were determined with the fluorogenic substrates l-leucine-7-amino-4-methylcoumarin and 4-methylumbelliferyl phosphate with the use of the multi-concentration kinetic method at concentrations ranging from 0.05 to 10 µM (ref. 24). Methane production rates were determined by measuring the production of methane with a gas chromatograph equipped with a flame ionization detector3. Sulphate reduction rates were determined by measuring the [35S]sulphide production from radiolabelled sulphate (1–2 µCi [35S]sulphate) with standard methods. ATP was measured on triplicate 10-ml samples filtered through 0.22-µm pore-size filters. ATP was extracted and measured directly on the filter with the luciferin–luciferase-based biomass test kit (Promicol), with a luminometer, following the instructions of the manufacturer. Relative luminescence units were converted to ATP concentrations with the use of a standard ATP curve.

DNA/RNA-based analysis and direct cell count

Samples were filtered onto sterile 0.22-µm pore-size filters that were stored at -20 °C in 2 ml of sterile lysis buffer (EDTA 40 mM, Tris-HCl 50 mM pH 5.8, sucrose 0.75 M). DNA was extracted from each filter as described previously3. Extracted DNA was quantified by comparative agarose-gel electrophoresis with known quantities of lambda phage DNA. 4,6-Diamidino-2-phenylindole (DAPI) staining was performed as described previously3. Real-time PCR experiments for the quantification of total Prokaryotes and Archaea were performed by Taqman assays as described elsewhere25.

From each DNA fraction along the salinity gradient, ARISAs were performed in accordance with an established protocol26. As a parameter for the structural diversity of a bacterial community a Shannon–Weaver index27 (H) was calculated for representative ARISA profiles by using the function H = -Σ(ni/N)ln(ni/N), where ni is the height of a peak and N is the sum of all peak heights in an ARISA profile. The construction and sequencing of 16S rRNA gene libraries from fractions with different salinities were performed with a previously described procedure3. Sequences were aligned with the use of ARB software28, and operational taxonomic unit (OTU) distribution at different salinities was calculated. Good coverage of the dominant OTU population was confirmed with rarefaction analysis of the clone libraries. Sequences of 16S rRNA genes showing more than 97% homology were considered to belong to the same OTU. For RNA extraction and reverse transcription the RNA/DNA Mini Isolation kit (Qiagen) and Superscript II reverse transcriptase (Invitrogen) were used.

Geochemical analyses and salinity-related measurements

The redox potential in each halocline fraction was measured on board ship with a portable Eh meter29. Major elements in brines were analysed by inductively coupled plasma atomic emission spectrometry at Geosciences Utrecht after dilution to 3.5% salinity; nitrate and ammonia were analysed by AutoAnalyser at NIOZ (Nederlands Instituut voor Onderzoek der Zee), Texel; Mn was measured at RIVM (Rijks Instituut voor Volksgezondheid en Milieu), Bilthoven, with the use of high-resolution inductively coupled plasma mass spectrometry.

Bacterial isolation and characterization

A collection of 84 isolates was obtained by inoculating anaerobic samples of the whole interface into a wide variety of media (Supplementary Information). Isolates were tentatively identified by sequencing the 16S rRNA gene, and physiological features tested included growth at different salinities and with different carbon and energy sources and terminal electron acceptors.