Lizards from warm and declining populations are born with extremely short telomeres.

Aging is the price to pay for acquiring and processing energy through cellular activity and life history productivity. Climate warming can exacerbate the inherent pace of aging, as illustrated by a faster erosion of protective telomere DNA sequences. This biomarker integrates individual pace of life and parental effects through the germline, but whether intra- and intergenerational telomere dynamics underlies population trends remains an open question. Here, we investigated the covariation between life history, telomere length (TL), and extinction risk among three age classes in a cold-adapted ectotherm (Zootoca vivipara) facing warming-induced extirpations in its distribution limits. TL followed the same threshold relationships with population extinction risk at birth, maturity, and adulthood, suggesting intergenerational accumulation of accelerated aging rate in declining populations. In dwindling populations, most neonates inherited already short telomeres, suggesting they were born physiologically old and unlikely to reach recruitment. At adulthood, TL further explained females' reproductive performance, switching from an index of individual quality in stable populations to a biomarker of reproductive costs in those close to extirpation. We compiled these results to propose the aging loop hypothesis and conceptualize how climate-driven telomere shortening in ectotherms may accumulate across generations and generate tipping points before local extirpation.

Population Genomics of Wall Lizards Reflects the Dynamic History of the Mediterranean Basin.

The Mediterranean Basin has experienced extensive change in geology and climate over the past six million years. Yet, the relative importance of key geological events for the distribution and genetic structure of the Mediterranean fauna remains poorly understood. Here, we use population genomic and phylogenomic analyses to establish the evolutionary history and genetic structure of common wall lizards (Podarcis muralis). This species is particularly informative because, in contrast to other Mediterranean lizards, it is widespread across the Iberian, Italian, and Balkan Peninsulas, and in extra-Mediterranean regions. We found strong support for six major lineages within P. muralis, which were largely discordant with the phylogenetic relationship of mitochondrial DNA. The most recent common ancestor of extant P. muralis was likely distributed in the Italian Peninsula, and experienced an "Out-of-Italy" expansion following the Messinian salinity crisis (∼5 Mya), resulting in the differentiation into the extant lineages on the Iberian, Italian, and Balkan Peninsulas. Introgression analysis revealed that both inter- and intraspecific gene flows have been pervasive throughout the evolutionary history of P. muralis. For example, the Southern Italy lineage has a hybrid origin, formed through admixture between the Central Italy lineage and an ancient lineage that was the sister to all other P. muralis. More recent genetic differentiation is associated with the onset of the Quaternary glaciations, which influenced population dynamics and genetic diversity of contemporary lineages. These results demonstrate the pervasive role of Mediterranean geology and climate for the evolutionary history and population genetic structure of extant species.

Impacts of temperature on O2 consumption of the Pyrenean brook newt (Calotriton asper) from populations along an elevational gradient.

Global warming impacts biodiversity worldwide, leading to species' adaptation, migration, or extinction. The population's persistence depends on the maintenance of essential activities, which is notably driven by phenotypic adaptation to local environments. Metabolic rate - that increases with temperature in ectotherms - is a key physiological proxy for the energy available to fuel individuals' activities. Cold-adapted ectotherms can exhibit a higher resting metabolism than warm-adapted ones to maintain functionality at higher elevations or latitudes, known as the metabolic cold-adaptation hypothesis. How climate change will affect metabolism in species inhabiting contrasting climates (cold or warm) is still a debate. Therefore, it is of high interest to assess the pace of metabolic responses to global warming among populations adapted to highly different baseline climatic conditions. Here, we conducted a physiological experiment in the endemic Pyrenean brook newt (Calotriton asper). We measured a proxy of standard metabolic rate (SMR) along a temperature gradient in individuals sampled among 6 populations located from 550 to 2189 m a.s.l. We demonstrated that SMR increased with temperature, but significantly diverged depending on populations' origins. The baseline and the slope of the relationship between SMR and temperature were both higher for high-elevation populations than for low-elevation populations. We discussed the stronger metabolic response observed in high-elevation populations suggesting a drop of performance in essential life activities for these individuals under current climate change. With the increase of metabolism as the climate warms, the metabolic-cold adaptation strategy selected in the past could compromise the sustainability of cold-adapted populations if short-term evolutionary responses do not allow to offset this evolutionary legacy.

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.

Vertical Transmission of a Nematode from Female Lizards to the Brains of Their Offspring.

Parasites have evolved a diversity of lifestyles that exploit the biology of their hosts. Some nematodes that parasitize mammals pass via the placenta or milk from one host to another. Similar cases of vertical transmission have never been reported in avian and nonavian reptiles, suggesting that egg laying may constrain the means of parasite transmission. However, here we report the first incidence of transovarial transmission of a previously undescribed nematode in an egg-laying amniote, the common wall lizard (Podarcis muralis). Nematodes enter the developing brain from the female ovary early in embryonic development. Infected lizard embryos develop normally and hatch with nematodes residing in their braincase. We present a morphological and molecular phylogenetic characterization of the nematode and suggest that particular features of lizard biology that are absent from birds and turtles facilitated the evolutionary origin of this novel life history.

Multiple glacial refugia and contemporary dispersal shape the genetic structure of an endemic amphibian from the Pyrenees.

Historical factors (colonization scenarios, demographic oscillations) and contemporary processes (population connectivity, current population size) largely contribute to shaping species' present-day genetic diversity and structure. In this study, we use a combination of mitochondrial and nuclear DNA markers to understand the role of Quaternary climatic oscillations and present-day gene flow dynamics in determining the genetic diversity and structure of the newt Calotriton asper (Al. Dugès, 1852), endemic to the Pyrenees. Mitochondrial DNA did not show a clear phylogeographic pattern and presented low levels of variation. In contrast, microsatellites revealed five major genetic lineages with admixture patterns at their boundaries. Approximate Bayesian computation analyses and linear models indicated that the five lineages likely underwent separate evolutionary histories and can be tracked back to distinct glacial refugia. Lineage differentiation started around the Last Glacial Maximum at three focal areas (western, central and eastern Pyrenees) and extended through the end of the Last Glacial Period in the central Pyrenees, where it led to the formation of two more lineages. Our data revealed no evidence of recent dispersal between lineages, whereas borders likely represent zones of secondary contact following expansion from multiple refugia. Finally, we did not find genetic evidence of sex-biased dispersal. This work highlights the importance of integrating past evolutionary processes and present-day gene flow and dispersal dynamics, together with multilocus approaches, to gain insights into what shaped the current genetic attributes of amphibians living in montane habitats.

Transplanting gravid lizards to high elevation alters maternal and embryonic oxygen physiology, but not reproductive success or hatchling phenotype.

Increased global temperatures have opened previously inhospitable habitats, such as at higher elevations. However, the reduction of oxygen partial pressure with increase in elevation represents an important physiological constraint that may limit colonization of such habitats, even if the thermal niche is appropriate. To test the mechanisms underlying the response to ecologically relevant levels of hypoxia, we performed a translocation experiment with the common wall lizard (Podarcis muralis), a widespread European lizard amenable to establishing populations outside its natural range. We investigated the impacts of hypoxia on the oxygen physiology and reproductive output of gravid common wall lizards and the subsequent development and morphology of their offspring. Lowland females transplanted to high elevations increased their haematocrit and haemoglobin concentration within days and maintained routine metabolism compared with lizards kept at native elevations. However, transplanted lizards suffered from increased reactive oxygen metabolite production near the oviposition date, suggesting a cost of reproduction at high elevation. Transplanted females and females native to different elevations did not differ in reproductive output (clutch size, egg mass, relative clutch mass or embryonic stage at oviposition) or in post-oviposition body condition. Developing embryos reduced heart rates and prolonged incubation times at high elevations within the native range and at extreme high elevations beyond the current range, but this reduced oxygen availability did not affect metabolic rate, hatching success or hatchling size. These results suggest that this opportunistic colonizer is capable of successfully responding to novel environmental constraints in these important life-history stages.

Variation of preferred body temperatures along an altitudinal gradient: A multi-species study.

Temperature affects the physiological functions of ectotherms. To maintain optimal body temperature and ensure physiological performance, these organisms can use behavioral adjustments to keep the body temperature in their specific temperature range, so-called preferred temperature (Tpref). It is therefore crucial to describe and understand how Tpref vary within and amongst populations to predict the effects of climate change of altitudinal range shifts in organisms. We aimed at determining the altitudinal variations in Tpref in three ectothermic species (the Pyrenean brook salamander - a semi-aquatic and thigmothermic amphibian - the European common lizard and the wall lizard - both heliothermic species). Using an experimental approach where Tpref were measured along a temperature gradient in laboratory conditions, we used a cross-sectional approach to compare the variation of Tpref measured in populations sampled along the altitudinal gradient in the Pyrenees. We hypothesized a complex and highly variable intra-specific response of Tpref along geographical clines, with a positive relationship between Tpref and altitude (as predicted by the countergradient variation), the reverse pattern (referring to the adaptation of local optima hypothesis), or no relationship at all. Our results corroborated partially the countergradient hypothesis in the salamander (middle to high elevation part). At high altitude level, individuals may compensate for lower opportunities of favorable conditions by choosing a high temperature which maximizes their activities. However, populations from low elevation level hence better supported the adaptation of local optima hypothesis, such as both lizard species, for which Tpref tended instead to decrease with altitude. Lizards from cold climates may be physiologically adapted to low temperature, with the ability to reach optimal functioning at lower temperature than individuals from low altitude. Our findings suggest that predicting future niche models should therefore take into account the variability of Tpref, where species might be differently affected by global warming.

Using GAMM to examine inter-individual heterogeneity in thermal performance curves for Natrix natrix indicates bet hedging strategy by mothers.

The thermal performance curve (TPC) illustrates the dependence on body- and therefore environmental- temperature of many fitness-related aspects of ectotherm ecology and biology including foraging, growth, predator avoidance, and reproduction. The typical thermal performance curve model is linear in its parameters despite the well-known, strong, non-linearity of the response of performance to temperature. In addition, it is usual to consider a single model based on few individuals as descriptive of a species-level response to temperature. To overcome these issues, we used generalized additive mixed modeling (GAMM) to estimate thermal performance curves for 73 individual hatchling Natrix natrix grass snakes from seven clutches, taking advantage of the structure of GAMM to demonstrate that almost 16% of the deviance in thermal performance curves is attributed to inter-individual variation, while only 1.3% is attributable to variation amongst clutches. GAMM allows precise estimation of curve characteristics, which we used to test hypotheses on tradeoffs thought to constrain the thermal performance curve: hotter is better, the specialist-generalist trade off, and resource allocation/acquisition. We observed a negative relationship between maximum performance and performance breadth, indicating a specialist-generalist tradeoff, and a positive relationship between thermal optimum and maximum performance, suggesting "hotter is better". There was a significant difference among matrilines in the relationship between Area Under the Curve and maximum performance - relationship that is an indicator of evenness in acquisition or allocation of resources. As we used unfed hatchlings, the observed matriline effect indicates divergent breeding strategies among mothers, with some mothers provisioning eggs unequally resulting in some offspring being better than others, while other mothers provisioned the eggs more evenly, resulting in even performance throughout the clutch. This observation is reminiscent of bet-hedging strategies, and implies the possibility for intra-clutch variability in the TPCs to buffer N. natrix against unpredictable environmental variability.

Loss of genetic diversity and increased embryonic mortality in non-native lizard populations.

Many populations are small and isolated with limited genetic variation and high risk of mating with close relatives. Inbreeding depression is suspected to contribute to extinction of wild populations, but the historical and demographic factors that contribute to reduced population viability are often difficult to tease apart. Replicated introduction events in non-native species can offer insights into this problem because they allow us to study how genetic variation and inbreeding depression are affected by demographic events (e.g. bottlenecks), genetic admixture and the extent and duration of isolation. Using detailed knowledge about the introduction history of 21 non-native populations of the wall lizard Podarcis muralis in England, we show greater loss of genetic diversity (estimated from microsatellite loci) in older populations and in populations from native regions of high diversity. Loss of genetic diversity was accompanied by higher embryonic mortality in non-native populations, suggesting that introduced populations are sufficiently inbred to jeopardize long-term viability. However, there was no statistical correlation between population-level genetic diversity and average embryonic mortality. Similarly, at the individual level, there was no correlation between female heterozygosity and clutch size, infertility or hatching success, or between embryo heterozygosity and mortality. We discuss these results in the context of human-mediated introductions and how the history of introductions can play a fundamental role in influencing individual and population fitness in non-native species.

Maternal exposure to predator scents: offspring phenotypic adjustment and dispersal.

Predation is a strong selective pressure generating morphological, physiological and behavioural responses in organisms. As predation risk is often higher during juvenile stages, antipredator defences expressed early in life are paramount to survival. Maternal effects are an efficient pathway to produce such defences. We investigated whether maternal exposure to predator cues during gestation affected juvenile morphology, behaviour and dispersal in common lizards (Zootoca vivipara). We exposed 21 gravid females to saurophagous snake cues for one month while 21 females remained unexposed (i.e. control). We measured body size, preferred temperature and activity level for each neonate, and released them into semi-natural enclosures connected to corridors in order to measure dispersal. Offspring from exposed mothers grew longer tails, selected lower temperatures and dispersed thrice more than offspring from unexposed mothers. Because both tail autotomy and altered thermoregulatory behaviour are common antipredator tactics in lizards, these results suggest that mothers adjusted offspring phenotype to risky natal environments (tail length) or increased risk avoidance (dispersal). Although maternal effects can be passive consequences of maternal stress, our results strongly militate for them to be an adaptive antipredator response that may increase offspring survival prospects.

Embryonic and juvenile snakes (Natrix maura, Linnaeus 1758) compensate for high elevation hypoxia via shifts in cardiovascular physiology and metabolism.

The colonization of novel environments requires a favorable response to conditions never, or rarely, encountered in recent evolutionary history. For example, populations colonizing upslope habitats must cope with lower atmospheric pressure at elevation, and thus reduced oxygen availability. The embryo stage in oviparous organisms is particularly susceptible, given its lack of mobility and limited gas exchange via diffusion through the eggshell and membranes. Especially little is known about responses of Lepidosaurian reptiles to reduced oxygen availability. To test the role of physiological plasticity during early development in response to high elevation hypoxia, we performed a transplant experiment with the viperine snake (Natrix maura, Linnaeus 1758). We maintained gravid females originating from low elevation populations (432 m above sea level [ASL]-normoxia) at both the elevation of origin and high elevation (2877 m ASL-extreme high elevation hypoxia; approximately 72% oxygen availability relative to sea level), then incubated egg clutches at both low and high elevation. Regardless of maternal exposure to hypoxia during gestation, embryos incubated at extreme high elevation exhibited altered developmental trajectories of cardiovascular function and metabolism across the incubation period, including a reduction in late-development egg mass. This physiological response may have contributed to the maintenance of similar incubation duration, hatching success, and hatchling body size compared to embryos incubated at low elevation. Nevertheless, after being maintained in hypoxia, juveniles exhibit reduced carbon dioxide production relative to oxygen consumption, suggesting altered energy pathways compared to juveniles maintained in normoxia. These findings highlight the role of physiological plasticity in maintaining rates of survival and fitness-relevant phenotypes in novel environments.

Body-size evolution on islands: are adult size variations in tiger snakes a nonadaptive consequence of selection on birth size?

Mean adult size has been used as the traditional measure of body size to explain trends of insular gigantism and dwarfism in a wide array of taxa. However, patterns of variation in body size at birth have received surprisingly little attention, leaving open the possibility that adult body-size differences are nonadaptive consequences of selection acting on neonate body size. Here I used an empirical and correlative approach to test this hypothesis in a mosaic of 12 island and mainland snake populations in Australia. Data collected on 597 adult and 1,084 neonate tiger snakes showed that (1) both adult and neonate mean body sizes varied strongly across populations; (2) prey diversity and size convincingly explained birth-size variations: birth size-notably, gape size-correlated with prey size; (3) neonate snout-vent length was significantly correlated with neonate gape size; and (4) neonate snout-vent length was significantly correlated with adult snout-vent length. Postnatal growth rates recorded under common-garden conditions differed across populations and were correlated with mean prey size. These data collectively suggest that (1) prey size is the main driver for the evolution of body size at birth in gape-limited predators, (2) adult size variations may reflect selective forces acting on earlier life stages, and (3) adult size variations may also reflect resource availability during ontogeny (notably, prey diversity).

Viability, behavior, and color expression in the offspring of matings between common wall lizard Podarcis muralis color morphs.

Color polymorphisms are widely studied to identify the mechanisms responsible for the origin and maintenance of phenotypic variability in nature. Two of the mechanisms of balancing selection currently thought to explain the long-term persistence of polymorphisms are the evolution of alternative phenotypic optima through correlational selection on suites of traits including color and heterosis. Both of these mechanisms can generate differences in offspring viability and fitness arising from different morph combinations. Here, we examined the effect of parental morph combination on fertilization success, embryonic viability, newborn quality, antipredator, and foraging behavior, as well as inter-annual survival by conducting controlled matings in a polymorphic lacertid Podarcis muralis, where color morphs are frequently assumed to reflect alternative phenotypic optima (e.g., alternative reproductive strategies). Juveniles were kept in outdoor tubs for a year in order to study inter-annual growth, survival, and morph inheritance. In agreement with a previous genome-wide association analysis, morph frequencies in the year-old juveniles matched the frequencies expected if orange and yellow expressions depended on recessive homozygosity at 2 separate loci. Our findings also agree with previous literature reporting higher reproductive output of heavy females and the higher overall viability of heavy newborn lizards, but we found no evidence for the existence of alternative breeding investment strategies in female morphs, or morph-combination effects on offspring viability and behavior. We conclude that inter-morph breeding remains entirely viable and genetic incompatibilities are of little significance for the maintenance of discrete color morphs in P. muralis from the Pyrenees.

Metal(loid) pollution, not urbanisation nor parasites predicts low body condition in a wetland bioindicator snake.

Urban ecosystems and remnant habitat 'islands' therein, provide important strongholds for many wildlife species including those of conservation significance. However, the persistence of these habitats can be undermined if their structure and function are too severely disrupted. Urban wetlands, specifically, are usually degraded by a monoculture of invasive vegetation, disrupted hydrology, and chronic-contamination from a suite of anthropogenic pollutants. Top predators-as bioindicators-can be used to assess and monitor the health of these ecosystems. We measured eight health parameters (e.g., parasites, wounds and scars, tail loss and body condition) in a wetland top predator, the western tiger snake, Notechis scutatus occidentalis. For three years, snakes were sampled across four wetlands along an urban gradient. For each site, we used GIS software to measure the area of different landscapes and calculate an urbanisation-landscape score. Previously published research on snake contamination informed our calculations of a metal-pollution index for each site. We used generalised linear mixed models to assess the relationship between all health parameters and site variables. We found the metal-pollution index to have the most significant association with poor body condition. Although parasitism, tail loss and wounds differed among sites, none of these parameters influenced body condition. Additionally, the suite of health parameters suggested differing health status among sites; however, our measure of contemporary landscape urbanisation was never a significant predictor variable. Our results suggest that the health of wetland predators surrounding a rapidly growing city may be offset by higher levels of environmental pollution.

Bioindicator snake shows genomic signatures of natural and anthropogenic barriers to gene flow.

Urbanisation alters landscapes, introduces wildlife to novel stressors, and fragments habitats into remnant 'islands'. Within these islands, isolated wildlife populations can experience genetic drift and subsequently suffer from inbreeding depression and reduced adaptive potential. The Western tiger snake (Notechis scutatus occidentalis) is a predator of wetlands in the Swan Coastal Plain, a unique bioregion that has suffered substantial degradation through the development of the city of Perth, Western Australia. Within the urban matrix, tiger snakes now only persist in a handful of wetlands where they are known to bioaccumulate a suite of contaminants, and have recently been suggested as a relevant bioindicator of ecosystem health. Here, we used genome-wide single nucleotide polymorphism (SNP) data to explore the contemporary population genomics of seven tiger snake populations across the urban matrix. Specifically, we used population genomic structure and diversity, effective population sizes (Ne), and heterozygosity-fitness correlations to assess fitness of each population with respect to urbanisation. We found that population genomic structure was strongest across the northern and southern sides of a major river system, with the northern cluster of populations exhibiting lower heterozygosities than the southern cluster, likely due to a lack of historical gene flow. We also observed an increasing signal of inbreeding and genetic drift with increasing geographic isolation due to urbanisation. Effective population sizes (Ne) at most sites were small (< 100), with Ne appearing to reflect the area of available habitat rather than the degree of adjacent urbanisation. This suggests that ecosystem management and restoration may be the best method to buffer the further loss of genetic diversity in urban wetlands. If tiger snake populations continue to decline in urban areas, our results provide a baseline measure of genomic diversity, as well as highlighting which 'islands' of habitat are most in need of management and protection.

PLASMA BIOCHEMISTRY PROFILES OF WILD WESTERN TIGER SNAKES (NOTECHIS SCUTATUS OCCIDENTALIS) BEFORE AND AFTER SIX MONTHS OF CAPTIVITY.

Urban wildlife often suffer poorer health than their counterparts living in more pristine environments due to exposure to anthropogenic stressors such as habitat degradation and environmental contamination. As a result, the health of urban versus nonurban snakes might be assessed by differences in their plasma biochemistries. We compared the plasma profiles of western tiger snakes (Notechis scutatus occidentalis) from a heavily urbanized wetland and a natural, nonurbanized wetland. Despite the urbanized snakes having lower body mass index, we found no significant difference between the plasma profiles of the two populations. We collected snakes from each population and kept them in captivity for 6 mo, providing them with stable conditions, uncontaminated (exempt from heavy metals and pesticides) food and water, and lowered parasite intensity in an attempt to promote better health through depuration. After captivity, snakes experienced a significant improvement in body mass index and significant changes in their plasma profiles. Snakes from the natural wetland initially had more variation of DNA damage; mean concentration of DNA damage in all snakes slightly decreased, but not significantly, after captivity. We present the plasma biochemistry profiles from western tiger snakes both before and after captivity and suggest a period of removal from natural stressors via captivity may offer a more reliable result of how plasma profiles of healthy animals might appear.

Evolutionary biology: adaptive developmental plasticity in snakes.

The morphology of organisms is generally well matched to their environment, presumably because expression of their genes is tailored either at the population or the individual level to suit local conditions: for example, snake populations that persistently encounter large prey may accumulate gene mutations that specify a large head size, or head growth may be increased in individual snakes to meet local demands (adaptive developmental plasticity). Here we test the relative contributions of genetics and environment to the jaw sizes of two tiger snake populations: one that consumes small prey on the mainland, and an island population that relies on larger prey and has a larger jaw size. Although the idea of adaptive plasticity in response to environmental pressures is controversial, we find that both factors influence the difference in jaw size between the two populations, and the influence of developmental plasticity is greater in the island population.

Comment on "Forest microclimate dynamics drive plant responses to warming".

Zellweger et al (Reports, 15 May 2020, p. 772) claimed that forest plant communities' response to global warming is primarily controlled by microclimate dynamics. We show that community thermophilization is poorly explained by the underlying components of microclimate, and that global warming primarily controls the climatic lag of plant communities. Deconstructing the underlying components of microclimate provides insights for managers.