Metadata und statistic analysis of archael and bacterial sequences originating from sediments of the Håkon Mosby mud volcano (Z1-Z3, all habitats)

DNA extraction was carried out as described on the MICROBIS project pages (http://icomm.mbl.edu/microbis ) using a commercially available extraction kit. We amplified the hypervariable regions V4-V6 of archaeal and bacterial 16S rRNA genes using PCR and several sets of forward and reverse primers (http://vamps.mbl.edu/resources/primers.php). Massively parallel tag sequencing of the PCR products was carried out on a 454 Life Sciences GS FLX sequencer at Marine Biological Laboratory, Woods Hole, MA, following the same experimental conditions for all samples. Sequence reads were submitted to a rigorous quality control procedure based on mothur v30 (doi:10.1128/AEM.01541-09) including denoising of the flow grams using an algorithm based on PyroNoise (doi:10.1038/nmeth.1361), removal of PCR errors and a chimera check using uchime (doi:10.1093/bioinformatics/btr381). The reads were taxonomically assigned according to the SILVA taxonomy (SSURef v119, 07-2014; doi:10.1093/nar/gks1219) implemented in mothur and clustered at 98% ribosomal RNA gene V4-V6 sequence identity. V4-V6 amplicon sequence abundance tables were standardized to account for unequal sampling effort using 1000 (Archaea) and 2300 (Bacteria) randomly chosen sequences without replacement using mothur and then used to calculate inverse Simpson diversity indices and Chao1 richness (doi:10.2307/4615964). Bray-Curtis dissimilarities (doi:10.2307/1942268) between all samples were calculated and used for 2-dimensional non metric multidimensional scaling (NMDS) ordinations with 20 random starts (doi:10.1007/BF02289694). Stress values below 0.2 indicated that the multidimensional dataset was well represented by the 2D ordination. NMDS ordinations were compared and tested using Procrustes correlation analysis (doi:10.1007/BF02291478). All analyses were carried out with the R statistical environment and the packages vegan (available at: http://cran.r-project.org/package=vegan), labdsv (available at: http://cran.r-project.org/package=labdsv), as well as with custom R scripts. Operational taxonomic units at 98% sequence identity (OTU0.03) that occurred only once in the whole dataset were termed absolute single sequence OTUs (SSOabs; doi:10.1038/ismej.2011.132). OTU0.03 sequences that occurred only once in at least one sample, but may occur more often in other samples were termed relative single sequence OTUs (SSOrel). SSOrel are particularly interesting for community ecology, since they comprise rare organisms that might become abundant when conditions change.16S rRNA amplicons and metagenomic reads have been stored in the sequence read archive under SRA project accession number SRP042162.

Microbial community, biomarker concentrations and their isotopic signature in sediment of Gullfaks and Tommeliten methane seeps in the Northern North Sea

Gullfaks is one of the four major Norwegian oil and gas fields, located in the northeastern edge of the North Sea Plateau. Tommeliten lies in the greater Ekofisk area in the central North Sea. During the cruises HE 208 and AL 267 several seep locations of the North Sea were visited. At the Heincke seep at Gullfaks, sediments were sampled in May 2004 (HE 208) using a video-guided multiple corer system (MUC; Octopus, Kiel). The samples were recovered from an area densely covered with bacterial mats where gas ebullition was observed. The coarse sands limited MUC penetration depth to maximal 30 centimeters and the highly permeable sands did not allow for a high-resolution, vertical subsampling because of pore water loss. The gas flare mapping and videographic observation at Tommeliten indicated an area of gas emission with a few small patches of bacterial mats with diameters <50 cm from most of which a single stream of gas bubbles emerged. The patches were spaced apart by 10-100 m. Sampling of sediments covered by bacterial mats was only possible with 3 small push cores (3.8 cm diameter) mounted to ROV Cherokee. These cores were sampled in 3 cm intervals. Lipid biomarker extraction from 10 -17 g wet sediment was carried out as described in detail elsewhere (Elvert et al., 2003; doi:10.1080/01490450303894). Briefly, defined concentrations of cholestane, nonadecanol and nonadecanolic acid with known delta 13C-values were added to the sediments prior to extraction as internal standards for the hydrocarbon, alcohol and fatty acid fraction, respectively. Total lipid extracts were obtained from the sediment by ultrasonification with organic solvents of decreasing polarity. Esterified fatty acids (FAs) were cleaved from the glycerol head group by saponification with methanolic KOH solution. From this mixture, the neutral fraction was extracted with hexane. After subsequent acidification, FAs were extracted with hexane. For analysis, FAs were methylated using BF3 in methanol yielding fatty acid methyl esters (FAMES). The fixation for total cell counts and CARD-FISH were performed on-board directly after sampling. For both methods, sediments were fixed in formaldehyde solution. After two hours, aliquots for CARD-FISH staining were washed with 1* PBS (10mmol/l sodium phosphate solution, 130mmol/l NaCl, adjusted to a pH of 7.2) and finally stored in a 1:1 PBS:ethanol solution at -20°C until further processing. Samples for total cell counts were stored in formalin at 4°C until analysis. For sandy samples, the total cell count/CARD-FISH protocol was optimized to separate sand particles from the cells. Cells were dislodged from sediment grains and brought into solution with the supernatant by sonicating each sample onice for 2 minutes at 50W. This procedure was repeated four times and supernatants were combined. The sediment samples were brought to a final dilution of 1:2000 to 1:4000 and filtered onto 0.2µm GTTP filters (Millipore, Eschbonn, Germany).