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.

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.

Enhancement of photosynthetic carbon assimilation efficiency by phytoplankton in the future coastal ocean

A mesocosm experiment was conducted to evaluate the effects of future climate conditions on photosynthesis and productivity of coastal phytoplankton. Natural phytoplankton assemblages were incubated in field mesocosms under the ambient condition (present condition: ca. 400 ppmv CO2 and ambient temp.), and two future climate conditions (acidification condition: ca. 900 ppmv CO2 and ambient temp.; greenhouse condition: ca. 900 ppmv CO2 and 3 °C warmer than ambient). Photosynthetic parameters of steady-state light responses curves (LCs; measured by PAM fluorometer) and photosynthesis-irradiance curves (P-I curves; estimated by in situ incorporation of 14C) were compared to three conditions during the experiment period. Under acidification, electron transport efficiency (alpha LC) and photosynthetic 14C assimilation efficiency (alpha) were 10% higher than those of the present condition, but maximum rates of relative electron transport (rETRm,LC) and photosynthetic 14C assimilation (PBmax) were lower than the present condition by about 19% and 7%, respectively. In addition, rETRm,LC and alpha LC were not significantly different between and greenhouse conditions, but PBmax and alpha of greenhouse conditions were higher than those of the present condition by about 9% and 30%, respectively. In particular, the greenhouse condition has drastically higher PBmax and alpha than the present condition more than 60% during the post-bloom period. According to these results, two future ocean conditions have major positive effects on the photosynthesis in terms of energy utilization efficiency for organic carbon fixation through the inorganic carbon assimilation. Despite phytoplankton taking an advantage on photosynthesis, primary production of phytoplankton was not stimulated by future conditions. In particular, biomass of phytoplankton was depressed under both acidification and greenhouse conditions after the the pre-bloom period, and more research is required to suggest that some factors such as grazing activity could be important for regulating phytoplankton bloom in the future ocean.

Light-modulated responses of growth and photosynthetic performance to ocean acidification in the model diatom Phaeodactylum tricornutum

Ocean acidification (OA) due to atmospheric CO2 rise is expected to influence marine primary productivity. In order to investigate the interactive effects of OA and light changes on diatoms, we grew Phaeodactylum tricornutum, under ambient (390 ppmv; LC) and elevated CO2 (1000 ppmv; HC) conditions for 80 generations, and measured its physiological performance under different light levels (60 µmol/m**2/s, LL; 200 µmol/m**2/s, ML; 460 µmol/m**2/s, HL) for another 25 generations. The specific growth rate of the HC-grown cells was higher (about 12-18%) than that of the LC-grown ones, with the highest under the ML level. With increasing light levels, the effective photochemical yield of PSII (Fv'/Fm') decreased, but was enhanced by the elevated CO2, especially under the HL level. The cells acclimated to the HC condition showed a higher recovery rate of their photochemical yield of PSII compared to the LC-grown cells. For the HC-grown cells, dissolved inorganic carbon or CO2 levels for half saturation of photosynthesis (K1/2 DIC or K1/2 CO2) increased by 11, 55 and 32%, under the LL, ML and HL levels, reflecting a light dependent down-regulation of carbon concentrating mechanisms (CCMs). The linkage between higher level of the CCMs down-regulation and higher growth rate at ML under OA supports the theory that the saved energy from CCMs down-regulation adds on to enhance the growth of the diatom.

Shotgun proteomics reveals physiological response to ocean acidification in Crassostrea gigas

Background. Ocean acidification as a result of increased anthropogenic CO2 emissions is occurring in marine and estuarine environments worldwide. The coastal ocean experiences additional daily and seasonal fluctuations in pH that can be lower than projected end of century open ocean pH reductions. Projected and current ocean acidification have wide-ranging effects on many aquatic organisms, however the exact mechanisms of the impacts of ocean acidification on many of these animals remains to be characterized. Methods. In order to assess the impact of ocean acidification on marine invertebrates, Pacific oysters (Crassostrea gigas) were exposed to one of four different pCO2 levels for four weeks: 400 µatm (pH 8.0), 800 µatm (pH 7.7), 1000 µatm (pH 7.6), or 2800 µatm (pH 7.3). At the end of 4 weeks a variety of physiological parameters were measured to assess the impacts of ocean acidification: tissue glycogen content and fatty acid profile, shell micromechanical properties, and response to acute heat shock. To determine the effects of ocean acidification on the underlying molecular physiology of oysters and their stress response, some of the oysters from 400 µatm and 2800 µatm were exposed to an additional mechanical stress and shotgun proteomics were done on oysters from high and low pCO2 and from with and without mechanical stress. Results. At the end of the four week exposure period, oysters in all four pCO2 environments deposited new shell, but growth rate was not different among the treatments. However, micromechanical properties of the new shell were compromised by elevated pCO2. Elevated pCO2 affected neither whole body fatty acid composition, nor glycogen content, nor mortality rate associated with acute heat shock. Shotgun proteomics revealed that several physiological pathways were significantly affected by ocean acidification, including antioxidant response, carbohydrate metabolism, and transcription and translation. Additionally, the proteomic response to a second stress differed with pCO2, with numerous processes significantly affected by mechanical stimulation at high versus low pCO2 (all proteomics data are available in the ProteomeXchange under the identifier PXD000835). Discussion. Oyster physiology is significantly altered by exposure to elevated pCO2, indicating changes in energy resource use. This is especially apparent in the assessment of the effects of pCO2 on the proteomic response to a second stress. The altered stress response illustrates that ocean acidification may impact how oysters respond to other changes in their environment. These data contribute to an integrative view of the effects of ocean acidification on oysters as well as physiological trade-offs during environmental stress.

Fatty acid and stable isotope composition of calanoid copepods in the open eastern Atlantic Ocean during POLARSTERN cruise ANT-XXIX/1

The majority of global ocean production and total export production is attributed to oligotrophic oceanic regions due to their vast regional expanse. However, energy transfers, food-web structures and trophic relationships in these areas remain largely unknown. Regional and vertical inter- and intra-specific differences in trophic interactions and dietary preferences of calanoid copepods were investigated in four different regions in the open eastern Atlantic Ocean (38°N to 21°S) in October/November 2012 using a combination of fatty acid (FA) and stable isotope (SI) analyses. Mean carnivory indices (CI) based on FA trophic markers generally agreed with trophic positions (TP) derived from d15N analysis. Most copepods were classified as omnivorous (CI ~0.5, TP 1.8 to ~2.5) or carnivorous (CI >=0.7, TP >=2.9). Herbivorous copepods showed typical CIs of <=0.3. Geographical differences in d15N values of epi- (200-0 m) to mesopelagic (1000-200 m) copepods reflected corresponding spatial differences in baseline d15N of particulate organic matter from the upper 100 m. In contrast, species restricted to lower meso- and bathypelagic (2000-1000 m) layers did not show this regional trend. FA compositions were species-specific without distinct intra-specific vertical or spatial variations. Differences were only observed in the southernmost region influenced by the highly productive Benguela Current. Apparently, food availability and dietary composition were widely homogeneous throughout the oligotrophic oceanic regions of the tropical and subtropical Atlantic. Four major species clusters were identified by principal component analysis based on FA compositions. Vertically migrating species clustered with epi- to mesopelagic, non-migrating species, of which only Neocalanus gracilis was moderately enriched in lipids with 16% of dry mass (DM) and stored wax esters (WE) with 37% of total lipid (TL). All other species of this cluster had low lipid contents (< 10% DM) without WE. Of these, the tropical epipelagic Undinula vulgaris showed highest portions of bacterial markers. Rhincalanus cornutus, R. nasutus and Calanoides carinatus formed three separate clusters with species-specific lipid profiles, high lipid contents (>=41% DM), mainly accumulated as WE (>=79% TL). C. carinatus and R. nasutus were primarily herbivorous with almost no bacterial input. Despite deviating feeding strategies, R. nasutus clustered with deep-dwelling, carnivorous species, which had high amounts of lipids (>=37% DM) and WE (>=54% TL). Tropical and subtropical calanoid copepods exhibited a wide variety of life strategies, characterized by specialized feeding. This allows them, together with vertical habitat partitioning, to maintain high abundance and diversity in tropical oligotrophic open oceans, where they play an essential role in the energy flux and carbon cycling.

Euechinoidea and Cidaroidea respond differently to ocean acidification

The impact of the chemical changes in the ocean waters due to the increasing atmospheric CO2 depends on the ability of an organism to control extracellular pH. Among sea urchins, this seems specific to the Euechinoidea, sea urchins except Cidaroidea. However, Cidaroidea survived two ocean acidification periods: the Permian-Trias and the Cretaceous-Tertiary crises. We investigated the response of these two sea urchin groups to reduced seawater pH with the tropical cidaroid Eucidaris tribuloides, the sympatric euechinoid Tripneustes ventricosus and the temperate euechinoid Paracentrotus lividus. Both euechinoid showed a compensation of the coelomic fluid pH due to increased buffer capacity. This was linked to an increased concentration of DIC in the coelomic fluid and thus of bicarbonate ions (most probably originating from the surrounding seawater as isotopic signature of the carbon -delta 13C- was similar). On the other hand, the cidaroid showed no changes within the coelomic fluid. Moreover, the delta 13C of the coelomic fluid did not match that of the seawater and was not significantly different between the urchins from the different treatments. Feeding rate was not affected in any species. While euechinoids are able to regulate their extracellular acid-base balance, many questions are still unanswered on the costs of this capacity. On the contrary, cidaroids do not seem affected by a reduced seawater pH. Further investigations need to be undertaken to cover more species and physiological and metabolic parameters in order to determine if energy trade-offs occur and how this mechanism of compensation is distributed among sea urchins.

Diversity of ocean acidification effects on marine N2 fixers

Considering the important role of N2 fixation for primary productivity and CO2 sequestration, it is crucial to assess the response of diazotrophs to ocean acidification. Previous studies on the genus Trichodesmium suggested a strong sensitivity towards ocean acidification. In view of the large functional diversity in N2 fixers, the objective of this study was to improve our knowledge of the CO2 responses of other diazotrophs. To this end, the single-celled Cyanothece sp. and two heterocystous species, Nodularia spumigena and the symbiotic Calothrix rhizosoleniae, were acclimated to two pCO2 levels (380 vs. 980 µatm). Growth rates, cellular composition (carbon, nitrogen and chlorophyll a) as well as carbon and N2 fixation rates (14C incorporation, acetylene reduction) were measured and compared to literature data on different N2 fixers. The three species investigated in this study responded differently to elevated pCO2, showing enhanced, decreased as well as unaltered growth and production rates. For instance, Cyanothece increased production rates with pCO2, which is in line with the general view that N2 fixers benefit from ocean acidification. Due to lowered growth and production of Nodularia, nitrogen input to the Baltic Sea might decrease in the future. In Calothrix, no significant changes in growth or production could be observed, even though N2 fixation was stimulated under elevated pCO2. Reviewing literature data confirmed a large variability in CO2 sensitivity across diazotrophs. The contrasting response patterns in our and previous studies were discussed with regard to the carbonate chemistry in the respective natural habitats, the mode of N2 fixation as well as differences in cellular energy limitation between the species. The group-specific CO2 sensitivities will impact differently on future biogeochemical cycles of open-ocean environments and systems like the Baltic Sea and should therefore be considered in models estimating climate feedback effects.

Experiment: Food availability outweighs ocean acidification effects in juvenile Mytilus edulis

Ocean acidification is expected to decrease calcification rates of bivalves. Nevertheless in many coastal areas high pCO2 variability is encountered already today. Kiel Fjord (Western Baltic Sea) is a brackish (12-20 g kg-1) and CO2 enriched habitat, but the blue mussel Mytilus edulis dominates the benthic community. In a coupled field and laboratory study we examined the annual pCO2 variability in this habitat and the combined effects of elevated pCO2 and food availability on juvenile M. edulis growth and calcification. In the laboratory experiment, mussel growth and calcification were found to chiefly depend on food supply, with only minor impacts of pCO2 up to 3350 µatm. Kiel Fjord was characterized by strong seasonal pCO2 variability. During summer, maximal pCO2 values of 2500 µatm were observed at the surface and >3000 µatm at the bottom. However, the field growth experiment revealed seven times higher growth and calcification rates of M. edulis at a high pCO2 inner fjord field station (mean pCO2 ca. 1000 µatm) in comparison to a low pCO2 outer fjord station (ca. 600 µatm). In addition, mussels were able to outcompete the barnacle Amphibalanus improvisus at the high pCO2 site. High mussel productivity at the inner fjord site was enabled by higher particulate organic carbon concentrations. Kiel Fjord is highly impacted by eutrophication, which causes bottom water hypoxia and consequently high seawater pCO2. At the same time, elevated nutrient concentrations increase the energy availability for filter feeding organisms such as mussels. Thus M. edulis can dominate over a seemingly more acidification resistant species such as A. improvisus. We conclude that benthic stages of M. edulis tolerate high ambient pCO2 when food supply is abundant and that important habitat characteristics such as species interactions and energy availability need to be considered to predict species vulnerability to ocean acidification.

Effects of Ocean Acidification on the Ballast of Surface Aggregates Sinking through the Twilight Zone

The dissolution of CaCO3 is one of the ways ocean acidification can, potentially, greatly affect the ballast of aggregates. A diminution of the ballast could reduce the settling speed of aggregates, resulting in a change in the carbon flux to the deep sea. This would mean lower amounts of more refractory organic matter reaching the ocean floor. This work aimed to determine the effect of ocean acidification on the ballast of sinking surface aggregates. Our hypothesis was that the decrease of pH will increase the dissolution of particulate inorganic carbon ballasting the aggregates, consequently reducing their settling velocity and increasing their residence time in the upper twilight zone. Using a new methodology for simulation of aggregate settling, our results suggest that future pCO2 conditions can significantly change the ballast composition of sinking aggregates. The change in aggregate composition had an effect on the size distribution of the aggregates, with a shift to smaller aggregates. A change also occurred in the settling velocity of the particles, which would lead to a higher residence time in the water column, where they could be continuously degraded. In the environment, such an effect would result in a reduction of the carbon flux to the deep-sea. This reduction would impact those benthic communities, which rely on the vertical flow of carbon as primary source of energy.

Field measurements of the atmosphere, ocean, sea ice and sub-ice platelet layer at Atka Bay in 2013

Ice shelves strongly interact with coastal Antarctic sea ice and the associated ecosystem by creating conditions favourable to the formation of a sub-ice platelet layer. The close investigation of this phenomenon and its seasonal evolution remain a challenge due to logistical constraints and a lack of suitable methodology. In this study, we characterize the seasonal cycle of Antarctic fast ice adjacent to the Ekström Ice Shelf in the eastern Weddell Sea. We used a thermistor chain with the additional ability to record the temperature response induced by cyclic heating of resistors embedded in the chain. Vertical sea-ice temperature and heating profiles obtained daily between November 2012 and February 2014 were analyzed to determine sea-ice and snow evolution, and to calculate the basal energy budget. The residual heat flux translated into an ice-volume fraction in the platelet layer of 0.18 ± 0.09, which we reproduced by a independent model simulation and agrees with earlier results. Manual drillings revealed an average annual platelet-layer thickness increase of at least 4m, and an annual maximum thickness of 10m beneath second-year sea ice. The oceanic contribution dominated the total sea-ice production during the study, effectively accounting for up to 70% of second-year sea-ice growth. In summer, an oceanic heat flux of 21 W/m**2 led to a partial thinning of the platelet layer. Our results further show that the active heating method, in contrast to the acoustic sounding approach, is well suited to derive the fast-ice mass balance in regions influenced by ocean/ice-shelf interaction, as it allows sub-diurnal monitoring of the platelet-layer thickness.

Photosynthetic activity buffers ocean acidification in seagrass meadows

Macrophytes growing in shallow coastal zones characterised by intense metabolic activity have the capacity to modify pH within their canopy and beyond. We observed diel pH changes in shallow (5-12 m) seagrass (Posidonia oceanica) meadows spanning 0.06 pH units in September to 0.24 units in June. The carbonate system (pH, DIC, and aragonite saturation state (omega Ar)) and O2 within the meadows displayed strong diel variability driven by primary productivity, and changes in chemistry were related to structural parameters of the meadow, in particular, the leaf surface area available for photosynthesis (LAI). LAI was positively correlated to mean, max and range pHNBS and max and range omega Ar. In June, vertical mixing (as Turbulent Kinetic Energy) influenced max and min omega Ar, while in September there was no effect of hydrodynamics on the carbonate system within the canopy. Max and range omega Ar within the meadow showed a positive trend with the calcium carbonate load of the leaves, pointing to a possible link between structural parameters, omega Ar and carbonate deposition.

Upper-ocean microstructure data from the tropical Atlantic

SST variability within the Atlantic cold tongue (ACT) region is of climatic relevance for the surrounding continents. A multi cruise data set of microstructure observations is used to infer regional as well as seasonal variability of upper ocean mixing and diapycnal heat flux within the ACT region. The variability in mixing intensity is related to the variability in large scale background conditions, which were additionally observed during the cruises. The observations indicate fundamental differences in background conditions in terms of shear and stratification below the mixed layer (ML) for the western and eastern equatorial ACT region causing critical Froude numbers (Fr) to be more frequently observed in the western equatorial ACT. The distribution of critical Fr occurrence below the ML reflects the regional and seasonal variability of mixing intensity. Turbulent dissipation rates (?) at the equator (2°N-2°S) are strongly increased in the upper thermocline compared to off-equatorial locations. In addition, ? is elevated in the western equatorial ACT compared to the east from May to November, whereas boreal summer appears as the season of highest mixing intensities throughout the equatorial ACT region, coinciding with ACT development. Diapycnal heat fluxes at the base of the ML in the western equatorial ACT region inferred from ? and stratification range from a maximum of 90 Wm-2 in boreal summer to 55 Wm-2 in September and 40 Wm-2 in November. In the eastern equatorial ACT region maximum values of about 25 Wm-2 were estimated during boreal summer reducing to about 5 Wm-2 towards the end of the year. Outside the equatorial region, inferred diapycnal heat fluxes are comparably low rarely exceeding 10 Wm-2. Integrating the obtained heat flux estimates in the ML heat budget at 10°W on the equator accentuates the diapycnal heat flux as the largest ML cooling term during boreal summer and early autumn. In the western equatorial ACT elevated meridional velocity shear in the upper thermocline contributes to the enhanced diapycnal heat flux within this region during boreal summer and autumn. The elevated meridional velocity shear appears to be associated with intra-seasonal wave activity.