Gene expression profiling in gills of the great spider crab Hyas araneus in respond to ocean acidification and warming, supplementary data

Hypercapnia and elevated temperatures resulting from climate change may have adverse consequences for many marine organisms. While diverse physiological and ecological effects have been identified, changes in those molecular mechanisms, which shape the physiological phenotype of a species and limit its capacity to compensate, remain poorly understood. Here, we use global gene expression profiling through RNA-Sequencing to study the transcriptional responses to ocean acidification and warming in gills of the boreal spider crab Hyas araneus exposed medium-term (10 weeks) to intermediate (1,120 µatm) and high (1,960 µatm) PCO2 at different temperatures (5°C and 10°C). The analyses reveal shifts in steady state gene expression from control to intermediate and from intermediate to high CO2 exposures. At 5°C acid-base, energy metabolism and stress response related genes were upregulated at intermediate PCO2, whereas high PCO2 induced a relative reduction in expression to levels closer to controls. A similar pattern was found at elevated temperature (10°C). There was a strong coordination between acid-base, metabolic and stress-related processes. Hemolymph parameters at intermediate PCO2 indicate enhanced capacity in acid-base compensation potentially supported by upregulation of a V-ATPase. The likely enhanced energy demand might be met by the upregulation of the electron transport system (ETS), but may lead to increased oxidative stress reflected in upregulated antioxidant defense transcripts. These mechanisms were attenuated by high PCO2, possibly as a result of limited acid-base compensation and metabolic down-regulation. Our findings indicate a PCO2 dependent threshold beyond which compensation by acclimation fails progressively. They also indicate a limited ability of this stenoecious crustacean to compensate for the effects of ocean acidification with and without concomitant warming.

Results of H2-consumption experiments during METEOR cruise M64/2

The discovery of deep-sea hydrothermal vents in 1977 revolutionized our understanding of the energy sources that fuel primary productivity on Earth. Hydrothermal vent ecosystems are dominated by animals that live in symbiosis with chemosynthetic bacteria. So far, only two energy sources have been shown to power chemosynthetic symbioses: reduced sulphur compounds and methane. Using metagenome sequencing, single-gene fluorescence in situ hybridization, immunohistochemistry, shipboard incubations and in situ mass spectrometry, we show here that the symbionts of the hydrothermal vent mussel Bathymodiolus from the Mid-Atlantic Ridge use hydrogen to power primary production. In addition, we show that the symbionts of Bathymodiolus mussels from Pacific vents have hupL, the key gene for hydrogen oxidation. Furthermore, the symbionts of other vent animals such as the tubeworm Riftia pachyptila and the shrimp Rimicaris exoculata also have hupL. We propose that the ability to use hydrogen as an energy source is widespread in hydrothermal vent symbioses, particularly at sites where hydrogen is abundant.

Hydrogen consumption in vent files of the equatorial Atlantic

The discovery of deep-sea hydrothermal vents in 1977 revolutionized our understanding of the energy sources that fuel primary productivity on Earth. Hydrothermal vent ecosystems are dominated by animals that live in symbiosis with chemosynthetic bacteria. So far, only two energy sources have been shown to power chemosynthetic symbioses: reduced sulphur compounds and methane. Using metagenome sequencing, single-gene fluorescence in situ hybridization, immunohistochemistry, shipboard incubations and in situ mass spectrometry, we show here that the symbionts of the hydrothermal vent mussel Bathymodiolus from the Mid-Atlantic Ridge use hydrogen to power primary production. In addition, we show that the symbionts of Bathymodiolus mussels from Pacific vents have hupL, the key gene for hydrogen oxidation. Furthermore, the symbionts of other vent animals such as the tubeworm Riftia pachyptila and the shrimp Rimicaris exoculata also have hupL. We propose that the ability to use hydrogen as an energy source is widespread in hydrothermal vent symbioses, particularly at sites where hydrogen is abundant.

Svalbard 2010 mesocosm experiment: Protistan diversity

Hierarchical clustering. Taxonomic assignment of reads was performed using a preexisting database of SSU rDNA sequences from including XXX reference sequences generated by Sanger sequencing. Experimental amplicons (reads), sorted by abundance, were then concatenated with the reference extracted sequences sorted by decreasing length. All sequences, experimental and referential, were then clustered to 85% identity using the global alignment clustering option of the uclust module from the usearch v4.0 software (Edgar, 2010). Each 85% cluster was then reclustered at a higher stringency level (86%) and so on (87%, 88%,.) in a hierarchical manner up to 100% similarity. Each experimental sequence was then identified by the list of clusters to which it belonged at 85% to 100% levels. This information can be viewed as a matrix with the lines corresponding to different sequences and the columns corresponding to the cluster membership at each clustering level. Taxonomic assignment for a given read was performed by first looking if reference sequences clustered with the experimental sequence at the 100% clustering level. If this was the case, the last common taxonomic name of the reference sequence(s) within the cluster was used to assign the environmental read. If not, the same procedure was applied to clusters from 99% to 85% similarity if necessary, until a cluster was found containing both the experimental read and reference sequence(s), in which case sequences were taxonomically assigned as described above.

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.

Dissolved organic matter composition derived from FT-ICR-MS analyses during HEINCKE cruise HE361

Dissolved organic matter (DOM) is the main substrate and energy source for heterotrophic bacterioplankton. To understand the interactions between DOM and the bacterial community (BC), it is important to identify the key factors on both sides in detail, chemically distinct moieties in DOM and the various bacterial taxa. Next-generation sequencing facilitates the classification of millions of reads of environmental DNA and RNA amplicons and ultrahigh-resolution mass spectrometry yields up to 10,000 DOM molecular formulae in a marine water sample. Linking this detailed biological and chemical information is a crucial first step toward a mechanistic understanding of the role of microorganisms in the marine carbon cycle. In this study, we interpreted the complex microbiological and molecular information via a novel combination of multivariate statistics. We were able to reveal distinct relationships between the key factors of organic matter cycling along a latitudinal transect across the North Sea. Total BC and DOM composition were mainly driven by mixing of distinct water masses and presumably retain their respective terrigenous imprint on similar timescales on their way through the North Sea. The active microbial community, however, was rather influenced by local events and correlated with specific DOM molecular formulae indicative of compounds that are easily degradable. These trends were most pronounced on the highest resolved level, that is, operationally defined 'species', reflecting the functional diversity of microorganisms at high taxonomic resolution.

Hydrogen consumption by gills of Bathymodiolus spp. collected at the Comfortless cove and Lilliput hydrothermal vent fields in May 2006

The discovery of deep-sea hydrothermal vents in 1977 revolutionized our understanding of the energy sources that fuel primary productivity on Earth. Hydrothermal vent ecosystems are dominated by animals that live in symbiosis with chemosynthetic bacteria. So far, only two energy sources have been shown to power chemosynthetic symbioses: reduced sulphur compounds and methane. Using metagenome sequencing, single-gene fluorescence in situ hybridization, immunohistochemistry, shipboard incubations and in situ mass spectrometry, we show here that the symbionts of the hydrothermal vent mussel Bathymodiolus from the Mid-Atlantic Ridge use hydrogen to power primary production. In addition, we show that the symbionts of Bathymodiolus mussels from Pacific vents have hupL, the key gene for hydrogen oxidation. Furthermore, the symbionts of other vent animals such as the tubeworm Riftia pachyptila and the shrimp Rimicaris exoculata also have hupL. We propose that the ability to use hydrogen as an energy source is widespread in hydrothermal vent symbioses, particularly at sites where hydrogen is abundant.

Photosynthetically Active Radiation (PAR) measurements from a Biospherical Instrument Inc. QCR-2150 surface PAR sensor mounted on a sensor mast during the Tara Oceans expedition 2009-2013

The Tara Oceans Expedition (2009-2013) sampled the world oceans on board a 36 m long schooner, collecting environmental data and organisms from viruses to planktonic metazoans for later analyses using modern sequencing and state-of-the-art imaging technologies. Tara Oceans Data are particularly suited to study the genetic, morphological and functional diversity of plankton. The present data set provides continuous measurements made with a Biospherical Instrument Inc. QCR-2150 surface PAR sensor mounted on a sensor mast at the stern of the ship (ca. 8m above deck) and time synchronized with the CTD recording unit. The sensor consists of a cosine collector and was also utilized to correct the CTD PAR sensor data. The dark was computed as the lowest 0.01% voltage of the signal that was found to be very stable (0.00965V) for all the legs except for the 2nd leg of the polar circle where there was no complete night (the manufacturer dark was 0.0097V). The manufacturer calibration slope from 12/ 2012 was used to transform the data to scientific units.

Biodiversity context of all samples from the Tara Oceans Expedition (2009-2013)

The Tara Oceans Expedition (2009-2013) sampled the world oceans on board a 36 m long schooner, collecting environmental data and organisms from viruses to planktonic metazoans for later analyses using modern sequencing and state-of-the-art imaging technologies. Tara Oceans Data are particularly suited to study the genetic, morphological and functional diversity of plankton. The present data set provides biodiversity context to all samples from the Tara Oceans Expedition (2009-2013), including various diversity indexes calculated for the sampling location using satellite and model climatologies (Darwin project, Physat) and results from the sequencing of Tara Oceans samples.

Properties of seawater and particulate matter (fluorescence) from a WETLabs Eco-FL sensor mounted on the continuous surface water sampling system during the Tara Oceans expedition 2009-2013

The Tara Oceans Expedition (2009-2013) sampled the world oceans on board a 36 m long schooner, collecting environmental data and organisms from viruses to planktonic metazoans for later analyses using modern sequencing and state-of-the-art imaging technologies. Tara Oceans Data are particularly suited to study the genetic, morphological and functional diversity of plankton. The present data set provides continuous measurements made with a WETLabs Eco-FL sensor mounted on the flowthrough system between June 4th, 2011 and March 30th, 2012. Data was recorded approximately every 10s. Two issues affected the data: 1. Periods when the water 0.2µm filtered water were used as blanks and 2. Periods where fluorescence was affected by non-photochemical quenching (NPQ, chlorophyll fluorescence is reduced when cells are exposed to light, e.g. Falkowski and Raven, 1997). Median data and their standard deviation were binned to 5min bins with period of light/dark indicated by an added variable (so that NPQ affected data could be neglected if the user so chooses). Data was first calibrated using HPLC data collected on the Tara (there were 36 data within 30min of each other). Fewer were available when there was no evident NPQ and the resulting scale factor was 0.0106 mg Chl m-3/count. To increase the calibration match-ups we used the AC-S data which provided a robust estimate of Chlorophyll (e.g. Boss et al., 2013). Scale factor computed over a much larger range of values than HPLC was 0.0088 mg Chl m-3/count (compared to 0.0079 mg Chl m-3/count based on manufacturer). In the archived data the fluorometer data is merged with the TSG, raw data is provided as well as manufacturer calibration constants, blank computed from filtered measurements and chlorophyll calibrated using the AC-S. For a full description of the processing of the Eco-FL please see Taillandier, 2015.

Properties of seawater (pH total scale) from a Satlantic SeaFET instrument mounted on the continuous surface water sampling system during the Tara Oceans expedition 2009-2013

The Tara Oceans Expedition (2009-2013) sampled the world oceans on board a 36 m long schooner, collecting environmental data and organisms from viruses to planktonic metazoans for later analyses using modern sequencing and state-of-the-art imaging technologies. Tara Oceans Data are particularly suited to study the genetic, morphological and functional diversity of plankton. The present data set provides continuous pH measurements made during 2013 expedition with a Satlantic SeaFET instrument that was connected to the flowthrough system. Data calibration was performed according to Bresnahan et al. (2014) (using spectrophotometric pH measurements on discrete samples (Clayton and Byrne 1993). pH_internal values were taken to calibrate the data (rather than pH_external) because of the better calibration coefficient (there was no trend associated with it). The equations of Clayton and Byrne (1993) was used to compute pH from the measured absorbance values at the temperature of measurement. The data was converted to in situ temperature using the "CO2-sys" program which can be downloaded from http://cdiac.ornl.gov/ftp/co2sys/.

Sedimentological investigations and age constraints of sediment cores from the outer Abbot trough, Amundsen Sea shelf

Knowing the extent of the West Antarctic Ice Sheet (WAIS) at the Last Glacial Maximum (LGM) is crucial for initiating and calibrating numerical ice sheet models that can predict future ice-sheet change and contributions to sea level. However, empirical data are lacking for key areas of outer continental shelves, where the LGM-WAIS must have terminated. We present detailed marine geophysical and geological data documenting an up to ~12 m-thick sequence of glaciomarine sediments within a relict glacial trough in the outer parts of the Amundsen Sea Embayment. Continuous deposition must have persisted here since at least >40 ka BP, pre-dating the established LGM by >13,000 years. Observations constrain the LGM grounding line to a distinct grounding-zone wedge ~100 km inland from the continental shelf edge. Thus, a substantial shelf area (~6000 km**2) remained ice free through the last glacial cycle.

Properties of seawater and particulate matter from a FRRF instrument, operating in a flow-through mode on the continuous surface water sampling system during the Tara Oceans expedition 2009-2013

The Tara Oceans Expedition (2009-2013) sampled the world oceans on board a 36 m long schooner, collecting environmental data and organisms from viruses to planktonic metazoans for later analyses using modern sequencing and state-of-the-art imaging technologies. Tara Oceans Data are particularly suited to study the genetic, morphological and functional diversity of plankton. The present data set provides continuous measurements made with a FRRF instrument, operating in a flow-through mode during the 2009-2012 part of the expedition. It operates by exciting chlorophyll fluorescence using a series of short flashes of controlled energy and time intervals (Kolber et al, 1998). The fluorescence transients produced by this excitation signal were analysed in real-time to provide estimates of abundance of photosynthetic pigments, the photosynthetic yields (Fv/Fm), the functional absorption cross section (a proxy for efficiency of photosynthetic energy acquisition), the kinetics of photosynthetic electron transport between Photosystem II and Photosystem I, and the size of the PQ pool. These parameters were measured at excitation wavelength of 445 nm, 470nm, 505 nm, and 535 nm, allowing to assess the presence and the photosynthetic performance of different phytoplankton taxa based on the spectral composition of their light harvesting pigments. The FRRF-derived photosynthetic characteristics were used to calculate the initial slope, the half saturation, and the maximum level of Photosynthesis vs Irradiance relationship. FRRF data were acquired continuously, at 1-minute time intervals.

Properties of seawater and particulate matter from a WETLabs ALFA hyperspectral laser spectrofluorometer mounted on the continuous surface water sampling system during the Tara Oceans expedition 2009-2013

The Tara Oceans Expedition (2009-2013) sampled the world oceans on board a 36 m long schooner, collecting environmental data and organisms from viruses to planktonic metazoans for later analyses using modern sequencing and state-of-the-art imaging technologies. Tara Oceans Data are particularly suited to study the genetic, morphological and functional diversity of plankton. The present data set provides continuous measurements made with an Aquatic Laser Fluorescence Analyzer (ALFA) (Chekalyuk et al., 2014), connected in-line to the TARA flow through system during 2013. The ALFA instrument provides dual-wavelength excitation (405 and 514 nm) of laser-stimulated emission (LSE) for spectral and temporal analysis. It offers in vivo fluorescence assessments of phytoplankton pigments, biomass, photosynthetic yield (Fv/Fm), phycobiliprotein (PBP)-containing phytoplankton groups, and chromophoric dissolved organic matter (CDOM) (Chekalyuk and Hafez, 2008; 2013A). Spectral deconvolution (SDC) is used to assess the overlapped spectral bands of aquatic fluorescence constituents and water Raman scattering (R). The Fv/Fm measurements are spectrally corrected for non-chlorophyll fluorescence background produced by CDOM and other constituents (Chekalyuk and Hafez, 2008). The sensor was cleaned weakly following the manufacturer recommended protocol.