Countergradient Variation in Reptiles: Thermal Sensitivity of Developmental and Metabolic Rates Across Locally Adapted Populations

Latitudinal Clines in an Ectothermic Vertebrate: Patterns in Body Size, Growth Rate, and Reproductive Effort Suggest Countergradient Responses in the Prairie Lizard

Although we have evidence that many organisms are exhibiting declines in body size in response to climate warming, we have little knowledge of underlying mechanisms or how associated phenotypic suites may coevolve. The better we understand coadaptations among physiology, morphology, and life history, the more accurate our predictions will be of organismal response to changing thermal environments. This is especially salient for ectotherms because they comprise 99% of species worldwide and are key to functioning ecosystems. Here, we assess body size, growth rates, and reproductive traits of a vertebrate ectotherm, the prairie lizard, Sceloporus consobrinus, for multiple populations along a latitudinal thermal gradient and compare body size clines between S. consobrinus and eastern fence lizard (S. undulatus) populations. We found that phenotypic values increased as environmental temperatures decreased for all traits examined, resulting in a pattern representative of countergradient variation. The positive covariation of phenotypes across the thermal gradient exemplifies the enigma of "master of all traits." This enigma was further illustrated by the energy expenditure toward growth and reproduction increasing as phenotypic values increased. The evolutionary responses appear to reveal overcompensation because annual energy expenditure toward growth and reproduction increased even as activity periods decreased. Overall, compensatory responses to cooler thermal environments were exhibited by prairie lizards in body size, growth rate, egg size, and clutch size, resulting in cold-adapted populations allocating more energy toward maintenance, growth, and reproduction than lower latitude, warm-adapted populations. If larger body size in ectotherms is a result of intrinsically faster physiological rates compensating for cooler temperatures and shorter activity periods, then smaller body sizes in warmer environments may be a result of greater reliance on available environmental temperatures for physiological rates and time for assimilating resources.

Vertical and diel niches modulate thermal selection by rainforest frogs

Thermoregulatory behaviour determines an organism's body temperature and therefore its physiological condition, and may differ for organisms situated across climate gradients. Species' preferred or selected temperatures may be higher in warmer locations-referred to as coadaptation-or lower in warmer temperatures-countergradient variation. Here, we tested if rainforest amphibians exhibited coadaptation or countergradient thermal selection across an underappreciated spatial climate gradient (vertical height from forest floor to canopy) and separating diel activity (diurnal versus nocturnal behaviour). We captured 2534 amphibians over 216 ground-to-canopy surveys, and conducted 282 thermal selection assays for 37 species while pairing microclimate measurements and mechanistic model predictions to understand vertical and daily thermal variation in the field. Amphibians exhibited countergradient thermal selection: species occupying cool nocturnal conditions in canopies selected warmer temperatures than species occupying hot diurnal conditions at the forest floor. Furthermore, amphibians selected warmer temperatures than the average conditions that they were exposed to when active, and this divergence was especially high for nocturnal arboreal species (8.68°C). This suggests that rainforest amphibians dramatically underfill the warm end of their thermal niches, a trend across local thermal gradients that reflects recent findings across elevational and latitudinal gradients. We show that considering multidimensional climate gradients is important to evaluate thermoregulatory behaviour, and its evolutionary underpinnings, for understanding species' niches and community assembly.

Metabolic Compensation Associated With Digestion in Response to the Latitudinal Thermal Environment Across Populations of the Prairie Lizard (Sceloporus consobrinus)

Environmental temperatures directly affect physiological rates in ectotherms by constraining the possible body temperatures they can achieve, with physiological processes slowing as temperatures decrease and accelerating as temperatures increase. As environmental constraints increase, as they do northward along the latitudinal thermal gradient, organisms must adapt to compensate for the slower physiological processes or decreased opportunity time. Evolving faster general metabolic rates is one adaptive response posited by the metabolic cold adaptation (MCA) hypothesis. Here we test the MCA hypothesis by examining metabolism of prairie lizard populations across the latitudinal thermal gradient. Our results show that populations from cooler environments have higher standard metabolic rates (SMRs), but these are explained by associated larger body sizes. However, metabolic rates of fed, postprandial individuals (MRFed) and metabolic energy allocated to digestion (MRΔ) were highest in the population from the coldest environment after accounting for the effect of body size. Our results suggest cold-adapted populations compensate for lower temperatures and shorter activity periods by increasing metabolic rates associated with physiological processes and thus support the MCA hypothesis. When examining energy expenditure, metabolic rates of individuals in a postprandial state (MRFed) may be more ecologically relevant than those in a postabsorptive state (SMR) and give a better picture of energy use in ectotherm populations.

Intraspecific variation in thermal performance curves for early development in Fundulus heteroclitus

Thermal performance curves (TPCs) provide a framework for understanding the effects of temperature on ectotherm performance and fitness. TPCs are often used to test hypotheses regarding local adaptation to temperature or to develop predictions for how organisms will respond to climate warming. However, for aquatic organisms such as fishes, most TPCs have been estimated for adult life stages, and little is known about the shape of TPCs or the potential for thermal adaptation at sensitive embryonic life stages. To examine how latitudinal gradients shape TPCs at early life stages in fishes, we used two populations of Fundulus heteroclitus that have been shown to exhibit latitudinal variation along the thermal cline as adults. We exposed embryos from both northern and southern populations and their reciprocal crosses to eight different temperatures (15°C, 18°C, 21°C, 24°C, 27°C, 30°C, 33°C, and 36°C) until hatch and examined the effects of developmental temperature on embryonic and larval traits (shape of TPCs, heart rate, and body size). We found that the pure southern embryos had a right-shifted TPC (higher thermal optimum (Topt) for developmental rate, survival, and embryonic growth rate) whereas pure northern embryos had a vertically shifted TPC (higher maximum performance (Pmax) for developmental rate). Differences across larval traits and cross-type were also found, such that northern crosses hatched faster and hatched at a smaller size compared to the pure southern population. Overall, these observed differences in embryonic and larval traits are consistent with patterns of both local adaptation and countergradient variation.

Morphological variations and demographic responses of the Mediterranean pond turtle Mauremys leprosa to heterogeneous aquatic habitats

Human activities affect terrestrial and aquatic habitats leading to changes at both individual and population levels in wild animal species. In this study, we investigated the phenotype and demographics of the Mediterranean pond turtle Mauremys leprosa (Schweigger, 1812) in contrasted environments of Southern France: two peri-urban rivers receiving effluents from wastewater treatment plants (WWTP), and another one without sewage treatment plant. Our findings revealed the presence of pesticides and pharmaceuticals in the three rivers of investigation, the highest diversities and concentrations of pollutants being found in the river subsections impacted by WWTP effluents. Principal component analysis and hierarchical clustering identified three levels of habitat quality, with different pollutant concentrations, thermal conditions, nutrient, and organic matter levels. The highest turtle densities, growth rates, and body sizes were estimated in the most disturbed habitats, suggesting potential adult benefits derived from harsh environmental conditions induced by pollution and eutrophication. Conversely, juveniles were the most abundant in the least polluted habitats, suggesting adverse effects of pollution on juvenile survival or adult reproduction. This study suggests that turtles living in polluted habitats may benefit from enhanced growth and body size, at the expense of reproductive success.

Geographic variation in incubation temperatures promoting viable offspring production in broadly co-distributed turtles

Organisms whose early life stages are environmentally sensitive produce offspring within a relatively narrow range of suitable abiotic conditions. In reptiles, development rate and survival are often maximized if incubation temperatures remain under 31°C, though this upper bound may vary within and among species. We addressed this expectation by comparing responses to egg incubation at 30°C versus 33°C in congeneric turtle species pairs with broad syntopic geographic distributions. In the two softshell turtles (Apalone spp.), the greatest changes in development rate and phenotypic variance were observed in the northernmost population, which had a low survival rate (40%) at 33°C. The presumably suboptimal temperature (33°C) for northern populations otherwise yielded 76%-93% survival rates and fast swimming speeds in more southern populations. Still, in one species, northern hatchlings incubated at 33°C matched the elevated speeds of their southern counterparts, revealing a countergradient response. In northern populations of the two map turtles (Graptemys spp.), survival was also reduced (28%-60%) at 33°C and the development rate (relative to 30°C) increased by up to 75%. Our experiments on divergent taxa with similar nesting ecologies substantiate that the optimal thermal range for offspring production is variable. These findings encourage further work on how population- and species-level differences relate to local adaptation in widely distributed oviparous species.

The influence of incubation temperature on offspring traits varies across northern and southern populations of the American alligator (Alligator mississippiensis)

Maternal provisioning and the developmental environment are fundamental determinants of offspring traits, particularly in oviparous species. However, the extent to which embryonic responses to these factors differ across populations to drive phenotypic variation is not well understood. Here, we examine the contributions of maternal provisioning and incubation temperature to hatchling morphological and metabolic traits across four populations of the American alligator (Alligator mississippiensis), encompassing a large portion of the species' latitudinal range. Our results show that whereas the influence of egg mass is generally consistent across populations, responses to incubation temperature show population-level variation in several traits, including mass, head length, head width, and residual yolk mass. Additionally, the influence of incubation temperature on developmental rate is greater at northern populations, while the allocation of maternal resources toward fat body mass is greater at southern populations. Overall, our results suggest that responses to incubation temperature, relative to maternal provisioning, are a larger source of interpopulation phenotypic variation and may contribute to the local adaptation of populations.

Population divergence in maternal investment and embryo energy use and allocation suggests adaptive responses to cool climates

The thermal sensitivity of early life stages can play a fundamental role in constraining species distributions. For egg-laying ectotherms, cool temperatures often extend development time and exacerbate developmental energy cost. Despite these costs, egg laying is still observed at high latitudes and altitudes. How embryos overcome the developmental constraints posed by cool climates is crucial knowledge for explaining the persistence of oviparous species in such environments and for understanding thermal adaptation more broadly. Here, we studied maternal investment and embryo energy use and allocation in wall lizards spanning altitudinal regions, as potential mechanisms that enable successful development to hatching in cool climates. Specifically, we compared population-level differences in (1) investment from mothers (egg mass, embryo retention and thyroid yolk hormone concentration), (2) embryo energy expenditure during development, and (3) embryo energy allocation from yolk towards tissue. We found evidence that energy expenditure was greater under cool compared with warm incubation temperatures. Females from relatively cool regions did not compensate for this energetic cost of development by producing larger eggs or increasing thyroid hormone concentration in yolk. Instead, embryos from the high-altitude region used less energy to complete development, that is, they developed faster without a concomitant increase in metabolic rate, compared with those from the low-altitude region. Embryos from high altitudes also allocated relatively more energy towards tissue production, hatching with lower residual yolk: tissue ratios than low-altitude region embryos. These results are consistent with local adaptation to cool climate and suggest that this is underpinned by mechanisms that regulate embryonic utilisation of yolk reserves and its allocation towards tissue, rather than shifts in maternal investment of yolk content or composition.

Climate aridity and habitat drive geographical variation in morphology and thermo-hydroregulation strategies of a widespread lizard species

Thermo-hydroregulation strategies involve concurrent changes in functional traits related to energy, water balance and thermoregulation and play a key role in determining life-history traits and population demography of terrestrial ectotherms. Local thermal and hydric conditions should be important drivers of the geographical variation of thermo-hydroregulation strategies, but we lack studies that examine these changes across climatic gradients in different habitat types. Here, we investigated intraspecific variation of morphology and thermo-hydroregulation traits in the widespread European common lizard (Zootoca vivipara louislantzi) across a multidimensional environmental gradient involving independent variation in air temperature and rainfall and differences in habitat features (access to free-standing water and forest cover). We sampled adult males for morphology, resting metabolic rate, total and cutaneous evaporative water loss and thermal preferences in 15 populations from the rear to the leading edge of the distribution across an elevational gradient ranging from sea level to 1750 m. Besides a decrease in adult body size with increasing environmental temperatures, we found little effect of thermal conditions on thermo-hydroregulation strategies. In particular, relict lowland populations from the warm rear edge showed no specific ecophysiological adaptations. Instead, body mass, body condition and resting metabolic rate were positively associated with a rainfall gradient, while forest cover and water access in the habitat throughout the season also influenced cutaneous evaporative water loss. Our study emphasizes the importance of rainfall and habitat features rather than thermal conditions for geographical variation in lizard morphology and physiology.

Population genetic differentiation and genomic signatures of adaptation to climate in an abundant lizard

Species distributed across climatic gradients will typically experience spatial variation in selection, but gene flow can prevent such selection from causing population genetic differentiation and local adaptation. Here, we studied genomic variation of 415 individuals across 34 populations of the common wall lizard (Podarcis muralis) in central Italy. This species is highly abundant throughout this region and populations belong to a single genetic lineage, yet there is extensive phenotypic variation across climatic regimes. We used redundancy analysis to, first, quantify the effect of climate and geography on population genomic variation in this region and, second, to test if climate consistently sorts specific alleles across the landscape. Climate explained 5% of the population genomic variation across the landscape, about half of which was collinear with geography. Linear models and redundancy analyses identified loci that were significantly differentiated across climatic regimes. These loci were distributed across the genome and physically associated with genes putatively involved in thermal tolerance, regulation of temperature-dependent metabolism and reproductive activity, and body colouration. Together, these findings suggest that climate can exercise sufficient selection in lizards to promote genetic differentiation across the landscape in spite of high gene flow.

Thermal-metabolic phenotypes of the lizard Podarcis muralis differ across elevation, but converge in high-elevation hypoxia

In response to a warming climate, many montane species are shifting upslope to track the emergence of preferred temperatures. Characterizing patterns of variation in metabolic, physiological and thermal traits along an elevational gradient, and the plastic potential of these traits, is necessary to understand current and future responses to abiotic constraints at high elevations, including limited oxygen availability. We performed a transplant experiment with the upslope-colonizing common wall lizard (Podarcis muralis) in which we measured nine aspects of thermal physiology and aerobic capacity in lizards from replicate low- (400 m above sea level, ASL) and high-elevation (1700 m ASL) populations. We first measured traits at their elevation of origin and then transplanted half of each group to extreme high elevation (2900 m ASL; above the current elevational range limit of this species), where oxygen availability is reduced by ∼25% relative to sea level. After 3 weeks of acclimation, we again measured these traits in both the transplanted and control groups. The multivariate thermal-metabolic phenotypes of lizards originating from different elevations differed clearly when measured at the elevation of origin. For example, high-elevation lizards are more heat tolerant than their low-elevation counterparts (counter-gradient variation). Yet, these phenotypes converged after exposure to reduced oxygen availability at extreme high elevation, suggesting limited plastic responses under this novel constraint. Our results suggest that high-elevation populations are well suited to their oxygen environments, but that plasticity in the thermal-metabolic phenotype does not pre-adapt these populations to colonize more hypoxic environments at higher elevations.

Latitudinal pattern of the thermal sensitivity of running speed in the endemic lizard Liolaemus multimaculatus

Physiological performance in lizards may be affected by climate across latitudinal or altitudinal gradients. In the coastal dune barriers in central-eastern Argentina, the annual maximum environmental temperature decreases up to 2°C from low to high latitudes, while the mean relative humidity of the air decreases from 50% to 25%. Liolaemus multimaculatus, a lizard in the family Liolaemidae, is restricted to these coastal dunes. We investigated the locomotor performance of the species at 6 different sites distributed throughout its range in these dune barriers. We inquired whether locomotor performance metrics were sensitive to the thermal regime attributable to latitude. The thermal performance breadth increased from 7% to 82% with latitude, due to a decrease in its critical thermal minimum of up to 5°C at higher latitudes. Lizards from high latitude sites showed a thermal optimum, that is, the body temperature at which maximum speed is achieved, up to 4°C lower than that of lizards from the low latitude. At relatively low temperatures, the maximum running speed of high-latitude individuals was faster than that of low-latitude ones. Thermal parameters of locomotor performance were labile, decreasing as a function of latitude. These results show populations of L. multimaculatus adjust thermal physiology to cope with local climatic variations. This suggests that thermal sensitivity responds to the magnitude of latitudinal fluctuations in environmental temperature.