Exploring the causes underlying the latitudinal variation in range sizes: Evidence for Rapoport's rule in spiny lizards (genus Sceloporus)

Species' range size is a fundamental unit of analysis in biodiversity research, given its association with extinction risk and species richness. One of its most notable patterns is its positive relationship with latitude, which has been considered an ecogeographical rule called Rapoport's rule. Despite this rule being confirmed for various taxonomic groups, its validity has been widely discussed and several taxa still lack a formal assessment. Different hypotheses have been proposed to explain their potential mechanisms, with those related to temperature and elevational being the most supported thus far. In this study, we employed two level of analyses (cross-species and assemblage) to investigate the validity of Rapoport's rule in spiny lizards (genus Sceloporus). Additionally, we evaluated four environmental-related hypotheses (minimum temperature, temperature variability, temperature stability since the last glacial maximum, and elevation) posed to explain such pattern, contrasting our results to those patterns expected under a null model of range position. Our results provided support for Rapoport's rule at both levels of analyses, contrasting with null expectations. Consistently, minimum temperature and elevation were the most relevant variables explaining the spatial variation in range size. At the cross-species level, our null simulations revealed that both variables deviated significantly from random expectations. Conversely, at the assemblage level, none of the variables were statistically different from the expected relationships. We discussed the implication of our findings in relation to the ecology and evolution of spiny lizards.

Seasonal and altitudinal variation in dorsal skin reflectance and thermic rates in a high-altitude montane lizard

Temperature is one of the most important factors in the life histories of ectotherms, as body temperature has an undeniable effect on growth, activity, and reproduction. Lizards have a wide variety of strategies to acquire and maintain body temperature in an optimal range. The "Thermal Melanism Hypothesis" proposes that individuals with lower skin reflectance can heat up faster as a result of absorbing more solar radiation compared to lighter conspecifics. Therefore, having a darker coloration might be advantageous in cold habitats. Dorsal skin reflectance has been found to change rapidly with body temperature in several lizard species, and it can also vary over longer, seasonal time scales. These variations may be important in thermoregulation, especially in lizards that inhabit areas with a large temperature variation during the year. Here, we study how dorsal reflectance fluctuates with body temperature and varies among seasons. We compared dorsal skin reflectance at three body temperature treatments, and measured thermal rates (i.e., heat and cool rate, thermic lapse, and net heat gain) by elevation (2500-4100 m) and seasons (spring, summer, and autumn) in the mesquite lizard, Sceloporus grammicus. Our results show that lizards were darker at high elevations and during the months with the lowest environmental temperatures. The rate of obtaining and retaining heat also varied during the year and was highest during the reproductive season. Our results indicate that the variation of dorsal skin reflectance and thermal rates follows a complex pattern in lizard populations and is affected by both elevation and season.

Flexibility in thermal requirements: a comparative analysis of the wide-spread lizard genus Sceloporus

Adaptation or acclimation of thermal requirements to environmental conditions can reduce thermoregulation costs and increase fitness, especially in ectotherms, which rely heavily on environmental temperatures for thermoregulation. Insight into how thermal niches have shaped thermal requirements across evolutionary history may help predict the survival of species during climate change. The lizard genus Sceloporus has a widespread distribution and inhabits an ample variety of habitats. We evaluated the effects of geographical gradients (i.e. elevation and latitude) and local environmental temperatures on thermal requirements (i.e. preferred body temperature, active body temperature in the field, and critical thermal limits) of Sceloporus species using published and field-collected data and performing phylogenetic comparative analyses. To contrast macro- and micro-evolutional patterns, we also performed intra-specific analyses when sufficient reports existed for a species. We found that preferred body temperature increased with elevation, whereas body temperature in the field decreased with elevation and increased with local environmental temperatures. Critical thermal limits were not related to the geographic gradient or environmental temperatures. The apparent lack of relation of thermal requirements to geographic gradient may increase vulnerability to extinction due to climate change. However, local and temporal variations in thermal landscape determine thermoregulation opportunities and may not be well represented by geographic gradient and mean environmental temperatures. Results showed that Sceloporus lizards are excellent thermoregulators, have wide thermal tolerance ranges, and the preferred temperature was labile. Our results suggest that Sceloporus lizards can adjust to different thermal landscapes, highlighting opportunities for continuous survival in changing thermal environments.

Maternal transmission of bacterial microbiota during embryonic development in a viviparous lizard

IMPORTANCE

We investigated the presence and diversity of bacteria in the embryos of the viviparous lizard Sceloporus grammicus and their amniotic environment. We compared this diversity to that found in the maternal intestine, mouth, and cloaca. We detected bacterial DNA in the embryos, albeit with a lower bacterial species diversity than found in maternal tissues. Most of the bacterial species detected in the embryos were also found in the mother, although not all of them. Interestingly, we detected a high similarity in the composition of bacterial species among embryos from different mothers. These findings suggest that there may be a mechanism controlling the transmission of bacteria from the mother to the embryo. Our results highlight the possibility that the interaction between maternal bacteria and the embryo may affect the development of the lizards.

Comparative analysis of two nonlethal methods for the study of the gut bacterial communities in wild lizards

Fecal samples or cloacal swabs are preferred over lethal dissections to study vertebrate gut microbiota for ethical reasons, but it remains unclear which nonlethal methods provide more accurate information about gut microbiota. We compared the bacterial communities of three gastrointestinal tract (GIT) segments, that is, stomach, small intestine (midgut), and rectum (hindgut) with the bacterial communities of the cloaca and feces in the mesquite lizard Sceloporus grammicus. The hindgut had the highest taxonomic and functional alpha diversity, followed by midgut and feces, whereas the stomach and cloaca showed the lowest diversities. The taxonomic assemblages of the GIT segments at the phylum level were strongly correlated with those retrieved from feces and cloacal swabs (rs > 0.84 in all cases). The turnover ratio of Amplicon Sequence Variants (ASVs) between midgut and hindgut and the feces was lower than the ratio between these segments and the cloaca. More than half of the core-ASVs in the midgut (24 of 32) and hindgut (58 of 97) were also found in feces, while less than 5 were found in the cloaca. At the ASVs level, however, the structure of the bacterial communities of the midgut and hindgut were similar to those detected in feces and cloaca. Our findings suggest that fecal samples and cloacal swabs of spiny lizards provide a good approximation of the taxonomic assemblages and beta diversity of midgut and hindgut microbiota, while feces better represent the bacterial communities of the intestinal segments at a single nucleotide variation level than cloacal swabs.

Distinct responses and range shifts of lizard populations across an elevational gradient under climate change

Ongoing climate change has profoundly affected global biodiversity, but its impacts on populations across elevations remain understudied. Using mechanistic niche models incorporating species traits, we predicted ecophysiological responses (activity times, oxygen consumption and evaporative water loss) for lizard populations at high-elevation (<3600 m asl) and extra-high-elevation (≥3600 m asl) under recent (1970-2000) and future (2081-2100) climates. Compared with their high-elevation counterparts, lizards from extra-high-elevation are predicted to experience a greater increase in activity time and oxygen consumption. By integrating these ecophysiological responses into hybrid species distribution models (HSDMs), we were able to make the following predictions under two warming scenarios (SSP1-2.6, SSP5-8.5). By 2081-2100, we predict that lizards at both high- and extra-high-elevation will shift upslope; lizards at extra-high-elevation will gain more and lose less habitat than will their high-elevation congeners. We therefore advocate the conservation of high-elevation species in the context of climate change, especially for those populations living close to their lower elevational range limits. In addition, by comparing the results from HSDMs and traditional species distribution models, we highlight the importance of considering intraspecific variation and local adaptation in physiological traits along elevational gradients when forecasting species' future distributions under climate change.

Elevation and blood traits in the mesquite lizard: Are patterns repeatable between mountains?

Ecogeographical patterns describe predictable variation in phenotypic traits between ecological communities. For example, high-altitude animals are expected to show elevated hematological values as an adaptation to the lower oxygen pressure. Mountains act like ecological islands and therefore are considered natural laboratories. However, the majority of ecophysiological studies on blood traits lack replication that would allow us to infer if the pattern reported is a local event or whether it is a widespread pattern resulting from larger-scale ecological processes. In lizards, in fact, the increase of hematological values at high altitudes has received mixed support. Here, for the first time, we compare blood traits in lizards along elevational gradients with replication. We tested the repeatability of blood traits in mesquite lizards between different elevations in three different mountains from the Trans-Mexican Volcanic Belt. We measured hematocrit, hemoglobin concentration, mean corpuscular hemoglobin concentration, and erythrocyte size in blood samples of low, medium, and high-elevation lizards. We obtained similar elevational patterns between mountains, but the blood traits differed among mountains. Middle-altitude populations had greater oxygen-carrying capacity than lizards from low and high altitudes. The differences found between mountain systems could be the result of phenotypic plasticity or genetic differentiation as a consequence of abiotic factors not considered.

Is Habitat More Important than Phylogenetic Relatedness for Elucidating the Gut Bacterial Composition in Sister Lizard Species?

The gut microbiota influences the phenotype and fitness of a host; however, limited information is currently available on the diversity and functions of the gut microbiota in wild animals. Therefore, we herein examined the diversity, composition, and potential functions of the gut microbiota in three Sceloporus lizards: Sceloporus aeneus, S. bicanthalis, and S. grammicus, inhabiting different habitats in a mountainous ecosystem. The gut bacterial community of S. bicanthalis from alpine grasslands at 4,150‍ ‍m a.s.l. exhibited greater taxonomic, phylogenetic, and functional alpha diversities than its sister species S. aeneus from cornfields and human-induced grasslands at 2,600‍ ‍m‍ ‍a.s.l. Bacteria of the genus Blautia and metabolic functions related to the degradation of aromatic compounds were more abundant in S. bicanthalis than in S. aeneus, whereas Oscillibacter and predicted functions related to amino acid metabolism and fermentation were more abundant in S. aeneus. The structure of the dominant and most prevalent bacteria, i.e., the core microbiota, was similar between the sister species from different habitats, but differed between S. grammicus and S. aeneus cohabiting at 2,600‍ ‍m‍ ‍a.s.l. and between S. grammicus and S. bicanthalis cohabiting at 4,150‍ ‍m a.s.l. These results suggest that phylogenetic relatedness defines the core microbiota, while the transient, i.e., non-core, microbiota is influenced by environmental differences in the habitats. Our comparisons between phylogenetically close species provide further evidence for the specialized and complex associations between hosts and the gut microbiota as well as insights into the roles of phylogeny and ecological factors as drivers of the gut microbiota in wild vertebrates.

Exceptional parallelisms characterize the evolutionary transition to live birth in phrynosomatid lizards

Viviparity, an innovation enhancing maternal control over developing embryos, has evolved >150 times in vertebrates, and has been proposed as an adaptation to inhabit cold habitats. Yet, the behavioral, physiological, morphological, and life history features associated with live-bearing remain unclear. Here, we capitalize on repeated origins of viviparity in phrynosomatid lizards to tease apart the phenotypic patterns associated with this innovation. Using data from 125 species and phylogenetic approaches, we find that viviparous phrynosomatids repeatedly evolved a more cool-adjusted thermal physiology than their oviparous relatives. Through precise thermoregulatory behavior viviparous phrynosomatids are cool-adjusted even in warm environments, and oviparous phrynosomatids warm-adjusted even in cool environments. Convergent behavioral shifts in viviparous species reduce energetic demand during activity, which may help offset the costs of protracted gestation. Whereas dam and offspring body size are similar among both parity modes, annual fecundity repeatedly decreases in viviparous lineages. Thus, viviparity is associated with a lower energetic allocation into production. Together, our results indicate that oviparity and viviparity are on opposing ends of the fast-slow life history continuum in both warm and cool environments. In this sense, the ‘cold climate hypothesis’ fits into a broader range of energetic/life history trade-offs that influence transitions to viviparity.

There have been five independent transitions from egg laying to live birth in the phrynosomatid lizards. Here, Domínguez-Guerrero et al. identify parallel changes in physiology, life history and behaviour that characterize these transitions to live birth.

Season-sex interaction induces changes in the ecophysiological traits of a lizard in a high altitude cold desert, Puna region

Functional traits are those characteristics of organisms that influence the ability of a species to develop in a habitat and persist in the face of environmental changes. The traits are often affected by a multiplicity of species-dependent and external factors. Our objective was to investigate thermal biology of Liolaemus ruibali in a high altitude cold desert at the arid Puna region, Argentina. We address the following question: do sex and seasonal variations in environmental temperature induce changes in the ecophysiological traits? We measured and compared the operative temperatures between fall and spring; and between sexes and seasons, we compared the ecophysiological traits of lizards, microenvironmental temperatures and thermoregulatory behavior. Air and operative temperatures were different between seasons. We found an effect of season-sex interaction on field body temperatures, preferred temperatures, panting threshold and thermal quality. The voluntary and critical temperatures presented seasonal variation in relation to changes in environmental temperatures, suggesting thermal acclimatization. We note behavioral changes between seasons, with the substrate being the main resource for gaining heat in spring. We conclude that Liolaemus ruibali is an efficient thermoregulator; it is a eurythermic lizard and presents phenotypic plasticity in different ecophysiological and behavioral traits induced by sex and seasonality. In addition, we predict that this population could buffer the effects of projected global warming scenarios.

Fast and dark: The case of Mezquite lizards at extreme altitude

Sprint speed is a major performance trait in animal fitness involved in escaping from predators, obtaining food, and defending territory. Biotic and abiotic factors may influence sprint speed in lizards. Temperature decreases at higher altitude. Therefore, lizards at high elevations may require longer basking times to reach optimal body temperatures, increasing their vulnerability to predation and decreasing their time for other activities such as foraging or reproduction. Here, we tested whether the maximum sprint speed of a lizard that shows conservative thermal ecology varied along an altitudinal gradient comprising low (2500 m), middle (3400 m) and high-altitude (4300 m) populations. We also tested whether sprint speed was related to dorsal reflectance at different ecologically relevant temperatures. Given that the lizard Sceloporus grammicus shows conservative thermal ecology with altitude, we expected that overall average sprint speed would not vary with altitude. However, given that darker lizards heat up quicker, we expected that darker lizards would be faster than lighter lizards. Our results suggest that S. grammicus at high altitude are faster and darker at 30 °C, while lizards from low and middle altitude are faster and lighter in color at 20 °C than high altitude lizards. Also, our results suggest a positive relationship between sprint speed and dorsal skin reflectance at 10 and 20 °C. Sprint speed was also affected by snout-vent length, leg length, and leg thickness at 10 °C. These results suggest that, even though predation pressure is lower at extreme altitudes, other factors such as vegetation cover or foraging mode have influenced sprint speed.

Thermoregulation of Liolaemus aparicioi (Iguania: Liolaemidae) along a 1000 m elevational gradient in La Paz Valley, La Paz, Bolivia

Lizard species have diverse behavioral and physiological responses to thermo-environmental conditions, which allow them to inhabit a broad range of latitudes and elevations. Because the availability of suitable thermal resources is limited and more variable at high-elevation environments than at lower elevations, we expect high-elevation lizards to be constrained in their thermoregulation relative to lizards at lower elevations by the fewer available thermal resources to reach optimal temperatures (colder environment). We studied the thermal biology of an endemic and Critically Endangered lizard, Liolaemus aparicioi, to assess its thermal responses along a 1000 m elevational gradient in La Paz Valley from May to August of 2015 (dry season). We took field body and microhabitat temperatures at capture sites (substrate and air above ground), and body size (snout-vent length and mass) of individuals at Taypichullo (3000 m asl), Gran Jardín de la Revolución Municipal Park (3500 m asl), and Taucachi (4000 m asl) localities. Operative temperatures were taken from calibrated models deployed in different available microhabitats. Preferred temperatures and thermal tolerance limits were determined in laboratory settings for lizards from each locality. Field body, microhabitat, and operative temperatures decreased with increasing elevation and differed between sexes. Lizards at the high elevation locality had the lowest thermoregulatory efficiency as compared with the mid and lower elevation localities. In laboratory measurements, while the preferred temperatures varied between sexes, pooled preferred temperatures and thermal tolerances were similar in all localities. Although thermal resources at high elevation can limit thermoregulatory possibilities in L. aparicioi, behavioral microhabitat use, time allocated to thermoregulation, and physiological adjustments seem to be possible strategies to counteract thermal costs along elevational gradients.

A nocturnally foraging gecko of the high-latitude alpine zone: Extreme tolerance of cold nights, with cryptic basking by day

Lizards that inhabit high-latitude alpine zones are exposed to extreme temperatures and long winters and most are diurnal heliotherms. Yet some poorly known nocturnal species exist in such locations, including several viviparous geckos from New Zealand. We studied the orange-spotted gecko (Mokopirirakau 'Roy's Peak'), a cryptic, nocturnal and viviparous lizard known only from the alpine zone (1150-1800 m a.s.l.) in the South Island (~44°S). Our field study investigated (1) the influence of female reproductive condition and sex on daytime body temperatures, including relationships with microhabitat rock temperatures, (2) the influence of temperature and other weather conditions on gecko emergence by night and day, and (3) the thermal microclimates available year-round to orange-spotted geckos. Building a better understanding of these lizards aids in species conservation efforts, for example in developing monitoring programmes, and provides insights into the evolution of thermal mechanisms in cold environments. Reproductive females maintained higher daytime body temperatures than non-reproductive females and males, suggesting pregnancy-related thermophily. On summer days, all reproductive groups reached similar body temperatures to New Zealand geckos from lower elevations, suggesting similar thermal preferences. Using trail cameras, we obtained evidence of geckos openly basking during the day (previously undocumented for this species) when temperatures of exposed lizard models (=T_exp) were 3.2-39.3 °C. We also observed emergence at night at low T_exp (-0.8-14.6 °C), when some T_bs were probably 0-6 °C. Diurnal activity increased as T_exp rose to peak at ~30 °C before dropping again at higher temperatures, whereas nocturnal activity unexpectedly decreased with increasing T_exp. Our study provides evidence of diurnal activity in a 'nocturnal' gecko that may be essential to squamate viviparity at high-latitude, high-elevation sites. It also suggests remarkable capacity for locomotor activity at extremely low T_b.

Living on the edge: Lower thermal quality but greater survival probability at a high altitude mountain for the mesquite lizard (Sceloporus grammicus)

A 20-month recapture analysis of 1001 individually marked mesquite lizards (Sceloporus grammicus) suggests that variation in thermal quality across three altitudes influences survival probability. Each additional unit of deviation from the temperature selected by these lizards in previous laboratory experiments (i.e. decreased thermal quality) meant an increase of roughly 1.01% in survival probability. Survival probabilities ranged from 0.80 to 0.90 at the lowest elevation site (2600 m), from 0.76 to 0.87 at the middle elevation site (3100 m) and from 0.90 to 0.94 at the highest elevation site (4150 m). These results suggest that in poor thermal quality environments mesquite lizards may employ thermoregulatory strategies (behavioral, physiological and/or morphological) to decrease their metabolic expenditure and their exposure to predators, maximizing survival. These findings highlight the relevance of thermal quality of the habitat in determining survival probability of ectotherms.

To be small and dark is advantageous for gaining heat in mezquite lizards, Sceloporus grammicus (Squamata: Phrynosomatidae)

Body temperature is important in determining individual performance in ectotherms such as lizards. Environmental temperature decreases with increasing altitude, but nevertheless many lizards inhabit high-altitude environments. The ‘thermal melanism hypothesis’ proposes that a dark dorsal coloration enables darker individuals to heat up faster because they absorb more solar radiation and thus being darker may be advantageous in cold habitats. The aim of the present study is to evaluate how heating rate, cooling rate and net heat gain vary with body size and dorsal skin coloration in Sceloporus grammicus lizards along an altitudinal gradient. We measured these traits multiple times in the same individuals with a radiation heat source and spectrophotometry under laboratory conditions. Our results showed that S. grammicus lizards are smaller and darker at high elevations than at low elevations. In addition, the smallest and darkest lizards showed the greatest heating rate and net heat gain. Therefore, in S. grammicus, we suggest that small body size and dark dorsal coloration provide thermoregulatory benefits in high-altitude environments. Hence, this study supports the thermal melanism hypothesis in a lizard species under varied thermal environments.

Comparative thermal ecology parameters of the mexican dusky rattlesnake (Crotalus triseriatus)

Montane habitats exhibit a high degree of thermal heterogeneity, and thus provide considerable thermoregulatory challenges for ectotherms. Comparative analyses provide an opportunity to understand how variation in abiotic factors (e.g., operative temperatures, thermal quality) can affect life history traits within species. We studied the thermal ecology of three populations of the rattlesnake Crotalus triseriatus inhabiting different volcanoes in the central region of Mexico using the Hertz et al. (1993) protocol. The average body temperature of dusky rattlesnakes from the three study sites was 22.4 °C; mean active body temperature was higher in site 2 than in sites 1 and 3, but no differences between females, males and juveniles nor an interaction among site and sex was found. The thermal quality was low in the three sites, particularly in sites 1 and 3. Thermoregulation accuracy statistically differed among populations: individuals from site 2 were more accurate thermoregulating, while individuals from site 1 were the least accurate. Compared to other snakes, dusky rattlesnakes can be considered as a eurythermic species, which can often be active at relatively low body temperatures.

Resting metabolic rates increase with elevation in a mountain-dwelling lizard

Individuals that inhabit broad elevational ranges may experience unique environmental challenges. Because temperature decreases with increased elevation, the ectotherms living at high elevations have to manage limited activity time and high thermoregulatory effort. The resting metabolic rate (RMR) of a postabsorptive animal is related to its total energy requirements as well as many other fitness traits. Mesquite lizards (Sceloporus grammicus) living on La Malinche Volcano, Mexico, inhabit a wide elevational range with some populations apparently thriving above the tree line. We measured the RMR of lizards from different elevations (i.e., 2,600, 3,200, and 4,100 m) at four ecologically relevant temperatures (i.e., 15, 25, 30, and 35 °C) and found that RMR of mesquite lizards increased with temperature and body mass. More importantly, lizards from the high-elevation population had mass specific RMR that was higher at all temperatures. While the higher RMRs of high-elevation populations imply higher metabolic costs at a given temperature these lizards were also smaller. Both of these traits may allow these high elevation populations to thrive in the face of the thermal challenges imposed by their environment.

Dietary effects on gut microbiota of the mesquite lizard Sceloporus grammicus (Wiegmann, 1828) across different altitudes

BACKGROUND

High-altitude ecosystems are extreme environments that generate specific physiological, morphological, and behavioral adaptations in ectotherms. The shifts in gut microbiota of the ectothermic hosts as an adaptation to environmental changes are still largely unknown. We investigated the food ingested and the bacterial, fungal, and protistan communities in feces of the lizard Sceloporus grammicus inhabiting an altitudinal range using metabarcoding approaches.

RESULTS

The bacterial phyla Bacteroidetes and Firmicutes, and the genera Bacteroides and Parabacteroides dominated the core fecal bacteriome, while Zygomycota and Ascomycota, and the species Basidiobolus ranarum and Basidiobolus magnus dominated the core fecal mycobiome. The diet of S. grammicus included 29 invertebrate families belonging to Arachnida, Chilopoda, and Insecta. The diversity and abundance of its diet decreased sharply at high altitudes, while the abundance of plant material and Agaricomycetes was significantly higher at the highest site. The composition of the fecal microbiota of S. grammicus was different at the three altitudes, but not between females and males. Dietary restriction in S. grammicus at 4150 m might explain the high fecal abundance of Akkermansia and Oscillopira, bacteria characteristic of long fasting periods, while low temperature favored B. magnus. A high proportion of bacterial functions were digestive in S. grammicus at 2600 and 3100, while metabolism of aminoacids, vitamins, and key intermediates of metabolic pathways were higher at 4150 m. Different assemblages of fungal species in the lizard reflect differences in the environments at different elevations. Pathogens were more prevalent at high elevations than at the low ones.

CONCLUSIONS

Limiting food resources at high elevations might oblige S. grammicus to exploit other food resources and its intestinal microbiota have degradative and detoxifying capacities. Sceloporus grammicus might have acquired B. ranarum from the insects infected by the fungus, but its commensal relationship might be established by the quitinolytic capacities of B. ranarum. The mycobiome participate mainly in digestive and degradative functions while the bacteriome in digestive and metabolic functions.