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.

Pollen analysis, particle size analysis, end-member modelling and x-ray fluorescence measurements from sediment cores off Cape Ghir, Morocco

Understanding past human-climate-environment interactions is essential for assessing the vulnerability of landscapes and ecosystems to future climate change. This is particularly important in southern Morocco where the current vegetation is impacted by pastoralism, and the region is highly sensitive to climate variability. Here, we present a 2000-year record of vegetation, sedimentation rate, XRF chemical element intensities, and particle size from two decadal-resolved, marine sediment cores, raised from offshore Cape Ghir, southern Morocco. The results show that between 650 and 850 AD the sedimentation rate increased dramatically from 100 cm/1000 years to 300 cm/1000 years, and the Fe/Ca and pollen flux doubled, together indicating higher inputs of terrestrial sediment. Particle size measurements and end-member modelling suggest increased fluvial transport of the sediment. Beginning at 650 AD pollen levels from Cichorioideae species show a sharp rise from 10% to 20%. Pollen from Atemisia and Plantago, also increase from this time. Deciduous oak pollen percentages show a decline, whereas those of evergreen oak barely change. The abrupt increase in terrestrial/fluvial input from 650 to 850 AD occurs, within the age uncertainty, of the arrival of Islam (Islamisation) in Morocco at around 700 AD. Historical evidence suggests Islamisation led to population increase and development of southern Morocco, including expanded pastoralism, deforestation and agriculture. Livestock pressure may have changed the vegetation structure, accounting for the increase in pollen from Cichorioideae, Plantago, and Artemisia, which include many weedy species. Goats in particular may have played a dominant role as agents of erosion, and intense browsing may have led to the decline in deciduous oak; evergreen oak is more likely to survive as it re-sprouts more vigorously after browsing. From 850 AD to present sedimentation rates, Fe/Ca ratios and fluvial discharge remain stable, whereas pollen results suggest continued degradation. Pollen results from the past 150 years suggest expanded cultivation of olives and the native argan tree, and the introduction of Australian eucalyptus trees. The rapidly increasing population in southern Morocco is causing continued pressure to expand pastoralism and agriculture. The history of land degradation presented here suggests that the vegetation in southern Morocco may have been degraded for a longer period than previously thought and may be particularly sensitive to further land use changes. These results should be included in land management strategies for southern Morocco.

Gasterosteus aculeatus egg traits and fish lengths vs. experiment setup

1) Our study addresses the role of non-genetic and genetic inheritance in shaping the adaptive potential of populations under a warming ocean scenario. We used a combined experimental approach (transgenerational plasticity and quantitative genetics) to partition the relative contribution of maternal vs. paternal (additive genetic) effects to offspring body size (a key component of fitness), and investigated a potential physiological mechanism (mitochondrial respiration capacities) underlying whole organism growth/size responses. 2) In very early stages of growth (up to 30 days), offspring body size of marine sticklebacks benefited from maternal transgenerational plasticity (TGP): offspring of mothers acclimated to17°C were larger when reared at 17°C, and offspring of mothers acclimated to 21°C were larger when reared at 21°C. The benefits of maternal TGP on body size were stronger and persisted longer (up to 60 days) for offspring reared in the warmer (21°C) environment, suggesting that maternal effects will be highly relevant for climate change scenarios in this system. 3) Mitochondrial respiration capacities measured on mature offspring (F1 adults) matched the pattern of TGP for juvenile body size, providing an intuitive mechanistic basis for the maternal acclimation persisting into adulthood. Size differences between temperatures seen at early growth stages remained in the F1 adults, linking offspring body size to maternal inheritance of mitochondria. 4) Lower maternal variance components in the warmer environment were mostly driven by mothers acclimated to ambient (colder) conditions, further supporting our tenet that maternal effects were stronger at elevated temperature. Importantly, all parent-offspring temperature combination groups showed genotype x environment (GxE) interactions, suggesting that reaction norms have the potential to evolve. 5) To summarise, transgenerational plasticity and genotype x environment interactions work in concert to mediate impacts of ocean warming on metabolic capacity and early growth of marine sticklebacks. TGP can buffer short-term detrimental effects of climate warming and may buy time for genetic adaptation to catch up, therefore markedly contributing to the evolutionary potential and persistence of populations under climate change.

Theoretical model simulations for variations in thermal reaction norms

Thermal reaction norms for growth rates of six Emiliania huxleyi isolates originating from the central Atlantic (Azores, Portugal) and five isolates from the coastal North Atlantic (Bergen, Norway) were assessed. We used the template mode of variation model to decompose variations in growth rates into modes of biological interest: vertical shift, horizontal shift, and generalist-specialist variation. In line with the actual habitat conditions, isolates from Bergen (Bergen population) grew well at lower temperatures, and isolates from the Azores (Azores population) performed better at higher temperatures. The optimum growth temperature of the Azores population was significantly higher than that of the Bergen population. Neutral genetic differentiation was found between populations by microsatellite analysis. These findings indicate that E. huxleyi populations are adapted to local temperature regimes. Next to between-population variation, we also found variation within populations. Genotype-by-environment interactions resulted in the most pronounced phenotypic differences when isolates were exposed to temperatures outside the range they naturally encounter. Variation in thermal reaction norms between and within populations emphasizes the importance of using more than one isolate when studying the consequences of global change on marine phytoplankton. Phenotypic plasticity and standing genetic variation will be important in determining the potential of natural E. huxleyi populations to cope with global climate change.

Molecular fingerprinting and amplicon sequencing of the bacterial biofilm on settlement tiles deployed along natural pH gradients in Papua New Guinea

Bacterial biofilms provide cues for the settlement of marine invertebrates such as coral larvae, and are therefore important for the resilience and recovery of coral reefs. This study aimed to better understand how ocean acidification may affect the community composition and diversity of bacterial biofilms on surfaces under naturally reduced pH conditions. Settlement tiles were deployed at coral reefs in Papua New Guinea along pH gradients created by two CO2 seeps, and upper and lower tiles surfaces were sampled 5 and 13 months after deployment. Automated Ribosomal Intergenic Spacer Analysis were used to characterize more than 200 separate bacterial communities, complemented by amplicon sequencing of the bacterial 16S rRNA gene of 16 samples. The bacterial biofilm consisted predominantly of Alpha-, Gamma- and Deltaproteobacteria, as well as Cyanobacteria, Flavobacteriia and Cytophaga, whereas putative settlement-inducing taxa only accounted for a small fraction of the community. Bacterial biofilm composition was heterogeneous with approximately 25% shared operational taxonomic units between samples. Among the observed environmental parameters, pH only had a weak effect on community composition (R² ~ 1%) and did not affect community richness and evenness. In contrast, there were strong differences between upper and lower surfaces (contrasting in light exposure and grazing intensity). There also appeared to be a strong interaction between bacterial biofilm composition and the macroscopic components of the tile community. Our results suggest that on mature settlement surfaces in situ, pH does not have a strong impact on the composition of bacterial biofilms. Other abiotic and biotic factors such as light exposure and interactions with other organisms may be more important in shaping bacterial biofilms than changes in seawater pH.

Transgenerational effects persist down the maternal line in marine sticklebacks: gene expression matches physiology in a warming ocean, supplementary material

Transgenerational effects can buffer populations against environmental change, yet little is known about underlying mechanisms, their persistence, or the influence of environmental cue timing. We investigated mitochondrial respiratory capacity (MRC) and gene expression of marine sticklebacks that experienced acute or developmental acclimation to simulated ocean warming (21°C) across three generations. Previous work showed that acute acclimation of grandmothers to 21°C led to lower (optimised) offspring MRCs. Here, developmental acclimation of mothers to 21°C led to higher, but more efficient offspring MRCs. Offspring with a 21°Cx17°C grandmother-mother environment mismatch showed metabolic compensation: their MRCs were as low as offspring with a 17°C thermal history across generations. Transcriptional analyses showed primarily maternal but also grandmaternal environment effects: genes involved in metabolism and mitochondrial protein biosynthesis were differentially expressed when mothers developed at 21°C, whereas 21°C grandmothers influenced genes involved in hemostasis and apoptosis. Genes involved in mitochondrial respiration all showed higher expression when mothers developed at 21° and lower expression in the 21°Cx17°C group, matching the phenotypic pattern for MRCs. Our study links transcriptomics to physiology under climate change, and demonstrates that mechanisms underlying transgenerational effects persist across multiple generations with specific outcomes depending on acclimation type and environmental mismatch between generations.

Data compilation on the biological response to ocean acidification: environmental and experimental context of data sets and related literature

The exponential growth of studies on the biological response to ocean acidification over the last few decades has generated a large amount of data. To facilitate data comparison, a data compilation hosted at the data publisher PANGAEA was initiated in 2008 and is updated on a regular basis (doi:10.1594/PANGAEA.149999). By January 2015, a total of 581 data sets (over 4 000 000 data points) from 539 papers had been archived. Here we present the developments of this data compilation five years since its first description by Nisumaa et al. (2010). Most of study sites from which data archived are still in the Northern Hemisphere and the number of archived data from studies from the Southern Hemisphere and polar oceans are still relatively low. Data from 60 studies that investigated the response of a mix of organisms or natural communities were all added after 2010, indicating a welcomed shift from the study of individual organisms to communities and ecosystems. The initial imbalance of considerably more data archived on calcification and primary production than on other processes has improved. There is also a clear tendency towards more data archived from multifactorial studies after 2010. For easier and more effective access to ocean acidification data, the ocean acidification community is strongly encouraged to contribute to the data archiving effort, and help develop standard vocabularies describing the variables and define best practices for archiving ocean acidification data.

Nitrogen fixation and nitrogen gene abundance of METEOR cruises M77/3 and M77/4

Nitrogen fixation, the biological reduction of dinitrogen gas (N2) to ammonium (NH4+), is quantitatively the most important external source of new nitrogen (N) to the open ocean. Classically, the ecological niche of oceanic N2 fixers (diazotrophs) is ascribed to tropical oligotrophic surface waters, often depleted in fixed N, with a diazotrophic community dominated by cyanobacteria. Although this applies for large areas of the ocean, biogeochemical models and phylogenetic studies suggest that the oceanic diazotrophic niche may be much broader than previously considered, resulting in major implications for the global N-budget. Here, we report on the composition, distribution and abundance of nifH, the functional gene marker for N2 fixation. Our results show the presence of eight clades of diazotrophs in the oxygen minimum zone (OMZ) off Peru. Although proteobacterial clades dominated overall, two clusters affiliated to spirochaeta and archaea were identified. N2 fixation was detected within OMZ waters and was stimulated by the addition of organic carbon sources supporting the view that non-phototrophic diazotrophs were actively fixing dinitrogen. The observed co-occurrence of key functional genes for N2 fixation, nitrification, anammox and denitrification suggests that a close spatial coupling of N-input and N-loss processes exists in the OMZ off Peru. The wide distribution of diazotrophs throughout the water column adds to the emerging view that the habitat of marine diazotrophs can be extended to low oxygen/high nitrate areas. Furthermore, our statistical analysis suggests that NO2- and PO43- are the major factors affecting diazotrophic distribution throughout the OMZ. In view of the predicted increase in ocean deoxygenation resulting from global warming, our findings indicate that the importance of OMZs as niches for N2 fixation may increase in the futur

Elevated CO2 affects predator-prey interactions through altered performance

Recent research has shown that exposure to elevated carbon dioxide (CO2) affects how fishes perceive their environment, affecting behavioral and cognitive processes leading to increased prey mortality. However, it is unclear if increased mortality results from changes in the dynamics of predator-prey interactions or due to prey increasing activity levels. Here we demonstrate that ocean pCO2 projected to occur by 2100 significantly effects the interactions of a predator-prey pair of common reef fish: the planktivorous damselfish Pomacentrus amboinensis and the piscivorous dottyback Pseudochromis fuscus. Prey exposed to elevated CO2 (880 µatm) or a present-day control (440 µatm) interacted with similarly exposed predators in a cross-factored design. Predators had the lowest capture success when exposed to elevated CO2 and interacting with prey exposed to present-day CO2. Prey exposed to elevated CO2 had reduced escape distances and longer reaction distances compared to prey exposed to present-day CO2 conditions, but this was dependent on whether the prey was paired with a CO2 exposed predator or not. This suggests that the dynamics of predator-prey interactions under future CO2 environments will depend on the extent to which the interacting species are affected and can adapt to the adverse effects of elevated CO2.