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.

(Table 1) Organic geochemistry of ODP Hole 134-832B and Site 134-833

Ocean Drilling Program (ODP) Leg 134 was located in the central part of the New Hebrides Island Arc, in the Southwest Pacific. Here the d'Entrecasteaux Zone of ridges, the North d'Entrecasteaux Ridge and South d'Entrecasteaux Chain, is colliding with the arc. The region has a Neogene history of subduction polarity reversal, ridge-arc collision, and back-arc spreading. The reasons for drilling in this region included the following: (1) to determine the differences in the style and time scale of deformation associated with the two ridge-like features (a fairly continuous ridge and an irregularly topographic seamount chain) that are colliding with the central New Hebrides Island Arc; (2) to document the evolution of the magmatic arc in relation to the collision process and possible Neogene reversal of subduction; and (3) to understand the process of dewatering of a small accretionary wedge associated with ridge collision and subduction. Seven sites were occupied during the leg, five (Sites 827-831) were located in the d'Entrecasteaux Zone where collision is active. Three sites (Sites 827, 828, and 829) were located where the North d'Entrecasteaux Ridge is colliding, whereas two sites (Sites 830 and 831) were located in the South d'Entrecasteaux Chain collision zone. Sites 828 (on North d'Entrecasteaux Ridge) and 831 (on Bougainville Guyot) were located on the Pacific Plate, whereas all other sites were located on a microplate of the North Fiji Basin. Two sites (Sites 832 and 831) were located in the intra-arc North Aoba Basin. Results of Leg 134 drilling showed that forearc deformation associated with the North d'Entrecasteaux Ridge and South d'Entrecasteaux Chain collision is distinct and different. The d'Entrecasteaux Zone is an Eocene subduction/obduction complex with a distinct submerged island arc. Collision and subduction of the North d'Entrecasteaux Ridge results in off scraping of ridge material and plating of the forearc with thrust sheets (flakes) as well as distinct forearc uplift. Some offscraped sedimentary rocks and surficial volcanic basement rocks of the North d'Entrecasteaux Ridge are being underplated to the New Hebrides Island forearc. In contrast, the South d'Entrecasteaux Chain is a serrated feature resulting in intermittent collision and subduction of seamounts. The collision of the Bougainville Guyot has indented the forearc and appears to be causing shortening through thrust faulting. In addition, we found that the Quaternary relative convergence rate between the New Hebrides Island Arc at the latitude of Espiritu Santo Island is as high as 14 to 16 cm/yr. The northward migration rate of the d'Entrecasteaux Zone was found the be ~2 to 4 cm/yr based on the newly determined Quaternary relative convergence rate. Using these rates we established the timing of initial d'Entrecasteaux Zone collision with the arc at ~3 Ma at the latitude of Epi Island and fixed the impact of the North d'Entrecasteaux Ridge upon Espiritu Santo Island at early Pleistocene (between 1.89 and 1.58 Ma). Dewatering is occurring in the North d'Entrecasteaux Ridge accretionary wedge, and the wedge is dryer than other previously studied accretionary wedges, such as Barbados. This could be the result of less sediment being subducted at the New Hebrides compared to the Barbados.

Dive depth and dive duration data of subadult and adult, male and female southern elephant seals from Marion Island between 2004 and 2008 with links to datasets

Sexual segregation in habitat use occurs in a number of animal species, including southern elephant seals, where differences in migration localities and dive behaviour between sexes have been recorded. Due to the extreme sexual size dimorphism exhibited by southern elephant seals, it is unclear whether observed differences in dive behaviour are due to increased physiological capacity of males, compared to females, or differences in activity budgets and foraging behaviour. Here we use a mixed-effects modelling approach to investigate the effects of sex, size, age and individual variation on a number of dive parameters measured on southern elephant seals from Marion Island. Although individual variation accounted for substantial portions of total model variance for many response variables, differences in maximum and targeted dive depths were always influenced by sex, and only partly by body length. Conversely, dive durations were always influenced by body length, while sex was not identified as a significant influence. These results support hypotheses that physiological capability associated with body size is a limiting factor on dive durations. However, differences in vertical depth use appear to be the result of differences in forage selection between sexes, rather than a by-product of the size dimorphism displayed by this species. This provides further support for resource partitioning and possible avoidance of inter-sexual competition in southern elephant seals.

(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.

(Appendix) Vertebral number of topsmelt, Atherinops affinis, adults and offspring from different latitudes along the North-American Pacific coast

Organisms that are distributed across spatial climate gradients often exhibit adaptive local variations in morphological and physiological traits, but to what extent such gradients shape evolutionary responses is still unclear. Given the strong natural contrast in latitudinal temperature gradients between the North-American Pacific and Atlantic coast, we asked how increases in vertebral number (VN, known as Jordan's Rule) with latitude would differ between Pacific (Atherinops affinis) and Atlantic Silversides (Menidia menidia), two ecologically equivalent and taxonomically similar fishes with similar latitudinal distributions. VN was determined from radiographs of wild-caught adults (genetic + environmental differences) and its genetic basis confirmed by rearing offspring in common garden experiments. Compared to published data on VN variation in M. menidia (a mean increase of 7.0 vertebrae from 32 to 46°N, VN slope = 0.42/lat), the latitudinal VN increase in Pacific Silversides was approximately half as strong (a mean increase of 3.3 vertebrae from 28 to 43°N, VN slope = 0.23/lat). This mimicked the strong Atlantic (1.11°C/lat) versus weak Pacific latitudinal gradient (0.40°C/lat) in median annual sea surface temperature (SST). Importantly, the relationship of VN to SST was not significantly different between the two species (average slope = -0.39 vertebrae/°C), thus suggesting a common thermal dependency of VN in silverside fishes. Our findings provide novel support for the hypothesis that temperature gradients are the ultimate cause of Jordan's Rule, even though its exact adaptive significance remains speculative. A second investigated trait, the mode of sex determination in Atlantic versus Pacific Silversides, revealed patterns that were inconsistent with our expectation: M. menidia displays temperature-dependent sex determination (TSD) at low latitudes, where growing seasons are long or unconstrained, but also a gradual shift to genetic sex determination (GSD) with increasing latitude due to more and more curtailed growing seasons. Sex ratios in A. affinis, on the other hand, were independent of latitude and rearing temperature (indicating GSD), even though growing seasons are thermally unconstrained across most of the geographical distribution of A. affinis. This suggests that additional factors (e.g., longevity) play an important role in shaping the mode of sex determination in silverside fishes.

Geochemistry of ODP Site 202-1233 and river sediments

Bulk sediment chemistry from three Chilean continental margin Ocean Drilling Program sites constrains regional continental erosion over the past 30,000 years. Sediments from thirteen rivers that drain the (mostly igneous) Andes and the (mostly metamorphic) Coast Range, along with existing rock chemistry datasets, define terrestrial provenance for the continental margin sediments. Andean river sediments have high Mg/Al relative to Coast-Range river sediments. Near 36°S, marine sediments have high-Mg/Al (i.e. more Andean) sources during the last glacial period, and lower-Mg/Al (less Andean) sources during the Holocene. Near 41°S a Ti-rich source, likely from coast-range igneous intrusions, is prevalent during Holocene time, whereas high-Mg/Al Andean sources are more prevalent during the last glacial period. We infer that there is a dominant ice-sheet control of sediment sources. At 36°S, Andean-sourced sediment decreased as Andean mountain glaciers retreated after ~17.6 ka, coincident with local oceanic warming and southward retreat of the Patagonian Forest and, by inference, westerly winds. At 41°S Andean sediment dominance peaks and then rapidly declines at ~19 ka, coincident with local oceanic warming and the earliest deglacial sea-level rise. We hypothesize that this decreased flux of Andean material in the south is related to rapid retreat of the marine-based portion of the Patagonian Ice Sheet in response to global sea-level rise, as the resulting flooding of the southern portion of the Central Valley created a sink for Andean sediments in this region. Reversal of the decreasing deglacial Mg/Al trend at 41°S from 14.5 to 13.0 ka is consistent with a brief re-advance of the Patagonian ice sheet coincident with the Antarctic Cold Reversal.

Chemistry of organic matter of glacial and Cretaceous sediments from the Prydz Bay

Organic matter in Miocene glacial sediments in Hole 739C on the Antarctic Shelf represents erosional recycled continental material. Various indications of maturity in bulk organic matter, kerogens, and extracts imply that an exposed section of mature organic carbon-rich material was present during the Miocene. Based on biomarker, n-alkane, and kerogen analysis, a massive diamictite of early Eocene/Oligocene age at Hole 739C contains immature organic matter. Visual and pyrolysis analyses of the kerogens suggest a predominance of terrestrial organic matter in all samples from Hole 739C. A reversal of thermal maturities, i.e., more-mature overlying less-mature sections, may be related to redeposition generated from glacial erosion. Siliciclastic fluviatile sediments of Lower Cretaceous age from Hole 741A were analyzed. The organic matter from this hole contains immature aliphatic and aromatic biomarkers as well as a suite of odd carbon number-dominated nalkanes. Visual examination and pyrolysis analysis of the kerogen suggests that predominantly immature terrestrial organic matter is present at Hole 741A. The similarities between Hole 739C Unit V and Hole 741A suggest that the source of the organic matter in the glacial sediments in Unit V at Hole 739C could be Cretaceous in age and similar to sediments sampled at Hole 741A in Prydz Bay.

Temperature and photoperiod effects on sex determination in Leuresthes tenuis (fish)

In some gonochoristic species, sex is influenced not only by genotype at conception but also by the environment that offspring experience during early ontogeny (termed environmental sex determination or ESD). ESD is thought to be adaptive when seasonal variations in environmental conditions provide a sex-specific fitness advantage. In vertebrates, temperature is the most common determinant of sex, and seasonal variation in temperature serves as a temporal cue of environmental quality such as length of the growing season. Some environments, however, lack strong seasonal temperature fluctuations and other cues, particularly photoperiod, may provide a more reliable indicator of the environment offspring enter. We tested this hypothesis by rearing the offspring of the California grunion (Leuresthes tenuis, Ayres), which experiences low seasonal temperature variation in nature, under common garden conditions at three temperature and two photoperiod treatments. Our experiments revealed that both temperature and photoperiod significantly affected sex ratios in L. tenuis. More females were produced at cooler temperatures and longer day lengths, which is consistent with female biased sex ratios early in the breeding season, and likely adaptive through increased female size and fecundity. To our knowledge, this is the first documented case of photoperiod-dependent sex determination in a gonochoristic vertebrate.

Sex-specific foraging during parental care in a size-monomorphic seabird

Sex differences in foraging behaviour are typically studied in size-dimorphic taxa. Data on sex-specific behavior in monomorphic taxa are needed to test theories of reproductive investment. It has been suggested that in seabirds foraging niche separation may be related to decreased intersexual competition for food between cooperating pair-bonded individuals. Alternatively, sex differences in foraging niches may be driven by different nutritional requirements of females associated with the reproductive costs of egg production and oviposition. To assess these possibilities, we studied a size-monomorphic colonial seabird, the Australasian Gannet (Morus serrator) at the Cape Kidnappers gannetry, New Zealand. We recorded maximum dive depths, and distinct diet composition of incubating females as indicated by stable isotopic signatures. Results suggested greater female foraging effort during early times of incubation, indicated by significantly deeper maximum dives. Sex-specific foraging patterns across other breeding stages were more variable. Nitrogen stable isotopic values showed that incubating females occupied a different trophic position compared to males at the same breeding stage, and also from those of gannets of both sexes at later stages of parental care. Overall, the data are consistent with cost-of-oviposition compensation in females necessitating male-bias in parental care in biparental breeders. Further research is needed to unravel the implications for the evolution of sex differences in behavior in this and other monomorphic taxa.

Stable isotope analysis of Porites coral core from the Andaman Islands

North Atlantic climate variations are reflected in sedimentary records from the northern Indian Ocean in which two basins, the Arabian Sea and the Bay of Bengal, are strongly affected by the monsoon. Contrary to the Bay of Bengal the Arabian Sea plays an important role in the global marine nitrogen cycle. In its mid-water oxygen minimum zone (OMZ) bioavailable fixed nitrogen is reduced to nitrogen gas (NO3- - > N2), whereas oxygen concentrations are slightly above the threshold of nitrate reduction in the OMZ of the Bay of Bengal. A coral colony (Porites lutea) growing south of Port Blair on the Andaman Islands in the Bay of Bengal was studied for its response to changes in the monsoon system and its link to temperature changes in the North Atlantic Ocean, between 1975 and 2006. Its linear extension rates, d13C and d18O values measured within the coral skeleton reveal a strong seasonality, which seems to be caused by the monsoon-driven reversal of the surface ocean circulation. The sampling site appears to be influenced by low salinity Bay of Bengal Water during the NE monsoon (boreal winter) and by the high salinity Arabian Sea Water during the SW monsoon in summer. The high salinity Arabian Sea Water circulates along with the Summer Monsoon Current (S-MC) from the Arabia Sea into the Bay of Bengal. Decreasing d18O and reconstructed salinity values correlate to the increasing SSTs in the North Atlantic Ocean indicating a reduced influence of the S-MC at the sampling site in the course of northern hemispheric warming. During such periods oxygen-depletion became stronger in the OMZ of the Arabian Sea as indicated by the sedimentary records. A reduced propagation of oxygen-depleted high salinity Arabian Sea Water into the Bay of Bengal could be a mechanism maintaining oxygen concentration above the threshold of nitrate reduction in the OMZ of the Bay of Bengal in times of global warming.

Lithium isotopic composition of pore fluids and sediments in the Costa Rica subduction zone

Pore fluid and sediment Li concentrations and isotopic ratios provide important insights on the hydrology, sediment contribution to the arc volcanoes and fluid-sediment reactions at the dominantly non-accretionary Costa Rica subduction zone. Ocean Drilling Program Site 1039 in the trench axis provides a reference section of 400 m of the incoming sediments, and Site 1040, situated arcward from the trench, consists of a deformed sedimentary wedge and apron sediments, the décollement, and the partially dewatered underthrust sediment section. At the reference site, pore fluids show important isotopic variations (delta6Li=-21.7 to -37.8 per mil), reflecting the interplay of in situ alteration of volcanic material and ion exchange with clay minerals. In the basal section, a reversal of Li concentration and delta6Li toward seawater values is observed, providing supporting evidence for a lateral seawater flow system in the upper oceanic basement underlying this sediment section. At Site 1040, pore fluid of the lower deformed wedge sediments and within the décollement is enriched in Li and the isotopic compositions are relatively light, suggesting infiltration of a deep-seated fluid. The delta6Li value of -22 per mil of this Li-enriched fluid (261 µM), when compared with the delta6Li value of the subducted sediment section (-11 per mil), suggests that the deep source fluid originates from mineral fluid dehydration and transformation reactions at temperatures of 100 to 150°C, consistent with the temperature range of the up-dip seismogenic zone and of transformation of smectite to illite. The distribution of Li and its isotopes in the underthrust section are similar to those at the reference site, indicating near complete subduction of the incoming sediments and that early dewatering of the underthrust sediments occurs predominantly by lateral flow into the ocean. The hemipelagic clay-rich sediment section of the subducting plate carries most of the Li into this subduction zone, and the pelagic diatomaceous and nannofossil calcareous oozes contain little Li. The Li isotopes of both the clay-rich hemipelagic sediments and of the pelagic oozes are, however, similar, with delta6Li values of -9 to -12 per mil. The observations that (1) the delta6Li values of the underthrust sediments are distinctly lower than that of the mantle, and (2) the lavas of the Costa Rican volcanoes are enriched in Li and 7Li, provide an approximation of the contribution of the subducted sediments to the arc volcanoes. A first order mass balance calculation suggests that approximately half of the Li flux delivered by subducted sediments and altered oceanic crust into the Middle American Trench is recycled to the Costa Rican arc and at most a quarter of sedimentary Li is returned into the ocean through thrust faults, primarily the décollement thrust.

Magnetic state of ocean basalts (Table 1)

The natural remanent magnetization (NRM) of ocean basalts, giving rise to the pattern of marine magnetic anomalies, is known to be of comparatively low intensity for about 20 Ma old oceanic crust. The aim of this study is to detect possible peculiarities in the rock magnetic properties of ocean basalts of this age, and to establish a link between magnetomineralogy, rock magnetic parameters, and the low NRM intensity. Ocean basalts covering ages from 0.7 to 135 Ma were selected for rock magnetic experiments and their room temperature hysteresis parameters, Curie temperature and temperature dependence of saturation magnetization MS(T) was determined and complemented by reflected light microscopy. The majority of samples is magnetically dominated by titanomagnetite and titanomaghemite with increasing oxidation state with age. For these, a strong dependence of hysteresis parameters on the age of the samples is found. The samples have a minimum in saturation magnetization and a maximum in magnetic stability in the age interval ranging from approximately 10 to 40 Ma, coinciding with the age interval of low NRM intensity. The observed change in saturation magnetization is in the same order as that for the NRM intensity. A further peculiarity of the titanomaghemites from this age interval is the shape of their MS(T) curves, which display a maximum above room temperature (Neel P-type) and, sometimes, a self-reversal of magnetization below room temperature (Neel N-type). These special rock magnetic properties can be explained by titanomagnetite low-temperature oxidation and highly oxidized titanomaghemites in the age interval 10-40 Ma. A corresponding measurement of the NRM at elevated temperature allows to identify a maximum in NRM intensity above room temperature for the samples in that age interval. This provides evidence that the NRM is equally carried by titanomaghemites and that the low NRM intensities for about 20 Ma old ocean basalts are caused consequently by the low saturation magnetization of these titanomaghemites.

(Table 1) Stable carbon and oxygen isotope ratios of benthic foraminifera of ODP Hole 184-1148A

Changes in intermediate and deep ocean circulation likely played a significant role in global carbon cycling and meridional heat/moisture transport during the middle Miocene climate transition (~14 Ma). High-resolution middle Miocene (16-13 Ma) benthic foraminifer stable isotope records from the South China Sea reveal a reorganization of regional bottom waters, which preceded the globally recognized middle Miocene ~1 per mil d18O increase (13.8 Ma) by 100,000 years. An observed reversal of the benthic foraminifera d13C gradient between ODP Sites 1146 (2092 m) and 1148 (3294 m; 13.9-13.5 Ma) is interpreted to reflect an increase in the southward flux of low d13C deep (> 2000 m) Pacific Ocean waters (Flower and Kennett, 1993, doi:10.1029/93PA02196; Shevenell and Kennett, 2004). Large-scale changes in Pacific intermediate and deep ocean circulation, coupled with enhanced global carbon cycling at the end of the Monterey Carbon Isotope excursion, likely acted as internal feedbacks to the Earth's climate system. These feedbacks reduced the sensitivity of Antarctica to lower latitude-derived heat/moisture and facilitated the transition of the Earth's climate system to a new, relatively stable glacial state.

Relative abundance of anaerobic methanotrophic archaea in bacterial mat Cascadia margin off Oregon

The microbially mediated anaerobic oxidation of methane (AOM) is the major biological sink of the greenhouse gas methane in marine sediments (doi:10.1007/978-94-009-0213-8_44) and serves as an important control for emission of methane into the hydrosphere. The AOM metabolic process is assumed to be a reversal of methanogenesis coupled to the reduction of sulfate to sulfide involving methanotrophic archaea (ANME) and sulfate-reducing bacteria (SRB) as syntrophic partners which were describes amongst others in Boetius et al. (2000; doi:10.1038/35036572). In this study, 16S rRNA-based methods were used to investigate the distribution and biomass of archaea in samples from sediments above outcropping methane hydrate at Hydrate Ridge (Cascadia margin off Oregon) and (ii) massive microbial mats enclosing carbonate reefs (Crimea area, Black Sea). Sediment samples from Hydrate Ridge were obtained during R/V SONNE cruises SO143-2 in August 1999 and SO148-1 in August 2000 at the crest of southern Hydrate Ridge at the Cascadia convergent margin off the coast of Oregon. The second study area is located in the Black Sea and represents a field in which there is active seepage of free gas on the slope of the northwestern Crimea area. Here, a field of conspicuous microbial reefs forming chimney-like structures was discovered at a water depth of 230 m in anoxic waters. The microbial mats were sampled by using the manned submersible JAGO during the R/V Prof. LOGACHEV cruise in July 2001. At Hydrate Ridge the surface sediments were dominated by aggregates consisting of ANME-2 and members of the Desulfosarcina-Desulfococcus branch (DSS) (ANME-2/DSS aggregates), which accounted for >90% of the total cell biomass. The numbers of ANME-1 cells increased strongly with depth; these cells accounted 1% of all single cells at the surface and more than 30% of all single cells (5% of the total cells) in 7- to 10-cm sediment horizons that were directly above layers of gas hydrate. In the Black Sea microbial mats ANME-1 accounted for about 50% of all cells. ANME-2/DSS aggregates occurred in microenvironments within the mat but accounted for only 1% of the total cells. FISH probes for the ANME-2a and ANME-2c subclusters were designed based on a comparative 16S rRNA analysis. In Hydrate Ridge sediments ANME-2a/DSS and ANME-2c/DSS aggregates differed significantly in morphology and abundance. The relative abundance values for these subgroups were remarkably different at Beggiatoa sites (80% ANME-2a, 20% ANME-2c) and Calyptogena sites (20% ANME-2a, 80% ANME-2c), indicating that there was preferential selection of the groups in the two habitats.