Dataset: North Atlantic subseafloor clay metagenomes and metatranscriptomes

Methodology: We optimized a DNA extraction protocol for ultra-organic lean sediments that provided increased yields of DNA, enabling sequencing of 16S rRNA and ammonia monooxygenase [amoA] genes, as well as metagenomes, from two deep oxic subseafloor sediment cores reaching up to ~15 million years old. In addition, we demonstrate the viability of uncultivated microbial populations via ¹⁸O-labeling in long-term (18-month) incubations in the presence of H₂¹⁸O, a method for identifying actively growing microbes.

Sampling and analytical procedures: All samples were taken by Cruise KN223 of the R/V Knorr, in the North Atlantic, from October 26th 2014 - December 2014 (Woods Hole, MA - Woods Hole, MA). At both Site 11 (22° 47.0’N, 56° 31.0’ W, water depth ~5,600 m) and Site 12 (29° 40.6’ N, 58° 19.7’ W, water depth ~5,400 m), successively deeper subseafloor cores were taken with a multicorer [to ~0.4 meters below seafloor (mbsf)], a gravity corer (to ~3 mbsf), and the 45-m WHOI "long corer" piston-coring device (https://www2.whoi.edu/site/longcore/) (to ~28 mbsf). Subsamples from the core sections for DNA extraction were sampled on board the ship immediately after retrieval with sterile 60 mL syringes with the Luer taper end cut off and frozen immediately at -80 °C prior to DNA extraction. Subsamples for the 18O-labeling experiment were sampled from the core sections in the same manner, but were stored at +4 °C prior to incubation set up.

DNA extraction: Subcores sampled aseptically with sterile syringes were sampled aseptically in a UV-sterilized DNA/RNA clean HEPA filtered laminar flow hood. To reduce contamination, the paraffin was removed and the outer 3 cm of sediment pushed out of the syringe, which was then sliced off with a red-hot, sterile spatula. A second unused, sterile, spatula was used to carefully sample the uncontaminated center of the remaining core sample inside the syringe. In brief, 10 g of sediment was transferred into 50 mL Lysing Matrix E tubes (MP Biomedicals) containing silica glass beads and homogenized for 40 sec at 6 m/s using a FastPrep 24 5-G homogenizer (MP Biomedicals) in the presence of 15 mL preheated (65 C) sterile filtered extraction buffer (76 vol% vol 1M NaPO₄ pH 8, 15 vol% 200 proof ethanol, 8 vol% MoBio's lysis buffer solution C1, and 1 vol% SDS). The samples were incubated at 99°C for 2 minutes and frozen overnight at -20°C, thawed, and frozen again at -20°C overnight, followed by additional incubation at 99°C for 2 minutes and a second homogenization using the settings described above. The additional freeze thaw steps, particularly freezing overnight, was found to increase DNA yield 2–10 fold. After the second homogenization, the samples were centrifuged for 15 minutes, and the supernatants concentrated to a volume of 100 uL using 50 KDa Amicon centrifugal filters (Millipore). Co-extracted PCR-inhibiting humic acids and other compounds were removed from the concentrated extract using the PowerClean Pro DNA Clean-up Kit (MoBio). Extraction blanks were performed alongside the samples to assess laboratory contamination during the extraction process.

qPCR: DNA was quantified fluorometrically using a Qubit with a dsDNA high-sensititivity kit (Life Technologies). Quantitative PCR (qPCR) was performed using the custom primer dual indexed approach that targets the V4 hypervariable region of the 16S rRNA gene using updated 16S rRNA gene primers 515F/806R (515F: 5' – GTGYCAGCMGCCGCGGTAA– 3', 806R: GGACTACNVGGGTWTCTAAT) that increase coverage of ammonia oxidizing Thaumarchaea and other marine strains(59). To measure the abundance of amoA genes from archaea, the primers Arch amoA-1F (STAATGGTCTGGCTTAGACG) and Arch amoA-2R (GCGGCCATCCATCTGTATGT) were used. qPCR reactions were prepared using an automated liquid handler (pipetting robot), the EpMotion 5070 (Eppendorf), was used to set up all qPCR reactions and standard curves. The efficiency values of the qPCR was <90% and R² values >0.95. qPCR was performed using white 96-well plates as this was found to increase the signal to noise in the SYBR green assay 2-fold compared to clear plates. The technical variability of 16S rRNA gene qPCR measurements was determined to be consistently <5% under the the EpMotion 5070.

16S rRNA and amoA gene sequencing: Barcoded V4 hypervariable regions of amplified 16S rRNA genes were sequenced on an Illumina MiniSeq. This yielded a total of >20,000,000 raw sequencing reads that were then subjected to quality control. In order to quality control the OTU picking algorithm for the data, we also sequenced a "mock community" alongside our environmental samples. The mock communities contained a defined number of species (n=18) all containing 16S rRNA genes >3% difference. USEARCH version 10.0.240 was used for quality control and OTU picking (61), OTUs were clustered at 97% sequence identity. The taxonomic relationship of OTU representative sequences were identified by BLASTn searches against SILVA database (www.arb-silva.de) version 128. To identify contaminants, 16S rRNA genes from extraction blanks and dust samples from the lab were also sequenced. These 16S rRNA gene sequences from contaminants were used to identify any contaminating bacteria in our oxic abyssal clay samples. All OTUs containing sequences from these 'contaminant' samples were removed prior to downstream analysis.

qPCR of 16S rRNA genes in DNA extraction blanks were consistently <10² copies per extraction, and thus we used 10² copies to define our detection limit for the abyssal clay samples. Consistent with this, high-throughput sequencing of amplicons with qPCR values <10² copies per g sediment had up to 50% sequence representation from contaminant taxa, whereas samples with values >10² copies per g sediment had <5% representation from contaminant taxa. This further supported our definition of <10² as a realistic detection limit. Using samples that had 16S rRNA gene copies >10² copies per g sediment, we were able to analyze microbial communities down to ca. 15 mbsf at Site 11 and ca. 8 mbsf at Site 12.

Thaumarchaeal amoA genes amplified via qPCR using the method described above were cloned using the TOPO TA cloning kit, and Sanger sequenced at the LMU Munich Sequencing Service at the Faculty of Biology (http://www.gi.bio.lmu.de/sequencing). Prior to phylogenetic analysis, the reads were quality trimmed in CLC Genomics using the default settings for quality control.

Experimental setups: DNA-SIP experiments with H₂¹⁸O were set up at two North Atlantic coring sites: Site 11 and Site 12, from depths 2.8 and 3.5 mbsf, respectively. Prior to setting up the incubations, the subcores were sampled with sterile syringes using the sample aseptic technique used for the DNA extraction. For each sample depth, seven grams of abyssal clay was placed into sterile 20-mL glass flasks and incubated with 4 mL of sterile artificial seawater composed of either H₂¹⁸O (97% atomic enrichment) or unlabeled artificial seawater. Vials were crimp sealed, with an oxygenated headspace of approximately 10 mL, and incubated at 8 C. The water content of the clay was measured to be approximately 40% (+/- 5%) of the total weight. This initial water content diluted the final concentration of added H₂¹⁸O to be ca. 60% of the total water content of the sample. The artificial seawater was different from the porewater at depth because there was no added nitrate, but there was also no added ammonia which should be similar to the in situ conditions where ammonia is generally below detection. Oxygen was measured continuously throughout the incubations using non-invasive fiberoptic measurements as described previously. Small fluctuations in the oxygen measurements in the killed control, and experimental incubations (Fig S3), were likely due to temperature fluctuations of the incubator itself (1°C), since the non-invasive fiber optic oxygen sensor spots are temperature sensitive.

To assess the preservation potential of extracellular DNA (eDNA), and its ability to bias our study that is based on DNA from living organisms, we monitored microbial growth in the presence and absence of added eDNA over a 210-day incubation experiment. eDNA extracted from a culture of Rhodococcus erythropolis was added to sediment slurries from 2.8 mbsf at Site 11 at a concentration of 5 ng g⁻¹. Microbial growth was measured over time with 16S rRNA gene qPCR, and also in a control that did not receive the eDNA. As a second control, we added eDNA to autoclaved (dead) sediment. DNA was extracted from each timepoint and measure with qPCR using the methods described above.

Density gradient fractionation, qPCR: We used Tag-SIP to measure the atom % ¹⁸O-enrichment of actively growing microbial taxa. In brief, after 7 and 18 months incubations DNA was extracted and subjected to Cesium Chloride (CsCl) density gradient centrifugation as described previously. The same 16S 515F/806R primers (described above) were used in qPCR (described above) to determine density shifts in the peak DNA of buoyant density (BD) for each incubation. 16S rRNA gene amplicons from each fractions resulting from the density gradient fractionation were Illumina sequenced as described above. To identify contaminants that may have entered during the fractionation process, we also included in the sequencing run extraction blanks from the SIP fractionation. OTUs containing sequences from extraction blanks were removed.

Metagenome library preparation, sequencing and bioinformatics analysis: Prior to library preparation, whole genome amplifications were performed on DNA extracts through a multiple displacement amplification (MDA) step of 6 to 7 hours, using the REPLI-g Midi Kit (QIAGEN) and following the manufacturer's instructions. We monitored amplification using SYBR green I (Invitrogen) on a CFX-Connect qPCR machine, stopping amplifications once the exponential phase was reached. Metagenomic libraries were prepared using the Nextera XT DNA Library Prep Kit (Illumina). Quality control and quantification of the libraries were obtained on an Agilent 2100 Bioanalyzer System, using the High Sensitivity DNA reagents and DNA chips (Agilent Genomics). was used to normalize library DNA concentrations. Metagenomic libraries were diluted to 1 nM using the Select-a-Size DNA Clean and Concentrator MagBead Kit (Zymo Research), and pooled for further sequencing on the Illumina MiniSeq platform.