Impact of thermal stress on evolutionary trajectories of pathogen resistance in three-spined stickleback (Gasterosteus aculeatus), supplementary material

Background: Pathogens are a major regulatory force for host populations, especially under stressful conditions. Elevated temperatures may enhance the development of pathogens, increase the number of transmission stages, and can negatively influence host susceptibility depending on host thermal tolerance. As a net result, this can lead to a higher prevalence of epidemics during summer months. These conditions also apply to marine ecosystems, where possible ecological impacts and the population-specific potential for evolutionary responses to changing environments and increasing disease prevalence are, however, less known. Therefore, we investigated the influence of thermal stress on the evolutionary trajectories of disease resistance in three marine populations of three-spined sticklebacks Gasterosteus aculeatus by combining the effects of elevated temperature and infection with a bacterial strain of Vibrio sp. using a common garden experiment. Results: We found that thermal stress had an impact on fish weight and especially on survival after infection after only short periods of thermal acclimation. Environmental stress reduced genetic differentiation (QST) between populations by releasing cryptic within-population variation. While life history traits displayed positive genetic correlations across environments with relatively weak genotype by environment interactions (GxE), environmental stress led to negative genetic correlations across environments in pathogen resistance. This reversal of genetic effects governing resistance is probably attributable to changing environment-dependent virulence mechanisms of the pathogen interacting differently with host genotypes, i.e. GPathogenxGHostxE or (GPathogenxE)x(GHostxE) interactions, rather than to pure host genetic effects, i.e. GHostxE interactions. Conclusion: To cope with climatic changes and the associated increase in pathogen virulence, host species require wide thermal tolerances and pathogen-resistant genotypes. The higher resistance we found for some families at elevated temperatures showed that there is evolutionary potential for resistance to Vibrio sp. in both thermal environments. The negative genetic correlation of pathogen resistance between thermal environments, on the other hand, indicates that adaptation to current conditions can be a weak predictor for performance in changing environments. The observed feedback on selective gradients exerted on life history traits may exacerbate this effect, as it can also modify the response to selection for other vital components of fitness.

Effects of seawater pCO2 and temperature on shell growth, shell stability, condition and cellular stress of Western Baltic Sea Mytilus edulis (L.) and Arctica islandica (L.)

Acidification of the World's oceans may directly impact reproduction, performance and shell formation of marine calcifying organisms. In addition, since shell production is costly and stress in general draws on an organism's energy budget, shell growth and stability of bivalves should indirectly be affected by environmental stress. The aim of this study was to investigate whether a combination of warming and acidification leads to increased physiological stress (lipofuscin accumulation and mortality) and affects the performance [shell growth, shell breaking force, condition index (Ci)] of young Mytilus edulis and Arctica islandica from the Baltic Sea. We cultured the bivalves in a fully-crossed 2-factorial experimental setup (seawater (sw) pCO2 levels "low", "medium" and "high" for both species, temperature levels 7.5, 10, 16, 20 and 25 °C for M. edulis and 7.5, 10 and 16 °C for A. islandica) for 13 weeks in summer. Mytilus edulis and A. islandica appeared to tolerate wide ranges of sw temperature and pCO2. Lipofuscin accumulation of M. edulis increased with temperature while the Ci decreased, but shell growth of the mussels only sharply decreased while its mortality increased between 20 and 25 °C. In A. islandica, lipofuscin accumulation increased with temperature, whereas the Ci, shell growth and shell breaking force decreased. The pCO2 treatment had only marginal effects on the measured parameters of both bivalve species. Shell growth of both bivalve species was not impaired by under-saturation of the sea water with respect to aragonite and calcite. Furthermore, independently of water temperatures shell breaking force of both species and shell growth of A. islandica remained unaffected by the applied elevated sw pCO2 for several months. Only at the highest temperature (25 °C), growth arrest of M. edulis was recorded at the high sw pCO2 treatment and the Ci of M. edulis was slightly higher at the medium sw pCO2 treatment than at the low and high sw pCO2 treatments. The only effect of elevated sw pCO2 on A. islandica was an increase in lipofuscin accumulation at the high sw pCO2 treatment compared to the medium sw pCO2 treatment. Our results show that, despite this robustness, growth of both M. edulis and A. islandica can be reduced if sw temperatures remain high for several weeks in summer. As large body size constitutes an escape from crab and sea star predation, this can make bivalves presumably more vulnerable to predation with possible negative consequences on population growth. In M. edulis, but not in A. islandica, this effect is amplified by elevated sw pCO2. We follow that combined effects of elevated sw pCO2 and ocean warming might cause shifts in future Western Baltic Sea community structures and ecosystem services; however, only if predators or other interacting species do not suffer as strong from these stressors.

Planktonic foraminifera in the latest Maastrichtian of Walvis Ridge, South Atlantic

An abrupt global warming of 3-4°C occurred near the end of the Maastrichtian at 65.45-65.10 Ma. The environmental effects of this warm event are here documented based on stable isotopes and quantitative analysis of planktonic foraminifera at the South Atlantic DSDP Site 525A. Stable isotopes of individual species mark a rapid increase in temperature and a reduction in the vertical water mass stratification that is accompanied by a decrease in niche habitats, reduced species diversity and/or abundance, smaller species morphologies or dwarfing, and reduced photosymbiotic activity. During the warm event, the relative abundance of a large number of species decreased, including tropical-subtropical affiliated species, whereas typical mid-latitude species retained high abundances. This indicates that climate warming did not create favorable conditions for all tropical-subtropical species at mid-latitudes and did not cause a massive retreat in the local mid-latitude population. A noticeable exception is the ecological generalist Heterohelix dentata Stenestad that dominated during the cool intervals, but significantly decreased during the warm event. However, dwarfing is the most striking response to the abrupt warming and occurred in various species of different morphologies and lineages (e.g. biserial, trochospiral, keeled globotruncanids). Dwarfing is a typical reaction to environmental stress conditions and was likely the result of increased reproduction rates. Similarly, photosymbiotic activity appears to have been reduced significantly during the maximum warming, as indicated by decreased delta13C values. The foraminiferal response to climate change is thus multifaceted resulting in decreased species diversity, decreased species populations, increased competition due to reduced niche habitats, dwarfing and reduced photosymbiotic activity.

Distribution of Recent benthic foraminifera in Arctic galciomarine sediments

Foraminifera were examined in recent (<100 years) fine-grained glaciomarine muds from surface sediments and cores from Nordensheld Bay, Novaja Zemlja, and Hornsund and Bellsund, Spitsbergen. This study presents the first data on modern foraminifera distribution for fjord environments in Novaja Zemlja, Russia. The data are interpreted with reference to the distribution of foraminiferal near Svalbard and the Barents Sea. In Nordensheld Bay, live and dead Nonionellina labradorica and Islandiella norcrossi are most abundant in the outer fjord. Cassidulina reniforme and Allogromiina spp. dominate in the middle and inner fjord. The dominant species are dissimilar to species occurring in other areas of the Barents Sea region, with the exception of Svalbard fjords. The number of live foraminifera (24 to 122 tests/10 cm1) in outer and middle Nordensheld Bay corresponds with values known from the open Barents Sea. However, the biomass (0.03 mg/10 cm**3) is two orders of magnitude less due to smaller foraminiferal test size, which in glaciomarine sediments reflects the absence of larger species, paucity of large specimens, and high occurrence of juvenile foraminifera. The smaller size indicates an opportunistic response to environmental stress due to glacier proximity. The presence of Quinqueloculina stalkeri is diagnostic of glaciomarine environments in fjords of Novaja Zemlja and Svalbard.

Maastrichtian planktic foraminifera of the South Atlantic Ocean

Stratigraphic, faunal and isotopic analyses of the Maastrichtian at DSDP sites 525A and 21 in the South Atlantic reveal a planktic foraminiferal fauna characterized by two major events, an early late Maastrichtian diversification and end-Maastrichtian mass extinction. Both events are accompanied by major changes in climate and productivity. The diversification event which occurred in two steps between 70.5 and 69.1 Ma increased species richness by a total of 43% and coincided with the onset of major cooling in surface and bottom waters and increased surface productivity. The onset of the terminal decline in Maastrichtian species richness began at 67.5 Ma and the first significant decline in surface productivity occurred at 66.2 Ma, coincident maximum cooling to 13°C in surface waters and the reduction of the surface-to-deep temperature gradient to less than 5°C. Major climatic and moderate productivity changes mark the mass extinction and the last 500 kyr of the Maastrichtian. Between 200 and 400 kyr before the K-T boundary surface and deep waters warmed rapidly by 3-4°C and cooled again during the last 100 kyr of the Maastrichtian. Surface productivity decreased only moderately across the K-T boundary. Species richness began to decline during the late Maastrichtian cooling and by K-T boundary time, the mass extinction had claimed 66% of the species. Viewed within the context of Maastrichtian climate and productivity changes, the K-T mass extinction could have resulted from extreme environmental stress even without the addition of an extraterrestrial impact.

Planktic foraminifera across the Cretaceous-Tertiary transition in the Antarctic Ocean

Three Antarctic Ocean K/T boundary sequences from ODP Site 738C on the Kerguelen Plateau, ODP Site, 752B on Broken Ridge and ODP Site 690C on Maud Rise, Weddell Sea, have been analyzed for stratigraphic completeness and faunal turnover based on quantitative planktic foraminiferal studies. Results show that Site 738C, which has a laminated clay layer spanning the K/T boundary, is biostratigraphically complete with the earliest Tertiary Zones P0 and P1a present, but with short intrazonal hiatuses. Site 752B may be biostratigraphically complete and Site 690C has a hiatus at the K/T boundary with Zones P0 and P1a missing. Latest Cretaceous to earliest Tertiary planktic foraminiferal faunas from the Antarctic Ocean are cosmopolitan and similar to coeval faunas dominating in low, middle and northern high latitudes, although a few endemic species are present. This allows application of the current low and middle latitude zonation to Antarctic K/T boundary sequences. The most abundant endemic species is Chiloguembelina waiparaensis, which was believed to have evolved in the early Tertiary, but which apparently evolved as early as Chron 30N at Site 738C. Since this species is only rare in sediments of Site 690C in the Weddell Sea, this suggests that a watermass oceanographic barner may have existed between the Indian and Atlantic Antarctic Oceans. The cosmopolitan nature of the dominant fauna began during the last 200,000 to 300,000 years of the Cretaceous and continued at least 300,000 years into the Tertiary. This indicates a long-term environmental crisis that led to gradual elimination of specialized forms and takeover by generalists tolerant of wide ranging temperature, oxygen, salinity and nutrient conditions. A few thousand years before the K/T boundary these generalists gradually declined in abundance and species became generally dwarfed due to increased environmental stress. There is no evidence of a sudden mass killing of the Cretaceous fauna associated with a bolide impact at the K/T boundary. Instead, the already declining Cretaceous taxa gradually disappear in the early Danian and the opportunistic survivor taxa (Ch. waiparaensis and Guembelitria cretacea) increase in relative abundance coincident with the evolution of the first new Tertiary species.

Ambulacral length with number of pore pairs/tubercles for both species of sea urchins and for six sampling sites to quantify fluctuating asymmetry

The surroundings of the Cortiou sewage are among the most polluted environments of the French Mediterranean Sea (Marseilles, France). So far, no studies have precisely quantified the impact of pollution on the development of organisms in this area.Methods: We used a fluctuating asymmetry (FA) measure of developmental instability (DI) to assess environmental stress in two species of radially symmetric sea urchins (Arbacia lixula and Paracentrotus lividus). For six sampling sites (Cortiou, Riou, Maire, East Maire, Mejean, and Niolon), levels of FA were calculated from continuous and discrete skeletal measures of ambulacral length, number of pore pairs and primary tubercles.Results: For both species, the most polluted sampling site, Cortiou, displayed the highest level of FA, while the Maire and East Maire sampling sites displayed the lowest levels. A. lixula revealed systematic differences in FA among sampling sites for all characters and P. lividus showed differences in FA for the number of primary tubercles.Conclusions: Statistical analyses of FA show a concordance between the spatial patterns of FA among sampling sites and the spatial distribution of sewage discharge pollutants in the Cortiou area. High developmental stress in these sampling sites is associated with exposure to high concentrations of heavy metals and many harmful organic substances contained in wastewater. FA estimated from structures with complex symmetry appears to be a fast and reliable tool to detect subtle differences in FA. Its use in biomonitoring programs for inferring anthropogenic and natural environmental stress is suggested.

(Table 1) Concentration of pollutants in hair of East Greenland polar bears (Ursus maritimus)

Hair sampled from 96 East Greenland polar bears (Ursus maritimus) over the periods 1892-1927 and 1988-2009 was analyzed for cortisol as a proxy to investigate temporal patterns of environmental stress. Cortisol concentration was independent of sex and age, and was found at significantly higher (p<0.001) concentrations in historical hair samples (1892-1927; n = 8) relative to recent ones (1988-2009; n = 88). In addition, there was a linear time trend in cortisol concentration of the recent samples (p< 0.01), with an annual decrease of 2.7%. The recent hair samples were also analyzed for major bioaccumulative, persistent organic pollutants (POPs). There were no obvious POP related time trends or correlations between hair cortisol and hair POP concentrations. Thus, polar bear hair appears to be a relatively poor indicator of the animal's general POP load in adipose tissue. However, further investigations are warranted to explore the reasons for the temporal decrease found in the bears' hair cortisol levels.

Warming and acidification experiment on early life-stage F. vesiculosus under natural fluctuation and seasonal variations (spring 2013 - 2014)

Genetic diversity of baltic F. vesiculosus is low compared to other populations which might jeopardize their potential for adaptation to climate change. Especially the early life-stage F. vesiculosus may be threaten by ocean warming and acidification. To test this, we exposed F. vesiculosus germlings to warming and acidification in the near-natural scenario in the "Kiel Outdoor Benthocosms" maintaining the natural variation of the Kiel Fjord, Germany (54°27 'N, 10°11 'W) in all seasons (spring 2013 - 2014). Warming was simulated by using a delta treatment adding 5 °C and by increasing pCO2 at 1000 µatm. Warming positively affected germlings' growth in spring and in summer but decreased non-photochemical quenching in spring and survival in summer. Acidified conditions showed much weaker effects than warming. The high genotypic variation in stress sensitivity as well as the enhanced survival at high diversity levels indicate higher potential for adaptation for genetically diverse populations. We conclude that the combination of stressors and season determines the sensitivity to environmental stress and that genetic variation is crucial for the adaptation to climate change stress.

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.

1H NMR Metabolomics Reveals Contrasting Response by Male and Female Mussels Exposed to Reduced Seawater pH, Increased Temperature, and a Pathogen

Human activities are fundamentally altering the chemistry of the world's oceans. Ocean acidification (OA) is occurring against a background of warming and an increasing occurrence of disease outbreaks, posing a significant threat to marine organisms, communities, and ecosystems. In the current study, 1H NMR spectroscopy was used to investigate the response of the blue mussel, Mytilus edulis, to a 90-day exposure to reduced seawater pH and increased temperature, followed by a subsequent pathogenic challenge. Analysis of the metabolome revealed significant differences between male and female organisms. Furthermore, males and females are shown to respond differently to environmental stress. While males were significantly affected by reduced seawater pH, increased temperature, and a bacterial challenge, it was only a reduction in seawater pH that impacted females. Despite impacting males and females differently, stressors seem to act via a generalized stress response impacting both energy metabolism and osmotic balance in both sexes. This study therefore has important implications for the interpretation of metabolomic data in mussels, as well as the impact of environmental stress in marine invertebrates in general.

Effect of ocean acidification and pH fluctuations on the growth and development of coralline algal recruits, and an associated benthic algal assemblage

Coralline algae are susceptible to the changes in the seawater carbonate system associated with ocean acidification (OA). However, the coastal environments in which corallines grow are subject to large daily pH fluctuations which may affect their responses to OA. Here, we followed the growth and development of the juvenile coralline alga Arthrocardia corymbosa, which had recruited into experimental conditions during a prior experiment, using a novel OA laboratory culture system to simulate the pH fluctuations observed within a kelp forest. Microscopic life history stages are considered more susceptible to environmental stress than adult stages; we compared the responses of newly recruited A. corymbosa to static and fluctuating seawater pH with those of their field-collected parents. Recruits were cultivated for 16 weeks under static pH 8.05 and 7.65, representing ambient and 4*preindustrial pCO2 concentrations, respectively, and two fluctuating pH treatments of daily (daytime pH = 8.45, night-time pH = 7.65) and daily (daytime pH = 8.05, night-time pH = 7.25). Positive growth rates of new recruits were recorded in all treatments, and were highest under static pH 8.05 and lowest under fluctuating pH 7.65. This pattern was similar to the adults' response, except that adults had zero growth under fluctuating pH 7.65. The % dry weight of MgCO3 in calcite of the juveniles was reduced from 10% at pH 8.05 to 8% at pH 7.65, but there was no effect of pH fluctuation. A wide range of fleshy macroalgae and at least 6 species of benthic diatoms recruited across all experimental treatments, from cryptic spores associated with the adult A. corymbosa. There was no effect of experimental treatment on the growth of the benthic diatoms. On the community level, pH-sensitive species may survive lower pH in the presence of diatoms and fleshy macroalgae, whose high metabolic activity may raise the pH of the local microhabitat.

Wave-induced changes in seaweed toughness entail plastic modifications in snail traits maintaining consumption efficacy, link to supplementary material

1. Environmental stress can influence species traits and performance considerably. Using a seaweed-snail system from NW (Nova Scotia) and NE (Helgoland) Atlantic rocky shores, we examined how physical stress (wave exposure) modulates traits in the seaweed Fucus vesiculosus and indirectly in its main consumer, the periwinkle Littorina obtusata. 2. In both regions, algal tissue toughness increased with wave exposure. Reciprocal-transplant experiments showed that tissue toughness adjusts plastically to the prevailing level of wave exposure. 3. Choice experiments tested the feeding preference of snails from sheltered, exposed, and very exposed habitats for algae from such wave exposures. Snails from exposed and very exposed habitats consumed algal tissues at similar rates irrespective of the exposure of origin of the algae. However, snails from sheltered habitats consumed less algal tissues from very exposed habitats than tissues from sheltered and exposed habitats. Choice assays using reconstituted algal food (triturated during preparation) identified high thallus toughness as the explanation for the low preference of snails from sheltered habitats for algae from very exposed habitats. 4. Ultrastructural analyses of radulae indicated that rachidian teeth were longest and the number of cusps in lateral teeth (grazing-relevant traits) was highest in snails from very exposed habitats, suggesting that radulae are best suited to rupture tough algal tissues in such snails. 5. No-choice feeding experiments revealed that these radular traits are also phenotypically plastic, as they adjust to the toughness of the algal food. 6. Synthesis. This study indicates that the observed plasticity in the feeding ability of snails is mediated by wave exposure through phenotypic plasticity in the tissue toughness of algae. Thus, plasticity in consumers and their resource species may reduce the potential effects of physical stress on their interaction.

Does Encapsulation Protect Embryos from the Effects of Ocean Acidification? The Example of Crepidula fornicata

Early life history stages of marine organisms are generally thought to be more sensitive to environmental stress than adults. Although most marine invertebrates are broadcast spawners, some species are brooders and/or protect their embryos in egg or capsules. Brooding and encapsulation strategies are typically assumed to confer greater safety and protection to embryos, although little is known about the physico-chemical conditions within egg capsules. In the context of ocean acidification, the protective role of encapsulation remains to be investigated. To address this issue, we conducted experiments on the gastropod Crepidula fornicata. This species broods its embryos within capsules located under the female and veliger larvae are released directly into the water column. C. fornicata adults were reared at the current level of CO2 partial pressure (pCO2) (390 µatm) and at elevated levels (750 and 1400 µatm) before and after fertilization and until larval release, such that larval development occurred entirely at a given pCO2. The pCO2 effects on shell morphology, the frequency of abnormalities and mineralization level were investigated on released larvae. Shell length decreased by 6% and shell surface area by 11% at elevated pCO2 (1400 µatm). The percentage of abnormalities was 1.5- to 4-fold higher at 750 µatm and 1400 µatm pCO2, respectively, than at 390 µatm. The intensity of birefringence, used as a proxy for the mineralization level of the larval shell, also decreased with increasing pCO2. These negative results are likely explained by increased intracapsular acidosis due to elevated pCO2 in extracapsular seawater. The encapsulation of C. fornicata embryos did not protect them against the deleterious effects of a predicted pCO2 increase. Nevertheless, C. fornicata larvae seemed less affected than other mollusk species. Further studies are needed to identify the critical points of the life cycle in this species in light of future ocean acidification.

Asian green mussels (Perna viridis) transplanted between Jakarta Bay and Lada Bay, West Java, Indonesia (April 2012 - November 2013)

The Asian green mussel Perna viridis is tolerant to environmental stress, but its robustness varies between populations from habitats that differ in quality. So far, it is unclear whether local adaptations through stressinduced selection or phenotypic plasticity are responsible for these inter-population differences. We tested for the relevance of both mechanisms by comparing survival under hypoxia in mussels that were transplanted from an anthropogenically impacted (Jakarta Bay, Indonesia) to a natural habitat (Lada Bay, Indonesia) and vice versa. Mussels were retrieved 8 weeks after transplantation and exposed to hypoxia in the laboratory. Additional hypoxia tests were conducted with juvenile mussels collected directly from both sites. To elucidate possible relationships between habitat quality and mussel tolerance, we monitored concentrations of inorganic nutrients, temperature, dissolved oxygen, salinity, phytoplankton density and the mussels' body condition index (BCI) for 20 months before, during and after the experiments. Survival under hypoxia depended mainly on the quality of the habitat where the mussels lived before the hypoxia tests and only to a small degree on their site of origin. Furthermore, stress tolerance was only higher in Jakarta than in Lada Bay mussels when the BCIs were substantially higher, which in turn correlated with the phytoplankton densities. We explain why phenotypic plasticity and high BCIs are more likely the causes of populationspecific differences in hypoxia tolerance in P. viridis than stress-induced selection for robust genotypes. This is relevant to understanding the role of P. viridis as mariculture organism in eutrophic ecosystems and invasive species in the (sub)tropical world.