Chemical characterization and quantification of the brown algal storage compound laminarin

In search of a meaningful stress indicator for Fucus vesiculosus we found that the often used quantitative determination procedures for the polysaccharide laminarin (beta-1,3-glucan) result in different kind of problems, uncertainties and limitations. This chemical long-term storage form of carbon enables perennial brown algae in seasonally fluctuating ecosystems to uncouple growth from photosynthesis. Because of this high ecological relevance a reliable and precise method for determination and quantification of laminarin is needed. Therefore, a simple, cold water extraction method coupled to a new quantitative liquid chromatography-mass spectrometrical method (LC-MS) was developed. Laminarin was determined in nine out of twelve brown algal species, and its expected typical molar mass distribution of 2000-7000 Da was confirmed. Furthermore, laminarin consisted of a complex mixture of different chemical forms, since fifteen chemical laminarin species with distinct molecular weights were measured in nine species of brown algae. Laminarin concentrations in the algal tissues ranged from 0.03 to 0.86% dry weight (DW). The direct chemical characterization and quantification of laminarin by LC-MS represents a powerful method to verify the biochemical and ecological importance of laminarin for brown algae. Single individuals of Laminaria hyperborea, L. digitata, Saccharina latissima, F. serratus, F. vesiculosus, F. spiralis, Himanthalia elongata, Cystoseira tamariscifolia, Pelvetia canaliculata, Ascophyllum nodosum, Halidrys siliquosa and Dictyota dichotoma were collected in fall (18.11.2013) during spring low tide from the shore of Finavarra, Co. Clare, west coast of Ireland (53° 09' 25'' N, 09° 06' 58'' W). After sampling, the different algae were immediately transported to the lab, lyophilized and sent to the University of Rostock. Laminarin was extracted with cold ultrapure water from the algal samples. Before extraction they were ground to < 1 mm grain size with an analytical mill (Ika MF 10 Basic). The algal material (approx. 1.5 g DW) was extracted in ultrapure water (8 mL) on a shaker (250 rpm) for 5 h. After the addition of surplus ultrapure water (4 mL) and shaking manually, 1 mL of the sample was filter centrifuged (45 µm) at 14,000 rpm (Hettich Mikro 22 R). The slightly viscous supernatant was free of suspended material and converted into a microvial (300 µL) for further analysis. The extracts were analyzed using liquid chromatography-mass spectrometry (LC-MS) analysis (LTQ Velos Pro ion trap spectrometer with Accela HPLC, Thermo Scientific). Laminarin species were separated on a KinetexTM column (2.6 µm C18, 150 x 3 mm). The mobile phase was 90 % ultrapure water and 10 % acetonitrile, run isocratically at a flow rate of 0.2 mL min-1. MS was working in ESI negative ion mode in a mass range of 100 - 4000 amu. Glucose contents were determined after extraction using high-performance liquid chromatography (HPLC). Extracted samples were analyzed in an HPLC (SmartLine, Knauer GmbH) equipped with a SUPELCOGELTM Ca column (30 x 7,8 mm without preColumn) and RI-detector (S2300 PDA S2800). Water was used as eluent at a flow rate of 0.8 mL min-1 at 75 °C. Glucose was quantified by comparison of the retention time and peak area with standard solutions using ChromGate software. Mannitol was extracted from three subsamples of 10-20 mg powdered alga material (L. hyperborea, L. digitata, S. latissima, F. serratus, F. vesiculosus, F. spiralis, H. elongata, P. canaliculata, A. nodosum, H. siliquosa) and quantified, following the HPLC method described by Karsten et al. (1991). For analyzing carbon and nitrogen contents, dried algal material was ground to powder and three subsamples of 2 mg from each alga thalli were loaded and packed into tin cartridges (6×6×12 mm). The packages were combusted at 950 °C and the absolute contents of C and N were automatically quantified in an elemental analyzer (Elementar Vario EL III, Germany) using acetanilide as standard according to Verardo et al. (1990).

Heterotrophic bacteria, Synechococcus, Prochlorococcus and Virus in the eastern Mediterranean Sea in March and April 2008 during SES_GR1

The dataset is based on samples collected in the framework of the project SESAME, in the Ionian, Libyan and Aegean Sea during March- April 2008. The objectives were to measure the standing stocks and calculate the production of the microbial compartment of the food web, describe the vertical distribution pattern and characterize its structure and function through the water column. Heterotrophic bacteria, Synechococcus, Prochlorococcus and Virus abundance: Subsamples for virus, heterotrophic bacteria and cyanobacteria (Synechococcus spp. and Prochlorococcus spp.) counting were analyzed using a FACSCalibur (Becton Dickinson) flow cytometer equipped with a standard laser (488 nm) and filter set and using deionized water as sheath fluid. Fluorescent beads with a diameter of 0.97 µm (Polysciences) were added to each sample as an internal standard, and all parameters were normalized to the beads and expressed as relative units. SYBRGreen I stain (Molecular Probe) was used to stain viral and heterotrophic bacterial DNA. Viruses were counted according to (Brussaard 1984). In order to avoid bulk consentrations of viruses samples we dilluted to Tris-EDTA (pH=8,0) buffer to a final sollution of 1/5 to 1/100. Total abundance and nucleid content classes were calculated using the Paint-A-Gate software (Becton Dickinson).

Pleistocene sedimentology of ODP Hole 133-820A from the northeastern Australian margin

A principal objective at Site 820, situated on the outer shelf, upper slope of the northeastern Australian continental margin, was to test the relationships between changes in Pleistocene sea level and sedimentary packages produced on a mixed carbonate-siliciclastic continental margin. To this end, we have examined the downcore distribution of grain size, magnetic susceptibility, and calcium-carbonate content throughout Hole 820A and, in particular, the top 35 meters below the seafloor (mbsf). These data are compared with variations in the oxygen-isotope signal defined for the same hole and are interpreted as indicating sea-level oscillations. The distribution of sand, mud, calcium carbonate of the mud fraction and total sample, and magnetic susceptibility during the last 20,000 yr defines the position of a sea-level regression (41,000-18,000 yr B.P.), a lowstand, early (18,000-9,400 yr B.P.) and late transgressions (9400-900 yr B.P.), and a highstand (4900 yr to the present). The regression is seen first in a high-carbonate content peak. Calcium carbonate constituents mainly comprise skeletal carbonate grains, with abundant planktonic and benthic foraminifers, and lime muds. The lowstand is characterized by a maximum abundance of the sand fraction, which contains dominantly skeletal carbonate grains and a minor abundance of lithoclasts. Sand-sized terrigenous sediments are proposed to have bypassed the continental shelf during a lowstand of sea level. Sedimentation rates throughout the regression and lowstand are low (3.0 cm/k.y.). The early transgression, marked by highest values in magnetic susceptibility, displays a rapid increase in sedimentation rate that coincided with an increase in terrigenous mud. Highest sedimentation rates of 82.3 cm/k.y. occurred during the late transgression, with increasing percentages of lime-mud. A decrease in noncarbonate constituents in the mud fraction during the late transgression and highstand of sea level is thought to be the result of restricted inner-shelf sedimentation of terrigenous sediments. The same relationship is also seen in the major sea-level oscillation, which is interpreted as isotope stage 6.

Geochemistry, clay mineralogy and grain size distribution in various sediment cores drilled during ODP Leg 119 on the Antarctic continental shelf off Prydz Bay

Drilling was undertaken at five sites (739-743) on ODP Leg 119 on a transect across the continental shelf of Prydz Bay, East Antarctica, to elucidate the long-term glacial history of the area and to examine the importance of the area with respect to the development of the East Antarctic ice sheet as a whole. In addition to providing a record of glaciation spanning 36 m.y. or more, Leg 119 has provided information concerning the development of a continental margin under the prolonged influence of a major ice sheet. This has allowed the development of a sedimentary model that may be applicable not only to other parts of the Antarctic continental margin, but also to northern high-latitude continental shelves. The cored glacial sedimentary record in Prydz Bay consists of three major sequences, dominated by diamictite: 1. An upper flat-lying sequence that ranges in thickness from a few meters in inner and western Prydz Bay to nearly 250 m in the outer or eastern parts of the bay. The uppermost few meters consist of Holocene diatom ooze and diatomaceous mud with a minor ice-rafted component overlying diamicton and diamictite of late Miocene to Quaternary age. The diamictite is mainly massive, but stratified varieties and minor mudstone and diatomite also occur. 2. An upper prograding sequence cored at Sites 739 and 743, unconformly below the flat-lying sequence. This consists of a relatively steep (4° inclination) prograding wedge with a number of discrete sedimentary packages. At Sites 739 and 743 the sequence is dominated by massive and stratified diamictite, some of which shows evidence of slumping and minor debris flowage. 3. A lower, more gently inclined, prograding sequence lies unconformably below the flat-lying sequence at Site 742 and the upper prograding sequence at Site 739. This extends to the base of both sites, to 316 and 487 mbsf, respectively. It is dominated by massive, relatively clast-poor diamictite which is kaolinite-rich, light in color, and contains sporadic carbonate-cemented layers. The lower part of Site 742 includes well-stratified diamictites and very poorly sorted mudstones. The base of this site has indications of large-scale soft-sediment deformation and probably represents proximity to the base of the glacial sequence. Facies analysis of the Prydz Bay glacial sequence indicates a range of depositional environments. Massive diamictite is interpreted largely as waterlain till, deposited close to the grounding line of a floating glacier margin, although basal till and debris flow facies are also present. Weakly stratified diamictite is interpreted as having formed close to or under the floating ice margin and influenced by the input of marine diatomaceous sediment (proximal glaciomarine setting). Well-stratified diamictite has a stronger marine input, being more diatom-rich, and probably represents a proximal-distal glaciomarine sediment with the glaciogenic component being supplied by icebergs. Other facies include a variety of mudstones and diatom-rich sediments of marine origin, in which an ice-rafted component is still significant. None of the recovered sediments are devoid of a glacial influence. The overall depositional setting of the prograding sequence is one in which the grounded ice margin is situated close to the shelf edge. Progradation was achieved primarily by deposition of waterlain till. The flat-lying sequence illustrates a complex sequence of advances and retreats across the outer part of the shelf, with intermittent phases of ice loading and erosion. The glacial chronology is based largely on diatom stratigraphy, which has limited resolution. It appears that ice reached the paleoshelf break by earliest Oligocene, suggesting full-scale development of the East Antarctic ice sheet by that time. The ice sheet probably dominated the continental margin for much of Oligocene to middle Miocene time. Retreat, but not total withdrawal of the ice sheet, took place in late Miocene to mid-Pliocene time. The late Pliocene to Pleistocene was characterized by further advances across, and progradation of, the continental shelf. Holocene time has been characterized by reduced glacial conditions and a limited influence of glacial processes on sedimentation.

Carbonate minerals, uranium and thorium concentration and U-Th ages of sediments from ODP Leg 166 sites

We have performed U-Th isotope analyses on pure aragonite samples from the upper sections of Leg 166 cores to assign each aragonite-rich sediment package to the correct sea-level highstand. The uppermost sediment package from each of the four sites investigated (Sites 1003, 1005, 1006, and 1007) yielded a Holocene U-Th age. Sediment packages from deeper in the cores have suffered diagenesis. This diagenesis consists of significant U loss (up to 40%) in the site nearest the platform (Site 1005), slight U gain in sites further from the platform, and continuous loss of pure 234U caused by alpha recoil at all sites. The difference in diagenesis between the sites can be explained by the different fluid-flow histories they have experienced. Site 1005 is sufficiently close to the platform to have probably experienced a change in flow direction whenever the banks have flooded or become exposed. Other sites have probably experienced continuous flow into the sediment. Although diagenesis prevents assignment of accurate ages, it is sufficiently systematic that it can be corrected for and each aragonite-rich package assigned to a unique highstand interval. Site 1005 has sediment packages from highstands associated with marine isotope Stages 1, 5, 7, 9, and 11. Site 1006 is similar, except that the Stage 7 highstand is missing, at least in Hole 1006A. Site 1003 has sediment only from Stage 1 and 11 highstands within the U-Th age range. And Site 1007 has sediment only from the stage 1 highstand. This information will allow the construction of better age models for these sites. No high-aragonite sediments are seen for Stage 3 or Substages 5a and 5c. Unless rather unusual erosion has occurred, this indicates that the banks did not flood during these periods. If true, this would require the sea level for Substages 5a and 5c to have remained at least ~10 m lower than today.

Bacterial abundance in the eastern Mediteranean Sea in August and September 2008 during SES_GR2

The SES_GR2_MICROBIAL PARAMETERS dataset is based on samples collected in the framework of the project SESAME, in the Ionian, Libyan and Aegean Sea during August-September 2008. The objectives were to measure the standing stocks and calculate the production of the microbial compartment of the food web, describe the vertical distribution pattern and characterize its structure and function through the water column. Subsamples for virus, heterotrophic bacteria and cyanobacteria (Synechococcus spp. and Prochlorococcus spp.) counting were analyzed using a FACSCalibur (Becton Dickinson) flow cytometer equipped with a standard laser (488 nm) and filter set and using deionized water as sheath fluid. Fluorescent beads with a diameter of 0.97 ?m (Polysciences) were added to each sample as an internal standard, and all parameters were normalized to the beads and expressed as relative units. SYBRGreen I stain (Molecular Probe) was used to stain viral and heterotrophic bacterial DNA. Viruses were counted according to (Brussaard 1984). In order to avoid bulk consentrations of viruses samples we dilluted to Tris-EDTA (pH=8,0) buffer to a final sollution of 1/5 to 1/100. Total abundance and nucleid content classes were calculated using the Paint-A-Gate software (Becton Dickinson).

Abundance and Biomass from heterotrophic bacteria, Synechococcus, Prochlorococcus and Virus in the eastern Mediteranean Sea in October 2008 during SES_GR2

The HCMR_SES_LAGRANGIAN_GR2_ MICROBIAL PARAMETERS dataset is based on samples collected in the framework of the project SESAME, in the North Aegean Sea during October 2008. The objectives were to measure the standing stocks and calculate the production of the microbial compartment of the food web, describe the vertical distribution pattern and characterize its structure and function through the water column as influenced by the BSW. Heterotrophic bacteria, Synechococcus, Prochlorococcus and Virus abundance: Subsamples for virus, heterotrophic bacteria and cyanobacteria (Synechococcus spp. and Prochlorococcus spp.) counting were analyzed using a FACSCalibur (Becton Dickinson) flow cytometer equipped with a standard laser (488 nm) and filter set and using deionized water as sheath fluid. Fluorescent beads with a diameter of 0.97 µm (Polysciences) were added to each sample as an internal standard, and all parameters were normalized to the beads and expressed as relative units. SYBRGreen I stain (Molecular Probe) was used to stain viral and heterotrophic bacterial DNA. Viruses were counted according to (Brussaard 1984). In order to avoid bulk consentrations of viruses samples we dilluted to Tris-EDTA (pH=8,0) buffer to a final sollution of 1/5 to 1/100. Total abundance and nucleid content classes were calculated using the Paint-A-Gate software (Becton Dickinson). Heterotrophic Nanoflagellate abundance: Subsamples (30-150 ml) were concentrated on 25mm black polycarbonate filters of porosity 0.6?m and stained with DAPI for 10 min (Porter and Feig 1980). Under epifluorescence microscopy heterotrophic nanoflagellates (HNAN) were distinguished using UV and blue excitation and enumerated. Nanoflagellates were classified in size categories and the biovolume was calculated. Ciliate abundance: For ciliate identification and enumeration, 100-3000 ml samples were left for 24h-4d for sedimentation and then observed under an inverted microscope. Ciliates were counted, distinguished into size-classes and major taxonomic groups and identified down to genus or species level where possible (Pitta et al. 2005). Heterotrophic bacteria, Synechococcus, Prochlorococcus bacteria: Subsamples for virus, heterotrophic bacteria and cyanobacteria (Synechococcus spp. and Prochlorococcus spp.) counting were analyzed using a FACSCalibur (Becton Dickinson) flow cytometer equipped with a standard laser (488 nm) and filter set and using deionized water as sheath fluid. Fluorescent beads with a diameter of 0.97 µm (Polysciences) were added to each sample as an internal standard, and all parameters were normalized to the beads and expressed as relative units. SYBRGreen I stain (Molecular Probe) was used to stain viral and heterotrophic bacterial DNA. Viruses were counted according to (Brussaard 1984). In order to avoid bulk consentrations of viruses samples we dilluted to Tris-EDTA (pH=8,0) buffer to a final sollution of 1/5 to 1/100. Total abundance and nucleid content classes were calculated using the Paint-A-Gate software (Becton Dickinson). Abundance data were converted into C biomass using 250 fgC cell-1 (Kana & Glibert 1987) for Synechococcus, 50 fgC cell-1 (Campbell et al. 1994) for Prochlorococcus and 20fgC cell-1 (Lee & Fuhrman 1987) for heterotrophic bacteria. Heterotrophic Nanoflagellate biomass: Subsamples (30-150 ml) were concentrated on 25mm black polycarbonate filters of porosity 0.6µm and stained with DAPI for 10 min (Porter and Feig 1980). Under epifluorescence microscopy heterotrophic nanoflagellates (HNAN) were distinguished using UV and blue excitation and enumerated. Nanoflagellates were classified in size categories and the biovolume was calculated. Abundance data were converted into C biomass using 183 fgC µm**3 (Caron et al. 1995). Ciliate biomass: For ciliate identification and enumeration, 100-3000 ml samples were left for 24h-4d for sedimentation and then observed under an inverted microscope. Ciliates were counted, distinguished into size-classes and major taxonomic groups and identified down to genus or species level where possible (Pitta et al. 2005). Ciliate cell sizes were measured and converted into cell volumes using appropriate geometric formulae using image analysis. For biomass estimation, the conversion factor 190 fgC µm**3 was used (Putt and Stoecker 1989).

PaleoTools - extended software package for quantitative paleoenvironmental reconstructions

The Imbrie and Kipp transfer function method (IKM) and the modern analog technique (MAT) are accepted tools for quantitative paleoenvironmental reconstructions. However, no uncomplicated, flexible software has been available to apply these methods on modern computer devices. For this reason the software packages PaleoToolBox, MacTransfer, WinTransfer, MacMAT, and PanPlot have been developed. The PaleoToolBox package provides a flexible tool for the preprocessing of microfossil reference and downcore data as well as hydrographic reference parameters. It includes procedures to randomize the raw databases; to switch specific species in or out of the total species list; to establish individual ranking systems and their application on the reference and downcore databasessemi; and to convert the prepared databases into the file formats of IKM and MAT software for estimation of paleohydrographic parameters.

Abundance and Biomass from heterotrophic bacteria, Synechococcus, Prochlorococcus and Virus in the Aegean Sea in April 2008 during SES_GR1

The HCMR_SES_LAGRANGIAN_GR1_ MICROBIAL PARAMETERS dataset is based on samples collected in the framework of the project SESAME, in the North Aegean Sea during April 2008. The objectives were to measure the standing stocks and calculate the production of the microbial compartment of the food web, describe the vertical distribution pattern and characterize its structure and function through the water column as influenced by the BSW. Heterotrophic bacteria, Synechococcus, Prochlorococcus and Virus abundance: Subsamples for virus, heterotrophic bacteria and cyanobacteria (Synechococcus spp. and Prochlorococcus spp.) counting were analyzed using a FACSCalibur (Becton Dickinson) flow cytometer equipped with a standard laser (488 nm) and filter set and using deionized water as sheath fluid. Fluorescent beads with a diameter of 0.97 µm (Polysciences) were added to each sample as an internal standard, and all parameters were normalized to the beads and expressed as relative units. SYBRGreen I stain (Molecular Probe) was used to stain viral and heterotrophic bacterial DNA. Viruses were counted according to (Brussaard 1984). In order to avoid bulk consentrations of viruses samples we dilluted to Tris-EDTA (pH=8,0) buffer to a final sollution of 1/5 to 1/100. Total abundance and nucleid content classes were calculated using the Paint-A-Gate software (Becton Dickinson). Heterotrophic Nanoflagellate abundance: Subsamples (30-150 ml) were concentrated on 25mm black polycarbonate filters of porosity 0.6µm and stained with DAPI for 10 min (Porter and Feig 1980). Under epifluorescence microscopy heterotrophic nanoflagellates (HNAN) were distinguished using UV and blue excitation and enumerated. Nanoflagellates were classified in size categories and the biovolume was calculated. Ciliate abundance: For ciliate identification and enumeration, 100-3000 ml samples were left for 24h-4d for sedimentation and then observed under an inverted microscope. Ciliates were counted, distinguished into size-classes and major taxonomic groups and identified down to genus or species level where possible (Pitta et al. 2005). Heterotrophic bacteria, Synechococcus, Prochlorococcus biomass: Subsamples for virus, heterotrophic bacteria and cyanobacteria (Synechococcus spp. and Prochlorococcus spp.) counting were analyzed using a FACSCalibur (Becton Dickinson) flow cytometer equipped with a standard laser (488 nm) and filter set and using deionized water as sheath fluid. Fluorescent beads with a diameter of 0.97 µm (Polysciences) were added to each sample as an internal standard, and all parameters were normalized to the beads and expressed as relative units. SYBRGreen I stain (Molecular Probe) was used to stain viral and heterotrophic bacterial DNA. Viruses were counted according to (Brussaard 1984). In order to avoid bulk consentrations of viruses samples we dilluted to Tris-EDTA (pH=8,0) buffer to a final sollution of 1/5 to 1/100. Total abundance and nucleid content classes were calculated using the Paint-A-Gate software (Becton Dickinson). Abundance data were converted into C biomass using 250 fgC cell-1 (Kana & Glibert 1987) for Synechococcus, 50 fgC cell-1 (Campbell et al. 1994) for Prochlorococcus and 20fgC cell-1 (Lee & Fuhrman 1987) for heterotrophic bacteria. Heterotrophic Nanoflagellate biomass: Subsamples (30-150 ml) were concentrated on 25mm black polycarbonate filters of porosity 0.6µm and stained with DAPI for 10 min (Porter and Feig 1980). Under epifluorescence microscopy heterotrophic nanoflagellates (HNAN) were distinguished using UV and blue excitation and enumerated. Nanoflagellates were classified in size categories and the biovolume was calculated. Abundance data were converted into C biomass using 183 fgC µm**3 (Caron et al. 1995). Ciliate biomass: For ciliate identification and enumeration, 100-3000 ml samples were left for 24h-4d for sedimentation and then observed under an inverted microscope. Ciliates were counted, distinguished into size-classes and major taxonomic groups and identified down to genus or species level where possible (Pitta et al. 2005). Ciliate cell sizes were measured and converted into cell volumes using appropriate geometric formulae using image analysis. For biomass estimation, the conversion factor 190 fgC µm**3 was used (Putt and Stoecker 1989).

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.