Nathalia Rossigalli-Costa
Tiana Kohlsdorf- Department of Biology, FFCLRP—University of São Paulo, Ribeirão Preto, São Paulo, Brazil
Environmental conditions often affect developmental processes and consequently influence the range of phenotypic variation expressed at population level. Expansion of urban sites poses new challenges for native species, as urbanization usually affects the intensity of solar exposure and shade availability, determining the thermal regimes organisms are exposed to. In this study, we evaluate the effects of different developmental conditions in a Tropidurus lizard commonly found in Brazilian urban sites. After incubating embryos of Tropidurus catalanensis in two different thermal regimes (Developmental Environments [DE]: cold 24°C and warm 30°C), we measured morphological traits in the neonates and quantified locomotor performance in horizontal and vertical surfaces at three temperatures [Test Temperatures (TT) = 24°C, 30°C and 36°C]. Results indicate effects of developmental temperatures on morphological features, expressing functional implications that might be decisive for the viability of T. catalanensis in urbanized areas. Lizards ran similarly on horizontal and vertical surfaces, and isolated analyses did not identify significant effects of DE or TT on the sprint speeds measured. Absolute Vmax (i.e., the maximum sprint speed reached among all TTs) positively correlated with body size (SVL), and neonates from the warm DE (30°C) were larger than those from the cold DE (24°C). Morpho-functional associations of absolute Vmax also involved pelvic girdle width and forelimb, hindlimb, trunk, and tail lengths. Emerging discussions aim to understand how animals cope with abrupt environmental shifts, a likely common challenge in urbanized sites. Our findings add a new dimension to the topic, providing evidence that temperature, an environmental parameter often affected by urbanization, influences the thermal sensitivity of locomotion and the morphological profile of T. catalanensis neonates. Thermal sensitivity of specific developmental processes may influence the ability of these lizards to remain in habitats that change fast, as those suffering rapid urbanization due to city growth.
Introduction
Embryo development in oviparous vertebrates often involves exposure to environmental conditions that, even when buffered by the presence of protective membranes in the eggs or by parental care, may affect homeostasis and deviate from optimal conditions for growth and differentiation (Refsnider and Janzen, 2010; Noble et al., 2018). In lizards, for example, the choice of oviposition sites by females can affect the offspring, as incubation conditions influence traits such as hatching time, hatchling size, growth rate, and may ultimately determine the phenotypic profile of neonates (Warner and Andrews, 2002; Warner and Shine, 2008; Tiatragul et al., 2020). Among environmental parameters that usually affect development, temperature is probably one of the most relevant variables shaping the phenotype in ectotherm vertebrates (Angilletta, 2009). Thermal responses during ontogeny may produce different phenotypes (e.g., large or small neonate lizards; see Rossigalli-Costa et al., 2021) and can be addressed in the context of developmental plasticity and reaction norms (see Angilletta, 2009; Refsnider et al., 2019; Singh et al., 2020). When the thermal regimes change, species that exhibit behavioral flexibility or quickly respond to new conditions through plasticity in morphological or physiological traits tend to prosper in the novel environmental conditions (Winchell et al., 2016; Lapiedra et al., 2017; Noble et al., 2021).
Thermal relationships in embryos and tadpoles are receiving rising attention in the past decades (e.g., Clusella-Trullas and Chown, 2014; Bodensteiner et al., 2021; Hall and Sun, 2021; Taylor et al., 2021), a reflection of concerns from the scientific community with possible impacts of global warming on biodiversity (e.g., Moritz and Agudo, 2013; Diele-Viegas and Rocha, 2018; Diele-Viegas et al., 2020). The set of anthropic factors that accelerate climatic changes includes urbanization, a process that probably imposes the fastest and more radical changes in thermal regimes experienced by native species of ectotherm vertebrates. Expansion of urban centers involves the substitution of vegetated areas by impermeable surfaces, structural changes in habitat components (e.g., presence of walls and sidewalks), and modification of temperature and hydrological regimes due to changes in rates of sun exposure, shade availability and pollution (Forman, 2014). Urbanization may challenge some species to remain in their original distribution, as differences in environmental parameters between urban and natural sites might represent distinct selection pressures (Forman, 2014; Winchell et al., 2016), and lizards seem to be a particularly good system to evaluate phenotypic patterns and functional relationships with special attention to species found in urban sites (see Winchell et al., 2016; Johnson and Munshi-South, 2017). Studies focusing on urban lizards often interpret specific patterns in morphology (e.g., long limbs, more lamellae, larger toepads and wider bodies, see Winchell et al., 2016; Winchell et al., 2018b) and performance (e.g., high sprint speeds, see Winchell et al., 2018b) as adaptations to urbanization. Although most studies using urban lizards do not focus on embryo development, this is a relevant ontogenetic window because the choice of favorable microhabitats for oviposition in lizard species that lack parental care may buffer adverse conditions for embryo development imposed by specific environmental conditions (Shine and Harlow, 1996; Li et al., 2018), such as those derived from urbanization. On the other hand, in several species, the embryos may be very susceptible to environmental parameters (Wang et al., 2016; Bonnet et al., 2017) and pollutants (Lazić et al., 2015; Lazić et al., 2017) and, if the chosen oviposition site is not suitable, developmental plasticity may still enable proper growth and differentiation despite challenging environmental conditions (see Du et al., 2019). Plastic responses may comprise morphological and physiological traits as well as thermal sensitivity of specific functions. Shifts in incubation regimes may differentially affect specific sets of phenotypic traits (e.g., morphology, physiology, behavior; see Braña and Ji, 2000; Braña and Ji, 2007), and stability of functional associations—despite thermal sensitivity of particular developmental processes - might contribute for perpetuating the occurrence of some lizard species in urban sites.
The increasing number of studies addressing urbanization effects in different animals (see Rivkin et al., 2019) reflects an effort for understanding how new habitats (composed by buildings, impermeable surfaces at different orientations and high-density human population) might have functional consequences (e.g., Winchell et al., 2018a) that ultimately determine the phenotypic profile of populations of native species found in the cities. Some lineages, as the Neotropical Tropidurus lizards, are widely distributed and occupy different habitats, including species present in urban sites. The genus Tropidurus comprises habitat specialists but also generalist species (see Werneck et al., 2015), and records suggest that some species are relatively prone to occupy new niches and persist despite the surrounding urban growth (Rand and Rand, 1966). Another positive aspect of Tropidurus lizards for studies focusing on phenotypic aspects of species found in urban sites is the extensive literature available for eco-morphology (e.g., Kohlsdorf et al., 2001; Kohlsdorf et al., 2008; Grizante et al., 2010; Tulli et al., 2016), behavior (e.g., Leirião et al., 2019; Rodrigues and Kohlsdorf, 2019), locomotion (Kohlsdorf and Navas, 2012; Brandt et al., 2015) and physiology (Longhini et al., 2019) in this group, including information for development (Py-Daniel et al., 2017; Rossigalli-Costa et al., 2021) and thermal relationships (Kohlsdorf and Navas, 2006; Leirião et al., 2019; Piantoni et al., 2019) of species found in urban areas. In this study, we isolated temperature effects from other environmental parameters likely affected by urbanization (e.g., hydrological regimes or pollutants) and evaluate developmental plasticity in offspring of a Tropidurus species commonly found in Brazilian cities. We developed T. catalanensis embryos in different thermal regimes in the laboratory [Developmental Environments (DE): cold = 24°C and warm = 30°C], to test the hypothesis that temperature during development affects the offspring phenotypic profile, evaluated from locomotor performance and morphological traits. Results provide evidence for phenotypic effects of developmental environments in morphological features among neonates of T. catalanensis incubated in different thermal regimes, with consequences for locomotion. The physiological and morpho-functional implications of different developmental conditions may be decisive for the permanence of T. catalanensis lizards in urbanized areas, a discussion that contributes for a better understanding of how native species cope with environmental shifts that might be accentuated in urbanized sites.
Materials and Methods
We evaluated the effects of developmental environments in the phenotypic profile of neonate lizards focusing on two sets of traits: 1) locomotor performance and 2) morphology. After incubating eggs of T. catalanensis in two different thermal regimes [Developmental Environments (DE): cold = 24°C and warm = 30°C], we measured sprint speeds at two orientations (horizontal and vertical surfaces) and three different temperatures [Test Temperatures (TT) = 24°C, 30oC and 36°C], assumed as a proxy of essential free-living activities, and quantified body size and limb and tail relative proportions. Data acquisition, experimental design and statistical analyses are described along this section.
Egg Collection and Incubation
A total of 20 gravid females of T. catalanensis were collected by lasso (Van Sluys, 1993) at the Brazilian city of São Simão, in the state of São Paulo (21° 28′ 44″ S 47° 33′ 03″ W), which was founded in 1824 and experienced intensified urbanization in the past two decades. Animals were transported in cloth bags to the Animal Facility of the Laboratory of Integrative Biology and Evolution (University of São Paulo at Ribeirão Preto, SP, Brazil) in the same day they were collected. Females were maintained in a maximum of three individuals per terrarium, which consisted of plastic boxes (40 cm × 40 cm × 60 cm) covered with a 15 cm layer of moist medium-grain vermiculite, until oviposition. Food was provided every two days and water was offered ad libitum; thermoregulation at a 12 h dark:light cycle was granted by incandescent UVA/UVB (for vitamin D production) and 60 W incandescent lamps (adapted from Sanger et al., 2008). We checked each terrarium daily for eggs, and carefully examined each female to identify which individual laid a given clutch in order to control for genetic and maternal effects. From the 20 clutches obtained (one from each female), we distributed 80 eggs into small plastic cylindric containers of 145 ml (7.5 cm × 5.0 cm) filled with moist vermiculite. We perforated the lids to enable air renovation and fixed the water potential at -150 cmkPa, replacing eventual evaporated water (based on Steele et al., 2018 and following; Rossigalli-Costa et al., 2021). To avoid genetic bias, we separated the clutches and equally distributed the eggs between two developmental environments: 24°C (cold, N = 40) and 30°C (warm, N = 40); we did not consider the thermal regime of 36°C (overheat) as a third developmental condition because this temperature seems unfeasible for the development of T. catalanensis embryos (see Rossigalli-Costa et al., 2021). Due to the lack of information for nest temperatures in T. catalanensis, we chose developmental temperatures within the range of air, substrate and body temperatures recorded for Tropidurus lizards at São Simão and other Brazilian cities where this species is found (see Kohlsdorf and Navas, 2006; Leirião et al., 2019; Maia-Carneiro and Navas, 2021), which also agree with the range of preferred temperatures recorded for Tropidurus lizards in thermal gradients (31°C to 37°C, see Kohlsdorf and Navas, 2006). Incubation temperatures were controlled by Bio-Oxygen Demand (BOD) incubators (Eletrolab®), and the precision of thermal regimes was confirmed using data loggers (NOVUS®). Eggs were daily inspected for monitoring the development of the embryos, being maintained in the incubators until hatching. These procedures followed Rossigalli-Costa et al. (2021).
Locomotor Performance
We started locomotor performance tests seven days after hatching. Until performing the locomotion trials, neonates were maintained in terraria at conditions similar to those described for gravid females. Specifically, they were placed in terraria maintained at a 25°C acclimatized room and 60 W incandescent lights were provided in each terrarium to grant thermoregulation. Photoperiod was established at a 12 h dark:light cycle. Tests were performed in an acclimatized room, and we ran each lizard individually. A test comprised one horizontal and one vertical race from each lizard in a given temperature [Test Temperatures (TT) = 24°C, 30°C, and 36°C], in a total of two races in a given day for each individual. Tests were intercalated every two days, and test temperatures were repeated at two different days. At the end of all locomotor performance tests, we gathered two sprint speeds measurements at each temperature (24°C, 30°C, and 36°C) and orientation (horizontal or vertical) for each individual, obtained in different days (total of 12 records for each lizard). The Test Temperatures chosen reflect the thermal variation experienced by Tropidurus lizards during daily hours of activity (Maia-Carneiro and Navas, 2021). We randomized the order of tests (i.e., in which temperature lizards ran) and also the orientation (horizontal or vertical) they started running each day. Before each race, lizards were maintained for at least one hour at the testing temperature until reaching body temperatures near those of the test (confirmed using an infrared thermometer). We recorded sprint speeds at horizontal and vertical directions of 20 neonates - ten from ‘cold DE’ (24°C) and ten from ‘‘warm DE’’ (30°C). From the total of 80 eggs incubated, 36 neonates from 30°C and 34 neonates from the cold developmental environment of 24°C actually hatched. We ran all the 70 neonates that hatched, but several individuals refused to run at one or more trials at a given temperature, and the time window we considered one lizard a ‘‘neonate’’ was 15 days. Therefore, and in order to include only individuals that performed satisfactorily at the three testing temperatures to infer running thermal sensitivity, our final dataset comprised races obtained for 10 lizards from each developmental environment. The horizontal sprint speed was measured at a pine board racetrack (90 cm length × 25 cm width × 35 cm high). The camera was placed 2 m above the center of the racetrack and captured its entire length. Vertical sprint speed was performed at 80 cm high × 25 cm width concrete blocks leaning on the wall, with the camera placed at a 2 m frontal distance from the blocks. For both tests, neonates were placed on the track and stimulated by to run by hand if they did not flee immediately; sprints speeds were recorded until lizards reached the end of the track. All races were recorded with a high-speed digital camera at 120fps (Canon® Sx50hs). Videos were analyzed using the Tracker® video analyzing software, in which we extracted the total time to complete the races (measured in seconds), maximum and average sprint speed (m/s) and maximum and average acceleration (m/s2) from each trial.
Morphological Measurements
We compiled a dataset of morphological measurements obtained in each lizard immediately after hatching to evaluate the effects of developmental conditions on body size and associated proportions. We avoided measuring animals along the days of performance tests to minimize stress from excessive manipulation. All data in this study were obtained in the first two weeks after hatching, and during this time window neonate lizards barely feed and seem to rely on remaining embryonic energetic reserves; in the lab, growth of T. catalanensis during the first days after hatching rarely exceeds 5% of their initial size. We focused on traits related to locomotor performance in Tropidurus lizards (see Kohlsdorf et al., 2001; Grizante et al., 2010; Kohlsdorf and Navas, 2012): Hand and Feet lengths (from the base of the wrist/ankle to the distal phalanx of Digit IV); Fore and Hindlimb lengths (the sum of humerus and radius/femur and tibia lengths); Height and Width of the Pectoral and Pelvic Girdles; Trunk length; Tail length and Snout-Vent length (SVL, a proxy of overall body size). Traits were measured externally and congregate osteological, muscular, and epithelial elements; all measurements were performed by the same researcher using a digital calliper (Mitutoyo Inc. ± 0.01 mm).
Statistical Analyses
Analyses were performed based on two parameters of sprint performance: 1) maximum sprint speed (Vmax; the maximum sprint speed the individual reached in each testing temperature); and 2) absolute maximum sprint speed (Absolute Vmax, the highest sprint speed value reached amongst all test temperatures). We also analyzed our data separately considering the thermal sensitivity of locomotor performance and functional associations with morphological traits, as further described. All analyses were performed using RStudio (1.4.1717) in R (version 4.0.5 - R Core Team, 2021).
We first performed descriptive analyses of the total duration of the race (in seconds), Vmax (m/s) at each test temperature and Absolute Vmax (m/s). A generalized linear mixed model (glmm; package glmmTMB, Brooks et al., 2017) was used to evaluate relationships of Vmax with surface orientation (vertical or horizontal), developmental environment (DE: cold or warm), and testing temperature (TT: 24°C, 30°C, and 36°C). Interactions between orientation, developmental environment (DE) and testing temperatures (TT) were considered in the models, and ‘‘individual’’ and ‘‘clutch’’ (maternal effects) were considered as random intercepts. When significant interaction effects were detected, we performed post-hoc pairwise analyses using the emmeans function (Russell, 2021). We corrected p-values for multiple comparisons using the Tukey method to compare a family of three estimates. Prior to these steps, we selected the statistical family distribution that best fitted our model through the Akaike information criterion and proceeded with the analyses using the gamma distribution. To visualize the thermal sensitivity of locomotor performance, we performed a linear regression, using ggplot2 (Wickham, 2016), between Vmax (m/s) at horizontal and vertical orientation and the body temperature (in °C) the animals reached during the locomotion trials (after being maintained for 1 h at the test temperature).
The morphological database was evaluated using a multivariate approach. First we extracted size-adjusted residuals of shape (see Shingleton, 2010), and then we performed a Principal Component Analysis (PCA) using the package vegan (Oksanen et al., 2019). We selected the Principal Components (PCs) having eigenvalues higher than 1 (PCs 1–4, which together explained 78.01% of the total variation) and plotted the corresponding scores with the predictor variable ‘‘development environment’’ (DE) for visual inspection of the morphospace (packages ggplot2 from Wickham, 2016 and ggord from Beck, 2021). We tested for morphological differences between individuals incubated at different developmental environments performing an analysis of variance (ANOVA) on the PC scores.
Associations between the morphological profile of neonates and locomotor performance were tested using linear models with the PCs scores as a function of Absolute Vmax, after selecting the best models for our data. In the models, Absolute Vmax was included as the response variable and different combinations of the four morphological PCs, the developmental environment (DE), and surface orientation (vertical or horizontal) were considered as predictor variables. We applied an automated model selection for each complex model using the package glmulti (Calcagno, 2020), starting from [glm (Absolute Vmax ∼ PC1 + PC2 + PC3 + PC4 + DE + Orientation)], and then sequentially removing variables that explained the least variation in each model tested (see Barros et al., 2021 for an example of similar analyses). The best models and their AICc support were used for model selection. For all models, ‘‘clutch’’ and ‘‘individual’’ were fixed as random effects. Integrated effects of the predictor variables on Absolute Vmax were tested using glmm. All analyses using PCs were size-corrected, so we complemented our statistical design using ANCOVAs to test Absolute Vmax (horizontal and vertical separately) as a function of integrated effects of body size (SVL) and DE, in order to address morpho-functional implications from the effect of DE on growth. All morphology-related models were tested using the Akaike Criterion Information for the best statistical family distribution and were implemented in the glmmTMB package (Brooks et al., 2017). Under request, the authors will provide detailed information on analyses and codes as an Rmarkdown file.
Results
Locomotor Performance of Tropidurus Neonates From Different Developmental Environments: Relationships With Test Temperatures
Races performed at the vertical wall (80 cm length) overall lasted longer than those at horizontal orientation (90 cm length), as detailed in Table 1 and Supplementary Table S1. In the locomotor trials, neonates reached body temperatures near 23.81 ± 0.62 (TT = 24°C), 28.38 ± 0.86 (TT = 30°C), and 31.80 ± 0.88 (TT = 36°C). We first evaluated thermal sensibility of locomotion by plotting Vmax at each orientation as a function of the body temperature reached at locomotor performance trials (Figure 1). Visual inspection of the curves suggested slight variation in the thermal sensitivity of Vmax between Tropidurus lizards from distinct developmental environments in both orientations (horizontal and vertical), despite the considerable variation among individuals observed in the graphs. When effects of Developmental Environment (DE), Test Temperature (TT), or the interaction among factors (Orientation*DE*TT) on maximum sprint speeds were tested separately, we did not identify any significant association (Supplementary Table S2). Optimal temperature for sprint speed was not consistent among individuals, even within the same developmental condition (Figure 1), so we used the highest value of maximum speed from all tests (Absolute Vmax) to evaluate variation in locomotor performance and morpho-functional relationships between neonates of Tropidurus lizards incubated at different thermal regimes.
TABLE 1. Locomotor performance (mean and standard error) of Tropidurus catalanensis neonates raised in different developmental environments (DE).
FIGURE 1. Thermal performance of maximum sprint speeds (Vmax) in lizards running at Horizontal and Vertical orientations in the body temperature reached at the Testing Temperature (TT) of 24°C, 30°C, and 36°C. Animals from different developmental environments (DE) coded as follows: 24°C, green dots; 30°C, purple triangles.
Morphological Profiles and Functional Relationships in Neonates From Different Developmental Environments
The morphological profiles of T. catalanensis neonates differed according to the thermal regime eggs were incubated (raw data presented at Supplementary Table S3). Lizards from the warm Developmental Environment (DE = 30°C) were bigger (i.e., larger SVL) than those from the cold DE of 24°C (p < 0.0001). The Principal Components Analyses performed with body proportions (size-adjusted residuals) retrieved four Principal Components (PCs) with eigenvalues higher than 1.0, which, together, explained 78.01% of total variation. The two first principal components explained together almost 50% of the variation (Table 2) and were mostly associated to variation in pelvic and pectoral girdles (PC1) and limbs (foot, hand and hindlimb, PC2). Among the four PCs retained, PC2 (hand, foot and hindlimb lengths) and PC3 (hindlimb, trunk and tail lengths) differed between individuals from distinct DEs (PC1: p = 0.381; PC2: p < 0.0001; PC3: p = 0.004, PC4: p = 0.960; see also Figure 2).
TABLE 2. Principal component analysis indicating the percentage of explained variance of each component retained.
FIGURE 2. Principal component analysis of morphological traits measured in 20 juveniles from the warm DE (30°C; purple triangles) and the cold DE (24°C; green dots). Components with eigenvalues higher than one were retained, and loadings of variables strongly associated to each PC are detailed in Table 2. The scores of PC2 and PC3 differed between neonates from the cold and warm DE (PC2: p = 0.001; PC3: p = 0.004). Codes for morphological traits are: H, Hand length; F, Foot length; FL, Frontlimb length; HL, Hindlimb length; PtGH, Pectoral Girdle Height; PtGW, Pectoral Girdle Width; PvGH, Pelvic Girdle Height; PvGW, Pelvic Girdle Width; Tr, Trunk length; T, Tail length.
Relationships between locomotor performance (Absolute Vmax) and body size (also considering effects of DE on SVL) were tested using ANCOVA, as lizards from the warm DE were larger (see Figure 3 and Supplementary Figure S1). We identified effects of body size on Absolute Vmax at vertical orientation (p = 0.0225), but no effects were detected from the interaction between SVL and DE (p > 0.05) in the speeds recorded at this orientation; when running at horizontal orientation, no effects on Absolute Vmax were detected from body size or its interaction with DE (Figure 3).
FIGURE 3. Values of Absolute Vmax (m/s) in relation to body size (SVL) at two orientations (Horizontal and Vertical). Lizards from the warm DE (30°C, represented by purple triangles) had larger SVL than those from the cold DE (represented by green dots), and body size positively affected Absolute Vmax at vertical races (p = 0.0225).
Functional relationships of locomotor performance in neonates from different developmental environments were better explained by the model comprising two morphological components (PC3 and PC4) and the Developmental Environment (DE) (Model 1: Absolute Vmax∼1 + PC3+PC4+DE; AICc = 9.647, weight = 0.161; Table 3). Morphological traits with higher loads in PC3 (18.59% of total variation) were hindlimb trunk and tail lengths, and those in PC4 (11.6% of total variation) were frontlimb and trunk lengths and pectoral girdle width (Table 2). We tested for specific associations between Absolute Vmax and the three factors of the best model selected (M1) - PC3, PC4, and DE - and all were statistically significant (Table 3). To summarize, our results provide evidence for associations between Absolute Vmax and body size (SVL) and relative proportions (limbs, trunk, and tail lengths, and pelvic girdle width) that are associated with the thermal environment where lizards developed, suggesting that effects of environmental conditions at early developmental stages on morphological traits also impact locomotor performance.
TABLE 3. Results of the three best generalized mixed models selected for morpho-functional relationships of Absolute Vmax.
Discussion
Thermal regimes are particularly important for embryos developing in environments susceptible to external conditions, including oviparous lizard species, because temperature may affect the establishment of features directly related to fitness. Environmental changes in specific settings due to urbanization processes received increasing attention in the past decades, due to expected impacts on biodiversity derived from differences in rates of sun exposure, shade availability, and structural components of the habitat during city growth (e.g., McDonnell and Hahs, 2015; Winchell et al., 2017; Rivkin et al., 2019; Miles et al., 2021), as well as variation in biotic interactions that include predation and structure of the community (Birnie-Gauvin et al., 2016; Andrade et al., 2018). Although urbanization processes may negatively affect native species and consequently demand conservation actions (see Bradley and Altizer, 2007; McKinney, 2008), some species benefit from new environmental contexts that represent increments in the availability of ecological niches (Eriksson, 2013).
In the present study, we evaluated effects of developmental temperatures in the phenotypic profiles of a lizard commonly found in Brazilian cities, inferring functional implications for locomotor performance in neonates of Tropidurus catalanensis incubated in different thermal regimes. Our results provide evidence that the thermal regime experienced during development influences the morphological profile of T. catalanensis neonates, which indirectly affects locomotor performance. Specifically, individuals from warm conditions were larger, and body size affected absolute maximum speeds of T. catalanensis especially at vertical surfaces. Moreover, morphological traits associated to PC2 (hand, foot and hindlimb lengths) and PC3 (hindlimb, trunk and tail lengths) that differed between individuals incubated at different temperatures correspond to variables often interpreted as adaptations to urbanization in other lizard species (see Winchell et al., 2016; Winchell et al., 2018b). Although the Tropidurus lizards overall ran at similar speeds in horizontal and vertical races and analyses performed on isolated traits did not indicate significant interactions of locomotor performance with Developmental Environment or Test Temperatures, neonates from the warm DE (30°C) were larger, and body size affected Absolute Vmax at vertical orientation. Moreover, some of the morphological traits that differed between neonates from different DEs (PC3: hindlimb, trunk and tail lengths) were associated with the best-fit model for morpho-functional associations of Absolute Vmax. Interestingly, effects of the DE on locomotor performance were only detected when the morphological profile of neonates was considered in the analyses, regardless of the test temperature. We performed locomotor trials at 24°C, 30°C, and 36°C, and most lizard species select preferred temperatures between 25°C and 35°C and exhibit sprint speeds usually optimized between 32°C and 36°C (Huey and Bennett, 1987). Lizards of T. catalanensis are diurnal and heliophilic, and experience temperatures between 25°C and 36°C during daily activities that include thermoregulation, foraging, and reproduction (Maia-Carneiro and Navas, 2021). Despite the lack of information about thermal regimes in the nests of T. catalanensis, literature records indicate an average temperature near 29°C in nests of other lizards found in urban areas, such as Anolis cristatellus (Tiatragul et al., 2020). Although the two developmental thermal regimes used here are within the range of field temperatures recorded for Tropidurus lizards (Kohlsdorf and Navas, 2006; Leirião et al., 2019; Maia-Carneiro and Navas, 2021), the developmental temperature of 24°C established in the lab for T. catalanensis might be interpreted as a cold environment, especially considering that preferred temperatures for tropidurine lizards are usually above 30°C (Kohlsdorf and Navas, 2006).
Morphological differences in Tropidurinae lizards are often associated to ecological divergence (Kohlsdorf et al., 2001; Grizante et al., 2010) and explain interspecific variation in locomotor performance (Kohlsdorf and Navas, 2012). In this study, most of the morphological variation among Tropidurus neonates from different developmental environments concentrates in two principal components that aggregate height and width of pelvic and pectoral girdles (PC1) and limb, hand and feet relative proportions (PC2). Locomotor performance was indirectly affected by the effects of incubation temperature in two other Principal Components (PC3 and PC4), which aggregate limb trunk and tail lengths and the pectoral girdle width. Interestingly, these parameters correspond to morphological traits interpreted in Anolis lizards as adaptive to urbanization (trunk width, limb lengths and autopodial features; see Winchell et al., 2016; Winchell et al., 2018). Our results therefore foster future studies addressing the role of developmental plasticity in adaptive processes of T. catalanensis populations in Brazilian urbanized sites.
Morphological and physiological changes mediate evolutionary processes in natural populations that may lead to local adaptation and rapid evolution in response to environmental changes (Winchell et al., 2016; Winchell et al., 2018; Noble et al., 2021). Despite expected negative effects of urbanization for native species (see McDonnell and Hahs, 2015), high densities of T. catalanensis lizards are observed in several Brazilian cities. Part of the success to remain in these areas despite rapid urbanization process may reside on how animals cope with the effects of developmental temperature on the overall morphological profile of neonates and their consequences for morpho-functional associations of sprint performance. We recognize that isolated effects of developmental environments in these lizards may influence other features related to fitness at different ages, but individual survival in the critical ontogenetic window immediately after birth might be granted by stability of sprint performance morpho-functional associations that likely enhances chances of escaping from predators (e.g., cats and birds) in Brazilian cities. Warmer developmental environments produce larger Tropidurus neonates, which might be favored in urban environments also due to the positive association between sprint speeds and body size. Studies evaluating how native species successfully remain in anthropized fields perfectly complement data available for the negative effects of urbanization that might lead to local extinction of specific lineages (see Winchell et al., 2018; Sparkman et al., 2018; Putman et al., 2019; Vaughn et al., 2021 for recent discussions). Our findings highlight the importance of a new dimension to this discussion (i.e., thermal effects in developmental processes), providing evidence that rapid shifts in specific environmental parameters, such as temperature, may influence the development of isolated traits in T. catalanensis and consequently affect morpho-functional relationships in the species, an aspect that must be considered by studies addressing how these lizards cope with habitats that change fast, as those suffering rapid urbanization due to city growth.
Data Availability Statement
The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.
Ethics Statement
The animal study was reviewed and approved by Ethical Committee for Animal Use from the University of São Paulo approved all experimental procedures involving live animals (CEUA-USP 2019.1.1509.59.4).
Author Contributions
NR-C and TK conceived and designed the study; NR-C conducted data acquisition and statistical analyses; both authors interpreted the results and wrote the manuscript.
Funding
This study was funded by the CAPES fellowship awarded to NR-C (88882.180050/2018-01) and FAPESP thematic grants awarded to TK (2015/07650-6 and 2020/14780-1).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We are thankful to Bianca Bonini-Campos, Mario Bronzati Filho, Stella Marco Kyomen and Gabriel Spanghero Vicente Ferreira for help during fieldwork and husbandry of gravid females and eggs clutches. We acknowledge Danilo Boscolo and Wilfried Klein (USP-Brazil) for lending equipment used in this study. Special thanks to Ana Paula Cipriano and Luísa Prado Prandini for help with data acquisition, to Melissa Bars Closel for lending high-speed camera, to Miranda Rew-Duffy and Natalie Freeman for training with Tracker and video analysis, and to Rodolfo Mei Pelinson and Gabriela Baroni Leite for statistical support in R Studio. Finally, we are especially thankful to the reviewers from Frontiers, which provided relevant suggestions that significantly contributed to this manuscript.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphys.2022.891545/full#supplementary-material
References
Andrade A. C., Franzini L. D., Mesquita D. O. (2018). Assessing the Effect of Urbanization on Tropical Forest Dwelling Teiid Lizards. Ecol. Indic. 99, 225–229. doi:10.1016/j.ecolind.2018.12.037
Angilletta M. J. (2009). Thermal Adaptation: A Theoretical and Empirical Synthesis. New York, NY: Oxford University Press Inc.
Barros F. C., Grizante M. B., Zampieri F. A. M., Kohlsdorf T. (2021). Peculiar Relationships Among Morphology, Burrowing Performance and Sand Type in Two Fossorial Microteiid Lizards. Zoology 144, 125880. doi:10.1016/j.zool.2020.125880
Birnie-Gauvin K., Peiman K. S., Gallagher A. J., De Bruijn R., Cooke S. J. (2016). Sublethal Consequences of Urban Life for Wild Vertebrates. Environ. Rev. 24, 416–425. doi:10.1139/er-2016-0029
Bodensteiner B. L., Agudelo‐Cantero G. A., Arietta A. Z. A., Gunderson A. R., Muñoz M. M., Refsnider J. M., et al. (2021). Thermal Adaptation Revisited: How Conserved Are Thermal Traits of Reptiles and Amphibians? J. Exp. Zool. 335, 173–194. doi:10.1002/jez.2414
Bonnet X., Naulleau G., Shine R. (2017). The Evolutionary Economics of Embryonic-Sac Fluids in Squamate Reptiles. Am. Nat. 189, 333–344. doi:10.1086/690119
Bradley C. A., Altizer S. (2007). Urbanization and the Ecology of Wildlife Diseases. Trends Ecol. Evol. 22, 95–102. doi:10.1016/j.tree.2006.11.001
Braña F., Ji X. (2000). Influence of Incubation Temperature on Morphology, Locomotor Performance, and Early Growth of Hatchling Wall Lizards (Podarcis Muralis). J. Exp. Zool. 286 (4), 422–433. doi:10.1002/(sici)1097-010x(20000301)286:4<422::aid-jez10>3.0.co;2-d
Braña F., Ji X. (2007). The Selective Basis for Intreased Egg Retention: Early Incubation Temperature Determines Hatchling Phenotype in Wall Lizards (Podarcis Muralis). Biol. J. Linn. Soc. 92, 441–447. doi:10.1111/j.1095-8312.2007.00871.x
Brandt R., Galvani F., Kohlsdorf T. (2015). Sprint Performance of a Generalist Lizard Running on Different Substrates: Grip Matters. J. Zool. 297, 15–21. doi:10.1111/jzo.12253
Brooks M. E., Kristensen K., Benthem K. J. v., Magnusson A., Berg C. W., Nielsen A., et al. (2017). glmmTMB Balances Speed and Flexibility Among Packages for Zero-Inflated Generalized Linear Mixed Modeling. R J. 9, 378–400. doi:10.32614/rj-2017-066
Calcagno V. (2020). Glmulti: Model Selection and Multimodel Inference Made Easy. R package version 1.0.8. Available at: https://CRAN.R-project.org/package=glmulti (Accessed January 2022).
Clusella-Trullas S., Chown S. L. (2014). Lizard Thermal Trait Variation at Multiple Scales: A Review. J. Comp. Physiol. B 184, 5–21. doi:10.1007/s00360-013-0776-x
Diele-Viegas L. M., Figueroa R. T., Vilela B., Rocha C. F. D. (2020). Are Reptiles Toast? A Worldwide Evaluation of Lepidosauria Vulnerability to Climate Change. Clim. Change 159, 581–599. doi:10.1007/s10584-020-02687-5
Diele-Viegas L. M., Rocha C. F. D. (2018). Unraveling the Influences of Climate Change in Lepidosauria (Reptilia). J. Therm. Biol. 78, 401–414. doi:10.1016/j.jtherbio.2018.11.005
Du W.-G., Shine R., Ma L., Sun B.-J. (2019). Adaptive Responses of the Embryos of Birds and Reptiles to Spatial and Temporal Variations in Nest Temperatures. Proc. R. Soc. B 286, 20192078. doi:10.1098/rspb.2019.2078
Eriksson O. (2013). Species Pools in Cultural Landscapes - Niche Construction, Ecological Opportunity and Niche Shifts. Ecography 36, 403–413. doi:10.1111/j.1600-0587.2012.07913.x
Grizante M. B., Navas C. A., Garland Jr T., Kohlsdorf T. (2010). Morphological Evolution in Tropidurinae Squamates: An Integrated View along a Continuum of Ecological Settings. J. Evol. Biol. 23, 98–111. doi:10.1111/j.1420-9101.2009.01868.x
Hall J. M., Sun B. j. (2021). Heat Tolerance of Reptile Embryos: Current Knowledge, Methodological Considerations, and Future Directions. J. Exp. Zool. 335, 45–58. doi:10.1002/jez.2402
Huey R. B., Bennett A. F. (1987). Phylogenetic Studies of Coadaptation: Preferred Temperatures versus Optimal Performance Temperatures of Lizards. Evolution 41 (5), 1098–1115. doi:10.1111/j.1558-5646.1987.tb05879.x
Johnson M. T. J., Munshi-South J. (2017). Evolution of Life in Urban Environments. Science 358, eaam8327. doi:10.1126/science.aam8327
Kohlsdorf T., Garland Jr. T., Navas C. A. (2001). Limb and Tail Lengths in Relation to Substrate Usage in Tropidurus Lizards. J. Morphol. 248, 151–164. doi:10.1002/jmor.1026
Kohlsdorf T., Grizante M. B., Navas C. A., Herrel A. (2008). Head Shape Evolution in Tropidurinae Lizards: Does Locomotion Constrain Diet? J. Evol. Biol. 21, 781–790. doi:10.1111/j.1420-9101.2008.01516.x
Kohlsdorf T., Navas C. A. (2006). Ecological Constraints on the Evolutionary Association between Field and Preferred Temperatures in Tropidurinae Lizards. Evol. Ecol. 20, 549–564. doi:10.1007/s10682-006-9116-x
Kohlsdorf T., Navas C. (2012). Evolution of Form and Function: Morphophysiological Relationships and Locomotor Performance in Tropidurine Lizards. J. Zool. 288, 41–49. doi:10.1111/j.1469-7998.2012.00918.x
Lapiedra O., Chejanovski Z., Kolbe J. J. (2017). Urbanization and Biological Invasion Shape Animal Personalities. Glob. Chang. Biol. 23, 592–603. doi:10.1111/gcb.13395
Lazić M. M., Carretero M. A., Zivković U., Crnobrnja-Isailović J. (2017). City Life Has Fitness Costs: Reduced Body Condition and Increased Parasite Load in Urban Common Wall Lizards, Podarcis Muralis. Salamandra 53 (1), 10–17.
Lazić M. M., Carretero M. A., Crnobrnja-Isailović J., Kaliontzopoulou A. (2015). Effects of Environmental Disturbance on Phenotypic Variation: an Integrated Assessment of Canalization, Developmental Stability, Modularity and Allometry in Lizard Head Shape. Am. Nat. 185 (1), 44–58. doi:10.1086/679011
Leirião L., Piantoni C., Ribeiro P. L., Navas C. A. (2019). Independent Influence of Thermoregulatory Cost on the Lower and Upper Set-Points of a Heliothermic Lizard. Behav. Process. 164, 17–24. doi:10.1016/j.beproc.2019.04.004
Li S. R., Hao X., Wang Y., Sun B. J., Bi J. H., Zhang Y. P., et al. (2018). Female Lizards Choose Warm, Moist Nests that Improve Embryonic Survivorship and Offspring Fitness. Funct. Ecol. 32, 416–423. doi:10.1111/1365-2435.12995
Longhini L. S., Porto L. S., Rocha A. C. G., Bícego K. C., Klein W., Gargaglioni L. H. (2019). Seasonal Variation of Hypoxic and Hypercarbic Ventilatory Responses in the Lizard Tropidurus Torquatus. Comp. Biochem. Physiol. A Mol. Integr. Physiol. 237, 110534. doi:10.1016/j.cbpa.2019.110534
Maia-Carneiro T., Navas C. A. (2021). Ecological Constraints to Match Field and Preferred Temperatures in Lizards Tropidurus Catalanensis (Squamata; Tropiduridae). J. Therm. Biol. 98, 102903. doi:10.1016/j.jtherbio.2021.102903
McDonnell M. J., Hahs A. K. (2015). Adaptation and Adaptedness of Organisms to Urban Environments. Annu. Rev. Ecol. Evol. Syst. 46, 261–280. doi:10.1146/annurev-ecolsys-112414-054258
McKinney M. L. (2008). Effects of Urbanization on Species Richness: A Review of Plants and Animals. Urban Ecosyst. 11, 161–176. doi:10.1007/s11252-007-0045-4
Miles L. S., Carlen E. J., Winchell K. M., Johnson M. T. J. (2021). Urban Evolution Comes into its Own: Emerging Themes and Future Directions of a Burgeoning Field. Evol. Appl. 14, 3–11. doi:10.1111/eva.13165
Moritz C., Agudo R. (2013). The Future of Species under Climate Change: Resilience or Decline? Science 341 (6145), 504–508. doi:10.1126/science.1237190
Noble D. W. A., Senior A. M., Uller T., Schwanz L. E. (2021). Heightened Among‐Individual Variation in Life History but Not Morphology is Related to Developmental Temperature in Reptiles. J. Evol. Biol. 34, 1793–1802. doi:10.1111/jeb.13938
Noble D. W. A., Stenhouse V., Schwanz L. E. (2018). Developmental Temperatures and Phenotypic Plasticity in Reptiles: a Systematic Review and Meta-Analysis. Biol. Rev. 93, 72–97. doi:10.1111/brv.12333
Oksanen J. F., Simpson G. L., Blanchet F. G., Kindt R., Legendre P., Minchin P. R., et al. (2019). Vegan: Community Ecology Package. R package version 2.5-6.
Piantoni C., Curcio F. F., Ibargüengoytía N. R., Navas C. A. (2019). Implications of Climate Change on the Habitat Shifts of Tropical Lizards. Austral Ecol. 44, 1174–1186. doi:10.1111/aec.12795
Putman B. J., Gasca M., Blumstein D. T., Pauly G. B. (2019). Downsizing for Downtown: Limb Lengths, Toe Lengths, and Scale Counts Decrease With Urbanization in Western Fence Lizards (Sceloporus occidentalis).. Urban Ecosyst. 22, 1071–1081. doi:10.1007/s11252-019-00889-z
Py-Daniel T. R., De-Lima A. K. S., Lima F. C., Pic-Taylor A., Pires Junior O. R., Sebben A. (2017). A Staging Table of Post-Ovipositional Development for the South American Collared Lizard Tropidurus Torquatus (Squamata: Tropiduridae). Anat. Rec. 300, 277–290. doi:10.1002/ar.23500
Rand A. S., Rand P. J. (1966). Aspects of the Ecology of the Iguanid Lizard Tropidurus Torquatus at Belém, Pará. Smithson. Misc. Collect. 151, 1–16. Available at: https://www.biodiversitylibrary.org/bibliography/8128.
Refsnider J. M., Clifton I. T., Vazquez T. K. (2019). Developmental Plasticity of Thermal Ecology Traits in Reptiles: Trends, Potential Benefits, and Research Needs. J. Therm. Biol. 84, 74–82. doi:10.1016/j.jtherbio.2019.06.005
Refsnider J. M., Janzen F. J. (2010). Putting Eggs in One Basket: Ecological and Evolutionary Hypotheses for Variation in Oviposition-Site Choice. Annu. Rev. Ecol. Evol. Syst. 41, 39–57. doi:10.1146/annurev-ecolsys-102209-144712
Rivkin L. R., Santangelo J. S., Alberti M., Aronson M. F. J., de Keyzer C. W., Diamond S. E., et al. (2019). A Roadmap for Urban Evolutionary Ecology. Evol. Appl. 12, 384–398. doi:10.1111/eva.12734
Rodrigues A. V., Kohlsdorf T. (2019). Learning Skills in Tropidurus Lizards Are Associated with Territory Harshness. J. Zool. 309, 250–258. doi:10.1111/jzo.12721
Rossigalli-Costa N., Barros F. C., Cipriano A. P., Prandini L. P., Medeiros T. A., Rothier P. S., et al. (2021). A Guide to Incubate Eggs of Tropidurus Lizards under Laboratory Conditions. J. Exp. Zool. B Mol. Dev. Evol. 336 (7), 576–584. doi:10.1002/jez.b.23095
Russell V. L. (2021). Emmeans: Estimated Marginal Means, Aka Least-Squares Means. R package version 1.5.5-1. Available at: https://CRAN.R-project.org/package=emmeans (Accessed January 2022).
Sanger T. J., Hime P. M., Johnson M. A., Diani J., Losos J. B. (2008). Laboratory Protocols for Husbandry and Embryo Collection of Anolis Lizards. Herpetol. Rev. 39, 58–63. Available at: https://lososlab.oeb.harvard.edu/files/lososlab/files/sanger_et_al_herp_review_.pdf.
Shine R., Harlow P. S. (1996). Maternal Manipulation of Offspring Phenotypes via Nest-Site Selection in an Oviparous Lizard. Ecology 77, 1808–1817. doi:10.2307/2265785
Shingleton A. (2010). Allometry: The Study of Biological Scaling. Nat. Educ. Knowl. 3 (10), 2. Available at: https://www.nature.com/scitable/knowledge/library/allometry-the-study-of-biological-scaling-13228439/.
Singh S. K., Das D., Rhen T. (2020). Embryonic Temperature Programs Phenotype in Reptiles. Front. Physiol. 11, 35. doi:10.3389/fphys.2020.00035
Sparkman A., Howe S., Hynes S., Hobbs B., Handal K. (2018). Parallel Behavioral and Morphological Divergence in Fence Lizards on Two College Campuses. PLoS ONE 13, e0191800. doi:10.1371/journal.pone.0191800
Steele A. L., Wibbels T., Warner D. A. (2018). Revisiting the First Report of Temperature‐Dependent Sex Determination in a Vertebrate, the African Redhead agama. J. Zool. 306, 16–22. doi:10.1111/jzo.12560
Taylor E. N., Diele‐Viegas L. M., Gangloff E. J., Hall J. M., Halpern B., Massey M. D., et al. (2021). The Thermal Ecology and Physiology of Reptiles and Amphibians: A User's Guide. J. Exp. Zool. 335, 13–44. doi:10.1002/jez.2396
Tiatragul S., Hall J. M., Warner D. A. (2020). Nestled in the City Heat: Urban Nesting Behavior Enhances Embryo Development of an Invasive Lizard. J. Urban Ecol. 6, juaa001. doi:10.1093/jue/juaa001
Tulli M. J., Cruz F. B., Kohlsdorf T., Abdala V. (2016). When a General Morphology Allows Many Habitat Uses. Integr. Zool. 11, 483–499. doi:10.1111/1749-4877.12193
Van Sluys M. (1993). Food Habits of the Lizard Tropidurus Itambere (Tropiduridae) in Southeastern Brazil. J. Herpetol. 27, 347–351. doi:10.2307/1565162
Vaughn P. L., Mcqueen W., Gangloff E. J. (2021). Moving to the City: Testing the Implications of Morphological Shifts on Locomotor Performance in Introduced Urban Lizards. Biol. J. Linn. Soc. 134, 141–153. doi:10.1093/biolinnean/blab076
Wang Y., Zeng Z.-G., Li S.-R., Bi J.-H., Du W.-G. (2016). Low Precipitation Aggravates the Impact of Extreme High Temperatures on Lizard Reproduction. Oecologia 182, 961–971. doi:10.1007/s00442-016-3727-x
Warner D. A., Andrews R. M. (2002). Nest-site Selection in Relation to Temperature and Moisture by the Lizard Sceloporus Undulatus. Herpetologica 58, 399–407. Available at: https://www.jstor.org/stable/pdf/3893295.pdf. doi:10.1655/0018-0831(2002)058[0399:nsirtt]2.0.co;2
Warner D. A., Shine R. (2008). Maternal Nest-Site Choice in a Lizard with Temperature-Dependent Sex Determination. Anim. Behav. 75, 861–870. doi:10.1016/j.anbehav.2007.07.007
Werneck F. P., Leite R. N., Geurgas S. R., Rodrigues M. T. (2015). Biogeographic History and Cryptic Diversity of Saxicolous Tropiduridae Lizards Endemic to the Semiarid Caatinga. BMC Evol. Biol. 15, 94. doi:10.1186/s12862-015-0368-3
Wickham H. (2016). ggplot2: Elegant Graphics for Data Analysis. New York: Springer-Verlag. ISBN 978-3-319-24277-4. Available at: https://ggplot2.tidyverse.org.
Winchell K. M., Carlen E. J., Puente‐Rolón A. R., Revell L. J. (2018a). Divergent Habitat Use of Two Urban Lizard Species. Ecol. Evol. 8, 25–35. doi:10.1002/ece3.3600
Winchell K. M., Maayan I., Fredette J. R., Revell L. J. (2018b). Linking Locomotor Performance to Morphological Shifts in Urban Lizards. Proc. R. Soc. B 285, 20180229. doi:10.1098/rspb.2018.0229
Keywords: development, eco-morphology, temperature, thermal physiology, Tropiduridae, urbanization
Citation: Rossigalli-Costa N and Kohlsdorf T (2022) Native Lizards Living in Brazilian Cities: Effects of Developmental Environments on Thermal Sensitivity and Morpho-Functional Associations of Locomotion. Front. Physiol. 13:891545. doi: 10.3389/fphys.2022.891545
Received: 07 March 2022; Accepted: 17 June 2022;
Published: 15 July 2022.
Edited by:
Turk Rhen, University of North Dakota, United StatesReviewed by:
Cayleih E. Robertson, McMaster University, CanadaHong-Liang Lu, Hangzhou Normal University, China
Miguel A. Carretero, Centro de Investigacao em Biodiversidade e Recursos Geneticos (CIBIO-InBIO), Portugal
Copyright © 2022 Rossigalli-Costa and Kohlsdorf. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Tiana Kohlsdorf, tiana@usp.br