Samples were collected with SCUBA divers at Cape Lookout Bight, North Carolina (34.6205°N, 76.5500°W), October 2, 2013. Thirty 20 cm PVC push cores were collected, capped, refrigerated, and then returned to the lab on ice in Tennessee within 48 hrs. Bubbles of methane were released from the sediments as each core was taken, indicating the presence of methane at the surface of the core. Using a plunger inserted from the bottom the first three centimeters of sediment taken from each core was placed in a 2L Erlenmeyer flask by way of a funnel. About ten core tubes were needed to fill each of the three flasks to 1.5 liters of sediment. About 100 ml of slurried sediment were autoclaved and incubated alongside the experiments under anoxic conditions as a negative control. Cell counts, hydrogen, sulfate, and methane were measured for the negative control on day 260. All measurements described below were performed for 18 weekly timepoints for all three incubations (0-122 days). An additional timepoint at 802 days was taken for 16S rRNA gene sequence analysis only.
Geochemistry
Each of the three flasks was fitted with a custom butyl rubber stopper with a hole drilled through the center to accommodate a wide bore (6 mm) glass and Teflon stopcock for the removal of samples. Two 18-gauge needles with stainless steel stopcocks were inserted into the stopper as well. Using the luer-lock fitting on the needles, ultra high purity nitrogen gas (99.999%) that had been scrubbed of oxygen using heated copper fillings was flowed through the bottles using the second needle for the outflow to make the headspace anoxic. Incubation 3 was then flushed with ≥ 99.0% methane gas (Sigma-Aldrich, St. Louis, MO). Then all the ports were closed and the flasks inverted so that sediment covered the gas ports, stopper, and stopcock, and placed in a ring stand at constant room temperature (21.4°C) in the dark.
The incubations were turned over once every seven days just before sampling. Prior to gas sampling, 2 ml of anoxic N2 gas (99.999%) was used to blow the needle clear of sediment. Hydrogen and methane gas samples were collected in glass gastight Hamilton syringes using the steel needle ports in the custom stopper. About 32 ml of sediment was removed through the glass and Teflon stopcock using a sterile 60 ml plastic catheter tip syringe. From this, two 15ml conical centrifuge tubes were filled and capped, one used for porewater analysis and the other frozen at -80°C for later molecular analysis. One ml of sediment was placed in a 2 ml screw cap tube to be fixed and used for cell counts as described below. After sampling, 30ml of oxygen- and hydrogen-scrubbed N2 was injected into the bottle to replace the lost volume. The 15 ml tube destined for porewater analysis was centrifuged at 5000 xG for 5 minutes. A syringe was used to remove the supernatant not in contact with the air. The porewater was then filtered using a 0.2 μm syringe filter into 100 μl of 10% HCl to a final volume of 1 ml. Porewater sulfate content was determined by ion chromatography (Dionex, Sunnyvale, CA), with the remaining porewater used for determining the pH.
500µl of headspace gas was injected into a Peak Performer 1 Reducing Compound Photometer (Peak Laboratories, Mountain View, CA). Premixed hydrogen ppm lab bottles (Sigma-Aldrich, St. Louis, MO) were used as standards. Hydrogen was assumed to be equilibrated between headspace and porewater, since the equilibration time for Cape Lookout Bight sediments is < 2 days (Hoehler et al. 1998), which is less than the 7 days between timepoints in our study. Therefore, gas phase partial pressures were converted to aqueous hydrogen concentrations using the solubility coefficient of hydrogen corrected for salinity of 35 ppt and temperature of 21.4°C (Crozier & Yamamoto 1980). Methane was determined by using injected 500 µl of gas from the headspace into an evacuated glass bottle to be later analyzed on a gas chromatograph with a flame ionized detector (Agilent, Santa Clara, CA). Methane concentrations were not assumed to be equilibrated with the aqueous phase, therefore concentrations are presented as partial pressures.
Cell quantification
Total cell counts were determined by direct epifluorescence microscopy SYBRGold DNA stain (Invitrogen, Carlsbad, CA). Sediments were sonicated at 20% power for 40 seconds to disaggregates cells from sediments and diluted 40-fold into PBS prior to filtration onto a 0.2 μm polycarbonate filter (Fisher Scientific, Waltham, MA) and mounted onto a slide. An autoclaved sediment sample was used as a negative control.
Quantitative PCR
DNA was extracted from sediment samples frozen at -80°C using the Fast DNA kit for Soil (MP Bio, Santa Ana, CA). Negative controls of sterilized sediment and a blank water extract were used as well. Quantitative PCR was used to determine 16S rRNA gene copy numbers of several taxa with Quantifast SYBRGreen kit (Qiagen) on a BioRad IQ5 machine. DNA standards were prepared from either existing stocks (Lloyd et al., 2011) or from TOPO plasmids (Invitrogen, Carlsbad, CA) containing PCR amplified 16S rRNA gene products of closely related relatives in the clades Methanomicrobiales and Methanosarcinales synthesized by Invitrogen (accession # AB236118 and AB679168 respectively). Standards were quantified using Hoechst dye in a flourimeter (Hoefer, Holliston, MA). 16s rRNA specific primers Methanomicrobiales and Methanosarcinales primers were selected to have good coverage of the taxa (Narihiro and Sekiguchi, 2011).
16S Ribosomal RNA Gene Amplicons
Extracted DNA was used for 16S rRNA gene amplicon analysis. The V4 region of each DNA extraction was amplified using primers 806r and 515f (Caporaso et al., 2012), as a universal primer pair for Bacteria and Archaea. Library preparations via Nexterra kit and sequencing using an Illumina MiSeq were performed at the Center for Environmental Biotechnology at the University of Tennessee in Knoxville. The Mothur MiSeq Standard Operating Procedure was used to make contigs of bidirectional sequences, cluster operational taxonomic units (OTUs) at 97% similarity, and classify them with the Silva reference set 119 (Schloss et al., 2009, Pruesse et al., 2007). 26.4% of unique sequences were removed as chimeric and then approximately 5% of total sequences were removed for failing to classify at the domain level. Reads were normalized against the sum of reads classified as bacteria and archaea.
Analyses of time course patterns of various microbial taxa were considered only for those with more than 20 reads when summed from the 18 timepoints from each of the three incubations, leaving between 593 and 669 genus-level clades of bacteria and archaea across the three incubations. Total reads ranged from 20,922 to 329,380 for the 54 libraries. 16S rRNA sequences were normalized to total classifiable bacterial and archaeal reads. ]