# Microbial diversity in a submarine carbonate edifice from the serpentinizing hydrothermal system of the Prony Bay (New Caledonia) over a 6-year period

Anne Postec, Marianne Quéméneur, Méline Bes, Nan Mei, Fatma Ben Aissa-Fennira, Claude Payri, Bernard Pelletier, Christophe Monnin, Linda Guentas-Dombrowsky, Bernard Ollivier, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Anne Postec, Marianne Quéméneur, Méline Bes, Nan Mei, Fatma Ben Aissa-Fennira, et al. Mi-crobial diversity in a submarine carbonate edifice from the serpentinizing hydrothermal system of the Prony Bay (New Caledonia) over a 6-year period. Frontiers in Microbiology, 2015, 6, pp.857. 10.3389/fmicb.2015.00857 . hal-01221183

HAL Id: hal-01221183

https://hal.sorbonne-universite.fr/hal-01221183

Submitted on 27 Oct 2015

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scien-tific research documents, whether they are pub-lished or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
OPEN ACCESS

**Edited by:**
Mark Alexander Lever,
Swiss Federal Institute of Technology
Zurich, Switzerland

**Reviewed by:**
Matthew Schrenk,
Michigan State University, USA
Julie A. Huber,
Marine Biological Laboratory, USA

**\*Correspondence:**
Anne Postec and Gaël Erauso,
Mediterranean Institute of
Oceanography, Aix-Marseille
Université, équipe Microbiologie
Environnementale et Biotechnologie,
Case 905, 163 Avenue de Luminy,
13288 Marseille, France
anne.postec@mio.osupytheas.fr;
gael.erauso@mio.osupytheas.fr

**Specialty section:**
This article was submitted to
Extreme Microbiology,
a section of the journal
Frontiers in Microbiology

**Received:** 26 May 2015

**Accepted:** 06 August 2015

**Published:** 27 August 2015

**Citation:**
Postec A, Quéméneur M, Bes M, Mei N, Benaïssa F, Payri C, Pelletier B, Monnin C, Guentas-Dombrowsky L, Ollivier B, Gérard E, Pisapia C, Gérard M, Ménez B and Erauso G (2015) Microbial diversity in a submarine carbonate edifice from the serpentinizing hydrothermal system of the Prony Bay (New Caledonia) over a 6-year period. *Front. Microbiol.* 6:857. doi: 10.3389/fmicb.2015.00857

# Microbial diversity in a submarine carbonate edifice from the serpentinizing hydrothermal system of the Prony Bay (New Caledonia) over a 6-year period

Anne Postec<fnref n="1" />*, Marianne Quéméneur<fnref n="1" />, Méline Bes<fnref n="1" />, Nan Mei<fnref n="1" />, Fatma Benaïssa<fnref n="1" />, Claude Payri², Bernard Pelletier², Christophe Monnin³, Linda Guentas-Dombrowsky<fnref n="1" />,², Bernard Ollivier<fnref n="1" />, Emmanuelle Gérard⁴, Céline Pisapia⁴, Martine Gérard⁵, Bénédicte Ménez⁴ and Gaël Erauso<fnref n="1" />*

Active carbonate chimneys from the shallow marine serpentinizing Prony Hydrothermal Field were sampled 3 times over a 6 years period at site ST09. Archaeal and bacterial communities composition was investigated using PCR-based methods (clone libraries, Denaturating Gel Gradient Electrophoresis, quantitative PCR) targeting 16S rRNA genes, methyl coenzyme M reductase A and dissimilatory sulfite reductase subunit B genes. *Methanosarcinales* (Euryarchaeota) and *Thaumarchaea* were the main archaeal members. The *Methanosarcinales*, also observed by epifluorescent microscopy and FISH, consisted of two phylotypes that were previously solely detected in two other serpentinizing ecosystems (The Cedars and Lost City Hydrothermal Field). Surprisingly, members of the hyperthermophilic order *Thermococcales* were also found which may indicate the presence of a hot subsurface biosphere. The bacterial community mainly consisted of *Firmicutes*, *Chloroflexi*, *Alpha*-, *Gamma*-, *Beta*-, and *Delta*-proteobacteria and of the candidate division NPL-UPA2. Members of these taxa were consistently found each year and may therefore represent a stable core of the indigenous bacterial community of the PHF chimneys. *Firmicutes* isolates representing new bacterial taxa were obtained by cultivation under anaerobic conditions. Our study revealed diverse microbial communities in PHF ST09 related to methane and sulfur compounds that share common populations with other terrestrial or submarine serpentinizing ecosystems.

**Keywords:** hydrothermal, serpentinization, alkaline, Prony, microbial communities

<footnote n="1">¹ Aix-Marseille Université, Centre National de la Recherche Scientifique/INSU, Université de Toulon, IRD, Mediterranean Institute of Oceanography, UM 110, Marseille, France, ² Institut pour la Recherche et le Développement Centre de Nouméa, Nouméa-Nouvelle-Calédonie, France, ³ Géosciences Environnement Toulouse, Université de Toulouse/Centre National de la Recherche Scientifique/IRD, Toulouse, France, ⁴ Institut de Physique du Globe de Paris, Sorbonne Paris Cité, Université Paris Diderot, Centre National de la Recherche Scientifique, UMR7154, Paris, France, ⁵ Institut de Minéralogie et de Physique des Milieux Condensés, Université Pierre et Marie Curie, Paris, France</footnote>
Introduction

Serpentinization is the process of aqueous alteration of ultramafic rocks that releases hydrogen (H₂), methane (CH₄), small organic compounds, and generates alkaline fluids with pH values commonly above pH 10, at reducing conditions (McCollom, 2007; Schrenk et al., 2013). This last decade, low temperature serpentinizing environments have attracted much scientific interest because they are considered as analogs of those from which life may have arisen on the early Earth. Indeed, their particular geochemical settings may have led, a long time ago, to the transition from the abiotic reactions to biologically mediated reactions probably similar to some present day microbial metabolisms (Schulte et al., 2006; Russell et al., 2010). However, there are still only few studies of the microbiology of modern serpentinizing systems and most concerns continental sites such as the alkaline groundwater (pH 11.4) of Cabeço de Vide Aquifer (CVA, Portugal) (Tiago and Veríssimo, 2013), the ophiolitic complexes of Leka (Norway) (Daae et al., 2013), Tablelands (Newfoundland) (Brazelton et al., 2012) or The Cedars springs (California) (Suzuki et al., 2013). All these sites are readily accessible for sampling and *in situ* analyses while exploration of active serpentinization in oceanic settings is more challenging (Schrenk et al., 2013).

The Lost City Hydrothermal Field (LCHF) located at a depth of about 800 m below sea level (mbls) on the top of the Atlantis massif at 30°N near the Mid-Atlantic Ridge is to date the most extensively studied submarine site driven by active serpentinization. There, high carbonate chimneys vent clear alkaline fluids (pH 9-11) at moderately high temperatures (40-90°C) (Kelley et al., 2001, 2005; Proskurowski et al., 2008). Molecular surveys have revealed a surprisingly low archaeal diversity inside these carbonate chimneys dominated by a single phylotype of *Methanosarcinales*, the LCMS (Lost City *MethanoSarcinales*), potentially involved in methane production or oxidation along with archaeal anaerobic methane oxidizers of the ANME-1 group (Schrenk et al., 2004; Brazelton et al., 2006, 2011). Bacterial sulfate reducers (namely *Desulfotomaculum*) were identified while microbial H₂ production by anaerobic *Clostridiales* was evidenced by metagenomics (Brazelton et al., 2012).

Until recently, LCHF was thought to be drastically different from any other ecosystem previously discovered. However, recent geochemical and microbiological studies of the hydrothermal system of the Prony Bay (New Caledonia, SW Pacific) have revealed important similarities with the LCHF ecosystem. PHF is a marine hydrothermal system located at shallow depths in the Prony bay, on a serpentinized peridotite nappe (Figure 1A). In 2004 and 2005, bathymetric mapping and scuba dives during scientific cruises on the *R/V Alis* revealed that, besides a impressive submarine edifice called “L’Aiguille de Prony,” the floor of the Prony Bay hosted several 2–10 m high domes and pinnacles at water depths between 30 and 50 m (Pelletier et al., 2006). A recent geochemical study has shown that the hydrothermal hyperalkaline (pH 11) end member has a meteoric origin (run off waters) and is enriched in H₂ and CH₄ (Monnin et al., 2014). The first study of the microbial community

of PHF showed that archaeal diversity was low and dominated by *Methanosarcinales* while members of *Chloroflexi*, *Firmicutes*, and *Proteobacteria* phyla were the most abundant bacteria. It also showed variation in microbial community structure, abundance, and diversity depending on the location of the sampling site in the Prony Bay, probably related to change in environmental factors such as water depth and maturity (age) of the carbonate chimneys (Quéméneur et al., 2014).

Very few culture-dependent approaches have been reported to characterize the bacterial populations from serpentinite environments and are so far restricted to alkaline groundwaters of Oman, CVA and the Cedars (Bath et al., 1987; Schrenk et al., 2013; Tiago and Veríssimo, 2013; Suzuki et al., 2014). In CVA, only the aerobic heterotrophic fraction of the cultivable microbial community was targeted and was found to be mainly constituted of Gram-positive bacteria belonging to *Actinobacteria*, *Bacillus* and *Staphylococcus* (Tiago et al., 2004). Recently, aerobic hydrogen-oxidizers growing optimally at pH 11 were isolated from the continental serpentinizing springs at The Cedars and proposed as representing a new genus, "Serpentinomonas," within the *Betaproteobacteria* (Suzuki et al., 2014).

Here, we describe the microbial communities inhabiting three active PHF chimneys collected over a 6-year period from one of the deepest hydrothermal sites not included in the previous study (Quéméneur et al., 2014) using a combination of PCR-based methods (DGGE, qPCR and sequencing of clone libraries of 16S rRNA gene and functional genes), microscopy (CSLM and FISH) and anaerobic cultivation methods.

Materials and Methods

Sampling Site and Sample Collection

The sampling site (22°21.653S, 166°52.777E, site ST09 on **Figure 1A**) is an active hydrothermal area located East of the Central Bay of Prony, in the south of New Caledonia. The general geological context and the composition of venting fluids and gasses are given by Monnin et al. (2014) (Supplementary Table 1).

The upper parts of three active chimneys, venting alkaline fluids, were collected by SCUBA diving on the western part of a hydrothermal carbonated edifice with a dome shape based in a depression at a water depth of about 50 m. The first chimney collected in June 2005 during the first exploratory campaign on *Alis* RV (Pelletier et al., 2006) at 49 mbls was 62 cm high with 22 cm maximal diameter; the second one collected in November 2010 at 47 mbls was about 30 cm high and 14 cm maximal diameter and the third one collected in November 2011 during the HYDROPRONY campaign on *Alis* RV (Monnin et al., 2014) at 45 mbls was 22 cm height with 15 cm maximal diameter (Figures 1B–H). The three sampled chimneys possessed very fragile crumbly white needles at their summit (Figure 1G), testifying to active fluid discharge. They were more or less encrusted with various invertebrates such as ascidia, coral and marine sponges (Figure 1D) attesting for their long-lasting existence and advanced stage. They were mainly composed of aragonite, brucite, calcite, and magnesium carbonates as detailed elsewhere (Pisapia et al., 2013). The chimneys collected in 2005
<figure><img src="image_3.png" /><figcaption>**FIGURE 1 |** Sampling site and chimney sample photographs.

(A) Map of the Bay of Prony in South New Caledonia indicating the hydrothermal sampling site ST09 located at 45–49 mbs in the East side of the Bay (after Maurizot and Vendé-Leclerc, 2009). (B) *In situ* photographs of the main edifice at ST09 around which chimneys were collected in 2011. (C) Chimney 2011 subsampling. (D) Top of the</figcaption></figure>

chimney collected in 2010 encrusted with coral ascidia and sponges.

(E) Longitudinal fracture of the 2010 top chimney showing a mucilaginous dark-pink biofilm-like texture around the conduit. (F) Chimney collected in 2005 (Pelletier et al., 2006), (G) white top of the 2005 chimney, (H) radial zonation in chimney 2005 with dark gray to green intermediary parts. Scale bars: 5 cm.

and 2011 were recovered aboard the *Alis* RV. In 2005, the chimney was quickly frozen aboard at –20°C, sent to the IRD research center in Nouméa, sawed in still frozen slices about 5 cm thick and stored at –80°C. The chimney collected in 2011 was directly sawed in sections and subsampled aboard the *Alis* RV, for molecular biology and culture analyses. In 2010, the chimney recovered aboard an outboard boat was kept on ice until back (within 4–5 h) at the IRD in Nouméa where it was sawed and subsampled as described below. Cross-sections of active chimneys showed a light beige highly porous material with sometimes friable beige to pink mucilaginous material visible around the central conduit, suggesting the presence of a bacterial biofilm (**Figures 1E,H**). Subsamples corresponding to internal parts of the chimneys were carefully collected trying to avoid contamination by the outer parts as much as possible and stored in sterile 50 mL Falcon tubes at –80°C for later DNA extraction. For fluorescent *in situ* hybridization (FISH) experiments, samples were fixed for 3 h at 4°C in a solution of phosphate-buffered saline pH 7.2 (PBS; 130 mM NaCl, 7 mM Na₂HPO₄, 30 mM NaH₂PO₄) containing 2% (vol/vol) formaldehyde (methanol free, Ultra Pure; Polysciences), washed twice with PBS and

resuspended in sterile Falcon tubes with 50% ethanol/PBS. Samples fixed in this buffer were then kept at –20°C.

Samples used for anaerobic cultivation experiments were ground aseptically with a pestle in a mortar under a continuous flux of N₂ gas to preserve anaerobiosis and distributed into sterile glass vials which were then tightly closed with butyl rubber stoppers and crimped with aluminum seals. The headspace in the vials was filled with 100% N₂ with a slight overpressure (100 kPa) to prevent oxygen entrance. They were stored at 4°C until use.

## DNA Extraction

Crushed chimney samples (0.3 g) were suspended in 1 mL TE-Na-1X lysis buffer (100 mM Tris, 50 mM EDTA, 100 mM NaCl, 1% SDS (w/v), 1% (w/v) Sarkosyl, pH 8.0). Cells were disrupted by bead beating using Lysing Matrix E (MP, Biomedicals). After digestion with proteinase K (1 mg.mL⁻¹ final) for 1 h at 55°C, DNA was extracted twice with a volume of phenol/chloroform/isoamyl alcohol (25/24/1; pH 8.0), then with a volume of chloroform. Aqueous and organic phases were separated by centrifugation at 13,000 g for 15 min at 4°C. DNA was precipitated with 0.7 volume of isopropanol, pelleted by
centrifugation at maximum speed for 15 min, then air-dried and finally resuspended in 30 µL of TE 1X (10 mM Tris, HCl, 2 mM EDTA, pH 8.0) buffer. The concentration of DNA extracts was measured using the Qubit® fluorometer (Invitrogen).

## PCR-amplification of 16S rRNA and Functional Genes

PCR set-up was carried out under aseptic conditions using autoclaved and/or UV-treated plasticware and pipettes. The PCR mixtures contained (total 50 µL) 0.4 pmol·µL⁻¹ of primers, 1 µL of chimney DNA template (1/10 dilution), 1X reaction buffer (GoTaq® Hot Start Green Master Mix, Promega). The primer sets used for 16S rRNA gene amplification (Supplementary Table 2) were (i) for cloning purposes: 27F/907R for *Bacteria* and 109F/958R for *Archaea*, (ii) for DGGE purpose: 341F-GC/907R for *Bacteria* and 344F-GC/958R for *Archaea*. PCR reactions were carried out using a T100 thermal cycler (BioRad) with the following conditions: initial denaturation at 95°C for 2 min followed by 30 cycles, each cycle consisting of denaturation at 94°C for 30 s, primers annealing at 50°C for 30 s and extension at 72°C for 90 s, followed by a final extension step of 5 min at 72°C. The *mcrA* and *dsrB* genes were amplified as described above using the primer sets MLF/MLR and DSRp2060F/DSR4R respectively, except that 35 cycles were performed and primer hybridization was carried out at 55°C.

## Cloning and Sequencing

Each gene library (archaeal and bacterial 16S rRNA, *mcrA*, and *dsrB*) was constructed from three to five independent PCR products that were pooled and purified using a Qiaquick PCR purification kit (Qiagen) before ligation in pGEM-T Easy vector (Promega) overnight at 4°C. *Escherichia coli* JM109 competent cells (Promega) were used for transformation by heat-shock using a standard protocol (Promega kit manual). White colonies were used to inoculate 96-wells microplates filled with LB plus ampicillin (100 µg·mL⁻¹) agar. Plasmid Miniprep and Sanger sequencing were operated by GATC Biotech AG, Germany.

## Phylogenetic and Diversity Analyses of Chimney Communities

Partial 16S rRNA gene sequences (average size around 900 bp) were checked with the online Bellerophon program version 3 at Greengenes website (http://greengenes.lbl.gov/) to exclude the chimeric artifacts. Each sequence was compared to 16S rRNA gene sequences of the NCBI GenBank nr database by using BLAST (http://www.ncbi.nlm.nih.gov/BLAST/) (Altschul et al., 1990). 16S rRNA sequences were also phylogenetically classified according to the Naive Bayesian rRNA Classifier (Version 1.0) of the Ribosomal Database Project II (http://rdp.cme.msu.edu/). Non-chimeric 16S rRNA gene sequences and translated *dsrB* and *mcrA* sequences were aligned using Muscle (Edgar, 2004) implemented in the MEGA6 software (Tamura et al., 2013). Alignments were edited manually and regions of ambiguous alignment removed. The MOTHUR program (Schloss et al., 2009) was used to group sequences into operational taxonomic units (OTUs) with a threshold of 97% sequences identity, commonly accepted as the minimum similarity level within a

prokaryote species (Gevers et al., 2005). The input files for MOTHUR were distance matrices generated by MEGA6 using a Maximum Composite Likelihood model for 16S rRNA genes. The coverage (C) of each clone library was calculated according to the equation: $C = 1-(n/N)$ where n is the number of OTU and N is the total number of clones in the library (Good, 1953). Coleman rarefaction curves, Chao1 values and Shannon indices were obtained from MOTHUR. At least one sequence representing each OTU, called phylotype, here was further aligned using Muscle with related sequences retrieved from NCBI databases using BLAST tools. The evolutionary history was inferred using both the Neighbor-Joining method with the Maximum Composite Likelihood model (Saitou and Nei, 1987; Tamura et al., 2004) and the Maximum Likelihood method based on the General Time Reversible model (Nei and Kumar, 2000) using MEGA6 software. The evolutionary distances were computed using the Kimura 2-parameters method (Kimura, 1980) for *dsrB* and *mcrA* translated sequences. The data were bootstrapped 1000 times to assess support for nodes. All positions containing gaps and missing data were eliminated.

## Quantitative PCR

Total bacterial and archaeal 16S rRNA genes were quantified using the primers sets 341F/518R and 344F/519R, respectively. The primer sets used for quantifying *mcrA* and *dsrB* genes were ME3MF:ME3MF-e (250:1)/ME2r' and DSRp2060F/DSR4R, respectively (see list of primers in Supplementary Table 2). The detailed procedure was described previously (Quéméneur et al., 2014).

## Epifluorescence Microscopy, FISH and CSLM

Direct microscopic observations were carried out on the sub-samples kept anaerobically at 4°C using a phase-contrast Nikon Eclipse 6500 epifluorescence microscope with an oil-immersion objective Nikon Plan Fluor 100X, by simply crushing a drop of chimney slurry between slide and cover slides. The blue autofluorescence of methanogen cells due to their co-factor F420 content was recorded after excitation at 420 nm. Fluorescence *in situ* hybridization was carried out at 46°C for 2–4 h using 5 ng·µL⁻¹ of probe in 0.9 M NaCl, 20 mM Tris-HCl pH 8, 0.01% SDS containing 35% (vol/vol) formamide. As experiments were targeting *Euryarchaeota*, the probe used was the FITC-labeled S-K-Eury-0498-a-A-14 (EURY498) (Supplementary Table 2). Samples were then washed for 15 min at 48°C in 20 mM Tris pH 8, 5 mM EDTA, 0.01% SDS, 0.084 mM NaCl, soaked in cold water for a few seconds and air dried. Samples were then stained with DAPI or Syto®9, (blue and green fluorescent nucleic acid stains at concentration of 1 µg·mL⁻¹ (Invitrogen) for 1 min, then washed for a few seconds in cold water and left to dry. Hybridized cells, covered by Citifluor AF3 (glycerol mountant solution from Emgrid), were examined at the Institut de Physique du Globe de Paris using an Olympus FluoView FV1000 confocal laser scanning microscope with a spectral resolution of 2 nm and a spatial resolution of 0.2 mm. The microscope is equipped with a 405-nm laser diode and multiline argon (458, 488, and 515 nm), helium–neon–green (543 nm) and helium–neon–red (633 nm) lasers. Observations were carried out
with an oil immersion UPLSAPO 60XO objective (Olympus; x60 magnification, numerical aperture = 1.35). Fluorescence image stacks were obtained with concomitant excitation at wavelengths of 405 and 488 nm by collecting the emitted fluorescence between 425-475 and 500-530 nm, respectively. The fluorescence spectra of the F420 cofactor, implied in the energy metabolism of methanogens (Doddema and Vogels, 1978) were recorded between 405 and 800 nm for a wavelength excitation of 405 nm, with band width analyses of 10 nm and steps of 5 nm. Three-dimensional images were acquired, visualized and processed using the F10-ASW FLUOVIEW software (Olympus).

## Denaturing Gradient Gel Electrophoresis (DGGE)

Total DNA of the 2005, 2010, and 2011 chimney samples were used as templates for the PCR amplification of the 16S rRNA gene V3 segments, which were profiled using DGGE. Aliquots of 25 µL of PCR products per lane were loaded on gels using a D-Code universal mutation detection apparatus (Bio-Rad, CA, USA). Electrophoresis was conducted in 6% (w/v) acrylamide gel with 40-70% denaturing gradient in 1 × TAE buffer (20 mM Tris acetate, 10 mM sodium acetate, and 0.5 mM EDTA) at a constant voltage of 200 V at 60°C for approximately 5 h. The gel was stained with SybRGold, de-stained with deionized water, and finally scanned using the Gel Doc image acquisition system (Bio-Rad, CA, USA). The predominant DGGE bands were then excised from the gels, eluted overnight in 50 µL sterile deionized water, and used as templates to re-amplify the 16S rRNA gene-V3 segments by using the same PCR primer pairs but without the GC clamp. The PCR products were then sequenced at GATC Biotech AG, Germany. Taxonomic affiliation of the PCR product sequence was based on the results of BLASTN against the nr database at NCBI. Local BLASTN was also used to compare these DGGE sequences (~600 pb) with our database containing all 16S rRNA sequences from both archaeal and bacterial clone libraries from PHF.

## Estimating the Growth Temperature of Archaea

The method is based on the strong correlation between the growth temperatures of cultured prokaryotes and the guanine-plus-cytosine contents ($P_{GC}$) of their 16S rRNA sequences, thermophilic and hyperthermophilic strains having higher $P_{GC}$ than those of psychrophilic and mesophilic isolates (Galtier et al., 1999; Khachane et al., 2005; Kimura et al., 2006, 2007; Wang et al., 2006). Kimura et al. (2010, 2013) proposed linear regression equations to infer minimum ($T_{min}$), optimum ($T_{opt}$) and maximum ($T_{max}$) growth temperatures of not-yet cultured prokaryotes and of microbial communities based on the $P_{GC}$ of their 16S rRNA sequences and they showed that these estimates correlate well with the *in situ* temperatures (Kimura et al., 2010, 2013). We therefore used their equations to estimate the growth temperatures of the main archaeal phylotypes encountered in the PHF chimneys. These calculations were possible since the sequences of our clones encompassed the same internal region of the 16S rRNA as that used by Kimura et al. (2013), which is between the *Archaea*-specific primers 109F and 915R. All the 16S rRNA sequences of our clones, of reference strains and from other studies, were aligned using MEGA 6 and manually trimmed

to be of the same length (800 bp) as those used in Kimura et al. (2013) linear regressions. Their GC content was calculated on line at http://www.endmemo.com/bio/gc.php and their respective $T_{min}$, $T_{opt}$, and $T_{max}$ were calculated using Kimura's equations in a spreadsheet.

## Anaerobic Cultivation Experiments

Various culture media were designed targeting enrichment of: (i) methanogens, using either acetate, formate, methanol, mono, and tri-methylamine (MMA and TMA), all at 20 mM, or $H_2/CO_2$ (80/20, 200 kPa) as substrates for methanogenesis, (ii) sulfate and sulfur-reducers using either yeast extract (5 g·L⁻¹), glucose (5 g·L⁻¹), acetate (20 mM), lactate (20 mM) or $H_2/CO_2$ (80/20, 200 kPa) as electron donors with either sulfate, thiosulfate (both 20 mM) or elemental sulfur (3 g·L⁻¹) as electron acceptors. The basal medium (BM) used for these cultures contained (g·L⁻¹): $NH_4Cl$ (1.0), $K_2HPO_4$ (0.3), $KH_2PO_4$ (0.3), $KCl$ (0.1), $CaCl_2·2H_2O$ (0.1), $NaCl$ (2, 10, 20, or 30), yeast extract (Difco) (0.1), cysteine hydrochloride (0.5) and 1 mL trace mineral element solution (Widdel et al., 1983). Media were boiled and distributed under strict anaerobic conditions in Hungate tubes or penicillin vials, then sterilized by autoclaving. Prior to inoculation, the pH was adjusted to 9 to 10.5, with a solution of 8% (w/v) $Na_2CO_3$. Hungate tubes (5 mL medium) were then supplemented with 0.1 mL of 3% (w/v) $Na_2S·9H_2O$, 0.1 mL of 15% (w/v) $MgCl_2·6H_2O$ and aliquots of stock solutions of the substrates (20 mM, otherwise 5 g·L⁻¹ for complex organic matter). Enrichment cultures were inoculated with about 5% (v/v) of anaerobic chimney slurry and incubated at 30, 37, 50, and 60°C. They were regularly checked for growth by visual observation of fresh aliquots with a light microscope or by measurement of the gas production ($H_2$, $CO_2$, $CH_4$, or $H_2S$) by gas chromatography as described previously (Mei et al., 2014). Strains were isolated by picking up colonies formed on 1.6% agar solid medium using Hungate's roll tubes method (Hungate, 1969). They were deemed pure after repeating the roll tubes step at least twice and finally transferred into liquid medium for stock cultures.

## Nucleotide Sequence Accession Numbers

All the sequences obtained in this study were deposited at Genbank under accession number from KJ149101 to KJ149254 and KJ159188 to KJ159206.

## Results and Discussion

### Microbial Community Abundance, Structure and Variation over the 6 Year Period

Real-time quantitative PCR was used to assess the relative abundance of archaea and bacteria within three PHF chimneys. The abundance of bacterial 16S rRNA genes was more than one order of magnitude larger than that of archaeal 16S rRNA genes each year (Table 1). The bacterial 16S rRNA gene abundance ranged from 4.0 ± 0.9 × 10⁷ (2010) to 9.1 ± 0.1 × 10⁷ (2005) whereas the archaeal 16S rRNA gene abundance ranged from 3.3 ± 0.8 × 10⁶ (2005) to 4.5 ± 0.9 × 10⁶ (2010) copies per g of wet weight (gww). No significant difference was observed
<table><caption>TABLE 1 | Gene abundance and cell number estimates from quantitative PCR on bacterial and archaeal 16S rRNA genes and on *dsrB* and *mcrA* genes and diversity indices for their respective clone libraries.</caption><thead><tr><th>Chimney<br>sampling year</th><th>Library</th><th>Number of gene<br>copies per g (ww)<br>of chimney</th><th>Number of cells<br>estimated per g (ww)<br>of chimney</th><th>OTUs <sup>a</sup>/Clones</th><th>Coverage</th><th>Chao1<sup>b</sup></th><th>Shannon<sup>b</sup></th></tr></thead><tbody><tr><td rowspan="4">2005</td><td>16S Bacteria</td><td>9.1 ± 0.1 × 10<sup>7</sup></td><td>2.5 ± 0.1 × 10<sup>7</sup></td><td>15/71</td><td>0.79</td><td>18.75 (15.63-37.02)</td><td>2.20 (2.47-1.97)</td></tr><tr><td>16S Archaea</td><td>3.3 ± 0.8 × 10<sup>6</sup></td><td>2.7 ± 0.8 × 10<sup>6</sup></td><td>3/96</td><td>0.97</td><td>3 (3-ND)</td><td>0.48 (0.33-0.63)</td></tr><tr><td>dsrB</td><td>1.1 ± 0.4 × 10<sup>6</sup></td><td>1.1 ± 0.4 × 10<sup>6</sup></td><td>6/28</td><td>0.79</td><td>5 (4.07-17.27)</td><td>0.92 (0.66-1.18)</td></tr><tr><td>mcrA</td><td>1.0 ± 0.1 × 10<sup>6</sup></td><td>1.0 ± 0.1 × 10<sup>6</sup></td><td>2/32</td><td>0.94</td><td>2 (2-2)</td><td>0.64 (0.52-0.76)</td></tr><tr><td rowspan="4">2010</td><td>16S Bacteria</td><td>4.0 ± 0.9 × 10<sup>7</sup></td><td>1.1 ± 0.9 × 10<sup>7</sup></td><td>55/109</td><td>0.54</td><td>385 (173.12-976.88)</td><td>3.46 (3.23-3.70)</td></tr><tr><td>16S Archaea</td><td>4.5 ± 0.9 × 10<sup>6</sup></td><td>3.7 ± 0.9 × 10<sup>6</sup></td><td>21/119</td><td>0.82</td><td>49 (28.41-126.69)</td><td>2.45 (2.26-2.64)</td></tr><tr><td>dsrB</td><td>3.4 ± 0.6 × 10<sup>5</sup></td><td>3.4 ± 0.6 × 10<sup>5</sup></td><td>1/29</td><td>0.97</td><td>1 (1-1)</td><td>0 (-0.047-0.047)</td></tr><tr><td>mcrA</td><td>1.3 ± 0.1 × 10<sup>6</sup></td><td>1.3 ± 0.1 × 10<sup>6</sup></td><td>1/27</td><td>0.96</td><td>1 (1-1)</td><td>0 (-0.081-0.081)</td></tr><tr><td rowspan="4">2011</td><td>16S Bacteria</td><td>6.7 ± 0.1 × 10<sup>7</sup></td><td>1.9 ± 0.1 × 10<sup>7</sup></td><td>33/80</td><td>0.64</td><td>64.66 (43.39-129.45)</td><td>3.03 (2.80-3.27)</td></tr><tr><td>16S Archaea</td><td>3.5 ± 0.6 × 10<sup>6</sup></td><td>2.9 ± 0.6 × 10<sup>6</sup></td><td>14/90</td><td>0.84</td><td>21.5 (15.32-56.53)</td><td>1.65 (1.37-1.93)</td></tr><tr><td>dsrB</td><td>1.0 ± 0.2 × 10<sup>6</sup></td><td>1.0 ± 0.2 × 10<sup>6</sup></td><td>6/33</td><td>0.82</td><td>6 (6-ND)</td><td>1.29 (0.99-1.58)</td></tr><tr><td>mcrA</td><td>1.8 ± 0.3 × 10<sup>6</sup></td><td>1.8 ± 0.3 × 10<sup>6</sup></td><td>4/29</td><td>0.86</td><td>5 (4.07-17.27)</td><td>0.92 (0.66-1.19)</td></tr></tbody></table>

between the sampling years (P > 0.05). Archaea accounted for 9.8 (2005) to 25.3% (2010) of the total prokaryotic cells, given an average 16S rDNA copy number in archaeal and bacterial genomes of 1.2 and 3.6, respectively (Fogel et al., 1999; Klappenbach et al., 2001). At LCHF, in comparison, the total number of DAPI stained cells counted by epifluorescence microscopy in the active venting chimneys varied from 0.02 to 8.4 × 10⁸ cells per gww, with archaea representing between 55 and 77% of the total microbial cells as determined by FISH experiments (Kelley et al., 2005). The *dsrB* gene abundance ranged from 3.4 ± 0.6 × 10⁵ (2010) to 1.1 ± 0.4 × 10⁶ (2005) copies per gww, while the *mcrA* abundance varied from 1.0 ± 0.1 × 10⁶ (2005) to 1.8 ± 0.3 × 10⁶ (2011) copies per gww. Assuming that archaeal and bacterial cells respectively contain only one copy of the *mcrA* or *dsrB* gene (Klein et al., 2001; Luton et al., 2002), the methanoarchaea accounted for 35.6% (2010) to 60.2% (2011) of the archaea, whereas the percentage of sulfate-reducers was less than 6% of the bacteria. While such quantitative PCR data are not available for other hyperalkaline systems, LCHF FISH experiments showed *Methanosarcinales* specific probe hybridized between 68 and 80% of the archaeal cells.

DGGE profiling of the bacterial and archaeal communities in the ST09 chimneys was performed for each sampling year. Bacterial profiles were very similar although for the 2010 sample, the major bands were less visible and included an additional band (n°5) (Supplementary Figure 1). The corresponding sequences were affiliated to *Firmicutes*, *Chloroflexi*, and *Bacteroidetes* (Table 2). Archaea profiles were also quite similar for 2005 and 2011 samples contrasting with the 2010 sample which contained only one distinct band (n°11) particularly intense, affiliated to the *Thaumarchaeota* (*Cenarchaeales*). All other archaea detected using DGGE belonged to the *Methanosarcinales*. Most of them (bands n° 7, 8, 9, 10, and 12) shared 97-99% sequence identity

with the Ced_A01 clone from The Cedars springs whilst one band (n°6) was affiliated to the LCMS cluster.

The diversity of the microbial community was assessed for each sample by clone library sequences of 16S rRNA and functional genes (*dsrB*, *mcrA*). The number of OTUs, the values of the richness estimators (Table 1) and the shape of the rarefaction curves (Supplementary Figure 2) obtained from 16S rRNA gene sequence analyses confirmed that archaea were much less diverse than bacteria in each sample. The coverage was the highest for the 2005 clone libraries and the lowest for the 2010 clone libraries, even if the latter counted a higher number of clones, reflecting a higher diversity in 2010 samples than in 2005.

Archaeal diversity was well covered with values ranging from 79 to 97% (Table 1), suggesting that dominant members of the archaeal community were represented in the clone libraries, as observed in recent studies on microbial communities from serpinization-driven subterrestrial ecosystems (Daae et al., 2013; Suzuki et al., 2013). Archaeal 16S rRNA sequences (n = 305) belonged to the phyla *Euryarchaeota* (69%) or *Thaumarchaeota* (31%).

Coverage values of the bacterial clone libraries ranged from 54% (2010) to 79% (2005), indicating that more clones may be needed for an adequate estimation of the bacterial diversity. Even so, we considered that the major bacterial groups were represented. Bacterial 16S rRNA gene sequences (n = 260) spanned many phyla (Figure 2). In average (over the three chimneys sampled), the *Firmicutes* phylum was the most represented (17.3% of total bacterial clones), followed by the *Chloroflexi* (16.9%), Candidate Divisions OP1 (13.8%), and NPL-UPA2 (13.8%), *Alphaproteobacteria* (12.7%), *Betaproteobacteria* (5.8%), *Gammaproteobacteria* (5.8%) and *Deltaproteobacteria* (5.8%). To a lesser extent (≤1.5% of bacterial clones), Candidate Division OD1, *Acidobacteria*, *Actinobacteria*, *Planctomycetes*, *Bacteroidetes*, *Deinococcus-Thermus*, *Epsilonproteobacteria*,

<footnote><sup>a</sup>OTUs (operational taxonomic units) defined at 3% sequence differences (nucleotide or amino acids) threshold.</footnote>

<footnote><sup>b</sup>Ranges of diversity indices include 95% confidence intervals as calculated by MOTHUR (Schloss et al., 2009). WW, Wet Weight; ND, Not Determined.</footnote>
<table><caption>TABLE 2 | Phylogenetic affiliations of the 16S rRNA sequences obtained from PCR-DGGE analyzes.</caption><thead><tr><th>DGGE band N°</th><th>Phylogenetic<br>affiliation</th><th>Closest PHF phylotype<br>(BLAST identity %)</th><th>Closest environmental sequence<br>(BLAST identity %)</th><th>Origin</th></tr></thead><tbody><tr><td>DG-Bac1-341F (KJ149114)</td><td>Firmicutes</td><td>PHF_13-B21_J06<br>(KJ149248) (99%)</td><td>clone B11-15_GoMY (AB476673)<br>(98%)</td><td>Sulfur containing freshwater source</td></tr><tr><td>DG-Bac2-341F (KJ149115)</td><td>Chloroflexi</td><td>PHFST12_B1 (KF886112)<br>(98%)</td><td>clone KM35B-155 (AB300104) (96%)</td><td>Holocene mud sediment</td></tr><tr><td>DG-Bac3-341F (KJ149116)</td><td>Firmicutes</td><td>PHF_15-B29_J20<br>(KJ149238) (92%)</td><td>clone ML635J-16 (AF507885) (91%)</td><td>Mono lake at a depth of 35 m</td></tr><tr><td>DG-Bac4-341F (KJ149117)</td><td>Firmicutes</td><td>PHF_15-B34_D22<br>(KJ159205) (100%)</td><td>clone 56S_1B_61 (DQ837275) (95%)</td><td>Pristine coastal aquifer, Spain</td></tr><tr><td>DG-Bac5-341F (KJ149118)</td><td>Bacteroidetes</td><td>PHF-2HY3BaC08<br>(KJ149173) (83%)</td><td>clone OTU199_Ref Clone01<br>(AB694464) (95%)</td><td>Deep-sea sediment</td></tr><tr><td>DG-Arc6-Saf (KJ149119)</td><td>Methanosarcinales</td><td>PHF_2AarcF09 (KJ149145)<br>(100%)</td><td>clone SGXU600 (FJ791432) (98%)</td><td>Carbonate chimney in Lost City</td></tr><tr><td>DG-Arc7-Saf (KJ149120)</td><td>Methanosarcinales</td><td>PHF_13-A45_G13<br>(KJ149156) (97%)</td><td>clone Ced_A01 (KC574884) (97%)</td><td>Serpentinizing springs in The Cedars</td></tr><tr><td>DG-Arc8-Saf (KJ149121)</td><td>Methanosarcinales</td><td>PHFST07_A2 (KF886032)<br>(99%)</td><td>clone Ced_A01 (KC574884) (99%)</td><td>Serpentinizing springs in The Cedars</td></tr><tr><td>DG-Arc9-Saf (KJ149122)</td><td>Methanosarcinales</td><td>PHF_13-A45_G13<br>(KJ149156) (98%)</td><td>clone Ced_A01 (KC574884) (98%)</td><td>Serpentinizing springs in The Cedars</td></tr><tr><td>DG-Arc10-Saf (KJ149123)</td><td>Methanosarcinales</td><td>PHF_13-A45_G13<br>(KJ149156) (99%)</td><td>clone Ced_A01 (KC574884) (98%)</td><td>Serpentinizing springs in The Cedars</td></tr><tr><td>DG-Arc11-Saf (KJ149124)</td><td>Thaumarchaeota</td><td>PHF_13-A02_M23<br>(KJ149153) (100%)</td><td>clone SPCiA-7 (KF036083) (99%)</td><td>Marine sponge, Cinachyra sp.</td></tr><tr><td>DG-Arc12-Saf (KJ149125)</td><td>Methanosarcinales</td><td>PHF_13-A45_G13<br>(KJ149156) (99%)</td><td>clone Ced_A01 (KC574884) (99%)</td><td>Serpentinizing springs in The Cedars</td></tr></tbody></table>

Gemmatimonadetes, candidate division OP9 and Nitrospira were also detected.

The predominance of *Euryarchaeota* (*Methanosarcinales*), *Firmicutes*, and *Chloroflexi* was also attested by DGGE. Both clone libraries composition and DGGE profiles suggest that the microbial communities were similar in 2005 and 2011, while *Thaumarchaeota*, *Gamma-*, and *Betaproteobacteria* corresponding to putative aerobic microorganisms were more abundant in 2010. The changes in population structure observed through the years may be explained by temporal variation in hydrothermal activity of the site in terms of discharge rate and composition of the fluid (Brazelton et al., 2006, 2010a). Unfortunately, we do not have such time series for ST09 fluid geochemistry because fluid was not sampled in 2005 and because of the difficulties to collect undiluted fluid samples by SCUBA diving (especially in deep sites). We obtained reliable chemical composition of the fluid in ST09 only since 2011 (Monnin et al., 2014). Distinct chimneys were sampled in 2005, 2010, and 2011. Although, they seemed to be at a roughly similar stage of formation, heterogeneity in the mineralogy or porosity of the subsamples as well as variation in the colonization of their outer wall by invertebrates may have influenced the permeability of the carbonate concretion and therefore the level of mixing with the surrounding seawater which can explain the marked difference of the microbial distribution in the chimney collected in 2010. Another possible explanation could lie in the fact that the 2010 chimney could not be processed as quickly as the two others and therefore that aerobic microorganisms from the outer parts of the chimney exposed to seawater may have diffused in its interior.

## Two Dominant *Methanosarcinales* Clusters (LCMS and TCMS) Coexist at PHF

One of the most salient results from the 16S rRNA diversity survey was the presence of two major clusters belonging to the *Methanosarcinales* order that largely dominated the whole archaeal population (up to 99% of the clones in 2005). Both clusters were detected in 2005, 2010, and 2011 samples (Figure 3) and their predominance was supported by DGGE results (Table 2, Supplementary Figure 1).

Three PHF phylotypes (PHF-2A-A13 I03, PHF-HY5ArA02, and PHF-13-A45 G13) totaling 142 clones closely matched with the unique (and minor) archaeal phylotype Ced-A01 detected so far only in The Cedars terrestrial serpentinizing springs (here named TCMS by analogy with the LCMS of Lost City), while the closest cultivated organisms were of the *Methanosaeta* genus (=86% identity) known as obligate acetoclastic methanogens (Suzuki et al., 2013). The second cluster contained five phylotypes (PHF-2AarcF09, PHF-2HY5ArE02, PHF-13-A03 K23, PHF-13-A49 I13, and PHF-15-A33 A09) totaling 34 clones that were clearly affiliated to the famous LCMS described at Lost City (93.4–95.8% identity with the LCMS clones LC1022a-1 and LC1149a-56) (Figure 3). Previous studies showed that LCMS cells represented up to 80% of the total cells in biofilms coating active carbonate chimneys of LCHF (Brazelton et al., 2010a,b). Accordingly, the seven *mcrA* OTUs obtained (over a total of 88 clones) fall into a single group within the *Methanosarcinales* order suggesting that all the methanoarchaea of PHF ST09 chimneys belonged to this order (Supplementary Figure 3). They were all closely affiliated with uncultured LCHF *mcrA* sequences
<figure><img src="image_4.png" /></figure>

**FIGURE 2** | Prokaryotic community structure in the interior parts of ST09 Prony chimneys based on phylum-level taxonomic distribution of the (A) bacterial and (B) archaeal 16S rRNA gene clones. *Phylum Proteobacteria* is divided into class-level taxons.

Sampling years are shown on the left of the bar charts (2005, 2010, and 2011) with the number of analyzed clones. "Cand. Div." stands for Candidate Division; SAGMEG, South Africa Gold Mine Euryarchaeotic Group; "Eury," Euryarchaeota.

(90.0–91.8% identity) (Kelley et al., 2005) and more distantly related to *mcrA* group f sequences of ANME-3 (75.6–83.7% identities) (Lösekann et al., 2007) (Supplementary Figure 3). There was no sequence of *mcrA* reported for TCMS (Suzuki et al., 2013) as reference, therefore, it was difficult to foresee where such sequences will branch in a *mcrA* phylogenetic tree. However, at the closest similarity level, two distinct clusters can be distinguished within PHF *mcrA* sequences, the first one (15-M25, 15-M26, and 15-M1) contains 16% of the total clones (all years) while the second (2-C-M1, 23-M22, 23-M1, 15-M3) represents 84% of the total clones, a ratio that roughly coincides with the ratio of LCMS-like vs. TCMS-like ratio of the 16S rRNA gene sequences, 19 vs. 81% respectively.

These two dominant *Methanosarcinales* phylotypes were up to now considered as unique and specific to two highly contrasted

serpentinizing ecosystems, one fully submarine (LCHF type) and the other fully continental (The Cedars). Their coexistence in PHF chimneys is remarkable and may reflect the singularity of PHF representing a transition between these two systems: the discharge of hydrothermal alkaline freshwater of continental origin in a shallow marine environment (Monnin et al., 2014). We propose that this *Methanosarcinales* group represents a hallmark of serpentinization-hosted ecosystems. The surprising higher abundance of TCMS over LCMS may be due to the low salinity of PHF fluids of meteoric origin, such as it is at The Cedars, but different from the fully marine LCHF. Thus, depending on the mixing degree of the low salinity fluid with surrounding seawater, both LCMS and TCMS may coexist.

Direct observations of carbonate samples by phase-contrast microscopy revealed dense microbial communities
<figure><img src="image_5.png" /></figure>

**FIGURE 3** | Maximum likelihood phylogenetic tree representing *Methanosarcinales* (among *Euryarchaeota*) 16S rRNA gene phylotypes from ST09 chimneys. Based on 489 aligned positions. The tree was constructed using the Maximum likelihood method. Bootstrap values <50% are not shown. Clone libraries are

distinguished by various symbols and colors: 2005 (■ in red), 2010 (▲ in green), and 2011 (● in blue). Purple font refers to phylotypes from marine hydrothermal serpentinized systems, pink font refers to phylotypes from other terrestrial alkaline sites. Scale bar, number of substitutions per site.

consisting mainly of thin long rods and cocci inhabiting
pore spaces within the chimneys of ST09 (not shown). The
cells seemed to be embedded in a matrix assumed to be
microbial exopolysaccharides that could be responsible for
the mucilaginous aspect around the central conduit of active
chimneys (dark pink to brown colored, **Figure 1E**). Only the
small cocci auto-fluoresced when excited at a wavelength
of 420 nm, a characteristic of archaea containing cofactor
420 (F₄₂₀) (Doddema and Vogels, 1978), a key coenzyme
involved in methanogenesis (Thauer, 1998). These small
fluorescent cocci occurred alone or typically in small spherical
aggregates (of around 5 µm in diameter) of 10 to more than
20 tiny cells (around 0.5 µm in diameter) (**Figure 4A**). This
organization resembles the structure known as cyst with a
common outer wall surrounding coccoid cell aggregates, typical
of the *Methanosarcina* genus and in particular of *M. acetivorans*
(Sowers et al., 1984). The coupling of FISH to CSLM showed that

the same aggregates of autofluorescent cocci hybridized with
an *Euryarchaeota* specific probe, supporting the hypothesis that
they belong to *Methanosarcinales* (**Figures 4B,C**).

LCMS were shown to be able to produce methane and to
oxidize it in anaerobic conditions at the same time in single-
species biofilms (Brazelton et al., 2012). The enhancement of
both metabolisms by dihydrogen (H₂) suggested syntrophic
interactions between the two metabolically differentiated cell
types (Brazelton et al., 2011). However, the cultivation of LCMS
in pure culture has not been reported so far.

Many attempts were carried out to cultivate PHF
*Methanosarcinales* by varying both the pH (from 8 to 10.5)
and temperature conditions (30–60°C) and the electron donors
(H₂, acetate, formate, methanol, MMA, TMA). Limited methane
production was detected after 5–12 days of incubation at
37°C during a few enrichment culture experiments using
TMA or H₂/CO₂. Scarce tiny coccoid cells and thin filaments
<figure><img src="image_6.png" /><figcaption>**FIGURE 4** | Confocal Laser Scanning Microscopy stacked images obtained with a sequential excitation at 405 and 488 nm and fluorescence detection in the ranges 425–475 and 500–530 nm.</figcaption></figure>

(A) Clusters of green Syto®9-stained cells interspersed in the chimney-forming minerals appearing as blue-violet. The archaeal cocci

appear in turquoise thanks to the superimposition of their autofluorescence in blue and DNA staining with Syto®9 in green.

(B) Archaeal cocci stained with DAPI (the self-fluorescence of F<sub>420</sub> cofactor is also collected) and (C) hybridized with the FITC-labeled specific probe for *Euryarchaeota* EURY498.

were observed in these mixed cultures, but the methanogenic populations could not be maintained in subcultures.

## Are there Thermophilic Archaea Living in PHF?

Besides the unexpected discovery of LCMS in PHF, another surprising finding was the identification of many sequences closely related to cultivated members of the hyperthermophilic order *Thermococcales* and of environmental sequences related to hot biotopes, such as deep sea black smokers (Figure 5). Remarkably, the *Thermococcales* 16S rRNA sequences of Prony formed a single phylotype, representing 13% of the 2011 archaeal clone library, and identical to *Thermococcus litoralis* which was isolated from a shallow submarine volcanic spring in Italy (Neuner et al., 1990). *Thermococcales* are mostly heterotrophic, anaerobic, sulfur-reducers commonly found in marine hydrothermal vents. A *priori*, their growth conditions (optimum pH 6.5, up to pH 8.5, and a temperature range of 55–98°C) are not compatible with the physical and chemical conditions found in the PHF venting chimney (pH of 10.6 at ST09, and temperature at max 40°C). Scarce *Thermococcales* sequences were also reported in Lost City but only in fluid samples (Brazelton et al., 2006), which were almost identical to an isolate (related to *T. celer*) from a deep-sea volcano of the East Pacific Rise (Huber et al., 2006). Intriguingly, previous DGGE and clone library analyses from a cold subsurface oceanic sediment sample (2–3°C, 1.5 mbsf) collected at 3019 m water depth in the New Caledonian basin revealed low archaeal diversity dominated by *Thermococcales* clones that all grouped in a unique cluster (Roussel et al., 2009). Since no local hydrothermal activity was detected, the authors suggested that these *Thermococcales* could have rafted from PHF located about 200 km North East. The minimum temperature required for the growth of a *Thermococcus* is 40°C (Miroshnichenko et al., 2001), which questions the thermophilic nature of these archaea found in PHF. Are they an indicator of a hot subseafloor biotope linked to the exothermic reactions of the serpentinization process (Allen and Seyfried Jr, 2004; Takai et al., 2004) as suggested by Brazelton et al. (2006) for LCHF or does this PHF phylo type represent a new moderately thermophilic ecotype of *T. litoralis*? At this stage we cannot conclude since our attempts to cultivate them at

various temperatures (from 37 to 75°C) failed. This question also relates to the LCMS-like cells of PHF since LCMS were described to colonize only venting chimneys of LCHF, i.e., at temperatures in the hyperthermophilic range (between 55 and 95°C). The presence of LCMS-like cells in PHF at moderate temperatures (max 40°C) could therefore be the indication of the existence of a hot biotope beneath the PHF seafloor or may simply indicate that the PHF LCMS-like cells are, in fact, mesophilic archaea. Nonetheless, as previously mentioned, our preliminary attempts to cultivate these elusive LCMS were also not conclusive in the temperature range of 30–60°C. This addresses more generally the question of estimating the temperature lifestyle of uncultured prokaryotes in their natural environments. This problem was recently treated by Kimura et al. (2010, 2013) who estimated the growth temperatures of archaeal communities from the guanine-plus-cytosine contents (P<sub>GC</sub>) based on 16S rRNA sequences. The growth temperatures of the archaeal phylotypes found in PHF were predicted using their method (Table 3). For the sake of comparison, clones from similar environments, i.e., serpentinized systems or hydrothermal vents, and of a few representative archaeal strains were also reported. The predicted growth temperatures of the PHF *Thermococcales* clearly qualify them as hyperthermophilic with a temperature growth range very similar to that predicted and experimentally verified for their closest cultivated relative, *T. litoralis*. If these predictions hold true, all the PHF phylotypes assigned to the South African Gold Mine Euryarchaeotic Group SAGMEG (4% of the total archaea), could also be considered as thermophiles with an average predicted optimal growth temperature of about 70°C. This group is frequently associated with subsurface sediments and methane-rich environments (Parkes et al., 2014). As they were discovered in gold mine alkaline groundwater linked to a deep hot aquifer, they were assumed to be thermophiles (Takai et al., 2001) (Figure 5) although their actual morphology, metabolism and lifestyle remain unknown. Both *Thermococcales* and SAGMEG phylotypes could therefore indicate the existence of a hot subsurface biosphere in PHF.

Using the same criterion, all the other archaea detected in PHF should be considered as mesophiles with a predicted optimal growth temperature (T<sub>opt</sub>) of 43°C at most. It notably
<figure><img src="image_7.png" /></figure>

includes all the *Methanosarcinales* (LCMS-like and TCMS-like) detected in PHF. Indeed, the average predicted growth temperatures of LCMS-like ($T_{opt}$ 19°C, $T_{max}$ 28°C) and those of TCMS-like ($T_{opt}$ 31°C, $T_{max}$ 40°C) were close to the temperatures recorded *in situ* at fluid vents. Based on this criterion, the predominance of *Methanosarcinales* sequences in

all our chimney samples and our microscopic observations of typical *Methanosarcinales* morphotypes in these samples, we conclude that the *Methanosarcinales* found in PHF are mesophilic. What is puzzling is the fact that the predicted growth temperatures for true LCMS sequences, that we included in our analyses (from Brazelton et al., 2006) for comparison
<table><caption>TABLE 3 | Predicted growth temperature of uncultivated archaea and reference strains as a function of their 16S rRNA gene G+C percent (PGC) calculated from (Kimura et al., 2013).</caption><thead><tr><th>Clone</th><th>%G+C</th><th>T<sub>min</sub></th><th>Topt</th><th>Tmax</th><th>Affiliation</th></tr></thead><tbody><tr><td>The Cedars ced_A01 (KC574884)</td><td>54.69</td><td>14</td><td>31</td><td>39</td><td>TCMS</td></tr><tr><td>PHF_A45_G13 (50 clones)</td><td>54.94</td><td>15</td><td>32</td><td>40</td><td>TCMS</td></tr><tr><td>PHF_HY5ArA02 (80 clones)</td><td>54.63</td><td>14</td><td>30</td><td>39</td><td>TCMS</td></tr><tr><td>PHF_2A-A13-I03 (12 clones)</td><td>54.93</td><td>15</td><td>32</td><td>40</td><td>TCMS</td></tr><tr><td>PHF_2HY5ArE02 (15 clones)</td><td>52.25</td><td>4</td><td>19</td><td>27</td><td>LCMS</td></tr><tr><td>PHF_2AarcF09 (10 clones)</td><td>52.75</td><td>6</td><td>21</td><td>30</td><td>LCMS</td></tr><tr><td>13-A03_K23 (7 clones)</td><td>52.50</td><td>5</td><td>20</td><td>29</td><td>LCMS</td></tr><tr><td>LC1022a1 (AY299516)</td><td>53.00</td><td>7</td><td>22</td><td>31</td><td>LCMS</td></tr><tr><td>LC1149a56 (AY299515)</td><td>53.25</td><td>8</td><td>24</td><td>32</td><td>LCMS</td></tr><tr><td>LC_AR37 (FJ640804)</td><td>51.88</td><td>2</td><td>17</td><td>25</td><td>LCMS</td></tr><tr><td>LC1231a51 (AY505046)</td><td>52.70</td><td>6</td><td>21</td><td>30</td><td>LCMS</td></tr><tr><td>Methanosota_harundinacea_6Ac (AY970347)</td><td>57.75</td><td>28 (25)</td><td>46 (37)</td><td>54 (45)</td><td>Methanosarcinales</td></tr><tr><td>Methanothrix_thermophila_(genome)</td><td>58.88</td><td>33</td><td>52 (55)</td><td>60</td><td>Methanosarcinales</td></tr><tr><td>PHF_2CarcG12 (4 clones)</td><td>54.12</td><td>12</td><td>28</td><td>36</td><td>MG-II</td></tr><tr><td>PHF_13A05_G23 (3 clones)</td><td>57.08</td><td>25</td><td>43</td><td>51</td><td>MG-II</td></tr><tr><td>LC_FS243A20_(DQ270599)</td><td>55.63</td><td>18</td><td>35</td><td>44</td><td>MG-II</td></tr><tr><td>PHF_2CA44_G23 (1 clone)</td><td>54.75</td><td>15</td><td>31</td><td>40</td><td>MG-II</td></tr><tr><td>PHF_13A48_A13 (1 clone)</td><td>55.98</td><td>20</td><td>37</td><td>46</td><td>Thermoplasmatales</td></tr><tr><td>PHF_2AarcE10 (1clone)</td><td>52.38</td><td>4</td><td>19</td><td>28</td><td>Unknown deep-sea vent</td></tr><tr><td>PHF_A16_C21 (3 clones)</td><td>61.25</td><td>43</td><td>63</td><td>71</td><td>SAGMEG</td></tr><tr><td>PHF_2CA32_019 (1 clone)</td><td>61.62</td><td>45</td><td>65</td><td>73</td><td>SAGMEG</td></tr><tr><td>PHF_2CarcH09_(4 clones)</td><td>62.75</td><td>50</td><td>71</td><td>79</td><td>SAGMEG</td></tr><tr><td>PHF_2HY7ArH02 (1 clone)</td><td>62.24</td><td>47</td><td>68</td><td>76</td><td>SAGMEG</td></tr><tr><td>3H3M_ARC52_(JN229767) (subseafloor_sediment)</td><td>62.25</td><td>47</td><td>68</td><td>76</td><td>SAGMEG</td></tr><tr><td>ODP1251A1.27_(AB177265) (deep_sediments,_Pacific)</td><td>61.38</td><td>44</td><td>64</td><td>72</td><td>SAGMEG</td></tr><tr><td>CK0606_Mud_MAS4A10_(AB369199) (drilling_mud_fluid,_Japan)</td><td>63.00</td><td>51</td><td>72</td><td>80</td><td>SAGMEG</td></tr><tr><td>PHF_A40_C15 (12 clones)</td><td>65.75</td><td>63</td><td>86</td><td>93</td><td>Thermococcales</td></tr><tr><td>Thermococcus_litoralis type strain: DSM5473 (FR749899)</td><td>65.87</td><td>63 (55)</td><td>86 (88)</td><td>94 (98)</td><td>Thermococcales</td></tr></tbody></table>

$T_{min}$, $T_{opt}$, and $T_{max}$: minimal, optimal, and maximal predicted temperatures for growth, respectively (in bracket are values experimentally determined for such growth parameters of described species when available in the literature). Color shading code: blue for mesophilic ($T_{opt}$ from 17 to 24° C); yellow for $T_{opt}$ from 28 to 43°C; orange: $T_{opt}$ from 46 to 72°C; red: $T_{opt}$ > 80°C. Text in bold: PHF, clones from this study; LCMS, Lost City Methanosarcinales; MG-II, Marine Group II; SAGMEG, South African Gold Mine Euryarchaeotic Group; TCMS, The Cedars Methanosarcinales.

purposes, were also in the same range as those of LCMS-like sequences from PHF and even slightly lower than that of TCMS and TCMS-like. The hypothesis that LCMS sequences constitute an exception in Kimura's prediction method is unlikely since it was recently shown to apply to putatively thermophilic anaerobic methanotrophic archaea (ANME) from hot deep-sea vents (Merkel et al., 2013). Therefore, we may consider that the actual temperature growth range of these enigmatic archaea is much cooler than expected based on the fluid temperature measured at venting chimney outlets (70-90°C), and is probably more within the mesophilic range.

## Cosmopolitan *Thaum-* and *Eury-archaeota*

Archaea at PHF also includes phylotypes assigned to *Thaumarchaeota* and to *Euryarchaeota* Marine Group II which is not surprising because these two groups are the most abundant and frequently detected planktonic archaea (Delong, 2003; Galand et al., 2009) in various oceanic provinces (see the example of a tropical atoll, Michotey et al., 2012). *Thaumarchaeota* detected in PHF were particularly abundant

in 2010 (70.6%) but less in 2011 (12.2%) clone libraries and showed high intra-phyllum diversity with representatives of both the *Cenarchaeales* and the *Nitrososphaerales* orders (Supplementary Figure 4). The closest relatives of the PHF phylotypes were environmental sequences from hydrothermal vents, marine sediments and marine sponges, including clones of LCHF chimneys. The first isolated *Thaumarchaeota* members showed chemoautotrophic growth by aerobic oxidation of ammonia to nitrite and played an important role in the biogeochemical cycling of nitrogen (Schleper and Nicol, 2010). Members of the *Euryarchaeota* MG II were less abundant than *Thaumarchaeota* in PHF but like them, they were related to clones from either deep or surface oceanic waters, including one clone from LCHF. There is yet no cultivated *Euryarchaeota* MG II but a partial genome representative of this group was reconstructed from metagenomic data of surface seawater. As stated by Iverson et al. (2012) *sic* "The genome describes a motile, photo-heterotrophic cell focused on degradation of protein and lipids and clarifies the origin of proteorhodopsin."
## Firmicutes Dominate the PHF Bacterial Community

Firmicutes was globally the most abundant bacterial phylum retrieved in this study although its proportion significantly varied between 2005, 2010, and 2011 (from 3.4 to 29.3%, **Figures 2, 6**). One particular cluster referred to here as “subsurface” was dominant and contained persistent phylotypes closely related to those found in a gold mine (Lin et al., 2006) and in CVA, Portugal (99.3% identity with CVCloAm2Ph102) (Tiago and Veríssimo, 2013). Another cluster of uncultivated organisms referred to here as “alkaline groundwater” encompasses two PHF phylotypes together with two OTUs of CVA and two OTU of The Cedars (**Figure 6**). The metabolism of these uncultivated members is so far unknown.

Assuming that the “true” indigenous fraction of the cultivable bacteria inhabiting the inside part of active PHF chimneys is suitably adapted to the *in situ* conditions, i.e., high pH fluids, low redox, and oxygen but high hydrogen and methane concentrations, we decided to apply anaerobic culture techniques to target these bacteria, guided with the results of the molecular survey. Only few enrichments cultures were successful. We tried to cultivate sulfate-reducing bacteria (SRB) from PHF chimney samples using various combinations of electron donors and acceptors at different pH and temperature conditions (see Materials and Methods). Unexpectedly, instead of SRB, we isolated from these cultures seven new strains of the order *Clostridiales*, all being strictly anaerobic fermentative heterotrophs. *Vallitalea pronyensis* sp. nov. was isolated from a sample of the ST09 chimney collected in 2010 (Ben Aissa et al., 2014). Since it is not alkaliphilic we assume that this bacterium develops in anaerobic niches of the chimney not directly in contact with pure hyperalkaline fluid. Four new strains of *Alkaliphilus* species were also isolated from PHF representing the first anaerobic alkaliphiles isolated from serpentinizing environments (Ben Aissa et al., 2015); they all are anaerobic and fermentative, able to grow by fermenting a variety of proteinaceous compounds and sugars in moderately alkaline conditions (optimally at pH 8.7–9.5). Organic matter supporting their growth in the field may be metabolites produced by chemosynthetic primary producers (i.e., hydrogenotrophs, methanogens) or by their decays. Alternatively, they might be fed by small organic molecules (short chain hydrocarbons) abiotically generated by Fisher-Tropsch-type reactions linked to serpentinization (Proskurowski et al., 2008). Moreover, thermogenic organic compounds generated by the thermal degradation of deep serpentinite-hosted ecosystems (microbes or metazoans) can significantly contribute to the organic carbon balance in peridotite-hosted hydrothermal fields (Pasini et al., 2013).

*Clostridia* detected in various serpentinizing environments (Lost City, CVA, and the Cedars) were assumed to be abundant in the most extreme zones, i.e., the deepest anoxic parts, due to their high metabolic versatility and their autotrophy or mixotrophy capabilities (**Figure 7**) (Schrenk et al., 2013). Metagenomic data indicated that *Clostridiales* are major hydrogen producers in such alkaline ecosystems (Brazelton et al., 2012) which was further confirmed by culture-based approaches (Mei et al., 2014).

## Sulfur-cycling Related Bacteria

Sulfate-reducing bacteria of the class *Deltaproteobacteria* represented 5.3% from total bacterial clones (6% according to quantitative PCR). Despite their significant abundance in the original samples, culture media designed to target sulfate-reducers failed to enrich any *Deltaproteobacteria* (see above). *Desulfonatronum* (*Desulfovibrioles*), highlighted by 16S rRNA and *dsrB* sequences, appeared dominant and persistent in ST09 PHF chimneys over the 6-year period, suggesting its ecological significance in this ecosystem. Most of the *dsrB* OTUs were affiliated to *D. lacustre* (88.4–97.0% identity, representing 100, 97, and 54% of the bacterial clones respectively in the 2010, 2011, and 2005 libraries) while 16S phylotypes were related to *D. cooperativum* (Supplementary Figures 5, 6). *Desulfonatronum* species were isolated from various soda lakes and described as lithotrophic haloalkaliphilic sulfate-reducing bacteria, able to use hydrogen and a few organic compounds as electron donors (Zhilina et al., 2005). *Desulfobacterales* were also well represented by both 16S rRNA and *dsrB* sequences. In particular, *Desulfobulbaceae*-related *dsrB* OTUs were affiliated to the hydrogenotrophic *Desulfurivibrio alkaliphilus* (Sorokin et al., 2008) (91.4–93.6% identity, representing 29% and 3% of the 2005 and 2011 clone libraries). *Desulfobacteraceae*-related *dsrB* OTUs (PHF-23-D10 and OTU PHF-23-D21) were detected only in the 2005 sample and were related to halophilic sulfate-reducing bacteria, *Desulfotignum balticum* (95.7% identity) and *Desulfosarcina variabilis* (91.4% identity) respectively, which have often been detected in both pristine and contaminated marine environments (Joulian et al., 2001; Kuever et al., 2001; Gillan et al., 2005). *Desulfobacterales* are known to be associated with anaerobic methanotrophic archaea (ANME) (Schreiber et al., 2010); they represent putative syntrophic partners of anaerobic oxidation of methane (AOM) in association with ANME-3 related *Methanosarcinales*, as previously reported in deep sea hydrothermal vents and methane seeps (Knittel et al., 2003). At Lost City, *Desulfovibrioales* and *Desulfobacterales* were only rarely detected by 16S rDNA tag-sequencing and no corresponding *dsrB* gene was found (Brazelton et al., 2010a). Instead, *Desulfotomaculum* represented the dominant SRB detected in LCHF carbonate chimneys (Brazelton et al., 2006; Gerasimchuk et al., 2010) while they were not found in PHF ST09 chimneys. The differences observed in the patterns of SRB populations between LCHF and PHF might be explained by differences in fluid chemical composition (Supplementary Table 1) and in particular the salinity which is about 10 times higher in LCHF fluids than in PHF, whilst sulfate is about the same (1–4 mM in LCHF and ~3 mM in PHF). Thus *Desulfovibrioales* and *Desulfobacterales*, which dominate in PHF chimneys, may be better adapted to low salinity environments.

Additionally, within *Firmicutes*, *Dethiobacter*-like sequences identified in PHF clustered with clones from the terrestrial serpentinizing sites The Cedars and CVA (**Figure 6**). The type species *D. alkaliphilus* is an alkaliphilic and facultative hydrogenotrophic autotroph isolated from soda lake sediments (Sorokin et al., 2008).
<figure><img src="image_8.png" /></figure>
<figure><img src="image_9.png" /></figure>

**FIGURE 7** | Geomicrobiological model of the serpentinization-driven submarine ecosystem of PHF. The physico-chemical conditions and mineralogy depicted in the zonation delimiting distinct microbial habitats are based on *in situ* measurements and lab analyses when available, or they have been hypothesized (such as in those in the deep subsurface settlement

of the carbonate edifice). The potential electron donors and acceptors in the fluid and in seawater are given along with the putative metabolisms assigned to all the stable prokaryotes populations detected over 6-year period covered by this study. OM: Organic Matter. The scheme was inspired from Figure 7 of Brazelton et al. (2006).
Putative sulfur-oxidizing chemolithoautotrophs were also identified, related to the genera *Sulfurimonas* (*Epsilonproteobateria*, Supplementary Figure 7), *Sulfuritalea* (*Betaproteobateria*, Supplementary Figure 8) and mainly *Thioprofundum* and *Thiomicrospira* (*Gammaproteobacteria* Supplementary Figure 9). *Thiomicrospira* phylotypes identified in PHF were closely related to sequences from LCHF which were shown to be dominant in the fluids (Zhang et al., 2012). These bacteria most likely thrive in the peripheral part of the chimney walls, where seawater mixing with hydrothermal fluid provides optimum concentrations of electron donors, H₂S, electron acceptor and O₂ (this scenario is depicted in **Figure 7**). Putative (moderately) thermophilic bacteria were also evidenced in PHF chimneys, which belong to the genera *Thiomicrospira*, *Thioprofundum*, *Schlegelella*, *Deinococcus-Thermus* (see below). As discussed above for thermophilic archaea, these bacteria were most likely carried out by the fluid, ascending from a hot subsurface environment.

## Chloroflexi and Candidate Divisions

*Chloroflexi* was the second most abundant bacterial phylum in the three studied samples (representing up to one third of the bacterial sequences in the 2010 library) (**Figure 2**). The same trend was observed in the other underwater sites of PHF (in ST08 and ST12) we previously studied, where *Chloroflexi* represented up to 47% of total bacteria (Quéméneur et al., 2014). Four phylotypes found in ST09 chimneys were closely related to The Cedars Ced_B01 sequence (up to 97.1% identity), belonging to the candidate order named MSBL5 (Mediterranean Sea Brine Lake group 5) (Supplementary Figure 10) and three other phylotypes clustered with *Dehalococcoidia* sequences from marine subseafloor and siliciclastic sediments. Altogether, these new data support the hypothesis that *Chloroflexi* play a major role in the functioning of the PHF chimneys ecosystem, although the exact metabolisms of its representatives in such a natural environment are still largely unknown (Wasmund et al., 2013).

Candidate divisions OP1 and NPL-UPA2 were especially well represented (each 13.8% of total clones) (Supplementary Figure 11). Interestingly, a distinct clade within the candidate division NPL-UPA2 contained six PHF phylotypes (found each year) together with several LCHF and The Cedars OTUs, these sequences sharing up to 95.6% identity.

## Conclusion

This study provides a detailed description of the microbial communities from carbonate samples collected at an underwater site of the Prony hydrothermal field between 2005 and 2011. In PHF, the discharge of low-salinity, high-pH waters (of meteoric origin) produced by serpentinization of ultrabasic basement into a coastal marine environment forms large carbonate concretions and needles similar to those of LCHF. Our study shows that the microbial diversity found in PHF reflects this intermediate situation with populations characteristic of serpentine environments similar to those found at on-land sites (i.e., *Chloroflexi* of the MSBL-5 cluster and TCMS from The Cedars, California) and to those characteristic of fully marine

environments (i.e., LCMS). For comparison purpose, reference clones from PHF ST09 chimneys were compared with those from other serpentinizing environments (Supplementary Table 3).

These microorganisms may occupy diverse microhabitats inside the highly porous structure of hydrothermal carbonate concretions, where the mixing of local oxygenated seawater with the highly reduced CH₄ and H₂ rich hydrothermal fluid ascending from the subsurface produces large ranges and very sharp gradients of pH, redox potential (Eh), salinity and dissolved element concentrations (**Figure 7**). For example, *Methanosarcinales* and *Firmicutes* dominate in anoxic microhabitats irrigated by hydrothermal end member fluid (warm, reducing, and alkaline) as confirmed by cultivation experiments. Minority groups representing (hyper)thermophilic lineages such as *Thermococcales* and *SAGMEG* in archaea or *Thermus* in bacteria were also consistently detected in the PHF chimneys but their predicted growth temperatures strongly suggest they were carried out by a fluid ascending from a hot subsurface environment heated by serpentinization reactions. Representatives of other metabolic groups are distributed along the physicochemical gradients depending mainly on their electron acceptor requirements and their tolerance to oxygen. Typically, SRB (*Deltaproteobacteria*), although most of them are oxygen tolerant (Muyzer and Stams, 2008), occupy anaerobic microhabitats where sulfate is provided by the mixing of the sulfate-free hydrothermal fluid with seawater. At the oxic pole of the gradient, microaerophilic or aerobic microorganisms such as methylotrophic bacteria (*Methylococcales*, *Methylophilales*), hydrogen-oxidizing (*Burkholderiales*) and sulfur-oxidizing bacteria (*Chromatiales*, *Thiotrichales*) thrives by oxidation of H₂, CH₄, and H₂S. Other marine aerobic bacteria and archaea found in PHF and LCHF do not seem to be specific of serpentine environments as they belong to cosmopolitan marine groups such as *Thaumarchaeota* or MG-II. Altogether, these microorganisms certainly play important roles in budgets of methane, hydrogen, and sulfur in PHF as depicted in **Figure 7**. However, these scenarios remain largely hypothetical since most of the microorganisms thriving in serpentinizing environments have not as yet been able to be cultivated. Yet, the first anaerobic alkaliphilic bacteria from a serpentinized ecosystem could be isolated from PHF ST09 chimneys. They belong to *Firmicutes*, a dominant phylum in PHF and in LCHF, generally known for their versatile metabolic capabilities but whose exact roles in the chimney ecosystem remain to be clarified.

## Acknowledgments

We thank the SCUBA divers: Jean-Louis Menou, Régis Hocdé, and Eric Folcher, John Butcher and Bertrand Bourgeois of the US IMAGO team at IRD research center in Nouméa for collecting the samples and for providing photos of the chimneys. We are grateful to all members of the crew of the R/V ALIS during the HYDROPRONY cruise for their assistance in sea operations. We are also indebted to the Scientific Direction of IRD who funded the scientific cruise HYDROPRONY in November 2011. The CNRS/INSU is acknowledged for its financial
support through the national program EC2CO-Biohefect/
Ecodyn/Dril/MicrobiEn to the MicroProny project devoted to
the study of microbial communities thriving inside chimneys of
the Prony Hydrothermal Field. This research was also supported
by The deepOASES ANR project (ANR-14-CE01-0008-01). CP
and BM gratefully acknowledge the support of the Deep Carbon
Observatory and the Deep Energy Community awarded by
Alfred P. Sloan Foundation (APS: 2011-06-01 (OSU) Grant;

DCO ID 11121/3515-2688-3796-9573-CC and APS: 2013-06-01
(OSU) Grant; DCO ID 11121/9902-5454-7041-1832-CC).

**Supplementary Material**

The Supplementary Material for this article can be found
online at: http://journal.frontiersin.org/article/10.3389/fmicb.
2015.00857

References

Allen, D. E., and Seyfried Jr, W. E. (2004). Serpentinization and heat generation: constraints from Lost City and Rainbow hydrothermal systems. *Geochim. Cosmochim. Acta* 68, 1347–1354. doi: 10.1016/j.gca.2003.09.003

Altschul, S. F., Gish, W., Miller, W., Myers, E. W., and Lipman, D. J. (1990). Basic local alignment search tool. *J. Mol. Biol.* 215, 403–410.

Bath, A. H., Christofi, N., Neal, C., Philip, J. C., Cave, M. R., Mckinley, I. G., et al. (1987). Trace element and microbiological studies of alkaline groundwaters in Oman, Arabian Gulf: a natural analogue for cement pore-waters. *British Geological Survey*, Fluid Processes Research Group, Report FLPU 87, 65.

Ben Aissa, F., Postec, A., Erauso, G., Payri, C., Pelletier, B., Hamdi, M., et al. (2014). *Vallitalea pronyensis* sp. nov., isolated from a marine alkaline hydrothermal chimney. *Int. J. Syst. Evol. Microbiol.* 64, 1160–1165. doi: 10.1099/ijs.0.055756-0

Ben Aissa, F., Postec, A., Erauso, G., Payri, C., Pelletier, B., Hamdi, M., et al. (2015). Characterization of *Alkaliphilus hydrothermalis* sp. nov., a novel alkaliphilic anaerobic bacterium, isolated from a carbonaceous chimney of the Prony hydrothermal field, New Caledonia. *Extremophiles* 19, 183–188. doi: 10.1007/s00792-014-0697-y

Brazelton, W. J., Ludwig, K. A., Sogin, M. L., Andreishcheva, E. N., Kelley, D. S., Shen, C. C., et al. (2010a). Archaea and bacteria with surprising microdiversity show shifts in dominance over 1,000-year time scales in hydrothermal chimneys. *Proc. Natl. Acad. Sci. U.S.A.* 107, 1612–1617. doi: 10.1073/pnas.0905369107

Brazelton, W. J., Mehta, M. P., Kelley, D. S., and Baross, J. A. (2011). Physiological differentiation within a single-species biofilm fueled by serpentinization. *mBio* 2, e00127-11. doi: 10.1128/mBio.00127-11

Brazelton, W. J., Nelson, B., and Schrenk, M. O. (2012). Metagenomic evidence for H₂ oxidation and H₂ production by serpentinite-hosted subsurface microbial communities. *Front. Microbiol.* 2:268. doi: 10.3389/fmicb.2011.00268

Brazelton, W. J., Schrenk, M. O., Kelley, D. S., and Baross, J. A. (2006). Methane- and sulfur-metabolizing microbial communities dominate the Lost City Hydrothermal Field Ecosystem. *Appl. Environ. Microb.* 72, 6257–6270. doi: 10.1128/aem.00574-06

Brazelton, W. J., Sogin, M. L., and Baross, J. A. (2010b). Multiple scales of diversification within natural populations of archaea in hydrothermal chimney biofilms. *Environ. Microbiol. Rep.* 2, 236–242. doi: 10.1111/j.1758-2229.2009.00097.x

Daae, F. L., Økland, I., Dahle, H., Jørgensen, S. L., Thorseth, I. H., and Pedersen, R. B. (2013). Microbial life associated with low-temperature alteration of ultramafic rocks in the Leka ophiolite complex. *Geobiology* 11, 318–339. doi: 10.1111/gbi.12035

Delong, E. F. (2003). Oceans of Archaea. *ASM NEWS* 69, 503–511.

Doddema, H. J., and Vogels, G. D. (1978). Improved identification of methanogenic bacteria by fluorescence microscopy. *Appl. Environ. Microb.* 36, 752–754.

Edgar, R. C. (2004). MUSCLE: multiple sequence alignment with high accuracy and high throughput. *Nucleic Acids Res.* 32, 1792–1797. doi: 10.1093/nar/gkh340

Fogel, G. B., Collins, C. R., Li, J., and Brunk, C. F. (1999). Prokaryotic genome size and SSU rDNA copy number, estimation of microbial relative abundance from a mixed population. *Microb. Ecol.* 38, 93–113.

Galand, P. E., Casamayor, E. O., Kirchman, D. L., and Lovejoy, C. (2009). Ecology of the rare microbial biosphere of the Arctic Ocean. *Proc. Natl. Acad. Sci. U.S.A.* 106, 22427–22432. doi: 10.1073/pnas.0908284106

Galtier, N., Tourasse, N., and Gouy, M. (1999). A nonhyperthermophilic common ancestor to extant life forms. *Science* 283, 220–221. doi: 10.1126/science.283.5399.220

Gerasimchuk, A., Shatalov, A., Novikov, A., Butorova, O., Pimenov, N., Lein, A., et al. (2010). The search for sulfate-reducing bacteria in mat samples from the Lost City hydrothermal field by molecular cloning. *Microbiology* 79, 96–105. doi: 10.1134/S0026261710010133

Gevers, D., Cohan, F. M., Lawrence, J. G., Spratt, B. G., Coenye, T., Feil, E. J., et al. (2005). Re-evaluating prokaryotic species. *Nat. Rev. Microbiol.* 3, 733–739. doi: 10.1038/nrmicro1236

Gillan, D. C., Danis, B., Pernet, P., Joly, G., and Dubois, P. (2005). Structure of sediment-associated microbial communities along a heavy-metal contamination gradient in the marine environment. *Appl. Environ. Microb* 71, 679–690. doi: 10.1128/AEM.71.2.679-690.2005

Good, I. J. (1953). The population frequencies of species and the estimation of population parameters. *Biometrika* 40, 237–264. doi: 10.1093/biomet/40.3.4.237

Huber, J. A., Butterfield, D. A., and Baross, J. A. (2006). Diversity and distribution of subseafloor *Thermococcales* populations in diffuse hydrothermal vents at an active deep-sea volcano in the northeast Pacific Ocean. *J. Geophys. Res.* 111, G04016. doi: 10.1029/2005JG000097

Hungate, R. E. (1969). A roll tube method for the cultivation of strict anaerobes. *Meth. Microbiol.* 3B, 117–132.

Iverson, V., Morris, R. M., Frazar, C. D., Berthiaume, C. T., Morales, R. L., and Armbrust, E. V. (2012). Untangling genomes from metagenomes: revealing an uncultured class of marine *Euryarchaeota*. *Science* 335, 587–590. doi: 10.1126/science.1212665

Joulian, C., Ramsing, N. B., and Ingvorsen, K. (2001). Congruent phylogenies of most common small-subunit rRNA and dissimilatory sulfite reductase gene sequences retrieved from estuarine sediments. *Appl. Environ. Microbiol.* 67, 3314–3318. doi: 10.1128/AEM.67.7.3314-3318.2001

Kelley, D. S., Karson, J. A., Blackman, D. K., Früh-Green, G. L., Butterfield, D. A., Lilley, M. D., et al. (2001). An off-axis hydrothermal vent field near the Mid-Atlantic Ridge at 30°N. *Nature* 412, 145–149. doi: 10.1038/35084000

Kelley, D. S., Karson, J. A., Fruh-Green, G. L., Yoerger, D. R., Shank, T. M., Butterfield, D. A., et al. (2005). A serpentinite-hosted ecosystem: the Lost City Hydrothermal Field. *Science* 307, 1428–1434. doi: 10.1126/science.1102556

Khachane, A. N., Timmis, K. N., and Dos Santos, V. A. P. M. (2005). Uracil content of 16S rRNA of thermophilic and psychrophilic prokaryotes correlates inversely with their optimal growth temperatures. *Nucleic Acids Res.* 33, 4016–4022. doi: 10.1093/nar/gki714

Kimura, H., Ishibashi, J. I., Masuda, H., Kato, K., and Hanada, S. (2007). Selective phylogenetic analysis targeting 16S rRNA genes of hyperthermophilic Archaea in the deep-subsurface hot biosphere. *Appl. Environ. Microbiol.* 73, 2110–2117. doi: 10.1128/aem.02800-06

Kimura, H., Mori, K., Tashiro, T., Kato, K., Yamanaka, T., Ishibashi, J., et al. (2010). Culture-independent estimation of optimal and maximum growth temperatures of archaea in subsurface habitats based on the G+C content in 16S rRNA gene sequences. *Geomicrobiol. J.* 27, 114–122. doi: 10.1080/01490450903546699

Kimura, H., Mori, K., Yamanaka, T., and Ishibashi, J. (2013). Growth temperatures of archaean communities can be estimated from the guanine-plus-cytosine contents of 16S rRNA gene fragments. *Environ. Microbiol. Rep.* 5, 468–474. doi: 10.1111/1758-2229.12035
Kimura, H., Sugihara, M., Kato, K., and Hanada, S. (2006). Selective phylogenetic analysis targeted at 16S rRNA genes of thermophiles and hyperthermophiles in deep-subsurface geothermal environments. *Appl. Environ. Microb.* 72, 21–27. doi: 10.1128/aem.72.1.21-27.2006

Kimura, M. (1980). A simple method for estimating evolutionary rates of base substitutions through comparative studies of nucleotide sequences. *J. Mol. Evol.* 16, 111–120.

Klappenbach, J. A., Saxman, P. R., Cole, J. R., and Schmidt, T. M. (2001). *rrndb*: the ribosomal RNA operon copy number database. *Nucleic Acids Res.* 29, 181–184. doi: 10.1093/nar/29.1.181

Klein, M., Friedrich, M., Roger, A. J., Hugenholtz, P., Fishbain, S., Abicht, H., et al. (2001). Multiple lateral transfers of dissimilatory sulfite reductase genes between major lineages of sulfate-reducing prokaryotes. *J. Bact.* 183, 6028–6035. doi: 10.1128/jb.183.20.6028-6035.2001

Knittel, K., Boetius, A., Lemke, A., Eilers, H., Lochte, K., Pfannkuche, O., et al. (2003). Activity, distribution, and diversity of sulfate reducers and other bacteria in sediments above gas hydrate (Cascadia Margin, OR). *Geomicrobiol. J.* 20, 269–294. doi: 10.1080/01490450308396

Kuever, J., Galushko, A., and Drzyzga, O. (2001). Reclassification of *Desulfobacterium phenolicum* as *Desulfobacula phenolica* comb. nov. and description of strain *SaxT* as *Desulfotignum balticum* gen. nov., sp. nov. *Int. J. Syst. Evol. Microbiol.* 51, 171–177. doi: 10.1099/00207713-51-1-171

Lin, L.-H., Hall, J., Onstott, T. C., Gihring, T., Sherwood Lollar, B., Boice, E., et al. (2006). Planktonic microbial communities associated with fracture-derived groundwater in a deep gold mine of South Africa. *Geomicrobiol. J.* 23, 475–497. doi: 10.1080/01490450600875829

Lösekann, T., Knittel, K., Nadalig, T., Fuchs, B., Niemann, H., Boetius, A., et al. (2007). Diversity and abundance of aerobic and anaerobic methane oxidizers at the Haakon Mosby Mud volcano, Barents Sea. *Appl. Environ. Microb.* 73, 3348–3362. doi: 10.1128/aem.000016-07

Luton, P. E., Wayne, J. M., Sharp, R. J., and Riley, P. W. (2002). The mcrA gene as an alternative to 16S rRNA in the phylogenetic analysis of methanogen populations in landfill. *Microbiology* 148, 3521–3530. doi: 10.1099/00221287-148-11-3521

Maurizot, P., and Vendé-Leclerc, M. (2009). *Carte Géologique* de la Nouvelle-Calédonie au 1/500 000. Nouméa: Direction de l'Industrie, des Mines et de l'Energie, Service de la Géologie de Nouvelle-Calédonie, Bureau de Recherches Géologiques et Minières.

McCollom, T. M. (2007). Geochemical constraints on sources of metabolic energy for chemolithoautotrophy in ultramafic-hosted deep-sea hydrothermal systems. *Astrobiology* 7, 933–950. doi: 10.1089/ast.2006.0119

Mei, N., Zergane, N., Postec, A., Erauso, G., Ollier, A., Payri, C., et al. (2014). Fermentative hydrogen production by a new alkaliphilic *Clostridium* sp. (strain PROH2) isolated from a shallow submarine hydrothermal chimney in Prony Bay, New Caledonia. *Int. J. Hydrogen Energ.* 39, 19465–19473. doi: 10.1016/j.ijhydene.2014.09.111

Merkel, A. Y., Huber, J. A., Chernyh, N. A., Bonch-Osmolovskaya, E. A., and Lebedinsky, A. V. (2013). Detection of putatively thermophilic anaerobic methanotrophs in diffuse hydrothermal vent fluids. *Appl. Environ. Microb.* 79, 915–923. doi: 10.1128/aem.03034-12

Michotey, V., Guasco, S., Boeuf, D., Morezzi, N., Durieux, B., Charpy, L., et al. (2012). Spatio-temporal diversity of free-living and particle-attached prokaryotes in the tropical lagoon of Ahe atoll (Tuamotu Archipelago) and its surrounding oceanic waters. *Mar. Poll. Bull.* 65, 525–537. doi: 10.1016/j.marpolbul.2012.01.009

Miroshnichenko, M., Hippe, H., Stackebrandt, E., Kostrikina, N., Chernyh, N., Jeanthon, C., et al. (2001). Isolation and characterization of *Thermococcus sibiricus* sp. nov. from a Western Siberia high-temperature oil reservoir. *Extremophiles* 5, 85–91. doi: 10.1007/s0079200100175

Monnin, C., Chavagnac, V., Boulart, C., Ménez, B., Gérard, M., Gérard, E., et al. (2014). Fluid chemistry of the low temperature hyperalkaline hydrothermal system of Prony Bay (New Caledonia). *Biogeosciences* 11, 5687–5706. doi: 10.5194/bg-11-5687-2014

Muyzer, G., and Stams, A. J. (2008). The ecology and biotechnology of sulphate-reducing bacteria. *Nat. Rev. Microbiol.* 6, 441–454. doi: 10.1038/nrmicro1892

Nei, M., and Kumar, S. (2000). *Molecular Evolution and Phylogenetics*. New York, NY: Oxford University Press.

Neuner, A., Jannasch, H., Belkin, S., and Stetter, K. (1990). *Thermococcus litoralis* sp. nov. a new species extremely thermophilic marine archae bacteria. *Arch. Microbiol.* 153, 205–207.

Parkes, R. J., Cragg, B., Roussel, E., Webster, G., Weightman, A., and Sass, H. (2014). A review of prokaryotic populations and processes in sub-seafloor sediments, including biosphere:geosphere interactions. *Mar. Geol.* 352, 409–425. doi: 10.1016/j.margeo.2014.02.009

Pasini, V., Brunelli, D., Dumas, P., Sandt, C., Frederick, J., Benzerara, K., et al. (2013). Low temperature hydrothermal oil and associated biological precursors in serpentinites from Mid-Ocean Ridge. *Lithos* 178, 84–95. doi: 10.1016/j.lithos.2013.06.014

Pelletier, B., Chevillon, C., Menou, J., Butscher, J., Folcher, E., Geoffray, C., et al. (2006). Plongées, forage et cartographie Baie du Prony et Banc Gail, lagon Sud de Nouvelle-Calédonie, campagne 2005-NC-PL du N.O. ALIS 13-17 Juin 2005 et cartographie baie du Prony et canal Woodin N.O. ALIS 25-26 September 2005: Nouméa IRD Sept. 2006. *Missions Sci. Terre, Géol-Géophys* 70, 44.

Pisapia, C., Gérard, E., Gérard, M., Ménez, B., and HYDROPRONY Scientific Party. (2013). Organomineralization drives early chimney edification at the hyperalkaline hydrothermal field of the Prony Bay (New Caledonia). *Mineralogical Magazine* 77, 1974.

Proskurowski, G., Lilley, M. D., Seewald, J. S., Früh-Green, G. L., Olson, E. J., Lupton, J. E., et al. (2008). Abiogenic hydrocarbon production at Lost City Hydrothermal Field. *Science* 319, 604–607. doi: 10.1126/science.1151194

Quéméneur, M., Bes, M., Postec, A., Mei, N., Hamelin, J., Monnin, C., et al. (2014). Spatial distribution of microbial communities in the shallow submarine alkaline hydrothermal field of the Prony Bay, New Caledonia. *Environ. Microbiol. Rep.* 6, 665–674. doi: 10.1111/1758-2229.12184

Roussel, E. G., Sauvadet, A. L., Chaduteau, C., Fouquet, Y., Charlou, J. L., Prieur, D., et al. (2009). Archaeal communities associated with shallow to deep subseafloor sediments of the New Caledonia Basin. *Environ. Microbiol.* 11, 2446–2462. doi: 10.1111/j.1462-2920.2009.01976.x

Russell, M. J., Hall, A. J., and Martin, W. (2010). Serpentinization as a source of energy at the origin of life. *Geobiology* 8, 355–371. doi: 10.1111/j.1472-4669.2010.00249.x

Saitou, N., and Nei, M. (1987). The neighbor-joining method: a new method for reconstructing phylogenetic trees. *Mol. Biol. Evol.* 4, 406–425.

Schleper, C., and Nicol, G. W. (2010). "Ammonia-oxidising archaea – Physiology, ecology and evolution," in *Advances Microbial Physiol*, Vol. 57, ed K. P. Robert (Oxford, UK: Academic Press), 1–41.

Schloss, P. D., Westcott, S. L., Ryabin, T., Hall, J. R., Hartmann, M., Hollister, E. B., et al. (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. *Appl. Environ. Microb* 75, 7537–7541. doi: 10.1128/aem.01541-09

Schreiber, L., Holler, T., Knittel, K., Meyerdierks, A., and Amann, R. (2010). Identification of the dominant sulfate-reducing bacterial partner of anaerobic methanotrophs of the ANME-2 clade. *Environ. Microbiol.* 12, 2327–2340. doi: 10.1111/j.1462-2920.2010.02275.x

Schrenk, M. O., Brazelton, W. J., and Lang, S. Q. (2013). Serpentinization, carbon, and deep life. *Rev. Mineral. Geochem.* 75, 575–606. doi: 10.2138/rmg.2013.75.18

Schrenk, M. O., Kelley, D. S., Bolton, S. A., and Baross, J. A. (2004). Low archaeal diversity linked to subseafloor geochemical processes at the Lost City Hydrothermal Field, Mid-Atlantic Ridge. *Environ. Microbiol.* 6, 1086–1095. doi: 10.1111/j.1462-2920.2004.00650.x

Schulte, M., Blake, D., Hoehler, T., and McCollom, T. (2006). Serpentinization and its implications for life on the early Earth and Mars. *Astrobiology* 6, 364–376. doi: 10.1089/ast.2006.6.364

Sorokin, D., Tourova, T., Mußmann, M., and Muyzer, G. (2008). Dethiobacter alkaliphilus; gen. nov. sp. nov., and Desulfurivibrio alkaliphilus; gen. nov. sp. nov.: two novel representatives of reductive sulfur cycle from soda lakes. *Extremophiles* 12, 431–439. doi: 10.1007/s00792-008-0148-8

Sowers, K. R., Baron, S. F., and Ferry, J. G. (1984). *Methanosarcina acetivorans* sp-nov, an acetotrophic methane-producing bacterium isolated from marine sediments. *Appl. Environ. Microbiol.* 47, 971–978.

Suzuki, S., Ishii, S. I., Wu, A., Cheung, A., Tenney, A., Wanger, G., et al. (2013). Microbial diversity in The Cedars, an ultrabasic, ultraredeucing, and low salinity serpentinizing ecosystem. *Proc. Natl. Acad. Sci. U.S.A.* 110, 15336–15341. doi: 10.1073/pnas.1302426110
Suzuki, S., Kuenen, J. G., Schipper, K., Van Der Velde, S., Ishii, S., Wu, A., et al. (2014). Physiological and genomic features of highly alkaliphilic hydrogen-utilizing *Betaproteobacteria* from a continental serpentinizing site. *Nat. Commun.* 5, 3900. doi: 10.1038/ncomms4900

Takai, K. H., Gamo, T., Tsunogai, U., Nakayama, N., Hirayama, H., Nealson, K. H., et al. (2004). Geochemical and microbiological evidence for a hydrogen-based, hyperthermophilic subsurface lithoautotrophic microbial ecosystem (HyperSLiME) beneath an active deep-sea hydrothermal field. *Extremophiles* 8, 269–282. doi: 10.1007/s00792-004-0386-3

Takai, K., Moser, D. P., DeFlaun, M., Onstott, T. C., and Fredrickson, J. K. (2001). Archaeal diversity in waters from deep South African gold mines. *Appl. Environ. Microbiol.* 67, 5750–5760. doi: 10.1128/aem.67.21.5750-5760.2001

Tamura, K., Nei, M., and Kumar, S. (2004). Prospects for inferring very large phylogenies by using the neighbor-joining method. *Proc. Natl. Acad. Sci. U.S.A.* 101, 11030–11035. doi: 10.1073/pnas.0404206101

Tamura, K., Stecher, G., Peterson, D., Filipski, A., and Kumar, S. (2013). MEGA6: molecular evolutionary genetics analysis version 6.0. *Mol. Biol. Evol.* 30, 2725–2729. doi: 10.1093/molbev/mst197

Thauer, R. K. (1998). Biochemistry of methanogenesis: a tribute to Marjory Stephenson. *Microbiology* 144, 2377–2406.

Tiago, I., Chung, A. P., and Veríssimo, A. (2004). Bacterial diversity in a nonsaline alkaline environment: heterotrophic aerobic populations. *Appl. Environ. Microbiol.* 70, 7378–7387. doi: 10.1128/aem.70.12.7378-7387.2004

Tiago, I., and Veríssimo, A. (2013). Microbial and functional diversity of a subterrestrial high pH groundwater associated to serpentinization. *Environ. Microbiol.* 15, 1687–1706. doi: 10.1111/1462-2920.12034

Wang, H. C., Xia, X., and Hickey, D. (2006). Thermal adaptation of the small subunit ribosomal RNA gene: a comparative study. *J. Mol. Evol.* 63, 120–126. doi: 10.1007/s00239-005-0255-4

Wasmund, K., Schreiber, L., Lloyd, K. G., Petersen, D. G., Schramm, A., Stepanauskas, R., et al. (2013). Genome sequencing of a single cell of the widely distributed marine subsurface *Dehalococcoidia*, phylum *Chloroflexi*. *ISME J.* 8, 383–397. doi: 10.1038/ismej.2013.143

Widdel, F., Kohring. G. W., and Mayer, F. (1983). Studies on dissimilatory sulfate-reducing bacteria that decompose fatty acids. III. Characterization of the filamentous gliding *Desulfonema limicola* gen. nov., and Desulfonema magnum sp. nov. *Arch. Microbiol.* 134, 286–294.

Zhang, Y., Zhao, Z., Chen, C. T., Tang, K., Su, J., and Jiao, N. (2012). Sulfur metabolizing microbes dominate microbial communities in andesite-hosted shallow-sea hydrothermal systems. *PLoS ONE* 7:e44593. doi: 10.1371/journal.pone.0044593

Zhilina, T. N., Zavarzina, D. G., Kuever, J., Lysenko, A. M., and Zavarzin, G. A. (2005). Desulfonatronum cooperativum sp. nov., a novel hydrogenotrophic, alkaliphilic, sulfate-reducing bacterium, from a syntrophic culture growing on acetate. *Int. J. Syst. Evol. Microbiol.* 55, 1001–1006. doi: 10.1099/ijs.0.63490-0

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2015 Postec, Quéméneur, Bes, Mei, Benaïssa, Payri, Pelletier, Monnin, Guañtas-Dombrowsky, Ollivier, Gérard, Pisapia, Gérard, Ménez and Erauso. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.