Influence of Igneous Basement on Deep Sediment Microbial Diversity on the Eastern Juan de Fuca Ridge Flank.

Microbial communities living in deeply buried sediment may be adapted to long-term energy limitation as they are removed from new detrital energy inputs for thousands to millions of years. However, sediment layers near the underlying oceanic crust may receive inputs from below that influence microbial community structure and/or activity. As part of the Census of Deep Life, we used 16S rRNA gene tag pyrosequencing on DNA extracted from a spectrum of deep sediment-basement interface samples from the subsurface of the Juan de Fuca Ridge flank (collected on IODP Expedition 327) to examine this possible basement influence on deep sediment communities. This area experiences rapid sedimentation, with an underlying basaltic crust that hosts a dynamic flux of hydrothermal fluids that diffuse into the sediment. Chloroflexi sequences dominated tag libraries in all sediment samples, with variation in the abundance of other bacterial groups (e.g., Actinobacteria, Aerophobetes, Atribacteria, Planctomycetes, and Nitrospirae). These variations occur in relation to the type of sediment (clays versus carbonate-rich) and the depth of sample origin, and show no clear connection to the distance from the discharge outcrop or to basement fluid microbial communities. Actinobacteria-related sequences dominated the basalt libraries, but these should be viewed cautiously due to possibilities for imprinting from contamination. Our results indicate that proximity to basement or areas of seawater recharge is not a primary driver of microbial community composition in basal sediment, even though fluids diffusing from basement into sediment may stimulate microbial activity.

Introduction

Over the past decade, several studies have expanded our understanding of the composition and function of “deep biosphere” microbial communities in deep marine sediments in both organic-rich (D’Hondt et al., 2004; Inagaki et al., 2006; Biddle et al., 2006, 2008; Lever et al., 2010) and organic-poor environments (D’Hondt et al., 2009; Kallmeyer et al., 2012; Røy et al., 2012). Studies of the form and function of microbial communities within or expelled from the deep oceanic crust are rarer and geographically more limited (Cowen et al., 2003; Thorseth et al., 2003; Huber et al., 2006; Santelli et al., 2008; Orcutt et al., 2011a; Lin et al., 2012; Jungbluth et al., 2013; Lever et al., 2013; Lee et al., 2015), even though this habitat type is spatially extensive and occupies a large volume of the crust (Bach and Edwards, 2003; Orcutt et al., 2011b). Recent studies have challenged traditional microbiological concepts of what it means for cells to be alive, active, or dormant (Hoehler and Jørgensen, 2013). Metabolically active cells (Orsi et al., 2013) that can respond to nutrient addition through growth (Morono et al., 2011) have been confirmed in organic-rich sediment, while investigations in organic-poor environments have confirmed the persistence of microbial activity and cells under extreme energy limitation all the way through the sediment column to the underlying basement (D’Hondt et al., 2009, 2015; Røy et al., 2012). Another study confirmed the persistence of microbial groups from initial deposition through burial down to 2.5 km below seafloor millions of years later, and documented strong decreases in microbial abundance at 40–60°C (Inagaki et al., 2015), supporting the notion that increased biomolecule damage combined with severe energy limitation may cause the upper temperature limit of life in subsurface sediments to be lower than in energy-rich surface environments (Lever et al., 2015).

While most of these studies have focused on what happens to microbial life as it is buried beneath increasing layers of sediment, several recent studies have also highlighted the impact that the underlying basement can have on stimulating microbial activity from below. Ocean bottom water is entrained into ocean crust as a result of hydrothermal circulation, undergoing fluid–rock reactions during transit through basement that can change the chemistry of the fluid (Elderfield et al., 1999; Wheat and Mottl, 2004). These altered fluids can then diffuse upward into the overlying sediment, potentially providing a source of electron acceptors, nutrients, and organic carbon compounds to deeply buried sediment microbial communities. For example, in organic-poor sediment environments, oxygen and nitrate diffuse upward from minimally altered oxic basement fluids into sediment to stimulate microbial activity (Ziebis et al., 2012; Orcutt et al., 2013; D’Hondt et al., 2015; Wankel et al., 2015), although it is currently unclear what the influence is on deep sediment microbial community structure at these sites. Likewise, sulfate diffusing upward from highly altered anoxic basement fluids into overlying sediment may stimulate anaerobic microbial activity (Expedition 301 Scientists, 2005; Engelen et al., 2008; Lever et al., 2010). However, the influence of these basement sources on the structure of the microbial community is not currently resolved.

The Juan de Fuca Ridge flank is a well-studied ridge flank environment where bottom seawater recharges into young basement (3.5 Ma) at crustal outcrops and undergoes thermal and chemical alteration until its discharge as hydrothermal (∼64°C) fluids >50 km away (Wheat and Mottl, 2000; Wheat et al., 2013) (Figure ). Sediment–basement chemical exchange in these locations is apparent from concentration profiles of sulfate, alkalinity, and calcium in the sediment (Supplementary Figure 1), as described elsewhere (Wheat et al., 2013). Increases in sulfate and calcium concentrations occur in bottom sediment toward the sediment-basement interface, indicating diffusion of sulfate and calcium from basement into overlying sediment. Bottom sediment also indicates the imprint of increasing chemical transformation of basement fluids during fluid–rock interaction: the increase in calcium concentration in bottom sediment along the flow path away from Grizzly Bare outcrop is an indication of the exchange reactions where magnesium is stripped from solution to replace calcium in crustal materials (Wheat et al., 2013; Lin et al., 2015).

Maps of sampling locations discussed in this study. (A) IODP Expedition 327 coring locations Sites U1362 and U1363 on the eastern flank of the Juan de Fuca Ridge, shown in proximity to the Grizzly Bare and Baby Bare outcrops. (B) Plan view of Grizzly Bare outcrop coring locations located along seismic line GeoB00-170, with bathymetric relief as shown in legend. (C) Cross-section view of Grizzly Bare outcrop coring locations plotted on seismic profile, with depth of the sediment-basement interface (SBI) at study Holes indicated in meters below seafloor (mbsf). In A, purple diamond represents Hole U1362A; in C, red circle represents U1363F, blue square represents U1363B, and green triangle represents U1363D. Images modified from Fisher et al. (2011) and Wheat et al. (2013) with permission.

In this study, we analyzed sediment microbial composition along a transect away from the edge of the Grizzly Bare seamount where seawater recharge occurs (Figure and Table ), to examine whether microbial diversity in deep subsurface sediment is influenced by basement fluid–rock–sediment interactions. We also aimed to analyze microbial composition in deeply buried basalts from this environment, to see how rock-hosted microbial communities compare to recent studies on the communities found in the altered hydrothermal fluids and rock biofilms in borehole observatories (Cowen et al., 2003; Orcutt et al., 2011a; Baquiran et al., 2016; Jungbluth et al., 2016; Smith et al., 2016) and within microfissures of basalt (Lever et al., 2013). As is common for low biomass studies from deep marine environments (Inagaki et al., 2015; Lee et al., 2015), sequencing contamination presented a challenge to fully resolve basalt microbial communities, despite advances in DNA extraction methods (Lever et al., 2015); however, trends in sediment-basement connectivity were still addressable.

Materials and Methods

Sampling Sites

Samples were collected from the eastern flank of the Juan de Fuca Ridge (JdF; Tables and Figure ) during Integrated Ocean Drilling Program (IODP) Expedition 327 in Summer 2010 in a transect away from the Grizzly Bare outcrop (Expedition 327 Scientists, 2010, 2011c). Samples from Hole U1363F (Figure ) were collected in relatively shallow sediment (<33 m total sediment depth to the basement interface) on the flank of the outcrop. Samples from Holes U1362B and U1363D were collected farther from the edge in sediment roughly 55 and 231 m thick, respectively. These sediments ranged from hemipelagic mud (clays and silts) to sand-silt-clay turbidite sequences to carbonate-rich claystones (Expedition 327 Scientists, 2011c). By contrast, Hole U1362A (Figure ) was located ∼7 km north of the Baby Bare outcrop, a location of hydrothermal fluid discharge, and within a few 100 m of Sites 1026 and U1301 (Expedition 327 Scientists, 2011b). Samples from Hole U1362A represent basaltic oceanic crust that has undergone hydrothermal alteration based on the presence of abundant secondary clays, iron oxides, and localized pyrite deposits (Expedition 327 Scientists, 2011b). For this study, we focused on the deepest available sediment samples from Hole U1363B, D, and F, and two deeper basalt samples from Hole U1362A.

Sample Collection and Handling

During IODP Expedition 327, hard rock samples for microbiological analysis were obtained by rotary core barrel (RCB) coring in Hole U1362A, and extended core barrel (XCB) coring at Site U1363 holes, as described in detail elsewhere (Expedition 327 Scientists, 2011a). Immediately following core delivery, rocks were exposed for subsampling in the core splitting room by shaking the recovered rocks into open split core liners. Rocks for microbiological sampling were identified immediately, photographed in place, and then collected using combusted aluminum foil for transport to the microbiology laboratory. There, rock pieces were transferred to a flame-sterilized stainless steel rock-processing box for subsampling. Following ethanol-rinsing and flame sterilization of the exterior, to remove potential contaminating organisms (Lever et al., 2006), rocks were broken into smaller pieces using flame-sterilized chisels and forceps. Portions of the interior and exterior were transferred into sterile vials containing distilled water for microsphere contamination analysis. The remaining rock pieces were transferred to sterile plastic Whirl-pak® bags and frozen at -80°C until analysis.

Sediment samples for microbiological analysis on IODP Expedition 327 were collected by advanced piston coring (APC), or by XCB coring, as described in detail elsewhere (Expedition 327 Scientists, 2011a). Whole round core (WRC) of 10 cm length were cut on the catwalk and processed using sterilized tools (autoclaved spatulas and ethanol-cleaned end caps). Sterile plastic cut-end syringes were also inserted into the interior and exterior of freshly cut core ends to collect sediment for microsphere contamination analysis and cell density measurements, as described elsewhere (Expedition 327 Scientists, 2011c). The remainder of the WRC was transferred to sterile Whirl-pak® bags, and frozen at -80°C. In the shore-based laboratory, frozen WRCs were subsampled with flame-sterilized tools, transferring the clean interior portions to sterile plastic vials for analysis.

Contamination Testing

To examine potential contamination of hard rock and sediment cores collected during IODP Expedition 327, yellow–green, 0.2-μm-diameter fluorescent microspheres (Fluoresbrite Carboxylate Microspheres, Polysciences, Inc. 15700) were deployed in every core barrel according to standard IODP protocol (Smith et al., 2000), as described elsewhere (Expedition 327 Scientists, 2011a). Perfluorocarbon tracer contamination checks recommended elsewhere (Lever et al., 2006) were not available during this expedition due to use of the tracer pumping apparatus for other purposes. From the sediment cores, sterile cut-end 3 ml plastic syringes were inserted into the interior (center) and exterior (next to core liner) portions of freshly cut core sections, and sediment of a known volume was immediately fixed on the catwalk with 10 ml of cold 3.7% paraformaldehyde in 1× phosphate-buffered saline (PBS) buffer. From the rock samples, rock fragments of known weight were vigorously mixed with a few milliliters of distilled water to liberate microspheres from the rock surface. Samples were stored cold until analysis in a shore-based laboratory. Aliquots of the homogenized samples were filtered onto 25-mm-diameter 0.2-μm-mesh polycarbonate filters. Filters were examined under 100–400× total magnification with epifluorescence microscopy to quantify the number of fluorescent microspheres per samples.

To further evaluate the degree of potential contamination of the samples, a plastic bag that had been used to deliver the fluorescent microsphere solution in a core catcher during one of the basalt samplings, and which was then recovered in the core liner during basement coring operations, was preserved for DNA analysis (referred to as “JdF control”). The bag was removed from the core liner, transferred to a sterile Whirl-pak® bags, and frozen at -80°C. This bag would have been exposed to contamination from human handling, seawater (drilling fluid matrix), drilling mud and drill pipe sources, as well as to the sample environment, and provides a background for evaluating potentially contaminating DNA.

As reported elsewhere (Expedition 327 Scientists, 2011c), the majority of samples showed no evidence of interior contamination, based on the absence of fluorescent microspheres (Table ). Some interior samples from the sediment-basement interface sediment samples (i.e., U1363F-4HCC, U1363D-5X1) or from the exterior portion of a basement core (i.e., U1362A-21R1A) were positive for contamination from drilling fluids based on the presence of fluorescent microspheres.

Cell Density Measurement

Following 6.5 years of storage at -20°C, cell densities in the non-core catcher sediment samples was determined on the same interior sediment samples used for fluorescent microsphere analysis described above, following published cell counting procedures (Kallmeyer et al., 2008). The formaldehyde-fixed sample slurry (1:6 sediment:buffer) was mixed with acid dissolution and detergent solutions to chemically release cells from the sediment matrix, and also gently sonicated to physically release cells from the matrix. Suspended cells were separated from sediment debris by density gradient centrifugation with Nycodenz. The cleaned cell pellet was concentrated on a 0.2-μm-mesh polycarbonate filter and stained with either acridine orange or propidium iodide DNA stains following methods described previously (Orcutt et al., 2011a), then visualized at 100x objective magnification on an Olympus BX60 epifluorescence microscope. Samples from the core catchers (i.e., U1363F-4HCC and U1363F-6HCC) were not preserved appropriately for direct cell counting.

DNA Extraction

Environmental DNA from the JdF sediment (0.2–5 g), basalts (∼0.01–0.02 g), and the “JdF control” plastic bag (small square removed with sterile razor blade) was extracted using a combined mechanical and chemical lysis-based method (Lever et al., 2015). Negative extraction blanks were processed with each batch of samples. For the basalt samples, care was taken to only use material from interior fractures with evidence of alteration, which resulted in the rather small sample size. Using solutions containing the protein denaturants guanidium hydrochloride and tris(2-carboxyethyl)phosphine (TCEP), the surfactant Triton X-100, and dNTP, which minimized sorptive losses of DNA, cells were lysed by a combination of bead-beating, freeze-thawing, and incubation at 50°C, as per the published protocol. DNA extracts were then washed with chloroform-isoamylalcohol (24:1), precipitated with ethanol and linear polyacrylamide, and purified using the CleanAll DNA/RNA Clean-Up and Concentration kit (Norgen, Biotek, Corp., Canada). Attempts to extract DNA using the commercially available MO BIO PowerSoil® DNA Isolation kit were unsuccessful. Quantification of DNA concentrations in the extracts was assessed using the Nanodrop PicoGreen quantification method according to manufacturer instructions.

Quantitative PCR

Quantitative PCR (qPCR) was performed on a Light Cycler 480 Real-Time PCR System and LightCycler 480 SYBR Green I Master mix (both Roche Applied Science, Indianapolis, IN, United States). Two qPCR replicates were run for each sample. Bacterial 16S rRNA genes were amplified using the Bac8Fmod (5′-AGA GTT TGAT YMT GGC TCA G; Juretschko et al., 1998) and Bac338Rabc (5′-GCW GCC WCC CGT AGG WGT; Daims et al., 1999) primer combination. Archaea were amplified using the Arc915Fmod (5′-AAT TGG CGG GGG AGC AC; Cadillo-Quiroz et al., 2006) and Arc1059R (5′-GCC ATG CAC CWC CTC T; Yu et al., 2005) primer combination. Both primer pairs had been subjected to extensive tests, including in silico, melting curve and DNA sequence analyses on cloned qPCR products, to confirm specificity and good phylogenetic coverage of their target domains (Chen, 2014).

Quantitative PCR protocols consisted of 1 × 5 min denaturation at 95°C, 50× cycles of 30 s denaturation at 95°C, 30 s annealing at 55°C, 30 s elongation at 72°C, and 5 s fluorescence acquisition at 80°C, followed by melting curve analyses. qPCR standards were prepared from stock solutions of pGEM-T (Promega Corporation, Madison, WI, United States) plasmids, containing 16S rRNA genes of Bacillus sp. for Bacteria and Methanosarcina sp. for Archaea, that had been quantified using a NanoDrop 1000 spectrophotometer (NanoDrop Products, Wilmington, DE, United States).

The limit of detection, set by the extraction negative control and PCR negative control, ranged from 35 to 285 16S rRNA gene copies μL-1 of extract for Bacteria, and was <1 16S rRNA gene copy μL-1 of extract for Archaea. In any given run, samples in which 16S rRNA copy numbers of both sample PCR replicates exceeded those of both extraction and PCR negative control replicates by a factor of >2 were considered above detection. The mean of the highest negative control (could be either extraction or PCR negative, depending on the run), was then subtracted from the mean of sample replicates to calculate 16S rRNA gene copy numbers. The variability in background contamination was largely driven by the volumes of extraction reagents used per sample, which depended on the amount of sample material used during each DNA extraction.

DNA Sequencing and Analysis

To examine the phylogenetic diversity of the microbial communities in these samples, the hypervariable V4–V6 region of the 16S rRNA gene in the environmental DNA extracts was targeted by 454 tag pyrosequencing. Tags were generated and sequenced at the Josephine Bay Paul Center at the Marine Biological Laboratory according to standard protocols using a 454 Life Sciences GS-FLX with Titanium chemistry (Sogin et al., 2006; Huse et al., 2007; Marteinsson et al., 2012) as part of the Deep Carbon Observatory’s Census of Deep Life. Sequencing included a Bacteria specific V4–V6 primer set (518F: 5′-CCAGCAGCYGCGGTAAN and 1064R: 5′-CGACRRCCATGCANCACCT), an Archaea specific V4–V6 primer set (517F: a mixture of 5′-GCCTAAAGCATCCGTAGC, 5′-GCCTAAARCGTYCGTAGC, 5′-GTCTAAAGGGTCYGTAGC, 5′-GCTTAAAGNGTYCGTAGC, and 5′-GTCTAAARCGYYCGTAGC and 1048R: 5′-CGRCRGCCATGYACCWC), and the Arch806F/Arch958R primer set (Delong, 1992), with 30 cycles of amplification with primer annealing at 60°C. The sequences were trimmed based on quality using the mothur software package as presented in the 454 Standard Operating Procedure (SOP) example analysis (Schloss et al., 2009, 2011). Following quality control, the sequences were trimmed to a minimum of 250 bp and a maximum of 350 bp to allow for a normalized sequence length for further processing.

To conservatively assess for the possibility of contamination of the samples from the drilling and/or handling processes, we followed an approach similar to recent reports (Inagaki et al., 2015; Lee et al., 2015). Sequences were clustered into operational taxonomic units (OTUs) at a >97% sequence similarity and then classified as recommended (Schloss et al., 2011) using the SILVA release version 119 full length and taxonomy reference alignment (Yilmaz et al., 2014) in mothur. The dataset was filtered using a conservative minimum OTU abundance cutoff threshold of 0.005% of total tags (five sequences), mitigating the generation of spurious OTUs (Bokulich et al., 2013); this filtering precludes calculating any richness estimates such as the Chao richness calculations that utilize number of singletons in calculations. The OTUs were then separated into four bins: (1) Control only – OTUs found only in the background contamination control (JdF control), (2) Samples only – OTUs found only in the basalt and sediment samples, (3) 10X in Sample – OTUs found in at least one sediment/basalt sample at >10× abundance as compared to the background sample, and (4) Control: Sample – OTUs found in at least one sediment/basalt sample at <10× abundance as compared to the background sample. All OTUs were also compared to the GenBank and RNA Database Project (RDP) databases to identify potential contaminants from DNA extraction and sequencing kits, as recommended elsewhere (Salter et al., 2014).

To examine the phylogeny of some OTUs, OTU sequences were compared with previously deposited sequences (SSU rRNA) using the RDP v10 Classifier (Cole et al., 2014) and National Center for Biotechnology Information (NCBI) BLASTnt nucleotide database (Altschul et al., 1997). The sequences were aligned with SINA (Pruesse et al., 2012). Phylogeny was performed using neighbor-joining implemented in Geneious (Biomatters, Auckland, New Zealand) with 1000 bootstrap replicated and the HKY model. Trees were constructed with FigTree[1].

Raw sequence data are publicly available through the VAMPS Archive from the Josephine Bay Paul Center under project “DCO_ORC” as well as through the NCBI SRA at accession number SRP051698. Fasta formatted files for OTU bins are available in the Supplementary Material of this manuscript, as well as a table indicating OTU abundance (Supplementary Table 1).

Results and Discussion

Microbial Abundance Based on Cell Counting and qPCR

Microbial abundance in the sediment and basalt samples was analyzed in a variety of ways: direct cell counting using two different staining approaches, and quantitative assessment of 16S rRNA gene copy number in DNA extracts.

Interior sediment core samples were analyzed for cell density measurements following procedures developed for low biomass deep sediment samples (Kallmeyer et al., 2008). Cell density measurements were not performed on basalt samples or core catcher samples, as these samples were not preserved in an appropriate manner for this analysis. Two different cell-staining approaches were used – based on acridine orange (AODC) and prodium iodide (PI) DNA-staining agents – to assess the influence of cell staining technique on cell quantification. In the six sediment samples, cell densities ranged from 2–6 × 106 cells cm-3 to 4–10 × 107 cells cm-3 using the AODC and PI staining methods, respectively (Table ). These values are within the ranges estimated previously for nearby sediment samples from IODP Hole U1301A using microscopic cell staining approaches (Engelen et al., 2008). The difference in cell densities estimated with the two staining approaches is unknown. As all cells were fixed with formaldehyde and rendered dead, difference in staining should not be due to exclusion of stains based on cell membrane integrity. Both stains suffer from non-specific binding to non-cell debris – despite attempts to separate and concentrate cells away from sediment particles using density centrifugation – which can make objective counting of small cell-sized particles challenging.

Following DNA extraction from 0.02 to 5 g sample, bulk DNA concentrations were quantifiable with the Nanodrop-3300 from only the four shallowest sediment samples (U1363F-4H1, U1363F-4H2, U1363F-4H3, and U1363B-8X2; Table ). Despite using optimized methods for DNA extraction (Lever et al., 2015), DNA concentrations in the extracts from basalts and deep sediments were below the limit of detection of ≤103 cells for a sample aliquot of 1 g. Despite the low DNA concentrations in the extracts, quantitative PCR analysis of bacterial and archaeal 16S rRNA genes abundance was robustly determined in comparison to template blanks and negative PCR extraction controls. Bacterial 16S rRNA genes were more abundant than archaeal 16S rRNA genes in all of the samples (Table ). In fact, in none of the three basalt sample extracts were archaeal gene copy numbers above the operational detection limit of >2 times higher gene copies compared to extraction and PCR negative controls.

Sediment samples contained between 4.2 × 103 and 8.4 × 105 bacterial 16S rRNA gene copies per gram, with the exception of one sediment-basement interface sample (U1363F-4HCC) that had 107 copies per gram. However, this one sample was recovered in the core catcher, visibly handled on the catwalk, and had abundant microspheres present; thus, this high gene copy estimate likely reflects external contamination during sampling or sample handling. The basalt samples had bacterial 16S rRNA gene copy densities of 106 copies per gram of vein scraping; however, note that the samples focused on the alteration veins in the samples that likely have higher biomass densities than the bulk rock. Assuming an average RNA operon copy number per bacterial cell of 4 (Stoddard et al., 2015) and a sediment density of 2.5 g cm-3, the 16S rRNA gene copies per gram (Table ) converts to a sediment cell density of 1 × 104–2 × 106 cells cm-3, which is within the low end of the range estimated for other deep sediments (Kallmeyer et al., 2012) and lower than previous direct cell counts in nearby sediment (Engelen et al., 2008). Archaeal rRNA gene copy numbers in sediments were one or more orders of magnitude lower (6.4 × 102–8.4 × 105) than bacterial rRNA gene copies and decreased toward the basaltic basement.

The low cell biomass predicted by qPCR for the sediments samples (Table ) in combination with the lack of amplification of genes in the negative controls provides confidence that contamination of the DNA extracts with exogenous DNA during DNA extraction was unlikely. The approximately two order of magnitude difference in biomass estimates between the microscopic cell counts and qPCR-based approach (Table ), which is a similar pattern to previous observations in nearby sediment (Engelen et al., 2008), may reflect possible overestimation of cell counts based on microscopy due to interference with other stained debris, as discussed above, or underestimation from the qPCR estimates due to inefficient DNA extraction or amplification inhibition. These differences in estimation method are important to consider when scaling such results to construct global inventories of deep biosphere biomass (Kallmeyer et al., 2012). Similarly, the estimates of biomass density in the basalt samples is a reflection of the density within the vein material, and not across the bulk basalt, so appropriate scaling of vein density would be required to convert these numbers to estimates of biomass within the basaltic crust for similar global estimates.

DNA Sequence Analysis, OTU Determination, and Evaluation of Potential Contamination

All samples generated sequences with the Bacteria-specific V4–V6 region 16S rRNA gene primer set (Table ), and two samples (U1363F-4H1 and U1363F-4H3) also generated sequence product with the Archaea specific V4–V6 primer set (marked with asterisk in Table ). Amplification of the archaeal 16S rRNA gene with the Arch806F/Arch958R primer set was not successful (data not shown). Samples yielded between 3,732 and 15,479 quality V4–V6 tags following quality control (Table ). A total of 1,394 OTUs with a >97% sequence similarity threshold were identified. After filtering for possibly spurious OTUs, 587 of these OTUs were comprised of >5 tags and were retained for further analysis (Supplementary Table 1).

To assess for possible contaminant sequences in the sample libraries – derived from either collection, handling, or sequence amplification sources – all of the 587 OTUs that passed quality control were divided into four sequence groups: (1) OTUs found only in the contamination JdF control (control only), (2) OTUs that were only present in samples (samples only), (3) OTUs that were at least 10 times more abundant in samples than in control (normalized to a common number of sequences per sample) (10X sample), and (4) OTUs that were less than 10 times more abundant in samples than in controls (control:sample). The third group (10X sample) represents cells that are relatively more abundant in the samples than in the control, or that preferentially over-amplified in the samples compared to the control. The fourth group (control:sample) represents cells that are found in both the samples and in the control. All sequences found in groups 3 and 4 could have been transferred from the subsurface environment to the control during collection or during laboratory handling.

The JdF control sample generated 3,732 bacterial sequence reads after quality control filtering, which clustered into 164 OTUs at 97% sequence similarity level, with 116 OTUs unique to the control and considered as sequence group 1 (Figure and Table ). No archaeal sequences were recovered from the JdF control sample from either primer set, even though Archaeal 16S rRNA genes were detected via qPCR (Table ). The majority of the unique-to-the-control bacterial OTUs grouped within the Gammaproteobacteria (25.6% of reads), Alphaproteobacteria (16.8%), Flavobacteriia (Bacteriodetes; 16.0%), and Cyanobacteria (12.8%) classes, with minor contributions from other classes. The majority of these OTUs are similar to bacteria commonly found in seawater, as demonstrated by the three most abundant OTUs belonging to SAR86 (Gammaproteobacteria; Dupont et al., 2012), SAR11 (Alphaproteobacteria; Rappé et al., 2002), and Flavobacteria NS4 group (Alonso et al., 2007). The dominant bacterial phyla in sequence group 1 were markedly different from the dominant phyla observed only in the sediment and basalt samples as sequence group 2 (Figure ). The remaining 48 OTUs from the control sample that comprised about 20% of the sequences (Figure ), and which were also observed in the samples as sequence groups 3 and 4, were predominantly phylogenetically related to known sequencing kit reagent contamination or to common seawater groups, as discussed below.

Taxonomic distribution at the class level of the bacterial tag sequences found in (A) the JdF control sample only (control only) and (B) in the sediment/basalt samples only (samples only). Colored symbols on right-hand side in B indicating sample location the same as in Figure . Number in parentheses in legend indicates how many OTUs are in each class. S, sediment; B, basalt.

Taxonomic distribution at the class level of the Bacteria tag sequences found in both the JdF control sample (a plastic bag used to deliver the fluorescent microsphere solution in a core catcher during basalt coring) and the sediment and basalt samples. Taxonomic groups that are present in both the control and samples, and highly abundant in the sediment/basalt samples (i.e., 10X more sequences than in JdF control, 10X sample), are represented in bold on the right side. Taxonomic groups that are in relatively the same abundance in both the sediment/basalt samples and in the JdF control are on the left side (control:sample). The sequences that were unique to samples are represented in white, indicating the fraction of sequences that were kept for comparison analyses (samples only). Number in parentheses in legend indicates how many OTUs are in each class. Colored symbols on right-hand side in B indicating sample location the same as in Figure .

Between 2 and 100% of the bacterial sequences from each sediment and basement sample (Figure ) had overlap with the JdF control sample [i.e., binning into sequence group 3 (10X sample) or 4 (control:sample); Figure ], with higher overlap in samples with lower DNA extract concentrations (i.e., the deepest sediment samples and the basalts). Sequence group 3 was much more abundant in the sediment and basalt samples than sequence group 4, which made up <0.5% of sequences (Figure ). The sequence group 4 OTUs (control:sample) were phylogenetically related to microbial groups that are commonly found in seawater, suggesting that these sequences could represent slight seawater contamination of the samples, or microbial groups that are common in seawater, sediment, and/or the basalt subsurface. Because of the possibility of these sequences potentially representing cross-contamination of seawater-derived groups during sample handling or sequencing, they were removed from further consideration of microbial community structure.

All of the sequence group 3 (10X sample) OTUs (shown in bold boxes in Figure ) had high sequence similarity with common cultured isolates known as laboratory contaminants of DNA extraction kits and laboratory reagents (Salter et al., 2014) (Supplementary Table 1). For example, the abundant Actinobacteria OTUs from this sequence group in the basalt were closely related to Arthrobacter and Rhodococcus isolates from laboratory clean rooms and humans, respectively (Supplementary Figure 2). As another example, Firmicutes-related OTUs were observed in all of the samples, although in relatively low abundance (Figure and Supplementary Table 1). The majority of these OTUs grouped with known laboratory contaminants from human metagenomes (Bacillales, Clostridiales, Lactobacillales, and Selenomonadales) and other suspect contaminant groups (Ruminococcaceae, Halanaerobiales 64K2). The sequence library from sample U1363F-4HCC, which tested positive for contamination with the fluorescent microspheres (Table ), was highly enriched (>99% of sequences) with Gammaproteobacteria Enterobacteriales OTUs also found in the background control (Figure ). Thus, the sequence group 3 OTUs were categorized as probable sequencing contaminants and excluded from further analyses of microbial community structure; all of the sequences from the suspect U1363F-4HCC sample were not considered further. In summary, only sequence group 2 (samples only) sequences were considered for further microbial community structure analysis of the sediment and basalt samples.

Sediment Microbial Community Structure

The sediment sample sequence libraries displayed little overlap with the background control sample that was collected during basalt coring operations (Figures ). The majority of sequences grouped within the Chloroflexi phylum (60.3–86.7% of all QC reads; Figure and Supplementary Table 1). The Chloroflexi-related OTUs grouped into seven different defined classes and two unclassified groups, with the majority of these grouping in the Dehalococcoidia class (Figure and Supplementary Table 1). Phylogenetic analysis of all 136 Dehalococcoidia-related OTUs showed that all were highly similar to environmental sequences from other sediment environments and distinct from cultured isolates (Supplementary Figure 2). In an earlier clone library-based study from sediment squeeze cakes and resulting porewater from Hole U1363B deep sediment samples, Chloroflexi sequences were also detected, as well as in deeper re-sequencing of some of those samples (Jungbluth et al., 2013, 2016). Chloroflexi from the terrestrial subsurface are known to be involved in dehalogenation as a metabolic process, and may include sulfur cycling (Wasmund et al., 2016), but the function of Chloroflexi in organic-rich deep sediment are not verified (Kaster et al., 2014; Wasmund et al., 2014; Fullerton and Moyer, 2016).

Diverse Proteobacteria were observed in all sediment samples at relatively low abundance (<5% of QC tags; Figure and Supplementary Table 1); however, the phylogenetic classification of some of these OTU sequences suggests that they should be viewed with caution. For example, Alphaproteobacteria-related sequences were the most abundant, grouping within the Rhizobiales and Rhodospirallales (Figure ). Rhizobiales are commonly found in soils, and the closest environmental relatives of the Rhodospirallales are from oxic/suboxic sediment from the South Pacific Gyre. Similarly, the Betaproteobacteria tags were related to Burkholderiales, which is known as a water filtration contaminant that can stymie environmental genomics surveys, although it is also found in sediment. The Gammaproteobacteria sequences were most closely related to Xanthomonadales, which is known to be enriched in drilling fluids (Masui et al., 2008). Only the Deltaproteobacteria sequences, grouped into five OTUs related to Desulfobacterales, were closely related to environmental sequences from anaerobic, sulfate-reducing sediment.

Other bacterial groups in relatively high abundance in some sediment samples included the Planctomycetes, Nitrospirae, and the recently described Aerophobetes and Atribacteria groups of uncultivated organisms (Figure and Supplementary Figure 3). Planctomycetes had higher abundance in the clay-rich samples, and the majority of the sequences grouped mostly closely to Phycisphaerae clade MSBL9 observed in other deep sediment samples (Inagaki et al., 2006; Durbin and Teske, 2011), where they can be associated with fermentation at the sulfate-methane transition zone (Harrison et al., 2009). This same Planctomycetes group was also observed in an earlier clone library from Hole U1362B deep sediment (Jungbluth et al., 2013). Two Nitrospirales order OTUs were present in all clay samples, but the function of Nitrospirae in marine sediment is unknown. Nitrospirae are indicator taxa of altered crustal fluid in this environment (Jungbluth et al., 2016). While they were not detected in previous studies of sediment from the Juan de Fuca Ridge flank (Jungbluth et al., 2013, 2016), the sequences in our dataset are quite divergent from the sequences identified in the crustal fluid (78% similarity). Similarly, Atribacteria (candidate division JS-1) OTUs were in high abundance in the shallower clay samples from Hole U1363F and U1363B (∼17% of all QC tags), although not in the deeper clays from Hole U1363D. This group was also observed in earlier clone libraries from Hole U1363B deep sediment (Jungbluth et al., 2013) as well as recent re-sequencing of those DNA extracts (Jungbluth et al., 2016). Atribacteria are common in anaerobic, organic-rich sediment, where they are thought to be involved in organotrophic processes including fermentation (Dodsworth et al., 2013; Carr et al., 2015). All sediment samples had Aerophobetes sequences (previously known as BH180-139), with the highest concentration in the foram-rich carbonate sample U1363B-8X5 (24.8% of all QC sequences). This group was also detected in an earlier clone library of a deep sediment sample from Hole U1363B, as well as in shallower samples (Jungbluth et al., 2013, 2016). Very little is known about Aerophobetes except the suggestion that they may be strict anaerobes (Rinke et al., 2013), which could explain their higher abundance in deeper samples. Other Bacteria phyla were observed in the sediment samples, but never in high abundance (<5%, at most, of all Bacteria tags in a sample, Supplementary Table 1). These rare phyla included Acidobacteria, Bacteriodetes, Gemmatinomadetes, Spirochaetes, and candidate divisions BRC1, KB1 (Calescamantes), OD1, OP11, OP3, OP8 (Aminicenantes), WS3, and TA06. Although sediment samples shared many of the same bacterial groups, there was sample-specific clustering of the OTUs found in these phyla (Supplementary Figure 2).

Six of the eight sediment samples yielded a few archaeal sequences from the sequencing runs with the Bacteria-specific V4–V6 primers, and two sediment samples (U1363F-4H1 and -4H3) generated thousands of QC sequences from the Archaea-specific primers (Table ). Note that no archaeal sequences were recovered from the JdF control sample, as discussed above, so it is not possible to evaluate the potential for cross-contamination in this dataset. Both samples sequenced with the Archaeal-specific primer set yielded sequences related to Thermoplasmata and Thaumarchaeota, although the Thaumarchaeota groups were different (Figure ). Previous clone libraries from deep sediment samples from Hole U1363B also detected Archaea, but from the Bathyarchaeota/Miscellaneous Crenarchaote Group (Jungbluth et al., 2013), which was not detected in the present study, while a more recent study of a “suite” of Site U1363 sediment samples also detected Marine Benthic Group E/Thermoplasmata (Jungbluth et al., 2016). The Archaea sequences picked up by the Bacteria-specific primer set also grouped within the Euryarchaeota and Thaumarchaeota (Figure ). No sequences were recovered that were related to methanogens or anaerobic methanotrophs, which have been sequenced in prior studies of basalts from the Juan de Fuca subsurface using assays targeting the alpha subunit of the gene for methyl coenzyme M reductase (Lever et al., 2013).

Taxonomic distribution at the class level of the Archaea tag sequences found in the sediment samples. No archaeal 16S rRNA genes could be sequenced from subseafloor basalt samples or the control. In bold are the classes that were found in higher abundance and are discussed in the main text. Number in parentheses in legend indicates how many OTUs are in each class. Due to the low number of reads in some samples, the graph is displayed as actual reads instead of being normalized to 100%. Colored symbols on at bottom indicating sample location the same as in Figure .

(A) Bacteria OTUs shared amongst more than one sample, ranked by abundance, and with color intensity reflecting abundance of OTU in the sample (top scale of legend in B indicating OTU abundance). (B) Heat map cross-comparison of sample community similarity using the Jaccard Index, with similarity levels scaled in bottom of legend.

Basalt Microbial Community Structure

The majority of sequences from the basalt samples grouped within sequence group 3 (i.e., they were also observed in the JdF control sample; Figure ), and these were predominantly related to known isolates of Actinobacteria from clean room and human origin (Figure and Supplementary Table 1). Of the remaining sequences unique to the basalt samples (i.e., sequence group 2), the majority (<95%) of sequences were also similar to Actinobacteria (Figure ). Less than 5% of the sequence group 2 basalt sequences grouped within other bacteria groups like the Firmicutes, Alphaproteobacteria, and Gammaproteobacteria (Figure ). No archaeal sequences were recovered from the basalt samples.

Phylogenetic analysis of the 16S rRNA gene from Actinobacteria OTUs from this study (Blue: samples only; Green: enriched in samples compared to control; Orange: similar abundance in samples and control; Red: control only), cultured microorganisms (black), and other closely related environmental sequences (gray). This tree shows the close phylogenetic relationship between the sequences from this study to common isolates as well as environmental sequences originating from human (mostly skin) and soil microbiome sequences, suggesting that the majority of the Actinobacteria-related sequences identified in our samples are laboratory contaminants. The tree was generated using neighbor-joining with 1000 bootstrap replicates, using the Jukes-Cantor model. Bootstrap replicates >90% are indicated by a black dot at the node. The first and second column represent the abundance of reads in this study (red; in log) and if the sequences were found in sediments only (blue), basalts only (blue), control only (orange), sediments and basalts (purple), or all three (brown). The remaining columns show the proportion of sequence found in each cluster that origin in our study (JdF), cultured microorganisms (C), human microbiome (H), soils (S), subsurface environments (Su), aquatic environments (A), or others (O). The majority of the sequences clustered into four supported groups where clades related to sediments are highlighted in purple and the others in gray.

Phylogenetic analysis of the Actinobacteria OTUs from this study compared to other studies provides a measure to assess the likelihood that an OTU represents a handling or sequencing contaminant or a true subsurface microbial group (Figure ). Only three OTUs grouped only with sequences from other subsurface environments (i.e., purple box in Figure ), and these three OTUs were only observed in sediment samples. The remaining Actinobacteria sequences in this study, including all of those from basalt samples, were most similar to either cultured microorganisms, or uncultured microorganisms found in various environments including human skin and soil/dust microbiomes (Figure ). The most abundant OTU found in basalts (OTU0001, 56.1% of all reads from basalts; 2.8% of the reads from sediments) was found in the control sample as well (sequence group 3) and has >99% sequence similarity to the type strain isolate Arthrobacter sp. TT4 (GenBank accession number JF706635; Figure ), a strain isolated from a clean air space laboratory (Li et al., 2004). Arthrobacter is a common laboratory contaminant of reagents and DNA extraction kits (Salter et al., 2014). It is thus highly likely that this OTU represents a sequencing contaminant and should not be considered representative of the subsurface. In contrast, the second most abundant OTU in basalt samples (OTU0002; 37.5% of the basalt reads; 0.7% of the sediment reads) was not observed in the control sample but has >99% sequence similarity with the cosmopolitan strain Rhodococcus erythropolis (GenBank accession number X79289), which is known to cause infection in humans. Notably, this was the only OTU that was shared with another sample library from the Census of Deep Life – in particular, with sequences from a subseafloor basalt (data not shown; Jason Sylvan, personal communication). As this OTU was absent from our control sample and sequencing blanks analyzed in the other study, this may indicate that this Rhodococcus-related OTU has a subsurface origin, even though it has cosmopolitan distribution in disparate environments (Figure ). All of the other Actinobacteria OTUs were far less abundant in the samples, and they also had phylogenetic similarity to isolates and uncultured species from various environments (Figure ).

Resolving the subsurface and/or basaltic origin (or not) of the Actinobacteria from the basalt samples is challenging with the current dataset. Considering that our study only examined the 350 bp V4–V6 region of the 16S rRNA gene, it is possible that the OTUs with close V4–V6 region sequence similarity to sequences from dust and skin microbiomes may reveal variation when comparing larger sequence regions. Actinobacteria have been observed under cool and oxic conditions on seafloor and subseafloor basalts (Lysnes et al., 2004; Mason et al., 2009; Lee et al., 2015; Jørgensen and Zhao, 2016), in cool crustal fluids (Meyer et al., 2016), and in weathered volcanic glass (Cockell et al., 2009, 2013). Many of the Actinobacteria groups found in those cool and oxic studies, as well as in the current study (e.g., Arthrobacter, Rhodococcus, Corynebacterium) are cosmopolitan, and occur widely in terrestrial, marine, aeolian, animalian, oxic, and anoxic environments (Figure ). For example, parsimony-based partial (non-overlapping) alignment of our V4–V6 tag sequences, the amplicon Actinobacteria sequences found in the recent analyses of the Mid-Atlantic Ridge flank basaltic subsurface (Jørgensen and Zhao, 2016; Meyer et al., 2016), and full length 16S rRNA gene data from environmental and cultured sequences reveals the cosmopolitan association of these groups, including the OTUs from this study that cannot be confirmed to be from the subsurface (data not shown). Actinobacteria have also been detected in the warm and anoxic crustal fluids from the Juan de Fuca (Jungbluth et al., 2012), but those specific lineages were not detected in the present study (data not shown).

Most known Actinobacteria have aerobic metabolisms, so finding them in the previous examinations of oxic environments is not surprising, whereas explaining their presence in the warm and anoxic JdF subsurface may be less straightforward. However, some Actinobacteria, including Arthrobacter, have been linked to anaerobic metabolisms including nitrate reduction and fermentation, both in the environment and in culture (Eschbach et al., 2003; Süß et al., 2006). Actinobacteria also have the ability to form spores, and could thus survive for long periods outside of suitable habitat without being metabolically active. Some Actinobacteria (Actinomycetes) form filaments, which could enable actinobacterial growth in rock environments in a manner previously suggested for Fungi (Ivarsson et al., 2016). Thus, Actinobacteria in anoxic deep biosphere environments could either be metabolically inactive survivors (e.g., remnants from aerobic surface sediment or seawater communities), or active in anaerobic metabolic processes. Further research is required to determine if subsurface Actinobacteria are metabolically active in subseafloor basalt, and the metabolic roles they play in this environment.

Microbial Community Transitions in Basement-Influenced Sediment

Beyond determining the overall structure of the microbial communities in deep sediment and basement of the Juan de Fuca Ridge flank – which had not been fully addressed in previous studies of this system (Engelen et al., 2008; Lever et al., 2010, 2013; Jungbluth et al., 2013) – the aim of this study was to evaluate the influence of basement on the microbial ecology of deep sediment. Does deep sediment close to areas of seawater recharge into basement exhibit signs indicating that this influences the structure of the microbial communities? Or does the complex nature of fluid flow at outcrops, with recharge into basement occurring at the same time as discharge (Wheat et al., 2013; Winslow and Fisher, 2015), overprint this influence?

In this study, samples from Holes U1363F, U1363B, and U1363D progress in their distance away from the site of recharge (Figure ), with samples U1363F-4H3, U1363B-8X5, and U1363D-5X1 being the sediment samples closest to the sediment-basement interface, thus comprising the set of samples that might display a trend. At a bulk level, these three samples are comprised of similar proportions of dominant Bacteria phyla – Chloroflexi are by the far most abundant sequence group, followed by Aerophobetes (Figure ). However, at the “species”/OTU level of taxonomy, the deeper and farthest U1363D-5X1 sample appears to be comprised of a different constellation of OTUs as compared to the shallower and more proximal samples (Figure ). While all of the sediment-basement interface samples share similar sulfate concentrations (26–27 mM, Table and Supplementary Figure 1), as well as similar concentrations of many other porewater constituents (Wheat et al., 2013), the shallower and more proximal samples have cooler temperatures (7–12°C) as compared to the distal and deeper location (33°C, Wheat et al., 2013), which may explain part of the variation in speciation. Notably, the deep sediments collected just a few 10s of meters away from the edge of net recharge outcrop do not indicate evidence of bottom seawater influence on microbial community structure following the conservative approach taken here (Figure ), nor does the chemistry of sediment porewaters (Wheat et al., 2013). This is a conservative assessment, given that a few percent of sediment sequences (0.9–3.7%) fell into the potential seawater contaminant category, so comparison to the background sample may have artificially removed these sequences if they were present in both sediment samples and the background sample (Figure ).

Because of our inability to confirm the origin of the sequences in the basalts, we cannot directly compare the sediment sequence libraries to the basalt sequence libraries to determine basement influence on sediment communities, and vice versa. However, recent deep sequencing work on highly altered crustal fluids from boreholes on the Juan de Fuca Ridge flank (Jungbluth et al., 2016) does provide an opportunity to assess the degree of basement fluid influence on basal sediment community structure. Although the fluid study was conducted on fluids in the vicinity of Hole U1362A, and not directly underneath the Site U1363 holes 52 km away, the porewater chemistry at Site U1363 is indicative of warm and anoxic fluids in basement similar to the Hole U1362A region. When comparing the sediment sample sequences in the present study to the sequences identified in the crustal fluids, the bacterial communities of the sediments are quite divergent from the indicator taxa identified in the crustal fluids, indicating a segregation between the microbial communities found in the sediments and in the fluids. For example, there was only 88–90% sequence similarity and <80% similarity to the indicator taxa Aminicenantes (OP8) and Calescamantes (BK1), respectively (data not shown).

The microbial communities living in deeply buried sediments have been cut off from new detrital energy inputs for thousands to millions of years and need to adapt to constant energy and nutrient limitation. The Juan de Fuca Ridge is influenced by hydrothermal activity and fast sedimentation rates, which influences energy availability (LaRowe and Amend, 2014). The Chloroflexi phylum dominated the sediment samples, suggesting that these microbes are well-adapted to organic rich anoxic sediments. However, there were differences amongst the less dominant groups (Figures ), and these differences in the microbial communities could be explained by the physical and chemical characteristics of the sediments. Nitrospirae was more prevalent in clays, indicating that the type of material is an important part of microbial adaptation. Actinobacteria OBP1 was more prevalent only in one sample, U1363F-4H1, which is the shallowest sample and with the “freshest” organic matter. Atribacteria and Planktomycetes were present only in the clay samples that were closer to the recharge area, while Aerophobetes were more abundant in the clay samples farther from the recharge area. The sediments display a gradient of temperature, which could be the main factor for the changes in microbial communities. Finally, there are still many unknown groups of bacteria that inhabit these sediments that remain uncharacterized, with a higher proportion of unclassified tags in the deepest samples.

Future Improvements for DNA-Based Assessment of Deep Sediment and Crust

Our DNA sequence data suggest DNA contamination in at least one of the sediment samples and all the basalt samples analyzed in this study. The sediment-basement interface sample of Hole U1363F (i.e., 4HCC), was contaminated during sample drilling and handling, as evidenced by the high concentration of fluorescently labeled microsphere at sampling, and the overwhelming dominance of Enterobacterales-related sequences (Table and Figure ). The basalt samples and low biomass sediment samples (U1363B-8X5 and U1363D-5X1) were likely contaminated during sample handling, DNA extraction or sequencing due to their very low biomass, as demonstrated by the enrichment of common laboratory reagent and airborne contaminants. These observations confirm that even with improved methods for core handling and sample processing (Expedition 327 Scientists, 2011a), DNA contamination of deep biosphere samples remains a challenge (Santelli et al., 2010; Lever, 2013; Inagaki et al., 2015).

Based on our results, the use of fluorescent microsphere proved to be a useful – though not quantitative (House et al., 2003) – tool to detect drilling fluid contamination, but more efforts need to be made to properly address contamination in low biomass samples. DNA contamination was also an issue in recent sampling of subsurface basalts from North Pond (Jørgensen and Zhao, 2016; Zhang et al., 2016) and from other low biomass environments where over 99% of raw sequences were suspect (Inagaki et al., 2015). For future expeditions, we suggest taking water column and drilling fluid samples as contamination controls and sequencing both separately, to enable more reliable distinctions between DNA sequences in drilling fluids that were introduced from the subsurface environment from ones that were introduced from surface seawater or during core handling. We encourage the systematic sequencing of a blank DNA extraction control to correctly remove contaminants that are present in the DNA extraction and sequencing reagents, as was done during IODP Expedition 337 (Inagaki et al., 2015). For extremely low biomass basalt samples, it will be important to collect larger rock samples for DNA extraction, as was done in several other studies, where larger samples were available for microbiological analyses (Lever et al., 2013; Jørgensen and Zhao, 2016; Zhang et al., 2016). This will increase the amount of sample material from rock interiors, which tend to be minimally contaminated (Lever et al., 2006). Increasing the amount of suitable sampling material will enable lower detection sensitivity of in situ microbial communities, lower the risk of significant contamination by laboratory handling, and maximize the sequencing of DNA coming from the sample. Finally, we caution all users of publically available sequence datasets, such as those from the Census of Deep Life, to critically evaluate for the possibility of sequence contamination in the samples.

Another strategy that will improve the reliability of nucleic acid data and help identify DNA sequences from indigenous microorganisms is to combine multiple interdisciplinary approaches. In the past, DNA sequence data have been combined with geochemical (e.g., electron donor and acceptor availability), isotopic (providing clues to ongoing microbial metabolic processes), and mineralogical data (e.g., redox/oxidation state, olivine content), as well as laboratory enrichments to demonstrate microbial carbon and sulfur cycling in subseafloor basalt (Lever et al., 2013). Microbial colonization experiments involving different rock types, and combining molecular biological and microbiological analyses with monitoring of surrounding fluid chemistry and rock alteration processes in CORKS have provided new insights into microbial colonization patterns and how rock biofilms form (Orcutt et al., 2011a; Smith et al., 2011, 2016; Baquiran et al., 2016). Future studies on unaltered rock samples and after cultivation experiments may include single-cell genome, metagenome, and metatranscriptome sequencing, and compare potential metabolisms inferred from genomic information to data on potential in situ microbial activities obtained with other methodologies to distinguish in situ microorganisms from microorganisms that were likely introduced during sampling, sample handling, or any other downstream laboratory-based research activities.

Conclusion

Our study provides a glimpse of microbial community dynamics and interactions across sediment-basement interfaces in the marine deep subsurface. Our results demonstrate that proximity of basal sediment to a site of seawater recharge into basement bears little correlation with sediment microbial community structure, which is more strongly influenced by depth below the seafloor and sediment type. Moreover, there appears to be little overlap in microbial community structure between basal sediments and basement fluids. The composition of basal sediment microbial communities is characteristic of anoxic, relatively organic rich subsurface sediment environments elsewhere, with dominance by subseafloor Chloroflexi, Aerophobetes, Atribacteria, and Nitrospirae. By contrast, the full phylogenetic composition of microbial biofilms within upper basement basalts from the Juan de Fuca subsurface remains elusive, unfortunately, due to typical contamination-related challenges of research on low biomass samples and due to strong phylogenetic overlaps not only with other basalt habitats but also with a wide range of other environments. Future studies involving deep sequencing of drilling fluids and all negative controls, using larger rock samples with a higher biomass of indigenous microorganisms inside, and complementing nucleic acid sequencing with other methodologies that provide insights into microbial energy sources and activities, will contribute to a more clear and accurate picture of the indigenous basalt microbiota.

Author Contributions

BO and KE designed the project, BO collected the samples, ML performed laboratory analyses, JL and ML performed data analysis, and JL and BO wrote the paper with input from coauthors.

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.

Table 1

Characteristics of the IODP Expedition 327 Juan de Fuca Ridge flank sampling locations discussed in this study.

ID	Location	Latitude (N)	Longitude (W)	Water depth (m)	SBI depth (mbsf)	SBI temperature (°C)
Hole U1362A	Second Rift	47.761047	127.761200	2661	236	64
Hole U1363B	Grizzly Bare Outcrop	47.289197	127.035100	2678	57	12
Hole U1363D	Grizzly Bare Outcrop	47.292873	128.029332	2678	231	33
Hole U1363F	Grizzly Bare Outcrop	47.288769	128.035623	2678	35	7

Table 2

Summary of IODP Expedition 327 sample depths, sediment sulfate concentration, sample descriptions, and potential for contamination (based on presence of fluorescent microspheres used during coring operations).

Core	Depth (mbsf)	SO4 (mM)	Description	μspheresa
U1363F-4H1	29.7	25.0	Clays	0
U1363F-4H2	31.2	25.6	Clays and Mn crusts	0
U1363F-4H3	31.9	26.2	Foram-rich and Mn crusts	0
U1363F-4HCC	32.3	25.2	Mn crusts, basalts, forams	1.3 × 104
U1363B-8X2	50.2	26.6	Clays	0
U1363B-8X5	53.6	27.5	Foram-rich carbonate	0
U1363D-5X1	222.3	27.3	Clays	5.8 × 103
U1363D-6XCC	231.2	n.a.	Altered massive basalt	n.m.
U1362A-16R2	456.7	n.a.	Altered massive basalt	0
U1362A-21R1A	491.4	n.a.	Basalt with chilled margin, exterior	7.9 × 101
U1362A-21R1B	491.4	n.a.	Basalt with chilled margin, interior	0
U1362A-JdF control	n.a.	n.a.	Contamination check	n.m.

Table 3

Summary of cells and sequence groups in IODP Expedition 327 samples, including density of cells per cubic centimeter of sample as measured using acridine orange (AODC) and propidium iodide (PI) staining methods, the concentration of bacterial and archaeal 16S rRNA genes per gram sample, the number of sequences that passed initial quality control (QC seqs), the number of predicted operational taxonomic units at the 97% or greater sequence similarity threshold (# OTUs), and the number of OTUs in each sample after screening for potential contamination (# OTUs-control).

Core	Cells cm-3 AODC	Cells cm-3 PI	Bacteria 16Sb (copies/g)	Archaea 16Sc (copies/g)	QC seqsd	# OTUse	# OTUs-control
U1363F-4H1∗	5 ± 3 × 106	5 ± 2 × 107	3.1 × 105	2.3 × 104	4,547	128	115
U1363F-4H2	4 ± 3 × 106	5 ± 1 × 107	2.4 × 105	8.6 × 103	15,582	257	239
U1363F-4H3∗	3 ± 3 × 106	5 ± 3 × 107	8.4 × 105	2.7 × 104	7,086	100	93
U1363F-4HCC	n.m.a	n.m.	1.4 × 107	6.4 × 102	15,479	10	3
U1363B-8X2	4 ± 3 × 106	1 ± 0.4 × 108	2.2 × 105	2.9 × 103	12,882	185	169
U1363B-8X5	6 ± 4 × 106	4 ± 2 × 107	4.9 × 103	1.3 × 102	10,634	89	65
U1363D-5X1	2 ± 2 × 106	9 ± 4 × 107	4.2 × 103	3.1 × 103	11,837	119	96
U1363D-6XCC	n.m.	n.m.	6.9 × 105	3.0 × 103	11,344	21	6
U1362A-16R2	n.m.	n.m.	1.1 × 106	2.3 × 103	9,081	24	10
U1362A-21R1A	n.m.	n.m.	1.1 × 106	ND	9,057	24	11
U1362A-21R1B	n.m.	n.m.	1.4 × 106	2.3 × 103	8,232	21	3
U1362A-JdF control	n.m.	n.m.	5.8 × 106	6.4 × 104	3,732	164	116f