(Table S2) Relative abundance of archaeal 16S rRNA genes in sediment cores

Deep drilling into the marine sea floor has uncovered a vast sedimentary ecosystem of microbial cells (Parkes et al., 1994, doi:10.1038/371410a0; D'Hondt et al., 2004, doi:10.1126/science.1101155). Extrapolation of direct counts of stained microbial cells to the total volume of habitable marine subsurface sediments suggests that between 56 Pg (Parkes et al., 1994, doi:10.1038/371410a0) and 303 Pg (Whitman et al., 1998) of cellular carbon could be stored in this largely unexplored habitat. From recent studies using various culture-independent techniques, no clear picture has yet emerged as to whether Archaea or Bacteria are more abundant in this extensive ecosystem (Schippers et al., doi:10.1038/nature03302; Inagaki et al., doi:10.1073/pnas.0511033103 ; Mauclaire et al., doi:10.1111/j.1472-4677.2004.00035.x; Biddle et al., doi:10.1073/pnas.0600035103). Here we show that in subsurface sediments buried deeper than 1 m in a wide range of oceanographic settings at least 87% of intact polar membrane lipids, biomarkers for the presence of live cells (Biddle et al., doi:10.1073/pnas.0600035103; Sturt et al., 2004, doi:10.1002/rcm.1378), are attributable to archaeal membranes, suggesting that Archaea constitute a major fraction of the biomass. Results obtained from modified quantitative polymerase chain reaction and slot-blot hybridization protocols support the lipid-based evidence and indicate that these techniques have previously underestimated archaeal biomass. The lipid concentrations are proportional to those of total organic carbon. On the basis of this relationship, we derived an independent estimate of amounts of cellular carbon in the global marine subsurface biosphere. Our estimate of 90 Pg of cellular carbon is consistent, within an order of magnitude, with previous estimates, and underscores the importance of marine subsurface habitats for global biomass budgets.

(Table 1) Individual codes, haplotype codes, GenBank accession numbers for DNA sequences and sample locality of the analysed Betamorpha fusiformis specimens from the ANDEEP III expedition

Based on our current knowledge about population genetics, phylogeography and speciation, we begin to understand that the deep sea harbours more species than suggested in the past. Deep-sea soft-sediment environment in particular hosts a diverse and highly endemic invertebrate fauna. Very little is known about evolutionary processes that generate this remarkable species richness, the genetic variability and spatial distribution of deep-sea animals. In this study, phylogeographic patterns and the genetic variability among eight populations of the abundant and widespread deep-sea isopod morphospecies Betamorpha fusiformis [Barnard, K.H., 1920. Contributions to the crustacean fauna of South Africa. 6. Further additions to the list of marine isopods. Annals of the South African Museum 17, 319-438] were examined. A fragment of the mitochondrial 16S rRNA gene of 50 specimens and the complete nuclear 18S rRNA gene of 7 specimens were sequenced. The molecular data reveal high levels of genetic variability of both genes between populations, giving evidence for distinct monophyletic groups of haplotypes with average p-distances ranging from 0.0470 to 0.1440 (d-distances: 0.0592-0.2850) of the 16S rDNA, and 18S rDNA p-distances ranging between 0.0032 and 0.0174 (d-distances: 0.0033-0.0195). Intermediate values are absent. Our results show that widely distributed benthic deep-sea organisms of a homogeneous phenotype can be differentiated into genetically highly divergent populations. Sympatry of some genotypes indicates the existence of cryptic speciation. Flocks of closely related but genetically distinct species probably exist in other widespread benthic deep-sea asellotes and other Peracarida. Based on existing data we hypothesize that many widespread morphospecies are complexes of cryptic biological species (patchwork hypothesis).

Soil characteristics and abundance of archaea in samples from the Ross Sea Region, Antarctica

The objective of this study was to examine the presence and diversity of Archaea within mineral and ornithogenic soils from 12 locations across the Ross Sea region. Archaea were not abundant but DNA sufficient for producing 16S rRNA gene clone libraries was extracted from 18 of 51 soil samples, from four locations. A total of 1452 clones were analysed by restriction fragment length polymorphism and assigned to 43 operational taxonomic units from which representatives were sequenced. Archaea were primarily restricted to coastal mineral soils which showed a predominance of Crenarchaeota belonging to group 1.1b (>99% of clones). These clones were assigned to six clusters (A through F), based on shared identity to sequences in the GenBank database. Ordination indicated that soil chemistry and water content determined archaeal community structure. This is the first comprehensive study of the archaeal community in Antarctic soils and as such provides a reference point for further investigation of microbial function in this environment.

Gene expression and physiological changes of the Antarctic krill Euphausia superba under different hypoxia intensities

Antarctic krill (Euphausia superba) from South Georgia comprise one of the most northern and abundant krill stocks. South Georgia waters are undergoing rapid warming, as a result of climate change, which in turn could alter the oxygen concentration of the water. We investigated gene expression in Antarctic krill related to aerobic metabolism, antioxidant defence, and heat-shock response under severe (2.5% O2 saturation or 0.6 kPa) and threshold (20% O2 saturation or 4 kPa) hypoxia exposure compared to in situ levels (normoxic; 100% O2 saturation or 21 kPa). Biochemical metabolic and oxidative stress indicators complemented the genic expression analysis to detect in vivo signs of stress during the hypoxia treatments. Expression levels of the genes citrate synthase (CS), mitochondrial manganese superoxide dismutase (SODMn-m) and one heat-shock protein isoform (E) were higher in euphausiids incubated 6 h at 20% O2 saturation than in animals exposed to control (normoxic) conditions. All biochemical antioxidant defence parameters remained unchanged among treatments. Levels of lipid peroxidation were raised after 6 h of severe hypoxia. Overall, short-term exposure to hypoxia altered mitochondrial metabolic and antioxidant capacity, but did not induce anaerobic metabolism. Antarctic krill are swarming organisms and may experience short periods of hypoxia when present in dense swarms. A future, warmer Southern ocean, where oxygen saturation levels are decreased, may result in smaller, less dense swarms as they act to avoid greater levels of hypoxia.

Spatial and seasonal distribution of Thaumarchaeota in the water column and sediment of the southern North Sea

We have examined the spatial and seasonal distribution of Thaumarchaeota in the water column and sediment of the southern North Sea using the specific intact polar lipid (IPL) hexose, phosphohexose (HPH) crenarchaeol, as well as thaumarchaeotal 16S rRNA gene abundances and expression. In the water column, a higher abundance of Thaumarchaeota was observed in the winter season than in the summer, which is in agreement with previous studies, but this was not the case in the sediment where Thaumarchaeota were most abundant in spring and summer. This observation corresponds well with the idea that ammonia availability is a key factor in thaumarchaeotal niche determination. In the surface waters of the southern North Sea, we observed a spatial variability in HPH crenarchaeol, thaumarchaeotal 16S rRNA gene abundance and transcriptional activity that corresponded well with the different water masses present. In bottom waters, a clear differentiation based on water masses was not observed; instead, we suggest that observed differences in thaumarchaeotal abundance with depth may be related to resuspension from the sediment. This could be due to suspension of benthic Thaumarchaeota to the water column or due to delivery of e.g. resuspended sediment or ammonium to the water column, which could be utilized by pelagic Thaumarchaeota. This study has shown that the seasonality of Thaumarchaeota in water and sediment is different and highlights the importance of water masses, currents and sedimentary processes in determining the spatial abundance of Thaumarchaeota in the southern North Sea.

Molecular composition in surface and subsurface sediments of the Helgoland mud area, North Sea

The role of microorganisms in the cycling of sedimentary organic carbon is a crucial one. To better understand relationships between molecular composition of a potentially bioavailable fraction of organic matter and microbial populations, bacterial and archaeal communities were characterized using pyrosequencing-based 16S rRNA gene analysis in surface (top 30 cm) and subsurface/deeper sediments (30-530 cm) of the Helgoland mud area, North Sea. Fourier Transform Ion Cyclotron Resonance Mass Spectrometry (FT-ICR MS) was used to characterize a potentially bioavailable organic matter fraction (hot-water extractable organic matter, WE-OM). Algal polymer-associated microbial populations such as members of the Gammaproteobacteria, Bacteroidetes, and Verrucomicrobia were dominant in surface sediments while members of the Chloroflexi (Dehalococcoidales and candidate order GIF9) and Miscellaneous Crenarchaeota Groups (MCG), both of which are linked to degradation of more recalcitrant, aromatic compounds and detrital proteins, were dominant in subsurface sediments. Microbial populations dominant in subsurface sediments (Chloroflexi, members of MCG, and Thermoplasmata) showed strong correlations to total organic carbon (TOC) content. Changes of WE-OM with sediment depth reveal molecular transformations from oxygen-rich [high oxygen to carbon (O/C), low hydrogen to carbon (H/C) ratios] aromatic compounds and highly unsaturated compounds toward compounds with lower O/C and higher H/C ratios. The observed molecular changes were most pronounced in organic compounds containing only CHO atoms. Our data thus, highlights classes of sedimentary organic compounds that may serve as microbial energy sources in methanic marine subsurface environments.

Environmental characteristics and gas composition of Lost Hammer, Canadian Arctic

We report the first microbiological characterization of a terrestrial methane seep in a cryo-environment in the form of an Arctic hypersaline (~24% salinity), subzero (-5 C), perennial spring, arising through thick permafrost in an area with an average annual air temperature of -15 C. Bacterial and archaeal 16S rRNA gene clone libraries indicated a relatively low diversity of phylotypes within the spring sediment (Shannon index values of 1.65 and 1.39, respectively). Bacterial phylotypes were related to microorganisms such as Loktanella, Gillisia, Halomonas and Marinobacter spp. previously recovered from cold, saline habitats. A proportion of the bacterial phylotypes were cultured, including Marinobacter and Halomonas, with all isolates capable of growth at the in situ temperature (-5 C). Archaeal phylotypes were related to signatures from hypersaline deep-sea methane-seep sediments and were dominated by the anaerobic methane group 1a (ANME-1a) clade of anaerobic methane oxidizing archaea. CARD-FISH analyses indicated that cells within the spring sediment consisted of ~84.0% bacterial and 3.8% archaeal cells with ANME-1 cells accounting for most of the archaeal cells. The major gas discharging from the spring was methane (~50%) with the low CH4/C2 + ratio and hydrogen and carbon isotope signatures consistent with a thermogenic origin of the methane. Overall, this hypersaline, subzero environment supports a viable microbial community capable of activity at in situ temperature and where methane may behave as an energy and carbon source for sustaining anaerobic oxidation of methane-based microbial metabolism. This site also provides a model of how a methane seep can form in a cryo-environment as well as a mechanism for the hypothesized Martian methane plumes.

Geochemical data from pore water and solids from cores GS08-153-91GC12 and GS08-153-53-GC6, Arctic Mid-Ocean Ridge

Microbial communities and their associated metabolic activity in marine sediments have a profound impact on global biogeochemical cycles. Their composition and structure are attributed to geochemical and physical factors, but finding direct correlations has remained a challenge. Here we show a significant statistical relationship between variation in geochemical composition and prokaryotic community structure within deep-sea sediments. We obtained comprehensive geochemical data from two gravity cores near the hydrothermal vent field Loki's Castle at the Arctic Mid-Ocean Ridge, in the Norwegian-Greenland Sea. Geochemical properties in the rift valley sediments exhibited strong centimeter-scale stratigraphic variability. Microbial populations were profiled by pyrosequencing from 15 sediment horizons (59,364 16S rRNA gene tags), quantitatively assessed by qPCR, and phylogenetically analyzed. Although the same taxa were generally present in all samples, their relative abundances varied substantially among horizons and fluctuated between Bacteria- and Archaea-dominated communities. By independently summarizing covariance structures of the relative abundance data and geochemical data, using principal components analysis, we found a significant correlation between changes in geochemical composition and changes in community structure. Differences in organic carbon and mineralogy shaped the relative abundance of microbial taxa. We used correlations to build hypotheses about energy metabolisms, particularly of the Deep Sea Archaeal Group, specific Deltaproteobacteria, and sediment lineages of potentially anaerobic Marine Group I Archaea. We demonstrate that total prokaryotic community structure can be directly correlated to geochemistry within these sediments, thus enhancing our understanding of biogeochemical cycling and our ability to predict metabolisms of uncultured microbes in deep-sea sediments.

Ladderane lipid analysis and pore water and sediment chemistry of sediment cores from the continental shelf and slope off northwest Africa

The presence and abundance of anaerobic ammonium-oxidizing (anammox) bacteria was investigated in continental shelf and slope sediments (300-3000 m water depth) off northwest Africa in a combined approach applying quantitative polymerase chain reaction (q-PCR) analysis of anammox-specific 16S rRNA genes and anammox-specific ladderane biomarker lipids. We used the presence of an intact ladderane monoether lipid with a phosphocholine (PC) headgroup as a direct indicator for living anammox bacteria and compared it with the abundance of ladderane core lipids derived from both living and dead bacterial biomass. All investigated sediments contained ladderane lipids, both intact and core lipids, in agreement with the presence of anammoxspecific 16S rRNA gene copies, indicating that anammox occurs at all sites. Concentrations of ladderane core lipids in core top sediments varied between 0.3 and 97 ng g**-1 sediment, with the highest concentrations detected at the sites located on the shelf at shallower water depths between 300 and 500 m. In contrast, the C20 [3]-ladderane monoether-PC lipid was most abundant in a core top sediment from 1500 m water depth. Both anammox-specific 16S rRNA gene copy numbers and the concentration of the C20 [3]-ladderane monoether-PC lipid increased downcore in sediments located at greater water depths, showing highest concentrations of 1.2 x 10**8 copies g**-1 sediment and 30 pg g**-1 sediment, respectively, at the deepest station of 3000 m water depth. This suggests that the relative abundance of anammox bacteria is higher in sediments at intermediate to deep water depths where carbon mineralization rates are lower but where anammox is probably more important than denitrification.

Seawater carbonate chemistry and phytoplankton and eubacterial community compositions in the northwest subarctic Pacific

On-deck CO2-Fe-manipulated incubation experiments were conducted using surface seawater collected from the Western Subarctic Gyre of the NW Pacific in the summer of 2008 to elucidate the impacts of ocean acidification and Fe enrichment on the abundance and community composition of phytoplankton and eubacteria in the study area. During the incubation, excluding the initial period, the mean partial pressures of CO2 in non-Fe-added bottles were 230, 419, 843, and 1124 µatm, whereas those in Fe-added treatments were 152, 394, 791, and 1008 µatm. Changes in the abundance and community composition of phytoplankton were estimated using HPLC pigment signatures with the program CHEMTAX and flow cytometry. A DGGE fingerprint technique targeting 16S rRNA gene fragments was also used to estimate changes in eubacterial phylotypes during incubation. The Fe addition induced diatom blooms, and subsequently stimulated the growth of heterotrophic bacteria such as Roseobacter, Phaeobacter, and Alteromonas in the post-bloom phase. In both the Fe-limited and Fe-replete treatments, concentrations of 19'-hexanoyloxyfucoxanthin, a haptophyte marker, and the cell abundance of coccolithophores decreased at higher CO2 levels (750 and 1000 ppm), whereas diatoms exhibited little response to the changes in CO2 availability. The abundances of Synechococcus and small eukaryotic phytoplankton (<10 µm) increased at the higher CO2 levels. DGGE band positions revealed that Methylobacterium of Alphaproteobacteria occurred solely at lower CO2 levels (180 and 380 ppm) during the post-bloom phase. These results suggest that increases in CO2 level could affect not only the community composition of phytoplankton but also that of eubacteria. As these microorganisms play critical roles in the biological carbon pump and microbial loop, our results indicate that the progression of ocean acidification can alter the biogeochemical processes in the study area.

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.