Over a decade of field physiology reveals life-history specific strategies to drought in garter snakes (Thamnophis legans).

Changing climates and severe weather events can affect population viability. Individuals need to buffer such negative fitness consequences through physiological plasticity. Whether certain life-history strategies are more conducive to surviving changing climates is unknown, but theory predicts that strategies prioritizing maintenance and survival over current reproduction should be better able to withstand such change. We tested this hypothesis in a meta-population of garter snakes having naturally occurring variation in life-history strategies. We tested whether slow pace-of-life (POL) animals, that prioritize survival over reproduction, are more resilient than fast POL animals as measured by several physiological biomarkers. From 2006 to 2019, which included two multi-year droughts, baseline and stress-induced reactivity of plasma corticosterone and glucose varied annually with directionalities consistent with life-history theory. Slow POL animals exhibited higher baseline corticosterone and lower baseline glucose, relative to fast POL animals. These patterns were also observed in stress-induced measures; thus, reactivity was equivalent between ecotypes. However, in drought years, measures of corticosterone did not differ between different life histories. Immune cell distribution showed annual variation independent of drought or life history. These persistent physiological patterns form a backdrop to several extirpations of fast POL populations, suggesting a limited physiological toolkit to surviving periods of extreme drought.

Multiple Paternity in Garter Snakes With Evolutionarily Divergent Life Histories.

Many animal species exhibit multiple paternity, defined as multiple males genetically contributing to a single female reproductive event, such as a clutch or litter. Although this phenomenon is well documented across a broad range of taxa, the underlying causes and consequences remain poorly understood. For example, it is unclear how multiple paternity correlates with life-history strategies. Furthermore, males and females may differ in mating strategies and these patterns may shift with ecological context and life-history variation. Here, we take advantage of natural life-history variation in garter snakes (Thamnophis elegans) to address these questions in a robust field setting where populations have diverged along a slow-to-fast life-history continuum. We determine both female (observed) and male (using molecular markers) reproductive success in replicate populations of 2 life-history strategies. We find that despite dramatic differences in annual female reproductive output: 1) females of both life-history ecotypes average 1.5 sires per litter and equivalent proportions of multiply-sired litters, whereas 2) males from the slow-living ecotype experience greater reproductive skew and greater variance in reproductive success relative to males from the fast-living ecotype males despite having equivalent average reproductive success. Together, these results indicate strong intrasexual competition among males, particularly in the fast-paced life-history ecotype. We discuss these results in the context of competing hypotheses for multiple paternity related to population density, resource variability, and life-history strategy.

Hormonal and metabolic responses to upper temperature extremes in divergent life-history ecotypes of a garter snake.

Extreme temperatures constrain organismal physiology and impose both acute and chronic effects. Additionally, temperature-induced hormone-mediated stress response pathways and energetic trade-offs are important drivers of life-history variation. This study employs an integrative approach to quantify acute physiological responses to high temperatures in divergent life-history ecotypes of the western terrestrial garter snake (Thamnophis elegans). Using wild-caught animals, we measured oxygen consumption rate and physiological markers of hormonal stress response, energy availability and anaerobic respiration in blood plasma across five ecologically relevant temperatures (24, 28, 32, 35 and 38°C; 3 h exposure). Corticosterone, insulin and glucose concentrations all increased with temperature, but with different thermal response curves, suggesting that high temperatures differently affect energy-regulation pathways. Additionally, oxygen consumption rate increased without plateau and lactate concentration did not increase with temperature, challenging the recent hypothesis that oxygen limitation sets upper thermal tolerance limits. Finally, animals had similar physiological thermal responses to high-temperature exposure regardless of genetic background, suggesting that local adaptation has not resulted in fixed differences between ecotypes. Together, these results identify some of the mechanisms by which higher temperatures alter hormonal-mediated energy balance in reptiles and potential limits to the flexibility of this response.

Insulin-like signaling (IIS) responses to temperature, genetic background, and growth variation in garter snakes with divergent life histories.

The insulin/insulin-like signaling pathway (IIS) has been shown to mediate life history trade-offs in mammalian model organisms, but the function of this pathway in wild and non-mammalian organisms is understudied. Populations of western terrestrial garter snakes (Thamnophis elegans) around Eagle Lake, California, have evolved variation in growth and maturation rates, mortality senescence rates, and annual reproductive output that partition into two ecotypes: "fast-living" and "slow-living". Thus, genes associated with the IIS network are good candidates for investigating the mechanisms underlying ecological divergence in this system. We reared neonates from each ecotype for 1.5years under two thermal treatments. We then used qPCR to compare mRNA expression levels in three tissue types (brain, liver, skeletal muscle) for four genes (igf1, igf2, igf1r, igf2r), and we used radioimmunoassay to measure plasma IGF-1 and IGF-2 protein levels. Our results show that, in contrast to most mammalian model systems, igf2 mRNA and protein levels exceed those of igf1 and suggest an important role for igf2 in postnatal growth in reptiles. Thermal rearing treatment and recent growth had greater impacts on IGF levels than genetic background (i.e., ecotype), and the two ecotypes responded similarly. This suggests that observed ecotypic differences in field measures of IGFs may more strongly reflect plastic responses in different environments than evolutionary divergence. Future analyses of additional components of the IIS pathway and sequence divergence between the ecotypes will further illuminate how environmental and genetic factors influence the endocrine system and its role in mediating life history trade-offs.

Physiological indices of stress in wild and captive garter snakes: correlations, repeatability, and ecological variation.

Glucocorticoids and leukocyte ratios have become the most widespread variables employed to test hypotheses regarding physiological stress in wild and captive vertebrates. Little is known, however, regarding how these two indices of stress covary in response to stressors, their repeatability within individuals, and differences in response time upon capture. Furthermore, few studies compare stress indices between captive and wild populations, to assess potential alteration of stress physiology in captivity. To address these issues, we examined corticosterone (CORT) and heterophil to lymphocyte (H:L) ratios in two ecotypes of the garter snake Thamnophis elegans. We found that CORT and H:L ratios were not correlated within individuals, and both variables showed little or no repeatability over a period of months. CORT levels, but not H:L ratios, were higher for individuals sampled after 10min from the time of capture. However, both variables showed similar patterns of ecotypic variation, and both increased over time in gravid females maintained in captivity for four months. We suggest that CORT and H:L ratios are both useful, but disparate indices of stress in this species, and may show complex relationships to each other and to ecological and anthropogenic variables.

Temperature-dependence of metabolism and fuel selection from cells to whole organisms.

Temperature affects nearly every aspect of how organisms interact with and are constrained by their environment. Measures of organismal energetics, such as metabolic rate, are highly temperature-dependent and governed through temperature effects on rates of biochemical reactions. Characterizing the relationships among levels of biological organization can lend insight into how temperature affects whole-organism function. We tested the temperature dependence of cellular oxygen consumption and its relationship to whole-animal metabolic rate in garter snakes (Thamnophis elegans). Additionally, we tested whether thermal responses were linked to shifts in the fuel source oxidized to support metabolism with the use of carbon stable isotopes. Our results demonstrate temperature dependence of metabolic rates across levels of biological organization. Cellular (basal, adenosine triphosphate-linked) and whole-animal rates of respiration increased with temperature but were not correlated within or among individuals, suggesting that variation in whole-animal metabolic rates is not due simply to variation at the cellular level, but rather other interacting factors across scales of biological organization. Counter to trends observed during fasting, elevated temperature did not alter fuel selection (i.e., natural-abundance stable carbon isotope composition in breath, δ13 Cbreath ). This consistency suggests the maintenance and oxidation of a single fuel source supporting metabolism across a broad range of metabolic demands.

Long-term monitoring of two snake species reveals immune-endocrine interactions and the importance of ecological context.

While there is huge promise in monitoring physiological parameters in free-living organisms, we also find high amounts of variability over time and space. This variation requires us to capitalize on long-term physiological monitoring to adequately address questions of population health, conservation status, or evolutionary trends as long-term sampling can examine ecoimmunological and endocrine interactions in wild populations while accounting for the variation that often makes ecophysiological field studies difficult to compare. In this study, we tested how immune efficacy and endocrinology interact while accounting for ecological context and environmental conditions in two snake species. Specifically, we measured bacterial killing ability, steroid hormones, and morphological characteristics in multiple populations of the Western Terrestrial Gartersnake (Thamnophis elegans) and Common Gartersnake (T. sirtalis) for multiple seasons over 6 years. Leveraging this long-term dataset, we tested how a broad immune measure and endocrine endpoints interact while accounting for individual traits, sampling date, and environmental conditions. Across both species, we found bacterial killing ability to be directly related to corticosterone (CORT) and temperature and greater overall in the spring compared to the fall. We found CORT and testosterone yielded relationships with individual sex, sampling temperature, and time of year. Wild populations can exhibit high amounts of variation in commonly collected physiological endpoints, highlighting the complexity and difficulty inherent in interpreting single endpoints without taking ecological and environmental conditions into account. Our study emphasizes the importance of reporting the environmental conditions under which the sampling occurred to allow for better contextualization and comparison between studies.

Mitochondria as central characters in a complex narrative: Linking genomics, energetics, pace-of-life, and aging in natural populations of garter snakes.

As a pacesetter for physiological processes, variation in metabolic rate can determine the shape of energetic trade-offs and thereby drive variation in life-history traits. In turn, such variation in metabolic performance and life-histories can have profound consequences for lifespan and lifetime fitness. Thus, the extent to which metabolic rate variation is due to phenotypic plasticity or fixed genetic differences among individuals or populations is likely to be shaped by natural selection. Here, we first present a generalized framework describing the central role of mitochondria in processes linking environmental, genomic, physiological, and aging variation. We then present a test of these relationships in an exemplary system: populations of garter snakes (Thamnophis elegans) exhibiting contrasting life-history strategies - fast-growing, early-reproducing, and fast-aging (FA) versus slow-growing, late-reproducing, and slow-aging (SA). Previous work has characterized divergences in mitochondrial function, reactive oxygen species processing, and whole-organism metabolic rate between these contrasting life-history ecotypes. Here, we report new data on cellular respiration and mitochondrial genomics and synthesize these results with previous work. We test hypotheses about the causes and implications of mitochondrial genome variation within this generalized framework. First, we demonstrate that snakes of the FA ecotype increase cellular metabolic rate across their lifespan, while the opposite pattern holds for SA snakes, implying that reduced energetic throughput is associated with a longer life. Second, we show that variants in mitochondrial genomes are segregating across the landscape in a manner suggesting selection on the physiological consequences of this variation in habitats varying in temperature, food availability, and rates of predation. Third, we demonstrate functional variation in whole-organism metabolic rate related to these mitochondrial genome sequence variants. With this synthesis of numerous datasets, we are able to further characterize how variation across levels of biological organization interact within this generalized framework and how this has resulted in the emergence of distinct life-history ecotypes that vary in their rates of aging and lifespan.

Integrating behaviour into the pace-of-life continuum: Divergent levels of activity and information gathering in fast- and slow-living snakes.

An animal's life history, physiology, and behaviour can be shaped by selection in a manner that favours strong associations among these aspects of an integrated phenotype. Recent work combining animal personality and life-history theory proposes that animals with faster life-history strategies (i.e., fast growth, high annual reproductive rate, short lifespan) should exhibit higher general activity levels relative to those with slower life-history strategies, but empirical tests of within-species variation in these traits are lacking. In garter snakes from ecotypes which are known to differ in ecology, life-history strategy, and physiology, we tested for differences in tongue-flick rate as a measure of information gathering and movement patterns as a measure of general activity. Tongue flicks and movement were strongly positively correlated and both behaviours were repeatable across trials. Snakes from the fast-living ecotype were more active and showed evidence of habituation. The slow-living ecotype maintained low levels of activity throughout the trials. We propose that environmental factors, such as high predation, experienced by the fast-living ecotype select for both increased information-gathering and activity levels to facilitate efficient responses to repeated challenges. Thus, we offer evidence that behaviour is an important component of co-evolved suites of traits forming a general pace-of-life continuum in this system.

Mitochondrial divergence between slow- and fast-aging garter snakes.

Mitochondrial function has long been hypothesized to be intimately involved in aging processes--either directly through declining efficiency of mitochondrial respiration and ATP production with advancing age, or indirectly, e.g., through increased mitochondrial production of damaging free radicals with age. Yet we lack a comprehensive understanding of the evolution of mitochondrial genotypes and phenotypes across diverse animal models, particularly in species that have extremely labile physiology. Here, we measure mitochondrial genome-types and transcription in ecotypes of garter snakes (Thamnophis elegans) that are adapted to disparate habitats and have diverged in aging rates and lifespans despite residing in close proximity. Using two RNA-seq datasets, we (1) reconstruct the garter snake mitochondrial genome sequence and bioinformatically identify regulatory elements, (2) test for divergence of mitochondrial gene expression between the ecotypes and in response to heat stress, and (3) test for sequence divergence in mitochondrial protein-coding regions in these slow-aging (SA) and fast-aging (FA) naturally occurring ecotypes. At the nucleotide sequence level, we confirmed two (duplicated) mitochondrial control regions one of which contains a glucocorticoid response element (GRE). Gene expression of protein-coding genes was higher in FA snakes relative to SA snakes for most genes, but was neither affected by heat stress nor an interaction between heat stress and ecotype. SA and FA ecotypes had unique mitochondrial haplotypes with amino acid substitutions in both CYTB and ND5. The CYTB amino acid change (Isoleucine → Threonine) was highly segregated between ecotypes. This divergence of mitochondrial haplotypes between SA and FA snakes contrasts with nuclear gene-flow estimates, but correlates with previously reported divergence in mitochondrial function (mitochondrial oxygen consumption, ATP production, and reactive oxygen species consequences).