Multiple models generate a geographical mosaic of resemblance in a Batesian mimicry complex.

Batesian mimics-benign species that receive protection from predation by resembling a dangerous species-often occur with multiple model species. Here, we examine whether geographical variation in the number of local models generates geographical variation in mimic-model resemblance. In areas with multiple models, selection might be relaxed or even favour imprecise mimicry relative to areas with only one model. We test the prediction that model-mimic match should vary with the number of other model species in a broadly distributed snake mimicry complex where a mimic and a model co-occur both with and without other model species. We found that the mimic resembled its model more closely when they were exclusively sympatric than when they were sympatric with other model species. Moreover, in regions with multiple models, mimic-model resemblance was positively correlated with the resemblance between the model and other model species. However, contrary to predictions, free-ranging natural predators did not attack artificial replicas of imprecise mimics more often when only a single model was present. Taken together, our results suggest that multiple models might generate a geographical mosaic in the degree of phenotype matching between Batesian mimics and their models.

Physiological thermal responses of three Mexican snakes with distinct lifestyles.

The impact of temperature on reptile physiology has been examined through two main parameters: locomotor performance and metabolic rates. Among reptiles, different species may respond to environmental temperatures in distinct ways, depending on their thermal sensitivity. Such variation can be linked to the ecological lifestyle of the species and needs to be taken into consideration when assessing the thermal influence on physiology. This is particularly relevant for snakes, which are a very functionally diverse group. In this study, our aim was to analyze the thermal sensitivity of locomotor performance and resting metabolic rate (RMR) in three snake species from central Mexico (Crotalus polystictus, Conopsis lineata, and Thamnophis melanogaster), highlighting how it is influenced by their distinctive behavioral and ecological traits. We tested both physiological parameters in five thermal treatments: 15 °C, 25 °C, 30 °C, 33 °C, and 36 °C. Using the performance data, we developed thermal performance curves (TPCs) for each species and analyzed the RMR data using generalized linear mixed models. The optimal temperature for locomotion of C. polystictus falls near its critical thermal maximum, suggesting that it can maintain performance at high temperatures but with a narrow thermal safety margin. T. melanogaster exhibited the fastest swimming speeds and the highest mass-adjusted RMR. This aligns with our expectations since it is an active forager, a high energy demand mode. The three species have a wide performance breadth, which suggests that they are thermal generalists that can maintain performance over a wide interval of temperatures. This can be beneficial to C. lineata in its cold habitat, since such a characteristic has been found to allow some species to maintain adequate performance levels in suboptimal temperatures. RMR increased along with temperature, but the proportional surge was not uniform since thermal sensitivity measured through Q10 increased at the low and high thermal treatments. High Q10 at low temperatures could be an adaptation to maintain favorable performance in suboptimal temperatures, whereas high Q10 at high temperatures could facilitate physiological responses to heat stress. Overall, our results show different physiological adaptations of the three species to the environments they inhabit. Their different activity patterns and foraging habits are closely linked to these adaptations. Further studies of other populations with different climatic conditions would provide valuable information to complement our current understanding of the effect of environmental properties on snake physiology.

Quantifying phenological diversity: a framework based on Hill numbers theory.

Background: Despite the great concern triggered by the environmental crisis worldwide, the loss of temporal key functions and processes involved in biodiversity maintenance has received little attention. Species are restricted in their life cycles by environmental variables because of their physiological and behavioral properties; thus, the timing and duration of species' presence and their activities vary greatly between species within a community. Despite the ecological relevance of such variation, there is currently no measure that summarizes the key temporal aspects of biological diversity and allows comparisons of community phenological patterns. Here, we propose a measure that synthesizes variability of phenological patterns using the Hill numbers-based attribute diversity framework. Methods: We constructed a new phenological diversity measure based on the aforementioned framework through pairwise overlapping distances, which was supplemented with wavelet analysis. The Hill numbers approach was chosen as an adequate way to define a set of diversity values of different order q, a parameter that determines the sensitivity of the diversity measure to abundance. Wavelet transform analysis was used to model continuous variables from incomplete data sets for different phenophases. The new measure, which we call Phenological Hill numbers (PD), considers the decouplings of phenophases through an overlapping area value between pairs of species within the community. PD was first tested through simulations with varying overlap in phenophase magnitude and intensity and varying number of species, and then by using one real data set. Results: PD maintains the diversity patterns of order q as in any other diversity measure encompassed by the Hill numbers framework. Minimum PD values in the simulated data sets reflect a lack of differentiation in the phenological curves of the community over time; by contrast, the maximum PD values reflected the most diverse simulations in which phenological curves were equally distributed over time. PD values were consistent with the homogeneous distribution of the intensity and concurrence of phenophases over time, both in the simulated and the real data set. Discussion: PD provides an efficient, readily interpretable and comparable measure that summarizes the variety of phenological patterns observed in ecological communities. PD retains the diversity patterns of order q characteristic of all diversity measures encompassed by the distance-based Hill numbers framework. In addition, wavelet transform analysis proved useful for constructing a continuous phenological curve. This methodological approach to quantify phenological diversity produces simple and intuitive values for the examination of phenological diversity and can be widely applied to any taxon or community's phenological traits.