Abstract
The East Lau Spreading Center (ELSC) and Valu Fa Ridge (VFR) comprise a ridge segment in the southwest Pacific Ocean where rapid transitions in the underlying mantle chemistry manifest themselves as gradients in seafloor rock geochemistry. We studied the geology and microbial diversity of three silicate rock samples and three inactive sulfide chimney samples collected, from north to south, at the vent fields Kilo Moana, ABE, Tui Malila, and Mariner. This is the first study of microbial populations on basaltic andesite, which was sampled at Mariner vent field. Silicate rock geochemistry exhibits clear latitudinal trends that are mirrored by changes in bacterial community composition. α-proteobacteria, ε-proteobacteria, and Bacteroidetes are most common on a silicate collected from Kilo Moana and their proportions decrease linearly on silicates collected further south. Conversely, a silicate from Mariner vent field hosts high proportions of a unique lineage of Chloroflexi unrelated (<90% sequence similarity) to previously recovered environmental clones or isolates, which decrease at ABE and are absent at Kilo Moana. The exteriors of inactive sulfide structures are dominated by lineages of sulfur oxidizing α-proteobacteria, γ-proteobacteria, and ε-proteobacteria, while the interior of one chimney is dominated by putative sulfur-reducing δ-proteobacteria. A comparison of bacterial communities on inactive sulfides from this and previous studies reveals the presence of a clade of uncultured Bacteroidetes exclusive to sulfidic environments, and a high degree of heterogeneity in bacterial community composition from one sulfide structure to another. In light of the heterogeneous nature of bacterial communities observed here and in previous studies of both active and inactive hydrothermal sulfide structures, the presence of numerous niches may be detected on these structures in the future by finer scale sampling and analysis.
Introduction
The Eastern Lau Spreading Center (ELSC) and Valu Fa Ridge (VFR) comprise the southern portion of the Lau Basin back-arc spreading center, located between the islands of Samoa and Tonga in the southwest Pacific Ocean. Hydrothermal venting was discovered in Lau Basin along the VFR in 1989 (Fouquet et al., ) and was subsequently discovered at various other sites along the ELSC and VFR (Langmuir et al., ; Ishibashi et al., ). The geochemistry of both host rock composition and hydrothermal fluids in Lau Basin changes along the north-south gradient from mid-ocean ridge-like basalt in the north to subduction influenced andesite in the south over the course of <600 km (Escrig et al., ; Dunn and Martinez, ; Mottl et al., ). This is a steeper gradient than seen anywhere else along the global mid-ocean ridge system.
Given the observed gradients in host rock composition and vent fluid geochemistry, Lau Basin is an ideal location to study how geochemistry influences the distribution and composition of biological communities. It has been shown that the gradient in chemistry influences distributions of megafauna (Podowski et al., 2010) and microbes on active sulfide chimneys (Flores et al., ), but no work currently exists examining microbiology on the host rock or inactive hydrothermal sulfides. Seafloor-exposed silicates (thus far basalts are the only seafloor silicates sampled) are known to host diverse microbial communities at mid-ocean ridges (Lysnes et al., ; Mason et al., ; Santelli et al., 2009), but no samples have been analyzed from back-arc systems, where geochemical controls on silicate microbiology can be tested explicitly. Recent work has also illustrated that microbial communities thrive on inactive hydrothermal sulfides long after venting ceases (Rogers et al., 2003; Suzuki et al., 2004; Kato et al., ; Sylvan et al., 2012a). A succession occurs on these sulfides whereby the microbial community present on inactive sulfides is different from that on active structures. This is likely due to a change in mineralogy brought on by the drastic decrease in temperature and disappearance of the reduced substrates in hydrothermal fluids once the vent dies. Similar mineralogical controls are present on other low temperature deep-sea deposits, including basalts, where the microbial community appears to be selected by substrate type (Toner et al., 2013); bacterial communities collected from seafloor basalts are more similar to each other than to communities on other substrates. The same is true for bacterial communities on inactive sulfides from the East Pacific Rise (EPR), the Okinawa Trough, and Indian Ocean Ridge (Suzuki et al., 2004; Toner et al., 2013). Both seafloor basalts and inactive hydrothermal sulfide structures can be considered extreme environments due to their high metal content and elevated concentrations of elements considered toxic to most life, such as copper.
We examined the geology and microbiology of low temperature deposits collected from four vent fields along the ELSC and VFR during summer 2009. Specifically, we seek to test the hypothesis that bacterial communities are selected by gradients in host rock composition on silicates [basalts and basaltic andesite, classified using total alkalis versus silica (Le Bas and Streckeisen, )] and to determine if bacterial communities on basaltic andesite, for which no data currently exists, differ from those on basalt. We also seek to further understand bacterial communities on inactive sulfides, for which limited data exists.
Materials and Methods
Sample collection
Five seafloor rock samples were collected from the ELSC and VFR (Figure 1) during cruise TN-235 on the R/V Thompson with ROV Jason II during 16 May–08 June 2009. Silicates located outside areas of diffuse flow and sulfides that appeared inactive on the seafloor were collected from four vent fields (Figure A1 in Appendix; Table 1). Samples were collected using ROV Jason II and then placed in sealed bioboxes on the sampling tray of the ROV for the remainder of the dive, isolating them in seawater from the collection site for the remainder of the dive. Once on deck, the samples were removed from the bioboxes and immediately processed in flame sterilized steel boxes, where they were separated from the seawater in the bioboxes, with a flame sterilized hammer and chisel. Aliquots of rock chips removed from the larger sample for DNA analysis were placed in 5 mL centrifuge tubes and immediately frozen at −80°C. Subsamples for optical mineralogy (via thin section analysis), X-ray diffraction (XRD), and elemental analysis were allowed to air dry on the ship and were subsequently stored at room temperature.
Figure 1
Table 1
| Sample | Vent field | Rock type | Rock mineralogy (thin section) | Rock mineralogy (XRD) | Bacterial cells g−1 rock* (% Bacteria versus Archaea) | Archaeal cells g−1 rock* |
|---|---|---|---|---|---|---|
| KiMba | Kilo Moana | Silicate; basalt | Augite (<5%), Plagioclase (5–10%), Aegirine-augite (<5%), amorphous glass, not porous | Major species: diopside, aegirine-augite Minor species: augite | 4.14 ± 0.06 × 105 (93.6%) | 2.82 ± 0.69 × 104 |
| ABEba | ABE | Silicate; transitional | Augite (35%), Plagioclase (20%), moderately altered, not porous | Major species: albite, augite Minor species: labradorite | 5.20 ± 0.12 × 104 (99.7%) | 1.49 ± 0.27 × 102 |
| Marba | Mariner | Silicate; basaltic andesite | Augite (5%), Plagioclase (5–10%), little to no alteration, 40% porous | Signal too amorphous | 9.97 ± 1.12 × 103 (15.4%) | 5.47 ± 1.86 × 104 |
| ABEsO** | ABE | Outside of an inactive sulfide | Poorly crystalline, pyrite, sphalerite | Sphalerite, pyrite, barite | 9.92 ± 0.34 × 107 (100%) | Below detection |
| ABEsIN | ABE | Inside of an inactive sulfide | nd | nd | 1.15 ± 0.04 × 107 (99.8%) | 2.87 ± 0.18 × 104 |
| TuiMs | Tui Malila | Inactive sulfide | Cinnabar, pyrite, sphalerite, anhydrite | Barite, sphalerite, pyrite | 3.65 ± 0.10 × 107 (100%) | Below detection |
Samples descriptions and bacterial biomass estimates.
*As determined by qPCR, ±Standard Error.
**Mineralogical analysis was completed for the bulk sample for ABEs, but DNA analysis was completed on two separate samples from the same chimney, one from the exterior and one from the interior conduit.
Thin section analysis
Standard thin sections (30 μm thick) were made from representative portions of each sample by Spectrum Petrographics, Inc. (Vancouver, WA, USA) and were analyzed using a polarizing petrographic microscope (Zeiss AxioImager.M2m with a Zeiss AxioCam HRc camera) with both reflected and transmitted light.
X-ray diffraction
Approximately 2 cm3 of each sample was ground to a powder using mortar and pestle and then mounted on a plastic holder using acetone and a disposable wooden stick to ensure random orientation of the particles. After air drying, samples were analyzed on a Siemens d-500 Diffractometer with a cobalt source. The software program JADE (Materials Data, Inc., v9.3.3) was used for phase identification by peak matching.
Elemental analysis
Major and minor elemental concentrations were determined from the same subsample used for XRD analysis at the University of Minnesota Analytical Geochemistry Lab. Powders were acid digested prior to analysis. Weight percent of major oxides were determined in triplicate on a Thermo Scientific iCAP 6500 dual view Inductively Coupled Plasma – Optical Emission Spectrometer (ICP-OES). Samples were diluted 40-fold prior to analysis with the addition of a Cs matrix modifier and Y as an internal standard and measured using free aspiration of the sample and integrations of five 10 s replicate readings per measurement. Trace elements were determined in duplicate on a Thermo Scientific XSERIES 2 ICP mass spectrometer with an electrospray ionization PC3 Peltier cooled spray chamber, SC-FAST injection loop, and SC-4 autosampler.
DNA extraction and analysis
DNA was extracted from ∼4 cm3 of sample using a CTAB phenol/chloroform extraction (Ausubel et al., ). qPCR for bacteria was carried out as described previously using primers 338f [5′-ACT CCT ACG GGA GGC AGC AG-3′) and 518r (5′-ATT ACC GCG GCT GCT GG-3′ (Einen et al., )]. qPCR for archaea was carried out using primers 806f [5′-ATT AGA TAC CCS BGT AGT-3′ (Takai and Horikoshi, 2000)] and 922r [5′-YCC GGC GTT GAN TCC AAT T-3′ (Delong, )]. For both bacterial and archaeal qPCR, 16S rRNA copy numbers g−1 rock were calculated by multiplying the mean copy number detected from triplicate reactions by the dilution factor (total DNA extraction volume divided by template volume in each qPCR reaction) and divided by the weight of the sample from which DNA was extracted, in grams. The cell number per gram rock was determined by assuming 3.9 16S rRNA gene copies per cell for bacteria and 1.8 16S rRNA gene copies per cell for archaea (Einen et al., ). The thermal program employed for both bacterial and archaeal qPCR primer sets was: 10 min at 95°C followed by 45 cycles of 30 s at 95°C, 30 s at 55°C and 25 s at 72°C. Negative controls [polymerase chain reaction (PCR) water included as a template] were included for all qPCR runs, and melt curves for all qPCR products were checked to ensure a single PCR product was generated. Reported values in Table 1 were all greater than the negative control and no samples yielded multiple PCR products. The reported qPCR reactions were run in triplicate.
Polymerase chain reaction of the 16S rRNA gene and subsequent cloning and sequencing of the product was carried out according to Sylvan et al. (2012b). Briefly, universal bacterial primers 27F (5′-GAG TTT GAT CCT GGC TCA G-3′) and 1492R (5′-RGY TAC CTT GTT ACG ACT T-3′) were used for PCR. Three reactions were combined and run out on an agarose gel and then excised and extracted using the Zymoclean Gel DNA Recovery Kit (Zymo Research, Irvine, CA, USA). DNA from the extracted PCR band was cloned into the pCR 4 TOPO vector using the TOPO TA Cloning Kit (Invitrogen, Grand Island, NY, USA) and transformants were plated on LB + 100 μg mL−1ampicillin according to the manufacturer’s instructions. Clones were sequenced at the Beckman Coulter Genomics center in Danvers, MA, USA.
16S rRNA contigs were generated using Geneious v5.6 (Drummond et al., ). Edited near full-length 16S rDNA sequences were classified and checked for chimeras using the Bellerophon tool of Greengenes (DeSantis et al., ). The resulting sequences were aligned using the Greengenes NAST server (DeSantis et al., ) and imported into ARB (Ludwig et al., ) for selection of sequences to include in phylogenetic trees. Closely related cultured strains to contigs were identified using the “Named Isolates” option with the BLAST function on the Greengenes website (http://greengenes.lbl.gov/cgi-bin/nph-blast_interface.cgi). If no closely related isolate to a contig existed, the nearest neighbor from the ARB database was identified. Closely related isolates and/or sequences were aligned with sequences from this study using MEGA 5 (Tamura et al., 2011). Phylogenetic trees were constructed following manual adjustment of this alignment using both the neighbor-joining method, based on the maximum composite likelihood model and gamma distribution, and maximum likelihood analysis, based on the Jukes–Cantor model with a Gamma distribution. Both types of phylogenies were tested using 1000 bootstrap replicates. Calculation of rarefaction curves and diversity estimates, as well as comparison between clone libraries and those of other studies, was carried out using the software Mothur (Schloss et al., 2009). Pre-clustering (Huse et al., ) and the average neighbor clustering algorithm were used to generate distance matrices from which rarefaction curves, estimates of shared richness (Jclass, Jest), and estimates of shared structure (theta Yue–Clayton, or ΘYC, and Bray–Curtis) were calculated. An operational taxonomic unit (OTU) cutoff of 97% was used for the rarefaction curves and an OTU cutoff of 95% was used for generating a cladogram for inter-sample comparison. In agreement with other prior work (Toner et al., 2013), we found that using a 95% cutoff to build cladograms yielded clearer results than a 97% cutoff. Using a 95% cutoff for OTUs, however, had little impact on the overall number of OTUs and we therefore maintained the commonly accepted value of 97% for OTUs. The cladograms were converted to a circular tree layout in Genious v5.6 (Drummond et al., ) and edited in Adobe Illustrator CS6.
DNA sequences generated for this project were deposited in the National Center for Biotechnology Information (NCBI) database under accession numbers KC682512–KC682862. For a few highly represented sequences, where multiple nearly identical clones were generated from a single sample, only one clone was deposited. They are listed here, followed by the number of additional clones they represent: ABEsO_A7 (25 additional clones), Marba_A2 (09 additional clones), ABEsIN_H1 (15 additional clones), ABEsIN_A1 (47 additional clones), and TuiMs_A6 (24 additional clones).
Results
Sample descriptions, silicates
Sample KiMba is a basalt collected from Kilo Moana vent field with a visible glassy rim underlain by ground mass. We sampled the glassy rim. The thin section revealed that KiMba is porphyritic with large phenocrysts in an amorphous glassy matrix that are encroached by spherules. KiMba is not porous, and plagioclase [(Na,Ca)(Si,Al)4O8], aegirine-augite [(Ca,Na)(Mg,Fe2+,Fe3+)(Si2O6)], and augite [(Ca,Na)(Mg,Fe,Al,Ti)(Si,Al)2O6] are represented in the thin section (Figure 2). XRD analysis revealed the presence of diopside [CaMg(Si2O6)] in this sample (Table 1). Sample ABEba is a basalt collected from ABE vent field with an oxidized rim. It was not porous and is rich in augite and plagioclase (Figure 2). The XRD pattern revealed the presence of albite (NaAlSi3O8), augite and labradorite [(Na,Ca)(Si,Al)4O8]. Sample Marba is basaltic andesite collected from Mariner vent field. It is very porous (Figure 2) with little to no alteration. Plagioclase and augite were detected by thin section analysis, but the XRD signal was too weak to identify specific minerals.
Figure 2
Inactive sulfides
Sample ABEs was broken off an inactive sulfide chimney located in ABE vent field. It is poorly crystalline and slightly oxidized. Light microscopy of samples in thin section revealed the presence of pyrite (FeS2) and sphalerite (ZnS; data not shown); and barite (BaSO4) was identified in the XRD pattern. A portion of this sample was collected from the outside wall of the chimney (ABEsO) and a separate section was sampled from the inside (ABEsIN) for analysis of the microbial community. Sample TuiMs is a zoned chimney that appeared inactive (no venting or shimmering fluid observed exiting the chimney) on the seafloor at the time of sampling from the Tui Malila vent field. Sulfide minerals are present on the outer rim of the thin section and sulfide and sulfate minerals are present away from the outer rim (Figure 3). On the outer tannish red to deep maroon rim, euhedral cinnabar (HgS) coats and encrusts anhedral pyrite. The interior edge of the thin section is characterized by anhydrite (CaSO4), sphalerite, and minor anhedral pyrite. XRD analysis reveals the presence of barite, sphalerite, and pyrite.
Figure 3
Geochemistry
Elemental analysis confirms transitions in host rock composition from north to south in the ELSC. Weight percent of Ba, K, P, Si, and Sr oxides and concentrations of V, Cu, Rb, Sr, Ba, La, Ce, Pr, Nd, and Th all increase from north to south in the silicates (Tables 2 and 3). Concentrations of Sc and Co follow the opposite trend. Based on analysis of total alkalis versus silica, KiMba, and ABEba are both sub-alkaline basalts while Marba is sub-alkaline basaltic andesite (data not shown).
Table 2
| Oxide | KiMba | ABEba | Marba | ABEs | TuiMs |
|---|---|---|---|---|---|
| Al2O3 | 14.76 | 14.5867 | 15.1733 | 0.1282 | 0.15713 |
| BaO | 0.0008 | 0.0148 | 0.0203 | 1.8943 | 3.175 |
| CaO | 12.0467 | 8.423 | 9.1697 | 0.2079 | 1.7117 |
| Fe2O3 | 11.22 | 12.9633 | 11.9967 | 11.9 | 2.8467 |
| K2O | 0.0366 | 0.2393 | 0.2588 | 0.0455 | 0.0514 |
| MgO | 7.7087 | 4.098 | 4.7983 | 0.0414 | 0.047 |
| MnO | 0.1839 | 0.3112 | 0.208 | 0.1575 | 0.0134 |
| Na2O | 2.106 | 3.085 | 2.9113 | 0.4025 | 0.1844 |
| P2O5 | 0.0787 | 0.1402 | 0.15 | 0.0179 | 0.0025 |
| SiO2 | 50.88 | 51.8933 | 52.5033 | 8.1507 | 8.1913 |
| SrO | 0.0088 | 0.0167 | 0.0183 | 0.3631 | 0.2826 |
| TiO2 | 1.0097 | 1.389 | 1.3247 | 0.0105 | 0.0086 |
| ZrO2 | 0.0115 | 0.0155 | 0.0139 | 0.0015 | 0.0017 |
Major oxide composition of the rocks sampled.
Note that totals for ABEs and TuiMs are not 100%; the remaining weight percent is composed primarily of barite for these two sulfides.
Table 3
| Element | KiMba | ABEba | Marba | ABEs | TuiMs |
|---|---|---|---|---|---|
| 7Li | 5.13 | 7.51 | 6.14 | 0.74 | 0.33 |
| 31P | 271.90 | 504.30 | 535.05 | 70.32 | 18.96 |
| 45Sc | 43.61 | 36.56 | 35.80 | 0.03 | 0.04 |
| 51V | 340.50 | 417.55 | 428.85 | 6.98 | 3.61 |
| 52Cr | 175.50 | 2.45 | 28.76 | 0.89 | 0 |
| 55Mn | 1500.50 | 2821.50 | 1685.50 | 1277.50 | 110.50 |
| 59Co | 42.91 | 36.95 | 33.30 | 0.05 | 0.01 |
| 60Ni | 69.31 | 23.89 | 23.84 | 3.04 | 0.48 |
| 65Cu | 77.42 | 77.74 | 80.19 | 2995.00 | 1950.00 |
| 66Zn | 86.48 | 133.50 | 84.49 | 102150.0 | 68325.0 |
| 71Ga | 15.66 | 18.13 | 17.66 | 2.82 | 11.93 |
| 72Ge | 2.49 | 2.73 | 2.63 | 33.32 | 16.25 |
| 85Rb | 0.74 | 3.96 | 5.19 | 1.70 | 1.08 |
| 86Sr | 75.41 | 143.00 | 155.65 | 3014.00 | 2267.00 |
| 89Y | 25.83 | 32.86 | 29.99 | 0 | 0 |
| 91Zr | 55.06 | 83.52 | 74.71 | 1.92 | 0.34 |
| 93Nb | 0.55 | 1.01 | 0.78 | 0 | 0 |
| 95Mo | 0 | 1.39 | 0.46 | 50.80 | 15.08 |
| 111Cd | 0.07 | 0.17 | 0.07 | 96.67 | 196.80 |
| 118Sn | 0.52 | 0.81 | 0.65 | 0 | 0 |
| 121Sb | 0 | 0 | 0 | 42.26 | 68.94 |
| 133Cs | 0 | 0 | 0 | 0.11 | 0 |
| 137Ba | 7.73 | 143.05 | 195.60 | 662.55 | 932.90 |
| 139La | 1.62 | 3.09 | 3.61 | 0.67 | 0.47 |
| 140Ce | 5.69 | 9.31 | 10.38 | 1.02 | 0.56 |
| 141Pr | 1.04 | 1.67 | 1.73 | 0.09 | 0.04 |
| 146Nd | 6.33 | 9.64 | 9.71 | 0.24 | 0.10 |
| 147Sm | 2.45 | 3.37 | 3.35 | 0.02 | 0.00 |
| 151Eu | 0.92 | 1.24 | 1.21 | 0.07 | 0.06 |
| 157Gd | 3.63 | 4.81 | 4.65 | 0.01 | 0.01 |
| 159Tb | 0.68 | 0.86 | 0.82 | 0.00 | 0.00 |
| 163Dy | 4.93 | 6.14 | 5.82 | 0 | 0 |
| 165Ho | 1.07 | 1.31 | 1.23 | 0.00 | 0.00 |
| 166Er | 3.30 | 4.04 | 3.75 | 0.00 | 0 |
| 169Tm | 0.48 | 0.59 | 0.55 | 0.00 | 0.00 |
| 172Yb | 3.23 | 3.94 | 3.65 | 0.00 | 0 |
| 175Lu | 0.50 | 0.61 | 0.56 | 0.00 | 0.00 |
| 178Hf | 1.76 | 2.49 | 2.40 | 0.02 | 0 |
| 181Ta | 0 | 0 | 0 | 0 | 0.01 |
| 182W | 0.07 | 0.12 | 0.07 | 0.05 | 0.00 |
| 208Pb | 0 | 3.07 | 1.41 | 3749.00 | 4972.50 |
| 232Th | 0.06 | 0.16 | 0.30 | 0.00 | 0 |
| 238U | 0 | 0 | 0 | 1.05 | 0.00 |
Trace element compositions of the rocks sampled.
All concentrations are in ppm.
Inactive sulfides ABEs and TuiMs were much more elevated in Ba, Cu, Zn, and Sr than the silicate samples (Tables 2 and 3). Sulfide TuiMs had nearly 10-fold more weight percent Ca than ABEs. A significant portion of ABEs and TuiMs was not acid digestible after 2 weeks in acid. This portion was collected and analyzed by XRD, which revealed that is composed of barite for both ABEs and TuiMs. With the exception of Sc and Ba, sample ABEs had higher trace element concentrations than sample TuiMs.
Microbiology
Bacterial and archaeal biomass were estimated via qPCR on all samples, including one sample each collected from the outside and inside of ABEs (ABEsO and ABEsIN, respectively). Bacteria accounted for >90% of the combined Bacteria + Archaea on all samples except for Marba and bacterial biomass was higher on the inactive sulfide samples (7 × 106–6 × 107 cells g−1) than on the silicates (∼1 × 104–6 × 105 cells g−1; Table 1). It must be noted, however, that because the primers used in qPCR are likely biased, this method provides an estimate, and the derived proportions are also estimates.
Given the high proportion of Bacteria versus Archaea, and their presence on all samples, we chose to analyze bacterial diversity. We PCR amplified and cloned near full-length bacterial 16S rRNA from all six samples; 469 total clones were recovered and sequenced, with 48–93 clones per sample (Figure 4). Rarefaction analysis using a 97% cutoff for OTU designation reveals that both basalt samples are more diverse than the other samples (Figure 4). Sample Marba exhibits the lowest diversity.
Figure 4
Further analyses of bacterial OTUs recovered at the phylum level reveals that some phyla exhibit north to south patterns in distribution on the silicates (Figure 4; Figure A2 in Appendix). α-Proteobacteria decrease in proportion from north to south, as do ε-proteobacteria and Bacteroidetes. In contrast, Planctomycetes and Chloroflexi increase in proportion from north to south. Because there are only two inactive sulfides, we cannot delineate any conclusions about north to south distributions of bacterial communities on hydrothermally inactive sulfides.
Bacterial populations on the inactive sulfide samples are markedly different from those on the silicates at the phylum level, and each of the three inactive sulfide samples are different from each other. No Chloroflexi were recovered from the sulfides, but given the smaller datasets generated by clone libraries compared to pyrosequencing (e.g., Flores et al.,
Recovered clones that group within the γ-proteobacteria include members of the orders Methylococcales, Thiotrichales, Chromatiales, and Pseudomonadales (Figure 5). Clones classified as Methylococcales are recovered only from sample KiMba. Those within the Thiotrichales are from the two basalt samples, while clones that fall within the Chromatiales order were recovered from all three silicate samples. Forty-seven clones recovered from ABEsO fall within the SUP05 clade of γ-proteobacteria, but these are exclusive to this sample. Nearly all of the clones that fall within the γ-proteobacteria are most similar to clones or isolates from other mid-ocean ridge or sedimentary settings.
Figure 5

Phylogenetic tree of representative clones from this study that fall within the γ-proteobacteria. Clones with more than one representative, or for which multiple members fell into a single clade, were chosen. Numbers in parentheses after some clones indicate total number of clones represented by the single clone depicted. The tree was generated in MEGA 5 using the maximum likelihood method with the Jukes–Cantor model and a Gamma distribution and 1000 bootstrap replicates. The same alignment was used to generate a neighbor-joining tree with the maximum composite likelihood method and 1000 bootstrap replicates. Nodes where both methods agree and bootstrap support was >50% are indicated with bootstrap values from the neighbor-joining tree on the left and the maximum likelihood tree on the right. Samples are color-coded and no representative clones are shared across multiple samples. D. psychrophila, S. deleyianum, and T. denitrificans were used as outgroups.
Some of the clones from samples KiMba and TuiMs are closely related to α-proteobacteria within the Roseobacter clade and the genus Hyphomicrobium (Figure 6). Clones that fall within the δ-proteobacteria are related to Nitrospina and Desulfocapsa and uncultured clones from other hydrothermal vent and sedimentary environments. Thirty-three clones recovered from TuiMs and three from ABEsIN are most similar to isolates in the genus Sulfurimonas [Thiomicrospira denitrificans was recently reclassified into this genus (Takai et al., 2006)].
Figure 6

Phylogenetic tree of representative clones from this study that fall within the α-, δ-, and ε-proteobacteria. Clones with more than one representative, or for which multiple members fell into a single clade, were chosen. Numbers in parentheses after some clones indicate total number of clones represented by the single clone depicted. The tree was generated in MEGA 5 using the maximum likelihood method with the Jukes–Cantor model and a Gamma distribution and 1000 bootstrap replicates. The same alignment was used to generate a neighbor-joining tree with the maximum composite likelihood method and 1000 bootstrap replicates. Nodes where both methods agree and bootstrap support was >50% are indicated with bootstrap values from the neighbor-joining tree on the left and the maximum likelihood tree on the right. Samples are color-coded and no representative clones are shared across multiple samples.
A diverse group of clones were recovered from our samples that group with phyla outside the Proteobacteria (Figure 7). Among these are a clade of Chloroflexi recovered from the two silicates ABEba and Marba. None of these clones are closely related to any cultured isolates, and the most similar environmental clone in the NCBI database is only 95% similar to clone ABEba_C8. Eighteen clones recovered from only inactive sulfides TuiMs and ABEsIN, represented by clones ABEsIN_H1 and TuiMs_LRBF3, group within a cluster of Bacteroidetes clones recovered from other inactive sulfides. This cluster falls within a larger clade, to which clone TuiMs_LR3E8 belongs, comprised of Bacteroidetes clones recovered exclusively from sulfidic environments.
Figure 7

Phylogenetic tree of representative clones from this study that fall within the non-proteobacterial phyla. Clones with more than one representative, or for which multiple members fell into a single clade, were chosen. Numbers in parentheses after some clones indicate total number of clones represented by the single clone depicted. The tree was generated in MEGA 5 using the maximum likelihood method with the Jukes-Cantor model and a Gamma distribution and 1000 bootstrap replicates. The same alignment was used to generate a neighbor-joining tree with the maximum composite likelihood method and 1000 bootstrap replicates. Nodes where both methods agree and bootstrap support was >50% are indicated with bootstrap values from the neighbor-joining tree on the left and the maximum likelihood tree on the right. Samples are color-coded and no representative clones are shared across multiple samples.
Nitrospirae were recovered from the three silicates, but not from the sulfides. These clones group most closely with Nitrospira marina (Figure 7). Four clones represented by clone ABEba_C10 are 85% similar to Ferrimicrobium acidiphilum str. T23 and Acidimicrobium ferrooxidans, both iron oxidizing acidophiles.
Correlations between geochemistry and microbiology
We used the elemental composition and phylogenetic data to determine correlations between geochemistry and microbiology. For this purpose, the same elemental data was used for ABEsO and ABEsIN because a bulk sample from this sulfide was analyzed for geochemistry. Ideally, one would compare the presence of individual OTUs in a sample with geochemical data to identify correlations between the two. However, using a 97% similarity cutoff for OTUs, only one OTU, comprised by clones TuiMs_LRB3A10 and ABEsIN_H6 (Figure 5), has membership from two different samples, and therefore no correlations exist at this cutoff between OTUs and geochemistry. Using a 10% cutoff only resulted in six OTUs with membership from multiple samples. Therefore, we measured correlations between the geochemical data and the abundance of bacterial phyla, as represented by their percentage on each sample. We do not report correlations between elemental variables alone because larger datasets are available for such purposes (e.g., Escrig et al.,
Nitrospirae and Chloroflexi are the only phyla with significant correlations to geochemistry. Both are strongly positively correlated to P, V, La, Ce, Pr, Nd, Th, P2O5, and SiO2 (Table 4). Nitrospirae are additionally positively correlated to the abundance of Y, Nb, Sn, Tb, Dy, Ho, Er, Tm, Yb, Lu, Li, La, Mn, Ga, Zr, Sm, Eu, Gd, Hf, W, MnO, Na2O, TiO2, and ZrO2, and negatively correlated to Sb. A negative correlation exists between abundance of γ-proteobacteria and δ-proteobacteria, Planctomycetes and Bacteroidetes, Planctomycetes and ε-proteobacteria, and Verrucomicrobia and δ-proteobacteria. There was a positive correlation between the proportion of Bacteroidetes and ε-proteobacteria, Nitrospirae and Chloroflexi, and Verrucomicrobia and γ-proteobacteria.
Table 4
| Variable 1 | By variable 2 | Kendall τ | Prob > | τ| |
|---|---|---|---|
| γ-proteobacteria | δ-Proteobacteria | −0.9429 | 0.0048 |
| Nitrospirae | Y, Nb, Sn, Tb, Dy, Ho, Er, Tm, Yb & Lu | 1 | <0.0001 |
| Nitrospirae | Li, Mn, Ga, Zr, Sm, Eu, Gd, Hf, W, MnO, Na2O, TiO2 & ZrO2 | 0.9549 | 0.0030 |
| Nitrospirae | P, V, La, Ce, Pr, Nd, Th, P2O5 & SiO2 | 0.8933 | 0.0165 |
| Nitrospirae | Sb | −0.8854 | 0.0190 |
| Nitrospirae | Chloroflexi | 0.8262 | 0.0427 |
| Bacteroidetes | ε-Proteobacteria | 0.9412 | 0.0051 |
| Planctomycetes | Bacteroidetes | −0.9276 | 0.0077 |
| Planctomycetes | ε-Proteobacteria | −0.8117 | 0.0499 |
| Chloroflexi | P, V, Rb, La, Ce, Pr, Nd, Th, K2O, P2O5, SiO2 | 0.8575 | 0.0291 |
| Verrucomicrobia | δ-Proteobacteria | −0.8452 | 0.0341 |
| Verrucomicrobia | γ-proteobacteria | 0.8452 | 0.0341 |
Kendall’s τ correlations between elemental data and bacterial distributions.
Correlations between pairs of variables that are both elemental data are not included here. Only correlations for which p < 0.05 are presented here. When multiple variables are listed in the “by Variable 2” column, this means that both Kendall τ and Prob > | τ| values are the same between Variable 1 and all variables listed in “by Variable 2” column.
Inactive sulfide biogeography
We compared the bacterial communities on inactive sulfides recovered during this study with those from previous studies (Suzuki et al., 2004; Kato et al.,
Figure 8

Cladogram of OTU-based comparison of bacterial communities sampled from inactive sulfides in this and prior studies using the Jclass algorithm. The three silicate samples from this study are included. Where the Jclass analysis agreed with community comparison using Jest, ΘYC, and/or Bray–Curtis algorithms, a dot for the agreeing methods is placed at the node in the cladogram. Overall agreement is strong between all four methods save for sample 3M33 from the EPR. Sampling location is indicated for samples from prior studies: EPR (Sylvan et al., 2012a), Mariana Trough (Kato et al.,
Discussion
Novelty of work
Here we investigated microbiology along a geochemical gradient in the ELSC. This back-arc spreading system is of particular interest because it displays clear trends in host rock geochemistry (Escrig et al.,
Silicates
Our geochemical and mineralogical analysis of the silicates collected reinforces north-south along axis patterns seen for Si, Sr, Ba, La, Ce, Nd, and Th in previous work that looked at fresh lava flows in the ELSC (Escrig et al.,
Bacterial communities on basalts from Lau Basin were more diverse than those on inactive sulfides from this study (Figure 4), in agreement with prior work indicating basalts host extremely diverse bacterial populations (Santelli et al., 2008). Within the sulfide samples, TuiMs is more diverse than both ABEsO and ABEsIN. Inactive sulfide ABEs is ∼150% more enriched in Zn than is TuiMs. A recent comparison of inactive sulfides from the EPR found that Zn-rich chimneys harbor fewer OTUs than Fe-rich sulfides (Toner et al., 2013). This is also true here.
North to south patterns in bacterial community composition occur at the phylum level in the ELSC. Kilo Moana (north) is analogous to a mid-ocean ridge vent field, whereas Mariner (south) is strongly influenced by subduction, and ABE exhibits transitional characteristics between the two end-members. As such, patterns of bacterial membership on the silicate rocks follow this pattern linearly for the α-proteobacteria, ε-proteobacteria, Bacteroidetes, Planctomycetes, and Chloroflexi, where the proportional representation of each phylum on sample ABEba falls between that on samples KiMba and Marba. This is supported in part by the high positive correlation between Bacteroidetes and ε-proteobacteria (both decreased from north to south) and the negative correlation between Planctomycetes and Bacteroidetes (opposite north to south patterns). These gradients are likely drivers for absence of shared OTUs between bacterial communities detected on each sample. This is reinforced by Libshuff analysis, which reveals that each silicate harbors a unique bacterial community. Further, when the inactive sulfides are included in this analysis, all six samples in this study are significantly different from each other.
The bacterial community on sample Marba is quite different from what has been observed on previously studied seafloor-exposed basalts (Lysnes et al.,
The high proportion on Marba of clones that fall within the phylum Chloroflexi is also unusual compared to previously analyzed basalts; while this phylum is often represented as a few percent of the total clones on seafloor-exposed basalt samples from mid-ocean ridge settings, the large community membership of a monophyletic clade of Chloroflexi, as seen here, has not previously been observed. Clones from the same clade were also recovered from sample ABEba and all members of this clade are <90% similar to the most closely related sequences in the NCBI database. As such, it is impossible to speculate on the ecological niche for the organisms represented by these clones, or whether their presence is or is not affected by substrate geochemistry. The Chloroflexi clones recovered from samples Marba and ABEba fall within the class Anaerolineae, which was also recovered from active hydrothermal vent sulfides at Mariner (Takai et al., 2008; Flores et al.,
Our finding of a microbial community on a silicate from Mariner vent field that is distinct from those collected at ABE and Kilo Moana mirrors a detailed study of bacterial and archaeal diversity on active sulfide chimneys in the ELSC, which also found the microbial communities at Mariner to be distinct (Flores et al.,
Clones classified within the phylum Nitrospirae and the genus Nitrospira and the γ-proteobacterial order Chromatiales were recovered from all three silicate samples. The clones related to Nitrospira fall within Nitrospira sublineage IV, which are known to oxidize nitrite to nitrate mixotrophically (Daims et al.,
Many of the bacterial clones recovered from silicate samples KiMba and ABEba are closely related to isolates involved in sulfur, methane, and hydrogen biogeochemical cycling. Within the γ-proteobacteria, this includes clones that fall within the order Methylococcales, aerobic methanotrophs, as well as OTUs allied to the orders Thiotrichales and Chromatiales (Figure 8). Known sulfur oxidizers populate both of these orders. Clones recovered from the basalts classified as ε-proteobacteria are most closely related to the genera Sulfurospirillum, Nitratifractor, and endosymbionts of vent macrofauna (Figure A3 in Appendix). Species within the genus Sulfurospirillum heterotrophically oxidize sulfur while those within Nitratifractor chemoautotrophically couple nitrate reduction with H2 oxidation (Nakagawa et al.,
Previous studies have also noted the co-occurrence of sulfur oxidizing and sulfur reducing bacteria on seafloor-exposed basalts (Santelli et al., 2009; Sudek et al., 2009). This indicates that pores within the rocks provide different niches for diverse microbial lifestyles. Fresh basalt is rich in reduced sulfur and therefore sulfur oxidizing lineages are expected (Bach and Edwards,
Inactive sulfides
We observed differences between the bacterial communities detected on silicate samples and those detected on the inactive sulfides. This is in agreement with recent work that shows geochemistry strongly influences bacterial community membership on a given substrate, even when temperature and location are similar (Toner et al., 2013), as is the case here for samples ABEba, ABEsO, and ABEsIN, which were all collected from ABE vent field. Further, we found significant differences between bacterial communities detected on the outer chimney wall and inner conduit of the same inactive sulfide chimney, ABEs. Similar zonation was noted in prior studies of both active and inactive sulfide chimneys (Schrenk et al., 2003; Suzuki et al., 2004; Kormas et al.,
Sample ABEsO is host to a population of γ-proteobacteria that belongs to the SUP05 clade (Figure 5). The SUP05 clade was originally detected in the hydrothermal plume of the Suiyo Seamount (Sunamura et al., 2004) and has since been detected in hydrothermal plumes in vent fields from Guaymas Basin (Dick and Tebo,
Clones TuiMs_LRB3C7, TuiMs_LRB3F6, TuiMs_C8, and TuiMs_B9 (Figure 6) all belong to the α-proteobacterial family Rhodobacteraceae and are most similar to cultured representatives belonging to the genera Sulfitobacter and Roseovarius. All the species from these genera whose genome has been sequenced contain the sox cluster of genes that imparts the ability to oxidize sulfur (Newton et al., 2010). Therefore, it is likely that inactive sulfide TuiMs harbors at least three different niches for sulfur oxidation to accommodate for S-oxidizing SUP05 bacteria, ε-proteobacteria, and Rhodobacteraceae.
The clones representative of δ-proteobacterial lineages recovered from ABEsIN are only distantly related to their closest cultured relatives; clone ABEsIN_D9 is 83% similar to Nitrospina griacilis, a nitrite oxidizer, and clone ABEsIN_A1 is 89% similar to both Mono Lake isolate strain MLMS-1 and Desulfovibrio alkaliphilus str. AHT2. This former isolate autotrophically couples sulfur oxidation with As-reduction (Hoeft et al.,
Clones TuiMs_LRB3F3 and ABEsIN_H1 (Figure 7) group with closely related clones recovered from inactive sulfides collected from the EPR (Sylvan et al., 2012b) and Indian Ocean (Suzuki et al., 2004). In the EPR study, these clones represented 17% of all the clones recovered, and were detected on six of seven samples. Similar clones detected via BlastN (Altschul et al.,
Figure 9

Phylogenetic tree of representative clones from this study that fall within the Sulfiphilic Bacteroidetes clade (red branches). One representative clone was chosen from each study and from each sample if multiple samples are represented from a study. The tree was generated in MEGA 5 (Tamura et al., 2011) using the maximum likelihood method with the Jukes–Cantor model, Gamma distribution and 1000 bootstrap replicates. The same alignment was used to generate a neighbor-joining tree with the maximum composite likelihood method and 1000 bootstrap replicates. Nodes where both methods agree and bootstrap support is >50% are indicated with bootstrap values from the neighbor-joining tree on the left and the maximum likelihood tree on the right. Stars indicate nodes supported by Bayesian analysis. Nitrospira marina was used as an outgroup.
Following cessation of venting on a sulfide chimney, a different bacterial assemblage is known to succeed the one present during active venting (Sylvan et al., 2012a). This assemblage is unique from that found on active chimneys (Kato et al.,
The mineralogy of inactive sulfides appears to influence the composition of the extant bacterial communities (Kato et al.,
Implications for weathering of seafloor rocks
The noted patterns in microbial communities, both from north to south and between different substrates, have implications for weathering of seafloor rocks. It is known that bacteria incubated with basalt enhance Si, Fe, and Mn release into the aqueous phase (Daughney et al.,
Conclusion
It is likely that the north-south gradients in bacterial community composition are driven by the differences in substrate chemistry between these fields. Many elements display linear transitions from Kilo Moana to Mariner, which makes it difficult to pinpoint which one (or ones) may be driving the observed differences in microbiology, but we were able to detect correlation between the elemental composition of the collected samples and bacterial phyla Nitrospirae and Chloroflexi (Table 4). Additionally, with the small sample size of this study, it is important to point out that the patterns observed here warrant further investigation and would benefit from additional and more intensive sampling. The samples collected for this study represent an initial foray into the microbiology of low temperature deposits in Lau Basin.
Bacterial communities detected on low temperature silicates and inactive sulfides along the ELSC and VFR in Lau Basin display distinct patterns that are driven by (1) gradients in rock geochemistry from north to south on the silicates, and (2) differences in substrate between the basalts, basaltic andesite, and inactive sulfides. The prevalence of Chloroflexi clones distantly related to any known isolates on basaltic andesite indicates that this substrate may host unique microbial populations. Clones related to organisms involved in multiple facets of sulfur and iron oxidation-reduction processes on the same rocks indicate that multiple micro-niches are present on both silicates and sulfides. Further, multiple niches for the same ecological function, like sulfur oxidation by three different clades of organisms on sulfide ABEs, indicate that fine scale differences in mineralogy likely support similarly nuanced microbial communities. Bacterial communities on inactive sulfides from Lau Basin are not unlike others studied at different sites, but it appears that each structure studied harbors a bacterial community most similar to itself and, therefore, further study of these ecosystems with finer resolution is warranted.
Statements
Acknowledgments
We thank the captain and crew of the RV Thomas G. Thompson and the pilots of ROV Jason II for their assistance and expertise in acquiring the samples studied here. Frank Corsetti (USC) kindly allowed the use of his microscopes for the thin section analysis. We thank Linda Sauer and Mike Manno (Characterization Facility, University of Minnesota) for help with XRD and Rick Knurr (Department of Earth Sciences, University of Minnesota) for his expertise running the samples for chemical composition. Parts of this work were carried out in the Characterization Facility, University of Minnesota, which receives partial support from NSF through the MRSEC program. This work was supported by National Science Foundation (NSF) Grant OCE-0732369 to PRG, the NSF funded Center for Dark Energy Biosphere Investigations (Katrina J. Edwards), and a C-DEBI postdoctoral fellowship to JBS. This is C-DEBI contribution number 150.
Conflict of interest
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.
References
1
AltschulS. F.GishW.MillerW.MyersE. W.LipmanD. J. (1990). Basic local alignment search tool. J. Mol. Biol.215, 403–410.10.1016/S0022-2836(05)80360-2
2
AusubelF. M.BrentR.KinsgstonR. E.MooreD. D.SeidmanJ. G.SmithJ. A.et al (1999). Short Protocols in Molecular Biology, 4th Edn.New York: Wiley.
3
BachW.EdwardsK. J. (2003). Iron and sulfide oxidation within the basaltic ocean crust: implications for chemolithoautotrophic microbial biomass production. Geochim. Cosmochim. Acta67, 3871–3887.10.1016/S0016-7037(03)00304-1
4
BaileyB.TempletonA.StaudigelH.TeboB. M. (2009). Utilization of substrate components during basaltic glass colonization by Pseudomonas and Shewanella isolates. Geomicrobiol. J.26, 648–656.10.1080/01490450903263376
5
BiddleJ. F.HouseC. H.BrenchleyJ. E. (2005). Microbial stratification in deeply buried marine sediment reflects changes in sulfate/methane profiles. Geobiology3, 287–295.10.1111/j.1472-4669.2006.00062.x
6
BrennerD. J.KriegN. R.StaleyJ. T.GarrityG. M.BooneD. R.De VosP.et al (2005). The proteobacteria, Part B, the Gammaproteobacteria, in Bergey’s Manual of Systematic Bacteriology, 2nd Edn, Ed. GarrityG. M. (East Lansing, MI: Springer).
7
CampbellB. J.EngelA. S.PorterM. L.TakaiK. (2006). The versatile epsilon-proteobacteria: key players in sulphidic habitats. Nat. Rev. Microbiol.4, 458–468.10.1038/nrmicro1414
8
DaimsH.NielsenJ. L.NielsenP. H.SchleiferK. H.WagnerM. (2001). In situ characterization of Nitrospira-like nitrite oxidizing bacteria active in wastewater treatment plants. Appl. Environ. Microbiol.67, 5273–5284.10.1128/AEM.67.11.5273-5284.2001
9
DaughneyC. J.RiouxJ. P.FortinD.PichlerT. (2004). Laboratory investigation of the role of bacteria in the weathering of basalt near deep sea hydrothermal vents. Geomicrobiol. J.21, 21–31.10.1080/01490450490253437
10
Dell’AnnoA.DanovaroR. (2005). Extracellular DNA plays a key role in deep-sea ecosystem functioning. Science309, 2179–2179.10.1126/science.1117475
11
DelongE. F. (1992). Archaea in coastal marine environments. Proc. Natl. Acad. Sci. U.S.A.89, 5685–5689.10.1073/pnas.89.12.5685
12
DeSantisT. Z.HugenholtzP.KellerK.BrodieE. L.LarsenN.PicenoY. M.et al (2006a). NAST: a multiple sequence alignment server for comparative analysis of 16S rRNA genes. Nucleic Acids Res.34, W394–W399.10.1093/nar/gkj156
13
DeSantisT. Z.HugenholtzP.LarsenN.RojasM.BrodieE. L.KellerK.et al (2006b). Greengenes, a chimera-checked 16S rRNA gene database and workbench compatible with ARB. Appl. Environ. Microbiol.72, 5069–5072.10.1128/AEM.03006-05
14
DickG. J.TeboB. M. (2010). Microbial diversity and biogeochemistry of the Guaymas Basin deep-sea hydrothermal plume. Environ. Microbiol.12, 1334–1347.10.1111/j.1462-2920.2010.02177.x
15
DrummondA. J.AshtonB.BuxtonS.CheungM.CooperA.DuranC.et al (2011). Geneious v5.4. Available at: http://wwwgeneiouscom
16
DunnR. A.MartinezF. (2011). Contrasting crustal production and rapid mantle transitions beneath back-arc ridges. Nature469, 198–202.10.1038/nature09690
17
EdwardsK. J.BachW.McCollomT. M.RogersD. R. (2004). Neutrophilic iron-oxidizing bacteria in the ocean: their habitats, diversity, and roles in mineral deposition, rock alteration, and biomass production in the deep-sea. Geomicrobiol. J.21, 393–404.10.1080/01490450490485863
18
EdwardsK. J.McCollomT. M.KonishiH.BuseckP. R. (2003). Seafloor bioalteration of sulfide minerals: results from in situ incubation studies. Geochim. Cosmochim. Acta67, 2843–2856.10.1016/S0016-7037(03)00089-9
19
EinenJ.ThorsethI. H.OvreasL. (2008). Enumeration of Archaea and Bacteria in seafloor basalt using real-time quantitative PCR and fluorescence microscopy. FEMS Microbiol. Lett.282, 182–187.10.1111/j.1574-6968.2008.01119.x
20
EscrigS.BezosA.GoldsteinS. L.LangmuirC. H.MichaelP. J. (2009). Mantle source variations beneath the Eastern Lau Spreading Center and the nature of subduction components in the Lau Basin-Tonga arc system. Geochem. Geophys. Geosyst.10:Q04014.10.1029/2008GC002281
21
FloresG. E.ShakyaM.MeneghinJ.YangZ. K.SeewaldJ. S.WheatC. G.et al (2012). Inter-field variability in the microbial communities of hydrothermal vent deposits from a back-arc basin. Geobiology10, 333–346.10.1111/j.1472-4669.2012.00325.x
22
FouquetY.VonstackelbergU.CharlouJ. L.DonvalJ. P.ErzingerJ.FoucherJ. P.et al (1991). Hydrothermal activity and metallogenesis in the Lau back-arc basin. Nature349, 778–781.10.1038/349778a0
23
GermanC. R.BowenA.ColemanM. L.HonigD. L.HuberJ. A.JakubaM. V.et al (2010). Diverse styles of submarine venting on the ultraslow spreading Mid-Cayman Rise. Proc. Natl. Acad. Sci. U.S.A.107, 14020–14025.10.1073/pnas.1009205107
24
HaymonR. M.KastnerM. (1981). Hot spring deposits on the East Pacific Rise at 21°N: preliminary description of mineralogy and genesis. Earth Planet. Sci. Lett.53, 363–381.10.1016/0012-821X(81)90041-8
25
HoeftS. E.KulpT. R.StolzJ. F.HollibaughJ. T.OremlandR. S. (2004). Dissimilatory arsenate reduction with sulfide as electron donor: experiments with mono lake water and isolation of strain MLMS-1, a chemoautotrophic arsenate respirer. Appl. Environ. Microbiol.70, 2741–2747.10.1128/AEM.70.5.2741-2747.2004
26
HuseS. M.WelchD. M.MorrisonH. G.SoginM. L. (2010). Ironing out the wrinkles in the rare biosphere through improved OTU clustering. Environ. Microbiol.12, 1889–1898.10.1111/j.1462-2920.2010.02193.x
27
IshibashiJ.LuptonJ. E.YamaguchiT.QuerellouJ.NunouraT.TakaiK. (2006). Expedition reveals changes in Lau Basin hydrothermal system. Eos (Washington DC)87, 13.
28
IvarssonM.BengstonS.BelivanovaV.StampanoniM.MaroneF.TehlerA. (2012). Fossilized fungi in subseafloor Eocene basalts. Geology40, 163–166.10.1130/G32590.1
29
KatoS.TakanoY.KakegawaT.ObaH.InoueK.KobayashiC.et al (2010). Biogeography and biodiversity in sulfide structures of active and inactive vents at deep-sea hydrothermal fields of the Southern Mariana Trough. Appl. Environ. Microbiol.76, 2968–2979.10.1128/AEM.00478-10
30
KormasK. A.TiveyM. K.Von DammK.TeskeA. (2006). Bacterial and archaeal phylotypes associated with distinct mineralogical layers of a white smoker spire from a deep-sea hydrothermal vent site (9 degrees N, East Pacific Rise). Environ. Microbiol.8, 909–920.10.1111/j.1462-2920.2005.00978.x
31
KristallB.NielsenD.HanningtonM. D.KelleyD. S.DelaneyJ. R. (2011). Chemical microenvironments within sulfide structures from the Mothra hydrothermal field: evidence from high-resolution zoning of trace elements. Chem. Geol.290, 12–30.10.1016/j.chemgeo.2011.08.008
32
LangmuirC. H.GermanC. R.MichaelP. J.YoergerD. R.FornariD. J.ShankG. C.et al (2004). Hydrothermal prospecting and petrological sampling in the Lau Basin: background data for the integrated study site. Eos (Washington DC)5, abstr.B13A–0189.
33
Le BasM. J.StreckeisenA. L. (1991). The IUGS systematics of igneous rocks. J. Geol. Soc. London148, 825–833.10.1144/gsjgs.148.5.0825
34
Lopez-GarciaP.VereshchakaA.MoreiraD. (2007). Eukaryotic diversity associated with carbonates and fluid-seawater interface in Lost City hydrothermal field. Environ. Microbiol.9, 546–554.10.1111/j.1462-2920.2006.01158.x
35
LudwigW.StrunkO.WestramR.RichterL.MeierH.Yadhukumaret al (2004). ARB: a software environment for sequence data. Nucleic Acids Res.32, 1363–1371.10.1093/nar/gkh293
36
LysnesK.ThorsethI. H.SteinsbuB. O.OvreasL.TorsvikT.PedersenR. B. (2004). Microbial community diversity in seafloor basalt from the Arctic spreading ridges. FEMS Microbiol. Ecol.50, 213–230.10.1016/j.femsec.2004.06.014
37
MasonO. U.Di Meo-SavoieC. A.Van NostrandJ. D.ZhouJ. Z.FiskM. R.GiovannoniS. J. (2009). Prokaryotic diversity, distribution, and insights into their role in biogeochemical cycling in marine basalts. ISME J.3, 231–242.10.1038/ismej.2008.92
38
MottlM. J.SeewaldJ. S.WheatC. G.TiveyM. K.MichaelP. J.ProskurowskiG.et al (2011). Chemistry of hot springs along the Eastern Lau Spreading Center. Geochim. Cosmochim. Acta75, 1013–1038.10.1016/j.gca.2010.12.008
39
NakagawaS.TakaiK.InagakiF.HorikoshiK.SakoY. (2005). Nitratiruptor tergarcus gen. nov., sp. nov. and Nitratifractor salsuginis gen. nov., sp. nov., nitrate-reducing chemolithoautotrophs of the epsilon-Proteobacteria isolated from a deep-sea hydrothermal system in the Mid-Okinawa Trough. Int. J. Syst. Evol. Microbiol.55, 925–933.10.1099/ijs.0.63351-0
40
NewtonR. J.GriffinL. E.BowlesK. M.MeileC.GiffordS.GivensC. E.et al (2010). Genome characteristics of a generalist marine bacterial lineage. ISME J.4, 784–798.10.1038/ismej.2009.150
41
PodowskiE. L.MaS.LutherG. W.IIIWardropD.FisherC. R. (2010). Biotic and abiotic factors affecting distributions of megafauna in diffuse flow on andesite and basalt along the Eastern Lau Spreading Center, Tonga. Mar. Ecol. Prog. Ser.418, 25–45.10.3354/meps08797
42
RogersD. R.SantelliC. M.EdwardsK. J. (2003). Geomicrobiology of deep-sea deposits: estimating community diversity from low-temperature seafloor rocks and minerals. Geobiology1, 109–117.10.1046/j.1472-4669.2003.00009.x
43
SantelliC. M.EdgcombV. P.BachW.EdwardsK. J. (2009). The diversity and abundance of bacteria inhabiting seafloor lavas positively correlate with rock alteration. Environ. Microbiol.11, 86–98.10.1111/j.1462-2920.2008.01743.x
44
SantelliC. M.OrcuttB. N.BanningE.BachW.MoyerC. L.SoginM. L.et al (2008). Abundance and diversity of microbial life in ocean crust. Nature453, 653–657.10.1038/nature06899
45
SchlossP. D.WestcottS. L.RyabinT.HallJ. R.HartmannM.HollisterE. B.et al (2009). Introducing mothur: open-source, platform-independent, community-supported software for describing and comparing microbial communities. Appl. Environ. Microbiol.75, 7537–7541.10.1128/AEM.01541-09
46
SchrenkM. O.KelleyD. S.DelaneyJ. R.BarossJ. A. (2003). Incidence and diversity of microorganisms within the walls of an active deep-sea sulfide chimney. Appl. Environ. Microbiol.69, 3580–3592.10.1128/AEM.69.6.3580-3592.2003
47
SmithA.PopaR.FiskM.NielsenM.WheatC. G.JannaschH. W.et al (2011). In situ enrichment of ocean crust microbes on igneous minerals and glasses using an osmotic flow-through device. Geochem. Geophys. Geosyst.12:Q06007.10.1029/2010GC003424
48
SorokinD. Y.TourovaT. P.MussmannM.MuyzerG. (2008). Dethiobacter alkaliphilus gen. nov. sp. nov., and Desulfurivibrio alkaliphilus gen. nov. sp. nov.: two novel representatives of reductive sulfur cycle from soda lakes. Extremophiles12, 431–439.10.1007/s00792-008-0148-8
49
SudekL. A.TempletonA. S.TeboB. M.StaudigelH. (2009). Microbial ecology of Fe (hydr)oxide mats and basaltic rock from Vailulu’u Seamount, American Samoa. Geomicrobiol. J.26, 581–596.10.1080/01490450903263400
50
SunamuraM.HigashiY.MiyakoC.IshibashiJ.MaruyamaA. (2004). Two bacteria phylotypes are predominant in the Suiyo Seamount hydrothermal plume. Appl. Environ. Microbiol.70, 1190–1198.10.1128/AEM.70.2.1190-1198.2004
51
SuzukiY.InagakiF.TakaiK.NealsonK. H.HorikoshiK. (2004). Microbial diversity in inactive chimney structures from deep-sea hydrothermal systems. Microb. Ecol.47, 186–196.10.1007/s00248-003-1014-y
52
SylvanJ. B.TonerB. M.EdwardsK. J. (2012a). Life and death of deep-sea vents: bacterial diversity and ecosystem succession on inactive hydrothermal sulfides. MBio3, e00279–e00211.10.1128/mBio.00279-11
53
SylvanJ. B.PyensonB. C.RouxelO.GermanC. R.EdwardsK. J. (2012b). Time series analysis of two hydrothermal plumes at 9°50’N East Pacific Rise reveals distinct, heterogeneous bacterial populations. Geobiology10, 178–192.10.1111/j.1472-4669.2011.00315.x
54
TakaiK.HorikoshiK. (2000). Rapid detection and quantification of members of the archaeal community by quantitative PCR using fluorogenic probes. Appl. Environ. Microbiol.66, 5066–5072.10.1128/AEM.66.11.5066-5072.2000
55
TakaiK.NunouraT.IshibashiJ. I.LuptonJ.SuzukiR.HamasakiH.et al (2008). Variability in the microbial communities and hydrothermal fluid chemistry at the newly discovered Mariner hydrothermal field, southern Lau Basin. J. Geophys. Res. Biogeosci.113:G02031.10.1029/2007JG000636
56
TakaiK.SuzukiM.NakagawaS.MiyazakiM.SuzukiY.InagakiF.et al (2006). Sulfurimonasparalvinellae sp. nov., a novel mesophilic, hydrogen- and sulfur-oxidizing chemolithoautotroph within the Epsilonproteobacteria isolated from a deep-sea hydrothermal vent polychaete nest, reclassification of Thiomicrospira denitrificans as Sulfurimonas denitrificans comb. nov. and emended description of the genus Sulfurimonas. Int. J. Syst. Evol. Microbiol.56, 1725–1733.10.1099/ijs.0.64255-0
57
TamuraK.PetersonD.PetersonN.StecherG.NeiM.KumarS. (2011). MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance and maximum parsimony methods. Mol. Biol. Evol.28, 2731–2739.10.1093/molbev/msr121
58
TonerB. M.LesniewskiR. A.MarlowJ. J.BriscoeL. J.SantelliC. M.BachW.et al (2013). Mineralogy drives bacterial biogeography of hydrothermally inactive seafloor sulfide deposits. Geomicrobiol. J.30, 313–326.10.1080/01490451.2012.688925
59
TonerB. M.SantelliC. M.MarcusM. A.WirthR.ChanC. S.McCollomT.et al (2009). Biogenic iron oxyhydroxide formation at mid-ocean ridge hydrothermal vents: Juan de Fuca Ridge. Geochim. Cosmochim. Acta73, 388–403.10.1016/j.gca.2008.09.035
60
WalshD. A.ZaikovaE.HowesC. G.SongY. C.WrightJ. J.TringeS. G.et al (2009). Metagenome of a versatile chemolithoautotroph from expanding oceanic dead zones. Science326, 578–582.10.1126/science.1174010
Appendix
Figure A1

Seafloor and shipboard photographs of samples collected during this study. Sample KiMba as it appeared on the seafloor (A) and in the sampling box (B). Sample ABEba on the seafloor (C) and onboard (D). Sample Marba as it appeared on the seafloor (E) and on the ship (F). Sample ABEs as it was collected by ROV Jason II (G) and on the ship (H). Sample TuiMs, photographed during collection (I) and on the ship (J).
Figure A2

Bacterial distributions of phyla/classes that either decrease (left column) or increase (right column) in proportionality from north (Kilo Moana) to south (Mariner) on silicate samples in this study. Key is the same as Figure 4, except that phyla/classes that decrease from north to south are all red, those that increase from north to south are all green, and those that do not show a north to south pattern are not included.
Figure A3

Phylogenetic tree of all representative clones from this study that fall within the ε-Proteobacteria. The tree was generated in MEGA 5 (Tamura et al., 2011) using the maximum likelihood method with the Jukes–Cantor model and a Gamma distribution and 500 bootstrap replicates. Nodes where bootstrap support is >50% are indicated. Samples are color-coded and no representative clones are shared across multiple samples.
Summary
Keywords
geomicrobiology, basalt, inactive sulfides, hydrothermal, bacteroidetes
Citation
Sylvan JB, Sia TY, Haddad AG, Briscoe LJ, Toner BM, Girguis PR and Edwards KJ (2013) Low Temperature Geomicrobiology Follows Host Rock Composition Along a Geochemical Gradient in Lau Basin. Front. Microbiol. 4:61. doi: 10.3389/fmicb.2013.00061
Received
02 December 2012
Accepted
04 March 2013
Published
27 March 2013
Volume
4 - 2013
Edited by
Andreas Teske, University of North Carolina at Chapel Hill, USA
Reviewed by
William D. Orsi, Woods Hole Oceanographic Institution, USA; Federico Lauro, University of New South Wales, Australia
Copyright
© 2013 Sylvan, Sia, Haddad, Briscoe, Toner, Girguis and Edwards.
This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in other forums, provided the original authors and source are credited and subject to any copyright notices concerning any third-party graphics etc.
*Correspondence: Katrina J. Edwards, Department of Biological Sciences, University of Southern California, 3616 Trousdale Parkway, Los Angeles, CA 90089, USA. e-mail: kje@usc.edu
This article was submitted to Frontiers in Extreme Microbiology, a specialty of Frontiers in Microbiology.
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.