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.

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.

Predicting changes in abundance and juvenile proportions of earthworms in the field using glmmADMB in R

The study aimed at investigating effects of three differently acting biocides; the insecticide esfenvalerate, the fungicide picoxystrobin and the bactericide triclosan, applied individually and as a mixture, on an earthworm community in the field. A concentration-response design was chosen and results were analyzed using univariate and multivariate approaches. Effects on juvenile proportions were less pronounced and more variable than effects on abundance, but effects in general were species- and chemical-specific, and temporal variations distinct. Esfenvalerate and picoxystrobin appeared to elicit stronger effects than triclosan at laboratory-based ECx values, which is in accordance with our previous laboratory study on Eisenia fetida. The mixture affected abundance and juvenile proportions, but the latter only at high mixture concentrations. Esfenvalerate and picoxystrobin appeared to be the main drivers for the mixture's toxicity. Species-specific toxicity patterns question the reliability of mixture toxicity predictions derived on E. fetida for field earthworms. Biocide concentrations equaling EC50s (reproduction) for E. fetida provoked effects on the field earthworms mainly exceeding 50%, indicating effect intensification from the laboratory to field as well as the influence of indirect effects produced by species interactions. The differing results of the present field study and the previous laboratory study imply that lower- and higher-tier studies may not be mutually exclusive, but to be used in complementary.

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.

Sedimentation and mineral analysis on sediment cores from the Lena Delta

Core and outcrop analysis from Lena mouth deposits have been used to reconstruct the Late Quaternary sedimentation history of the Lena Delta. Sediment properties (heavy mineral composition, grain size characteristics, organic carbon content) and age determinations (14C AMS and IR-OSL) are applied to discriminate the main sedimentary units of the three major geomorphic terraces, which form the delta. The development of the terraces is controlled by complex interactions among the following four factors: (1) Channel migration. According to the distribution of 14C and IR-OSL age determinations of Lena mouth sediments, the major river runoff direction shifted from the west during marine isotope stages 5-3 (third terrace deposits) towards the northwest during marine isotope stage 2 and transition to stage 1 (second terrace), to the northeast and east during the Holocene (first terrace deposits). (2) Eustasy. Sea level rise from Last Glacial lowstand to the modern sea level position, reached at 6-5 ka BP, resulted in back-filling and flooding of the palaeovalleys. (3) Neotectonics. The extension of the Arctic Mid-Ocean Ridge into the Laptev Sea shelf acted as a halfgraben, showing dilatation movements with different subsidence rates. From the continent side, differential neotectonics with uplift and transpression in the Siberian coast ridges are active. Both likely have influenced river behavior by providing sites for preservation, with uplift, in particular, allowing accumulation of deposits in the second terrace in the western sector. The actual delta setting comprises only the eastern sector of the Lena Delta. (4) Peat formation. Polygenetic formation of ice-rich peaty sand (''Ice Complex'') was most extensive (7-11 m in thickness) in the southern part of the delta area between 43 and 14 ka BP (third terrace deposits). In recent times, alluvial peat (5-6 m in thickness) is accumulated on top of the deltaic sequences in the eastern sector (first terrace).

C and H geochemistry of altered oceanic crus of ODP/IODP Hole 1256D

Carbon and hydrogen concentrations and isotopic compositions were measured in 19 samples from altered oceanic crust cored in ODP/IODP Hole 1256D through lavas, dikes down to the gabbroic rocks. Bulk water content varies from 0.32 to 2.14 wt% with dD values from -64per mil to -25per mil. All samples are enriched in water relative to fresh basalts. The dD values are interpreted in terms of mixing between magmatic water and another source that can be either secondary hydrous minerals and/or H contained in organic compounds such as hydrocarbons. Total CO2, extracted by step-heating technique, ranges between 564 and 2823 ppm with d13C values from -14.9per mil to -26.6per mil. As for water, these altered samples are enriched in carbon relative to fresh basalts. The carbon isotope compositions are interpreted in terms of a mixing between two components: (1) a carbonate with d13C = -4.5per mil and (2) an organic compound with d13C = -26.6per mil. A mixing model calculation indicates that, for most samples (17 of 19), more than 75% of the total C occurs as organic compounds while carbonates represent less than 25%. This result is also supported by independent estimates of carbonate content from CO2 yield after H3PO4 attack. A comparison between the carbon concentration in our samples, seawater DIC (Dissolved Inorganic Carbon) and DOC (Dissolved Organic Carbon), and hydrothermal fluids suggests that CO2 degassed from magmatic reservoirs is the main source of organic C addition to the crust during the alteration process. A reduction step of dissolved CO2 is thus required, and can be either biologically mediated or not. Abiotic processes are necessary for the deeper part of the crust (>1000 mbsf) because alteration temperatures are greater than any hyperthermophilic living organism (i.e. T > 110 °C). Even if not required, we cannot rule out the contribution of microbial activity in the low-temperature alteration zones. We propose a two-step model for carbon cycling during crustal alteration: (1) when "fresh" oceanic crust forms at or close to ridge axis, alteration starts with hot hydrothermal fluids enriched in magmatic CO2, leading to the formation of organic compounds during Fischer-Tropsch-type reactions; (2) when the crust moves away from the ridge axis, these interactions with hot hydrothermal fluids decrease and are replaced by seawater interactions with carbonate precipitation in fractures. Taking into account this organic carbon, we estimate C isotope composition of mean altered oceanic crust at ? -4.7per mil, similar to the d13C of the C degassed from the mantle at ridge axis, and discuss the global carbon budget. The total flux of C stored in the altered oceanic crust, as carbonate and organic compound, is 2.9 ± 0.4 * 10**12 molC/yr.

Diatom and hydrocarbon analysis of sediment core NBP00-03-JPC38

In this study, we present a unique high-resolution Holocene record of oceanographic and climatic change based on analyses of diatom assemblages combined with biomarker data from a sediment core collected from the Vega Drift, eastern Antarctic Peninsula (EAP). These data add to the climate framework already established by high-resolution marine sedimentary records from the Palmer Deep, western Antarctic Peninsula (WAP). Heavy sea ice conditions and reduced primary productivity were observed prior to 7.4 ka B.P. in relation with the proximity of the glacial ice melt and calving. Subsequent Holocene oceanographic conditions were controlled by the interactions between the Westerlies-Antarctic Circumpolar Current (ACC)-Weddell Gyre dynamics. A warm period characterized by short seasonal sea ice duration associated with a southern shift of both ACC and Westerlies field persisted until 5 ka B.P. This warm episode was then followed by climate deterioration during the middle-to-late Holocene (5 to 1.9 ka B.P.) with a gradual increase in annual sea ice duration triggered by the expansion of the Weddell Gyre and a strong oceanic connection from the EAP to the WAP. Increase of benthic diatom species during this period was indicative of more summer/autumn storms, which was consistent with changes in synoptic atmospheric circulation and the establishment of low- to high-latitude teleconnections. Finally, the multicentennial scale variability of the Weddell Gyre intensity and storm frequency during the late Holocene appeared to be associated with the increased El Niño-Southern Oscillation frequency.

Dissolution index of Globigerina bulloides in recent and Last Glacial Maximum sediments

The modern Atlantic Ocean, dominated by the interactions of North Atlantic Deep Water (NADW) and Antarctic Bottom Water (AABW), plays a key role in redistributing heat from the Southern to the Northern Hemisphere. In order to reconstruct the evolution of the relative importance of these two water masses, the NADW/AABW transition, reflected by the calcite lysocline, was investigated by the Globigerina bulloides dissolution index (BDX?). The depth level of the Late Glacial Maximum (LGM) calcite lysocline was elevated by several hundred metres, indicating a more corrosive water mass present at modern NADW level. Overall, the small range of BDX? data and the gradual decrease in preservation below the calcite lysocline point to a less stratified Atlantic Ocean during the LGM. Similar preservation patterns in the West and East Atlantic demonstrate that the modern west-east asymmetry did not exist due to an expansion of southern deep waters compensating for the decrease in NADW formation.

Abundance of organic-walled dinoflagellate cysts of sediment core GeoB3603-2

Here, we reconstruct the varying influence of the Agulhas Current (AgC), the South Atlantic Current and the Antarctic Circumpolar Current on the Agulhas Retroflection (AgR) in the eastern South Atlantic Ocean for the last 160,000 years on the basis of the dinoflagellate cysts, pollen and spores present in a sediment core (GeoB 3603-2) from the southeastern Cape Basin offshore South Africa, where the Agulhas Current enters the Atlantic Ocean. Our analyses reveal strong orbital forcing on the heat exchange between the Indian Ocean and the South Atlantic Ocean during the Late Quaternary. Maxima in local productivity appear to be primarily related to a strengthening of the ocean circulation as a result of the high seasonal contrast during precession maxima. During precession minima, seasonal contrast was low and stratified, oligo- to mesotrophic conditions prevailed, notably when these minima coincided with the glacial terminations. The clear presence of periodicities on a sub-Milankovitch scale as well as modulations of the primary frequencies demonstrate that the Agulhas Retroflection furthermore is modulated substantially by complex interactions between the subtropical, 'precession-driven', climate and ocean circulation systems, and the southern, 'obliquity-driven', high latitudes.