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.

Lake Baikal amphipods under climate change: Thermal constraints and ecological consequences, links to supplementary material

Lake Baikal, the world's most voluminous freshwater lake, has experienced unprecedented warming during the last decades. A uniquely diverse amphipod fauna inhabits the littoral zone and can serve as a model system to identify the role of thermal tolerance under climate change. This study aimed to identify sublethal thermal constraints in two of the most abundant endemic Baikal amphipods, Eulimnogammarus verrucosus and Eulimnogammarus cyaneus, and Gammarus lacustris, a ubiquitous gammarid of the Holarctic. As the latter is only found in some shallow isolated bays of the lake, we further addressed the question whether rising temperatures could promote the widespread invasion of this non-endemic species into the littoral zone. Animals were exposed to gradual temperature increases (4 week, 0.8 °C/d; 24 h, 1 °C/h) starting from the reported annual mean temperature of the Baikal littoral (6 °C). Within the framework of oxygen- and capacity-limited thermal tolerance (OCLTT), we used a nonlinear regression approach to determine the points at which the changing temperature-dependence of relevant physiological processes indicates the onset of limitation. Limitations in ventilation representing the first limits of thermal tolerance (pejus (= "getting worse") temperatures (Tp)) were recorded at 10.6 (95% confidence interval; 9.5, 11.7), 19.1 (17.9, 20.2), and 21.1 (19.8, 22.4) °C in E. verrucosus, E. cyaneus, and G. lacustris, respectively. Field observations revealed that E. verrucosus retreated from the upper littoral to deeper and cooler waters once its Tp was surpassed, identifying Tp as the ecological thermal boundary. Constraints in oxygen consumption at higher than critical temperatures (Tc) led to an exponential increase in mortality in all species. Exposure to short-term warming resulted in higher threshold values, consistent with a time dependence of thermal tolerance. In conclusion, species-specific limits to oxygen supply capacity are likely key in the onset of constraining (beyond pejus) and then life-threatening (beyond critical) conditions. Ecological consequences of these limits are mediated through behavioral plasticity in E. verrucosus. However, similar upper thermal limits in E. cyaneus (endemic, Baikal) and G. lacustris (ubiquitous, Holarctic) indicate that the potential invader G. lacustris would not necessarily benefit from rising temperatures. Secondary effects of increasing temperatures remain to be investigated.

Temperature tolerance of different larval stages of the spider crab Hyas araneus exposed to elevated seawater PCO2

Exposure to elevated seawater PCO2 limits the thermal tolerance of crustaceans but the underlying mechanisms have not been comprehensively explored. Larval stages of crustaceans are even more sensitive to environmental hypercapnia and possess narrower thermal windows than adults. In a mechanistic approach, we analysed the impact of high seawater CO2 on parameters at different levels of biological organization, from the molecular to the whole animal level. At the whole animal level we measured oxygen consumption, heart rate and activity during acute warming in zoea and megalopa larvae of the spider crab Hyas araneus exposed to different levels of seawater PCO2. Furthermore, the expression of genes responsible for acid-base regulation and mitochondrial energy metabolism, and cellular responses to thermal stress (e.g. the heat shock response) was analysed before and after larvae were heat shocked byrapidly raising the seawater temperature from 10°C rearing temperature to 20°C. Zoea larvae showed a high heat tolerance, which decreased at elevated seawater PCO2, while the already low heat tolerance of megalopa larvae was not limited further by hypercapnic exposure. There was a combined effect of elevated seawater CO2 and heat shock in zoea larvae causing elevated transcript levels of heat shock proteins. In all three larval stages, hypercapnic exposure elicited an up-regulation of genes involved in oxidative phosphorylation, which was, however, not accompanied by increased energetic demands. The combined effect of seawater CO2 and heat shock on the gene expression of heat shock proteins reflects the downward shift in thermal limits seen on the whole animal level and indicates an associated capacity to elicit passive thermal tolerance. The up-regulation of genes involved in oxidative phosphorylation might compensate for enzyme activities being lowered through bicarbonate inhibition and maintain larval standard metabolic rates at high seawater CO2 levels. The present study underlines the necessity to align transcriptomic data with physiological responses when addressing mechanisms affected by an interaction of elevated seawater PCO2 and temperature extremes.

Gene expression under thermal stress varies across a geographical range expansion front

Acknowledgements: We would like to thank Hanna Bensch and Hannes Weise for assistance with the collection of samples in the field. This work was supported by the Biodiversity and Ecosystem Services in a Changing Climate (BECC; a joint Lund-Gothenburg University initiative), the Swedish Research Council (EIS, BH), the Crafoord Foundation (EIS, BH), the Swedish Royal Society (EIS), ‘Gyllenstiernska Krapperupstiftelsen (EIS), the Wenner-Gren Foundations (postdoctoral stipend to RYD), EU FP7 (Marie Curie International Incoming Fellowship to RYD), the Kungliga Fysiografiska Sällskapet i Lund (MW) and the Helge Ax:son Johnson Stiftelse (MW). B.H. and E.I.S. conceived of the study. L.L. developed the hypotheses to be tested. L.L. and R.D. collected the field data and samples. All six authors contributed to planning RNA-seq analyses. P.C. and L.L. analysed the data. L.L. wrote the manuscript, which all six authors edited.

A Predictive Thermal Habitat Model for Harbor Seals in the Northwest Atlantic

We analyzed projections of current and future ambient temperatures along the eastern United States in relationship to the thermal tolerance of harbor seals in air. Using the earth systems model (HadGEM2-ES) and representative concentration pathways (RCPs) 4.5 and 8.5, which are indicative of two different atmospheric CO2 concentrations, we were able to examine possible shifts in distribution based on three metrics: current preferences, the thermal limit of juveniles, and the thermal limits of adults. Our analysis focused on average ambient temperatures because harbor seals are least effective at regulating their body temperature in air, making them most susceptible to rising air temperatures in the coming years. Our study focused on the months of May, June, and August from 2041-2060 (2050) and 2061-2080 (2070) as these are the historic months in which harbor seals are known to annually come ashore to pup, breed, and molt. May, June, and August are also some of the warmest months of the year. We found that breeding colonies along the eastern United States will be limited by the thermal tolerance of juvenile harbor seals in air, while their foraging range will extend as far south as the thermal tolerance of adult harbor seals in air. Our analysis revealed that in 2070, harbor seal pups should be absent from the United States coastline nearing the end of the summer due to exceptionally high air temperatures.

Intraspecific individual variation of temperature tolerance associated with oxygen demand in the European sea bass (Dicentrarchus labrax)

The European sea bass (Dicentrarchus labrax) is an economically important fish native to the Mediterranean and Northern Atlantic. Its complex life cycle involves many migrations through temperature gradients that affect the energetic demands of swimming. Previous studies have shown large intraspecific variation in swimming performance and temperature tolerance, which could include deleterious and advantageous traits under the evolutionary pressure of climate change. However, little is known of the underlying determinants of this individual variation. We investigated individual variation in temperature tolerance in 30 sea bass by exposing them to a warm temperature challenge test. The eight most temperature-tolerant and eight most temperature-sensitive fish were then studied further to determine maximal swimming speed (U-CAT), aerobic scope and post-exercise oxygen consumption. Finally, ventricular contractility in each group was determined using isometric muscle preparations. The temperature-tolerant fish showed lower resting oxygen consumption rates, possessed larger hearts and initially recovered from exhaustive exercise faster than the temperature-sensitive fish. Thus, whole-animal temperature tolerance was associated with important performance traits. However, the temperature-tolerant fish also demonstrated poorer maximal swimming capacity (i.e. lower UCAT) than their temperature-sensitive counterparts, which may indicate a trade-off between temperature tolerance and swimming performance. Interestingly, the larger relative ventricular mass of the temperature-tolerant fish did not equate to greater ventricular contractility, suggesting that larger stroke volumes, rather than greater contractile strength, may be associated with thermal tolerance in this species.

James's rule and causes and consequences of a latitudinal cline in the demography of John's Snapper (Lutjanus johnii) in coastal waters of Australia

Demographic parameters were derived from sectioned otoliths of John’s Snapper (Lutjanus johnii) from 4 regions across 9° of latitude and 23° of longitude in northern Australia. Latitudinal variation in size and growth rates of this species greatly exceeded longitudinal variation. Populations of John’s Snapper farthest from the equator had the largest body sizes, in line with James’s rule, and the fastest growth rates, contrary to the temperature-size rule for ectotherms. A maximum age of 28.6 years, nearly 3 times previous estimates, was recorded and the largest individual was 990 mm in fork length. Females grew to a larger mean asymptotic fork length (L∞) than did males, a finding consistent with functional gonochorism. Otolith weight at age and gonad weight at length followed the same latitudinal trends seen in length at age. Length at maturity was ~72–87% of L∞ and varied by ~23% across the full latitudinal gradient, but age at first maturity was consistently in the range of 6–10 years, indicating that basic growth trajectories were similar across vastly different environments. We discuss both the need for complementary reproductive data in age-based studies and the insights gained from experiments where the concept of oxygen- and capacity-limited thermal tolerance is applied to explain the mechanistic causes of James’s rule in tropical fish species.

Nutritional Control of L1 Arrest and Recovery in Caenorhabditis elegans by Insulin-like Peptides and Signaling

Animals must coordinate development with fluctuating nutrient availability. Nutrient availability governs post-embryonic development in Caenorhabditis elegans: larvae that hatch in the absence of food do not initiate post-embryonic development but enter "L1 arrest" (or "L1 diapause") and can survive starvation for weeks, while rapidly resume normal development once get fed. Insulin-like signaling (IIS) has been shown to be a key regulator of L1 arrest and recovery. However, the C. elegans genome encodes 40 insulin-like peptides (ILPs), and it is unknown which peptides participate in nutritional control of L1 arrest and recovery. Work in other contexts has identified putative receptor agonists and antagonists, but the extent of specificity versus redundancy is unclear beyond this distinction.

We measured mRNA expression dynamics with high temporal resolution for all 40 insulin-like genes during entry into and recovery from L1 arrest. Nutrient availability influences expression of the majority of insulin-like genes, with variable dynamics suggesting complex regulation. We identified 13 candidate agonists and 8 candidate antagonists based on expression in response to nutrient availability. We selected ten candidate agonists (daf-28, ins-3, ins-4, ins-5, ins-6, ins-7, ins-9, ins-26, ins-33 and ins-35) for further characterization in L1 stage larvae. We used destabilized reporter genes to determine spatial expression patterns. Expression of candidate agonists was largely overlapping in L1 stage larvae, suggesting a role of the intestine, chemosensory neurons ASI and ASJ, and the interneuron PVT in systemic control of L1 development. Transcriptional regulation of candidate agonists was most significant in the intestine, as if nutrient uptake was a more important influence on transcription than sensory perception. Scanning in the 5' upstream promoter region of these 40 ILPs, We found that transcription factor PQM-1 and GATA putative binding sites are depleted in the promoter region of antagonists. A novel motif was also found to be over-represented in ILPs.

Phenotypic analysis of single and compound deletion mutants did not reveal effects on L1 recovery/developmental dynamics, though simultaneous disruption of ins-4 and daf-28 extended survival of L1 arrest without enhancing thermal tolerance, while overexpression of ins-4, ins-6 or daf-28 shortened L1 survival. Simultaneous disruption of several ILPs showed a temperature independent, transient dauer phenotype. These results revealed the relative redundancy and specificity among agonistic ILPs.

TGF- β and steroid hormone (SH) signaling have been reported to control the dauer formation along with IIS. Our preliminary results suggest they may also mediate the IIS control of L1 arrest and recovery, as the expression of several key components of TGF-β and SH signaling pathway genes are negatively regulated by DAF-16, and loss-of-function of these genes partially represses daf-16 null phenotype in L1 arrest, and causes a retardation in L1 development.

In summary, my dissertation study focused on the IIS, characterized the dynamics and sites of ILPs expression in response to nutrient availability, revealed the function of specific agonistic ILPs in L1 arrest, and suggested potential cross-regulation among IIS, TGF-β signaling and SH signaling in controlling L1 arrest and recovery. These findings provide insights into how post-embryonic development is governed by insulin-like signaling and nutrient availability.

Prey landscapes help identify potential foraging habitats for leatherback turtles in the NE Atlantic

Identifying key marine megavertebrate habitats has become ever more important as concern increases regarding global fisheries bycatch and accelerated climate change. This will be aided by a greater understanding of the patterns and processes determining the spatiotemporal distribution of species of conservation concern. We identify probable foraging grounds for leatherback turtles in the NE Atlantic using monthly landscapes of gelatinous organism distribution constructed from Continuous Plankton Recorder Survey data. Using sightings data (n = 2013 records, 1954 to 2003) from 9 countries (UK, Ireland, France, Belgium, The Netherlands, Denmark, Germany, Norway and Sweden), we show sea surface temperatures of approximately 10 to 12 degree C most likely indicate the lower thermal threshold for accessible habitats during seasonal foraging migrations to high latitudes. Integrating maps of gelatinous plankton as a possible indicator of prey distribution with thermal tolerance parameters demonstrates the dynamic (spatial and temporal) nature of NE Atlantic foraging habitats. We highlight the importance of body size- related thermal constraints in structuring leatherback foraging populations and demonstrate a latitudinal gradient in body size (Bergmann's rule) where smaller animals are excluded from higher latitude foraging areas. We highlight the marine area of the European continental shelf edge as being both thermally accessible and prey rich, and therefore potentially supporting appreciable densities of foraging leatherbacks, with some suitable areas not yet extensively surveyed.

Disentangling the impacts of heat wave magnitude, duration and timing on the structure and diversity of sessile marine assemblages

Extreme climatic events, including heat waves (HWs) and severe storms, influence the structure of marine and terrestrial ecosystems. Despite growing consensus that anthropogenic climate change will increase the frequency, duration and magnitude of extreme events, current understanding of their impact on communities and ecosystems is limited. Here, we used sessile invertebrates on settlement panels as model assemblages to examine the influence of HW magnitude, duration and timing on marine biodiversity patterns. Settlement panels were deployed in a marina in southwest UK for ≥5 weeks, to allow sufficient time for colonisation and development of sessile fauna, before being subjected to simulated HWs in a mesocosm facility. Replicate panel assemblages were held at ambient sea temperature (∼17 °C), or +3 °C or +5 °C for a period of 1 or 2 weeks, before being returned to the marina for a recovery phase of 2–3 weeks. The 10-week experiment was repeated 3 times, staggered throughout summer, to examine the influence of HW timing on community impacts. Contrary to our expectations, the warming events had no clear, consistent impacts on the abundance of species or the structure of sessile assemblages. With the exception of 1 high-magnitude long-duration HW event, warming did not alter not assemblage structure, favour non-native species, nor lead to changes in richness, abundance or biomass of sessile faunal assemblages. The observed lack of effect may have been caused by a combination of (1) the use of relatively low magnitude, realistic heat wave treatments compared to previous studies (2), the greater resilience of mature adult sessile fauna compared to recruits and juveniles, and (3) the high thermal tolerance of the model organisms (i.e., temperate fouling species, principally bryozoans and ascidians). Our study demonstrates the importance of using realistic treatments when manipulating climate change variables, and also suggests that biogeographical context may influence community-level responses to short-term warming events, which are predicted to increase in severity in the future.

One size fits all: stability of metabolic scaling under warming and ocean acidification in echinoderms

Responses by marine species to ocean acidification (OA) have recently been shown to be modulated by external factors including temperature, food supply and salinity. However the role of a fundamental biological parameter relevant to all organisms, that of body size, in governing responses to multiple stressors has been almost entirely overlooked. Recent consensus suggests allometric scaling of metabolism with body size differs between species, the commonly cited 'universal' mass scaling exponent (b) of A3/4 representing an average of exponents that naturally vary. One model, the Metabolic-Level Boundaries hypothesis, provides a testable prediction: that b will decrease within species under increasing temperature. However, no previous studies have examined how metabolic scaling may be directly affected by OA. We acclimated a wide body-mass range of three common NE Atlantic echinoderms (the sea star Asterias rubens, the brittlestars Ophiothrix fragilis and Amphiura filiformis) to two levels of pCO(2) and three temperatures, and metabolic rates were determined using closed-chamber respirometry. The results show that contrary to some models these echinoderm species possess a notable degree of stability in metabolic scaling under different abiotic conditions; the mass scaling exponent (b) varied in value between species, but not within species under different conditions. Additionally, we found no effect of OA on metabolic rates in any species. These data suggest responses to abiotic stressors are not modulated by body size in these species, as reflected in the stability of the metabolic scaling relationship. Such equivalence in response across ontogenetic size ranges has important implications for the stability of ecological food webs.

Clams from Tagus estuary : microbiological, physiological and chemical responses to depuration, transpot and environmental stress

Tese de doutoramento, Biologia (Biologia Marinha e Aquacultuta), Universidade de Lisboa, Faculdade de Ciências, 2014

Thermoaclimatory variations in the microenvironment of hemoglobin in the rainbow trout, Salmo gairdneri and the carp, Cyprinus carpio

Several inorganic substances (e.g., C£ , Mg , Ca , H ) are potent negative modulators of hemoglobin-oxygen affinity. To evaluate the possibility that potentially adaptive changes in the red cell ionic environment of hemoglobin may take place during acclimation of fishes to increased environmental temperature, hematological status (hemoglobin, hematocrit, red cell numbers, mean erythrocytic volume and hemoglobin content), plasma + + 2+ 2+ and packed red cell electrolyte levels (Na , K , Ca , Mg , C£ ) were evaluated in summer and winter populations of the stenothermal rainbow trout, Salmo gairdneri, following acclimation to 2°, 10°, 18°C, and in a spring population of eurythermal carp, Cyprinus carpio, held at 2°, 16° and 30°C. From these data cell ion concentrations and ion:hemoglobin ratios were estimated. In view of the role of red cell carbonic anhydrase in the reductions of blood C02 tensions and the recruitment of Na and C£~ lost by fishes, a preliminary investigation of thermoacclimatory changes in the activity of this system in rainbow trout erythrocytes was conducted. Few changes in hematological status were encountered following acclimation. There was, however, some evidence of weight-specific differential hematological response in carp. This lead to markedly greater increases in hemoglobin, hematocrit and red cell numbers in smaller rather than in larger specimens at higher temperatures; variations which were 2+ well correlated with changes in plasma Ca . Plasma composition in summer trout was not altered by acclimation. In winter trout plasma Na and K increased at higher temperatures. Carp were characterized by increases in plasma calcium, and reductions in sodium and magnesium under these conditions. Several significant seasonal differences in plasma ion levels were observed in the trout. (n) In trout, only erythrocytic K and K :Hb were altered by acclimation, rising at higher temperatures. In carp Na , Na :Hb, C£~ and C£~:Hb in- 2+ 2+ creased with temperature, while Mg and Mg :Hb declined. Changes in overall ionic composition in carp red cells were consistent with increases in H content. In both species significant reciprocal variations in C£~ 2+ - + and Mg were found. In mammalian systems increases in C£ and H reduce hemoglobin-oxygen affinity by interaction with hemoglobin. Reduction in 2+ 2+ Mg maximizes organophosphate modulator availability by decreasing ATP»Mg complex formation. Thus, the changes observed may be of adaptive value in reducing hemoglobin-oxygen affinity, and facilitating oxygen release to cells at higher temperatures. Trout appear to maintain a high chloridelow magnesium state over the entire thermal tolerance zone. Carp, however, achieved this state only at higher temperatures. In both species mean erythrocytic volume was decreased at higher temperatures and this may facilitate branchial oxygen loading. Since mean erythrocytic volume was inversely related to red cell ion content, it is hypothesized that reductions in cell volume are achieved by export of some unidentified solute or solutes. Variations in the carbonic anhydrase activity that could be attributed to the thermoacclimatory process were quite modest. On the other hand, assays performed at the temperature of acclimation showed a large temperature effect where under in vivo conditions of temperature fish acclimated to higher temperatures might be expected to have higher activities. Furthermore, since hematocrit increased with temperature in these fish, while carbonic anhydrase is present only in the erythrocyte, the whole blood levels of this enzyme are expected to increase and further augment the temperature effect. This, in turn, could aid in the reduction of C02 (111) tension and increase the production of H and HC0~~ used in the active uptake of Na and C£ at higher temperatures.

ORGANOSILICON BIOTECHNOLOGY: A BIO-INSPIRED APPROACH TO THE HYDROLYSIS OF ALKOXYSILANES and THE LIPASE-CATALYZED SYNTHESIS OF SILOXANE-CONTAINING POLYESTERS AND POLYAMIDES

The first part of this thesis studied the capacity of amino acids and enzymes to catalyze the hydrolysis and condensation of tetraethoxysilane and phenyltrimethoxysilane. Selected amino acids were shown to accelerate the hydrolysis and condensation of tetraethoxysilane under ambient temperature, pressure and at neutral pH (pH 7±0.02). The nature of the side chain of the amino acid was important in promoting hydrolysis and condensation. Several proteases were shown to have a capacity to hydrolyze tri- and tet-ra- alkoxysilanes under the same mild reaction conditions. The second part of this thesis employed an immobilized Candida antarctica lipase B (Novozym-435, N435) to produce siloxane-containing polyesters, polyamides, and polyester amides under solvent-free conditions. Enzymatic activity was shown to be temperature dependent, increasing until enzyme denaturation became the dominant pro-cess, which typically occurred between 120-130ᵒC. The residual activity of N435 was, on average, greater than 90%, when used in the synthesis of disiloxane-containing polyesters, regardless of the polymerization temperature except at the very highest temperatures, 140-150ᵒC. A study of the thermal tolerance of N435 determined that, over ten reaction cycles, there was a decrease in the initial rate of polymerization with each consecutive use of the catalyst. No change in the degree of monomer conversion after a 24 hour reaction cycle was found.