Environmental conditions and habitats of GOS sites and bottle data of HOT-ALOHA station

The Global Ocean Sampling (GOS) expedition is currently the largest and geographically most comprehensive metagenomic dataset, including samples from the Atlantic, Pacific, and Indian Oceans. This study makes use of the wide range of environmental conditions and habitats encompassed within the GOS sites in order to investigate the ecological structuring of bacterial and archaeal taxon ranks. Community structures based on taxonomically classified 16S ribosomal RNA (rRNA) gene fragments at phylum, class, order, family, and genus rank levels were examined using multivariate statistical analysis, and the results were inspected in the context of oceanographic environmental variables and structured habitat classifications. At all taxon rank levels, community structures of neritic, oceanic, estuarine biomes, as well as other exotic biomes (salt marsh, lake, mangrove), were readily distinguishable from each other. A strong structuring of the communities with chlorophyll a concentration and a weaker yet significant structuring with temperature and salinity were observed. Furthermore, there were significant correlations between community structures and habitat classification. These results were used for further investigation of one-to-one relationships between taxa and environment and provided indications for ecological preferences shaped by primary production for both cultured and uncultured bacterial and archaeal clades.

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.

EPOCA Svalbard 2009 benthic experiment: Serripes study, 2009

Ocean acidification influences sediment/water nitrogen fluxes, possibly by impacting on the microbial process of ammonia oxidation. To investigate this further, undisturbed sediment cores collected from Ny Alesund harbour (Svalbard) were incubated with seawater adjusted to CO2 concentrations of 380, 540, 760, 1,120 and 3,000 µatm. DNA and RNA were extracted from the sediment surface after 14 days' exposure and the abundance of bacterial and archaeal ammonia oxidising (amoA) genes and transcripts quantified using quantitative polymerase chain reaction. While there was no change to the abundance of bacterial amoA genes, an increase to 760 µatm pCO2 reduced the abundance of bacterial amoA transcripts by 65 %, and this was accompanied by a shift in the composition of the active community. In contrast, archaeal amoA gene and transcript abundance both doubled at 3,000 µatm, with an increase in species richness also apparent. This suggests that ammonia oxidising bacteria and archaea in marine sediments have different pH optima, and the impact of elevated CO2 on N cycling may be dependent on the relative abundances of these two major microbial groups. Further evidence of a shift in the balance of key N cycling groups was also evident: the abundance of nirS-type denitrifier transcripts decreased alongside bacterial amoA transcripts, indicating that NO3 ? produced by bacterial nitrification fuelled denitrification. An increase in the abundance of Planctomycete-specific 16S rRNA, the vast majority of which grouped with known anammox bacteria, was also apparent at 3,000 µatm pCO2. This could indicate a possible shift from coupled nitrification-denitrification to anammox activity at elevated CO2.