We carried out the study from April 2021 to May 2022 in peninsular India (Fig. 1). We initially targeted 10 species but only eight species could be sufficiently sampled for active and sleep sites (excluding Sarada darwini and Sitana laticeps). These species were chosen to cover a range of habitats: arboreal (Salea horsfieldii, Salea anamalayana), semi-arboreal (Calotes versicolor, Monilesaurus rouxii), ground (Sitana visiri, Sitana marudhamneydhal, Sarada superba), and rupicolous (Psammophilus dorsalis). All these species belong to the subfamily Draconinae and are closely related to each other overall (Grismer et al., 2016; Pal et al., 2018). At a finer phylogenetic scale, Monilesaurus and Psammophilus are sister genera and together, are sister to the genus Calotes (Pal et al., 2018); Sitana and Sarada are also sister genera to each other (Deepak et al., 2016). Habitat types of the sampled species are tightly linked to the fine-scaled pattern of phylogeny (Fig. 1), with both arboreal species coming from the same genus, both semi-arboreal species from sister genera and all three ground-dwelling species from sister genera. Apart from differences in habitat, all species share broad similarities in ecology, as they are diurnal, heliothermic, and insectivorous.
We chose sampling sites to be representative of typical habitats of the species and devoid of direct human disturbance (Online Resource 1 Table 1). The only exceptions were the sampling sites for Salea spp., which were transformed from native forests to Eucalyptus and Acacia plantations adjoining Shola-grasslands. Day sampling spanned from 0700h to 1100h, and night sampling from 1900h to 2300h, the timing of which was informed by previous studies on these and other species of diurnal agamid lizards (Radder et al., 2005; Mohanty et al., 2016, 2021). Two to three researchers searched with equal attention to all habitat strata, from the ground up to 4m; this approach was particularly important in the absence of a priori expectations for most species and to avoid bias. The sampling protocol broadly followed previous studies on sleep sites of semi-arboreal species (Mohanty et al. 2016; Bors et al. 2020) and for rupicolous species (Mohanty et al. 2021). We equated nocturnal perch sites to sleep sites, as individuals fulfilled some criteria for behavioural sleep at night (Piéron 1912; Tobler 1995): prolonged inactivity, relatively low responsiveness (e.g., easy to approach and capture by hand, without evoking movement), and typical sleep postures (head resting on perch, eyes closed).
For lizards of all habitat types, we categorised the substrate type (e.g., shrub, tree, rock), perch type (e.g., leaf, branch), perch angle (e.g., horizontal, angular, or vertical), and recorded perch height (perpendicular distance from lizard’s head to ground). For lizards observed on plants, we also measured the distance to main trunk (non-linear distance along branches), perch diameter (branch circumference at or closest to the location of lizard’s head), and trunk girth (maximum girth of trunk). Compliance of a perch is generally calculated using standard weights and measuring the resulting displacement (Samson and Hunt 2014), but diameter is a good proxy of compliance for branches and trunks (Gilman and Irschick 2013). We also recorded ‘head direction’ with respect to the potential approach path of a climbing predator as ‘inward’, ‘outward’, or ‘perpendicular’ (following Mohanty et al. 2016). For lizards on rocks, we defined head direction as ‘upward’, ‘angular-upward’, ‘perpendicular’, ‘downward’, or ‘angular-downward’ (the latter case was never observed). Individuals observed on the ground were not assigned any head direction.
As we were interested in the site use contrast between active and sleep phase in the same habitat, we did not estimate site selection by quantifying availability. Further, we did not sample the thermal-hydric conditions as an axis of habitat, as we considered it meaningful only in the presence of large spatial-scale capturing of availability, which we were unable to collect. We marked individuals using a felt pen on capture and did not re-collect habitat information for any individuals (i.e., different individuals were sampled for active sites and for sleep sites).
Data analyses
We retained observations of only adults with complete information in the dataset, excluding juveniles as they had limited sample sizes. We followed the protocol for data exploration by Zuur et al. (2010). We identified and excluded outliers (n = 21) by visualizing the raw data and subsequently excluded them by Tukey’s two fence method (values beyond the interval between Quartile 1 − 1.5x Interquartile range and Quartile 3 + 1.5x Interquartile range). Data was screened for multicollinearity by plotting pairwise scatterplots and examining correlation coefficients. We did not pool data at the level of habitat as preliminary results showed distinct patterns between species within habitat types (e.g., Calotes versicolor and Monilesaurus rouxii). Although we sampled in a range of habitats, species within each habitat were closely related, precluding analysis of phylogenetic signals. All statistical analyses were performed using R (version 4.4.1).
We compared microhabitat use across active and sleep phases on two scales: at a broad scale of parameter space composed of all variables and at a fine scale of each perch characteristic. To examine the overall parameter space across diel phase, we ran one principal component analysis (PCA) for each species. Continuous variables (perch height, trunk girth, perch diameter, and distance to trunk) were standardized and categorical variables (substrate type, perch angle, and head direction) were dummy coded (i.e., each category of a variable was defined as a new variable and scored as 0 (no) or 1 (yes); Bertolo et al., 2012; Miró et al., 2017). We did not include perch type in the PCAs, as we considered substrate type to be an adequate categorical descriptor at this scale and to avoid including a large number of perch types that are possible across substrate types. The PCAs for arboreal and semi-arboreal species included trunk girth, perch diameter, and distance to trunk, in addition to substrate type, perch height, perch angle, and head direction. For ground-dwelling and rupicolous species, PCAs included substrate type, perch height, perch angle, and head direction. However, for all three ground-dwelling species these variables had limited to no variation within or across diel phases, i.e., substrate type was predominantly ‘ground’, perch height was 0mm, perch angle was ‘horizontal’. Thus, a PCA could not be run for ground-dwelling Sitana marudhamneydhal; for Sarada superba and Sitana visiri, most observations coincided on one or two points on the PCA in the sleep phase, precluding any demarcation of parameter space by ellipses (see below).
We overlaid observations on the first two components of the PCA (PC1 and PC2) for each species and tested if they differed in their relative positions and variance using Permutational Multivariate Analysis of Variance (PERMANOVA; with 999 permutations) and multivariate homogeneity of groups dispersions, respectively. We then constructed multivariate normal ellipses (at 95% confidence interval) on these observations for each diel phase in all species, except the three ground-dwelling species which had limited variability. For these pairs of ellipses, we computed area (breadth), overlap, and uniqueness in absolute and in percentage of the total parameter space (union of ellipses of both diel phases).
At the fine scale, we compared each perch characteristic across diel phase (active vs sleep) at the species-level. Observation of an individual at a specific perch location is the outcome of a series of decisions, starting from choosing a substrate among several available substrates of varying girth, to climbing a certain distance from the ground and moving away from the trunk, and finally to orienting the head in a particular direction. Therefore, we examined the differences across diel phase (‘active’, ‘sleep’) in each perch characteristic individually. We used chi-square test of proportions for categorical variables (substrate type, perch type, perch angle, and head direction) and Ordinary Least Squares regressions (OLS) for continuous response variables (i.e., perch structure measurements). Substrate types and perch types were reclassified prior to analysis to reduce the number of categories; for example, saplings were merged along with shrubs and observations with a combination of perch types were assigned to one type (e.g., branch-leaf to leaf). For OLS regressions, all perch measurements were log transformed (1 + log10); normality and homogeneity of variance of residuals, as well as model validity were checked using diagnostic plots (Zuur and Ieno, 2016) along with Levene’s test. As trunk girth and perch diameter were highly correlated across species (r > 0.5), we considered only trunk girth and not perch diameter in subsequent models. Models with trunk girth as the response variable had phase as the only predictor. For models of perch height and distance to trunk in semi-arboreal and arboreal species, we included both phase and trunk girth as predictors because trunk girth can influence other perch measures (following Mohanty et al., 2016). In the case of rupiculous and ground-dwelling species, models for perch height had only phase as the predictor. We did not include sex as a predictor as we considered our sample sizes too low for testing sex effects by phase in each species. Body size was similar for samples of both phases in all species (Online Resource 1 Table 1) and therefore was not included in the analyses as a covariate.