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.

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.

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.

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.