key: cord-0072817-9n7g4dy2 authors: Mühlenhaupt, Max; Baxter-Gilbert, James; Makhubo, Buyisile G.; Riley, Julia L.; Measey, John title: No evidence for innate differences in tadpole behavior between natural, urbanized, and invasive populations date: 2022-01-05 journal: Behav Ecol Sociobiol DOI: 10.1007/s00265-021-03121-1 sha: a3c159f4b8421a20cde79b064d71af1f1832f4d4 doc_id: 72817 cord_uid: 9n7g4dy2 ABSTRACT: Animals are increasingly challenged to respond to novel or rapidly changing habitats due to urbanization and/or displacement outside their native range by humans. Behavioral differences, such as increased boldness (i.e., propensity for risk-taking), are often observed in animals persisting in novel environments; however, in many cases, it is unclear how these differences arise (e.g., through developmental plasticity or evolution) or when they arise (i.e., at what age or developmental stage). In the Guttural Toad (Sclerophrys gutturalis), adult urban toads from both native and invasive ranges are bolder than conspecifics in natural habitats. Here, we reared Guttural Toad tadpoles in a common garden experiment, and tested for innate differences in boldness across their development and between individuals whose parents and lineage came from rural-native, urban-native, and urban-invasive localities (i.e., origin populations). Tadpoles did not differ in their boldness or in how their boldness changed over ontogeny based on their origin populations. In general, tadpoles typically became less bold as they aged, irrespective of origin population. Our findings indicate that differences in boldness in free-living adult Guttural Toads are not innate in the tadpole stage and we discuss three possible mechanisms driving phenotypic divergence in adult boldness for the focus of future research: habitat-dependent developmental effects on tadpole behavior, decoupled evolution between the tadpole and adult stage, and/or behavioral flexibility, learning, or acclimatization during the adult stage. SIGNIFICANCE STATEMENT: To determine if animals can persist in urban areas or become invasive outside their native ranges, it is important to understand how they adapt to life in the city. Our study investigates if differences in boldness that have been found in adult Guttural Toads (Sclerophrys gutturalis) represent heritable differences that can also be found in early life stages by rearing tadpoles from eggs in a common garden experiment. We did not find any differences in boldness among tadpoles from rural-native, urban-native, and urban-invasive origin populations. Our findings suggest that differences in boldness are not innate and/or that boldness is a behavioral trait that is decoupled between the tadpole and the adult stage. Anthropogenic global change poses a host of novel challenges to species worldwide (Grimm et al. 2008; Martin et al. 2010; Pautasso et al. 2010) . Two important components of direct or indirect human impact on animals are the urbanization of natural environments and the displacement of animals into non-native ranges (Vitousek et al. 1997; Lambin and Meyfroidt 2011) . Interestingly, invasive species distributions are often centered around cities and other anthropogenically modified habitats (Pyšek et al. 2020 ). This might simply be a result of the more frequent transportation to, and release within, human-modified habitats due to human activity and global inter-connectedness that is inherent to the concept of invasive populations (Hulme 2009; Gaertner et al. 2017; Padayachee et al. 2017) . However, recent evidence suggests that the invasiveness of species may be enhanced by prior adaptation to urban environments before they are subsequently released into novel urban ranges, bolstering a population's ability to establish and spread (Hufbauer et al. 2012; Borden and Flory 2021) . Behavioral flexibility is a key mechanism that can enable animals to cope with rapid human-induced environmental changes and can even allow them to thrive in heavily modified habitats like cities (Candolin and Wong 2012; Sol et al. 2013; van Baaren and Candolin 2018) . Commonly reported behavioral shifts in response to novel anthropogenically driven conditions, such as urban landscapes or extra-limital introductions, include increases in activity (e.g., Pintor et al. 2008; Pintor and Sih 2009; Philpott et al. 2019) , aggression (e.g., Pintor et al. 2008; Evans et al. 2010; Kralj-Fišer and Schneider 2012) , boldness (e.g., Atwell et al. 2012; Batabyal et al. 2017; Hewes and Chaves-Campos 2018; Damas-Moreira et al. 2020; Mazza et al. 2020) , and tendency to explore (e.g., Thompson et al. 2018; Damas-Moreira et al. 2019; Oliveira et al. 2020) , but can vary between taxa and locations. Research on why these behavioral shifts occur indicates that urban-adapted and/or invasive populations can achieve higher growth rates and reproductive output in anthropogenically modified habitats (Peach et al. 2008; Pintor and Sih 2009 ; Sargent and Lodge 2014; Putman and Tippie 2020; Thawley and Kolbe 2020) by using novel ecological opportunities (e.g., González-Bernal et al. 2012; González-Bernal et al. 2016; Thawley and Kolbe 2020) and outcompeting native species (e.g., Pintor and Sih 2009; Damas-Moreira et al. 2020) while benefitting from reduced predation pressure (Brownscombe and Fox 2013; Eötvös et al. 2018; Westrick et al. 2019 ) and higher resource availability (Lowry et al. 2013; Tryjanowski et al. 2015; Hradsky et al. 2017) ; however, many questions remain regarding how these behavioral adaptations arise. Similar behaviors appear to be beneficial for animals in urbanized or invasive populations (e.g., increased boldness; Atwell et al. 2012; Damas-Moreira et al. 2020; Mazza et al. 2020 ), yet the mechanism driving these behavioral tendencies is often unclear. Anthropogenically derived phenotypic change may be driven by prior adaptation in the donor habitat (Hufbauer et al. 2012; Borden and Flory 2021) , rapid evolution in the novel environment (Whitney and Gabler 2008; Baxter-Gilbert et al. 2019) , or realization through developmental plasticity and behavioral flexibility (e.g., Hudson et al. 2017; Hewes and Chaves-Campos 2018; Mazza et al. 2020) . Furthermore, exactly when trait divergence occurs during ontogeny is largely unknown (e.g., Brodin et al. 2013; Gruber et al. 2017 ). This question is especially relevant when different ontogenetic stages of an animal exhibit different morphologies and use diverse ecological niches and habitats (Werner 1988; Kelleher et al. 2018 ). This is more commonplace than one may expect-approximately 80% of all animals, including insects and amphibians, have these types of complex life cycles (Werner 1988) . We aimed to study evolutionary and ontogenetic divergence of a behavioral trait (i.e., boldness; herein defined as an animal's propensity to exhibit risk-taking behavior) known to differ among populations of an urban-invasive amphibian, the Guttural Toad (Sclerophrys gutturalis). A recent study by Baxter-Gilbert et al. (2021) found that adult toads are bolder in urban habitats than they are in natural areas across native and invasive ranges. This study gives rise to the questions of how (i.e., through evolution or plasticity) and when (i.e., at what age or developmental stage) these differences arise. In a common garden experiment, we raised tadpoles produced from wild adults captured from rural-native, urban-native, and urban-invasive populations (i.e., origin populations). We then quantified the tadpoles' boldness (i.e., latency to leave a hide after a standardized simulated predator attack; Baxter-Gilbert et al. 2021) every 10 days across development to investigate whether differences in boldness between habitat types are genetically fixed, and if so, when do these differences arise ontogenetically. We expected that (1) all tadpoles increase in boldness as they age, due to an increased energy requirement with larger body size, as well as the well-established link between boldness and foraging activity (Wilson and Krause 2012b; Urszán et al. 2015a, b; Kelleher et al. 2018) , and (2) urban-invasive tadpoles will be consistently bolder across ontogeny than rural-native tadpoles, with urban-native tadpoles being intermediary, given the reported coupling of boldness across life stages in anurans (Wilson and Krause 2012b; Kelleher et al. 2018 ). If prediction (2) is supported, then this would also suggest that prior adaptation to invade urban habitats is already present in the tadpole stage of this species (Hufbauer et al. 2012; Borden and Flory 2021) . The Guttural Toad is a large bufonid with a maximum reported snout-vent length (SVL) of 140 mm, that has a wide distribution across much of sub-Saharan Africa (du Preez et al. 2004; Telford et al. 2019) . Tadpoles are aquatic and develop over the course of 6 to 10 weeks, after which metamorphosis is complete and they leave the water with a terrestrial body form (du Preez et al. 2004; Channing et al. 2012 ). This toad is a habitat generalist that is encountered across a range of grassland and forest ecosystems commonly associated with natural wetlands and waterbodies (du Preez et al. 2004; Channing et al. 2012 ). Yet, it also has a proclivity for synanthropic behavior and can be found in anthropogenic landscapes (e.g., urban/suburban) where they will exploit anthropogenic water bodies such as garden ponds for breeding habitat (du Preez et al. 2004; Vimercati et al. 2017) . The species has successfully established invasive populations on the islands of Mauritius and Réunion and near the city of Cape Town, South Africa (Telford et al. 2019; Baxter-Gilbert et al. 2020) . A study by Telford et al. (2019) showed that all invasive populations came from the port city of Durban, South Africa. In the case of the invasive population near Cape Town, eggs and tadpoles are presumed to have been unintentionally introduced via a consignment of aquatic ornamental plants in 1998 (De Villiers 2006 Measey et al. 2017 ). This urban-invasive Cape Town population has since shown rapid adaptation to the drier climate in the Cape region in the adult stage, altering their behavior, physiology, and reproductive strategy (Vimercati et al. 2018 (Vimercati et al. , 2019 Madelaire et al. 2020; Barsotti et al. 2021) . For the purposes of our study, we collected reproductively sized adult toads from November 2019 to February 2020. We chose two locations to collect toads from natural-native origin populations, hereafter referred to as "Durban Rural" (from sites at 29° 51′ 31″ S, 30° 43′ 18″ E and 29° 28′ 17″ S, 31° 13′ 25″ E), which were areas with little human modification that are 20-60 km from the city of Durban. We also collected toads from two urbannative origin populations, hereafter referred to as "Durban Urban" (29° 50′ 55″ S, 31° 00′ 30″ E and 29° 49′ 08″ S, 30° 56′ 37″ E), which were from suburban areas within Durban. Lastly, we collected toads from a single urbaninvasive origin population, hereafter referred to as "Cape Town" (34° 01′ 29″ S, 18° 26′ 03″ E) from the suburb of Constantia, located near Cape Town, where an invasive population of Guttural Toads had become established. The Durban Rural sites consisted of wetlands with a lake or a series of ponds that showed little anthropogenic modification or signs of low levels of use by humans (e.g., nature reserves and infrequent recreational fishing) and had mostly native plant communities. The Durban Urban sites had artificial water bodies (e.g., ornamental ponds with non-native fish species, ditches, or fountains), frequent anthropogenic structures (e.g., buildings and roads), and high levels of human presence. Similarly, the Cape Town site also had many garden ponds, buildings, roads, and high levels of human presence. After collection, adults were brought to a greenhouse located at the University of KwaZulu-Natal Westville Campus. In this facility, toads were housed by sampling location and sex in large plastic tubs (1.1 m [L] × 1.3 m [W] × 0.50 m [H]; n = 10) filled with ~50 mm of soil and another ~50 mm of leaf litter over the soil. Water was provided in plastic bowls ad libitum. We ensured the soil was moist, but not saturated, and fed toads 2-5 adult crickets (Acheta domesticus) every other day. Toad breeding was carried out from December 2019 to March 2020 using individuals that were available in the greenhouse at the time. Thus, clutches were not deposited synchronously, but rather in clumped batches over the course of this ~3-month period. Every toad was used only once for breeding. To initiate mating behavior, and subsequent oviposition and fertilization, we used a 1:20 solution of the synthetic gonadotrophin leuroprorelin acetate-"Lucrin" (Lucrin Depot, Abbott) and Ringer's solution. We injected females with 0.666 mL and males with 0.333 mL (Hamilton et al. 2005) at ~17:00 h. After injection, pairs of males and females were placed inside separate plastic containers (260 × 410 mm) filled with 40 mm of aged tap water and left undisturbed overnight. The next morning at ~9:00 h, we checked the containers for fertilized clutches, removed the adults from the containers, and replaced the water inside these containers with clean, aged tap water to a height of approximately 150 mm. We regularly replaced the water inside these containers and made sure clutches were well aerated. Over the next few days, we followed embryonic development closely. Typically, 1 day after the clutch was laid, most embryos had already hatched (i.e., Gosner stage 20; Gosner 1960) . At this point, individuals were assigned a tadpole age of 0 days, allowing us to track age in days for the entirety of our common garden experiment. We housed tadpoles in large mesocosms (large plastic tubs; 1.1 m [L] × 1.3 m [W] × 0.5 m [H]; n = 30) filled with 600 L of tap water placed halfway below workbenches inside the facility, so that one half of a mesocosm was shaded and the other half received sunlight through the roof panels. We let the water inside the mesocosms age for a week, before adding ~1 L of water from a standing water tank to each mesocosm to induce zoo-and phyto-plankton establishment. This water tank was located inside the facility (i.e., preventing access from toads and fish) and had a well-established community of aquatic vegetation and algae within it. After another week, we added 50 g of rabbit chow (Rabbit Pellets, Westerman's Premium) to supply the plankton with additional nutrients (Semlitsch and Boone 2009) . At an age of 3 days, we separated 1000 enumerated tadpoles from one clutch (i.e., full siblings) and translocated them to a mesocosm that had been prepared for their arrival approximately 1 month before and contained a well-established and abundant zoo-and phyto-plankton community. At the age of 10 days, and then every subsequent 10 days (Table 1) , we haphazardly selected 20 tadpoles from each mesocosm for testing. If we were not able to conduct the assay on the specific day the tadpoles were scheduled for, we conducted the assay on the next possible day (see Table 1 ). We did not assess Gosner stage (Gosner 1960) of each tadpole as this would have required additional handling. Generally, developmental stages of individual tadpoles varied greatly within a clutch (e.g., between Gosner stages 25 and 30 for younger clutches and Gosner stages 30 and 41 for older clutches; Gosner 1960) . However, for older tadpoles, we ensured that only tadpoles that had not developed past Gosner stage 41 (Gosner 1960) were used, as strong morphological and physiological changes might affect behavior in subsequent developmental stages (Dodd and Dodd 1976) . We used 190 [L] × 300 [W] mm plastic containers filled with 20 mm of aged tap water as boldness arenas. Above the boldness arenas, we installed digital video cameras (recording at a frame rate of 25 fps and a resolution of 101.38 kP), that allowed us to remotely record the tadpole behavior (i.e., without a researcher present). Before the start of each boldness assay, tadpoles were tested for swimming performance (Mühlenhaupt et al. 2021) and were introduced to the boldness arenas immediately afterwards. Before starting the simulated predator attack, we documented the water temperature in the boldness arenas using a digital thermometer (Control Company, product code: 4052; ± 0.1°C). Using the tip of a red silicon spatula, a single researcher (MM) chased the tadpoles for 10 s before placing a circular hide (i.e., a white opaque plastic cup with a diameter of 116 mm with a 50 [W] × 30 [H] mm exit hole cut into the side) over them. The hide, while it was surrounding the tadpole, was then centered within the arena and the exit hole aligned to be visible by the camera. For the next 20 min, we recorded tadpole behavior. From these video recordings, we scored whether tadpoles left the hide during the duration of recording (binary variable) and documented the time (s) it took the tadpoles to leave the hide (i.e., latency to leave the hide). For the latter behavior, we considered a tadpole had left the hide when it had passed through the exit hole with the vent exposed and the tip of the snout facing perpendicular to the hide. Video scoring was carried out by one observer (MM) who was blind to the origin populations of tadpoles. Tadpoles that did not leave the hide were assigned the maximum assay time (1200 s; as per Riley et al. 2017; Damas-Moreira et al. 2019; Baxter-Gilbert et al. 2021 ). This protocol closely followed the boldness assay used by Baxter-Gilbert et al. (2021) . We measured the same variables and began each boldness assay with a standardized "mock predation" event (i.e., chasing the tadpoles for 10 s with a rubber spatula) to ensure comparability of our results with those of Baxter-Gilbert et al. (2021). Furthermore, we ensured that boldness, as quantified here, was not conflated with general activity by simulating a risky situation through the "mock predation" event at the beginning of each assay (see Réale et al. 2007; Kelleher et al. 2018) . After the boldness assay, we returned the tadpoles to their respective mesocosms. Thus, we cannot rule out the possibility that any individual tadpole was also sampled for a boldness assay at a different age, and we acknowledge that recurrent behavioral testing can influence tadpole behavior in later trials (Urszán et al. 2015a ). However, given the high initial number of tadpoles per mesocosm (n = 1000), there is Table 1 The number of tadpoles whose boldness was assayed in our study, separated by age and origin population. Numbers in parentheses denote the total number of clutches the tadpoles belonged to. For Durban Rural, 3 and 7 clutches were produced from toads collected in sampling location 1 or 2, respectively, and for Durban Urban, 6 and 6 clutches from sampling location 1 or 2, respectively (see "Methods"). Please note that due to the staggered breeding of adults over a period of 3 months, as well as the onset of COVID-19 pandemic restrictions in South Africa in March 2020, sample sizes are lower for higher ages as we were not able to continue the data collection after this time a low likelihood that many tadpoles were repeatedly used for testing and thus, we conducted statistical analysis assuming non-repeated measures. All statistical analyses were carried out using R (Version 4.0.1; R Core Team 2020). Before analyzing our data, we explored the data set following Zuur et al. (2010) and did not find any unexplainable outliers. We verified the assumptions of normality and homoscedasticity of residuals for all models prior to interpretation and present data as predicted means ± standard error (SE) in the text. The significance level (α) was set at 0.05 for all models. The water temperature during the boldness trials and the age of tadpoles were significantly correlated (confirmed by using the "lm" and "Anova" function in the R "stats" package; R Core Team 2020: F 1, 2493 = 11.12, p < 0.01). Therefore, we decided to exclude temperature from our models because these two factors are confounded. However, due to our staggered breeding protocol (see "Husbandry of adults and breeding protocols"), younger and older tadpoles were tested at the end of the experiment, so seasonal variation in temperature likely did not affect our results strongly. First, we investigated the proportion of tadpoles that did not leave the hide during the 20-min trial. The variability in this behavioral trait was very low (overall less than 1% of tadpoles did not exit the hide; Table 2 ). Thus, we present a summary of this behavioral trait, but it was not analyzed statistically. Instead, we used the tadpole's latency to leave the hide (s) as our boldness metric and, before analyses, we rank-transformed these values because they were strongly positively skewed (Riley et al. 2017; Damas-Moreira et al. 2019) . We used a linear mixed effect model (LMM; using the "lmer" function in the R package "lme4"; Bates et al. 2015) with our boldness metric as the response variable that included the fixed effects of tadpole origin population (categorical: Cape Town, Durban Rural, Durban Urban), age (continuous), and initially included the interaction between origin population and age. This LMM also included a random intercept and slope of mesocosm number to control for dependency in our data due to relatedness and shared housing for individuals from the same tub, and a random intercept of parentage location to control for dependency in our data due to similarities between specific study sites where the parents were collected. Post hoc, we tested for an interaction effect using "emtrends" in the R package "emmeans" (Lenth et al. 2018 ) and did not find any significant effects (all p corr > 0.16; p-values were corrected using a "Tukey" adjustment). Thus, we removed the interaction term and reran the LMM to allow interpretation of main effects. We also ran a separate LMM investigating differences in the boldness of tadpoles that were 40/41 days old among the different origin populations. This LMM had our boldness metric as the response variable, included tadpole origin population (categorical: Cape Town, Durban Rural, Durban Urban) as a fixed factor, and a random intercept of mesocosm number and parentage site. We chose this age because tadpoles had sufficient time to develop and possibly diverge in boldness, sample size was still reasonably high (Table 1) , and, around this age, Guttural Toad tadpoles are known to begin emergence as toadlets (du Preez et al. 2004; Channing et al. 2012) . Furthermore, our previous work has indicated that tadpoles from different origin populations show differences in developmental rate (Mühlenhaupt et al. 2021 ). Thus, with higher age (e.g., 50 or 60 days), our results might be biased towards slow-developing Durban Urban tadpoles because fast-developing Durban Urban tadpoles might have already exceeded Gosner stage 41 (Gosner 1960) . Following analysis of both LMMs, we conducted post hoc multiple comparisons between tadpole origin populations using the R package "emmeans" and p-values generated for these comparisons were "Tukey" adjustment-corrected using a Monte Carlo method to produce "exact" Tukey corrections (Lenth et al. 2018) . The proportion of tadpoles that did not leave the hide during our boldness assay was very low, with less than 1% of tadpoles remaining in the hide during the trials (Table 2) . Tadpoles became less bold (i.e., shyer) as they aged, but boldness did not differ between origin populations (Table 3 , Fig. 1) . Similarly, at the age of 40/41 days, tadpole boldness did not differ among origin populations (Table 4 , Fig. 2 ). Table 2 A summary of the number of tadpoles that did not leave the hide during the 20-min boldness assays separated by age and origin population. Numbers in parentheses denote the specific percentage (%) of tadpoles that did not leave the hide as calculated from the number of tadpoles that did not leave the hide divided by the total number of tadpoles assayed at that age and specific origin population (see Table 1 Contrary to our predictions, tadpole boldness did not differ among rural-native, urban-native, and urban-invasive origin populations, and decreased with tadpole age in our common garden experiment. Given the clear divergence in boldness among natural and urban adult Guttural Toads (Baxter- , this study gives valuable insight in how urban-derived boldness in Guttural Toads may or may not arise. Notably, this would suggest that the increased boldness observed in adult urban toads is either a result of experience or plasticity at either the tadpole or adult stage (i.e., requiring exposure to actual urban/rural ecosystems), or if it is heritable and fixed, then the phenotype is not expressed until after metamorphosis. Within the scope of this study, since bolder phenotypes are not uniquely arising in the larval stage of urbanized Guttural Toads, we would assert that prior-adaptation is not enhancing the species' invasion success (Hufbauer et al. 2012; Borden and Flory 2021) , with respect to tadpole boldness. In the context of the recent surge in studies examining rapid evolution in urban and invasive populations, our study results demonstrate a break with many significant findings on larval stages, suggesting substantial gaps in our understanding of the mechanisms and timing by which these phenotypic differences arise, and support the recommendation for more robust methodologies (see Lambert et al. 2020 ). Our common garden experimental design gives us the ability to demonstrate strong evidence that Guttural Toad tadpoles are not expressing heritable or innate differences in their boldness between populations (Lambert et al. 2020) . These findings are contrary to previous research which showed heritable differences in tadpole boldness between mainland and island populations of Common Frogs (Rana temporaria) that were carried over to the adult stage (Brodin et al. 2013) . In Common Frogs, the bolder phenotype of island populations was presumed to be the result of a disperser-dependent founder effect through bolder, more exploratory adults colonizing the new islands that rose above sea level approximately 350-800 years ago (Brodin et al. 2013) , which is similar to the increased levels of boldness seen in invasive range-front (i.e., expanding) adult Cane Toads (Rhinella marina) in Australia after less than 100 years (Gruber et al. 2017 ). In the context of Guttural Toads, however, it is more likely that a population around the area of Durban already existed before the city was founded and became urbanized during the city's founding and growth . Furthermore, the Cape Town population of Guttural Toads was likely established through tadpoles or eggs, both of which are non-dispersive life stages (De Villiers 2006; Measey et al. 2017 ). Thus, a disperserdependent founder effect is unlikely in the urban and urbaninvasive populations of the Guttural Toad, emphasizing the context-dependent nature of behavioral divergence between populations. Similarly, a previous study showed that urbannative and urban-invasive tadpoles do not differ in swimming performance and underlying morphological traits compared to rural-native tadpoles in the Guttural Toad (Mühlenhaupt et al. 2021) , suggesting that predation pressure might not be reduced in urban and/or invasive tadpole habitats. Although our findings suggest that urbanized or invasive Guttural Toad tadpoles are not exhibiting heritable increases in boldness, the drivers of the observed differences in adult boldness between urban and rural adults (Baxter-Gilbert et al. 2021) remain unknown. One possibility is that phenotypic differences could be the result of developmental effects experienced during the tadpole stage, which may be habitat-dependent for their expression in tadpoles and are carried over and expressed during post-metamorphic stages. Tadpoles are well known to show high levels of plasticity in many traits as a response to different ecological factors (Richardson 2001; Relyea 2002; Van Buskirk 2002; Relyea 2004) and differences in stress response between urban, agricultural, and rural Common Toad (Bufo bufo) tadpoles are seen to be a result of phenotypic plasticity (Bókony et al. 2021) . Furthermore, more active, explorative, and bolder phenotypes in Lake Frog (Pelophylax ridibundus) tadpoles are seen to be carried over to the adult stage (Wilson and Krause 2012b) and if these traits arise through plasticity, then variable environments would be required to trigger their expression. As such, due to the controlled nature of our common garden experimental design, if the bolder phenotypes are a result of plasticity, then it is not surprising that we were unable to detect them in this study. A reciprocal transplant experiment might help to elucidate what role phenotypic plasticity of boldness and the evolution of plasticity in this trait play in shaping the bolder urban phenotype of adult toads . For example, plasticity in Fig. 2 Guttural Toad tadpole boldness (i.e., rank-transformed latency to leave the hide after a standardized simulated predation attempt) at the age of 40/41 days. Depicted are boxplots of predicted data from a linear mixed effect model testing for differences between Durban Rural (green), Durban Urban (black), and Cape Town (blue) tadpoles. The thick horizontal line in the boxplots corresponds with the predicted median, the boxes represent the second and third quartile ranges, and the whiskers encompass the minimum and maximum of the data. We did not find any significant differences in boldness among tadpoles from different origin populations (Table 2 ). The blue triangle on the left depicts the directionality of the trait (i.e., low values correspond with bold tadpoles and vice versa) boldness might have evolved only in urban populations, and then may only be observable if the correct stimuli are present during development to trigger a shift in boldness (e.g., Iraeta et al. 2006; Kaiser et al. 2019 ). Conversely, differential expression of boldness between adults but not tadpoles might be the result of decoupled evolution between life stages in Guttural Toads (Collet and Fellous 2019). Decoupling of traits has recently been shown in the invasive African Clawed Frog (Xenopus laevis; Kruger et al. 2021) . Adults from the periphery of the invasive range in France show a distinct dispersal-mediating phenotype compared to adults from the core (Louppe et al. 2018; Padilla et al. 2019 ). However, when reared from the egg to the post-metamorphic frog, those differences were not apparent, indicating that differences in adult traits emerge well after metamorphosis (Kruger et al. 2021) . Similarly, in a previous study, we presented data indicating that adult Guttural Toads are larger in urban areas but did not find differences in growth rates of tadpoles when raised in a common garden environment (Mühlenhaupt et al. 2021 ). The amphibian larval stage is typically more sedentary and adapted for foraging and growth, whereas the adult stage is more mobile and adapted for dispersal and reproduction (Wilson and Krause 2012a) . Combined with the far-reaching changes in the environment from the aquatic to the terrestrial habitat (Dodd and Dodd 1976) , differential selective forces are likely acting on both stages and decoupling of traits across metamorphosis could enhance phenotypic match in the environment experienced by the specific ontogenetic stage (Wilbur 1980; Collet and Fellous 2019; Goedert and Calsbeek 2019) . Therefore, if selection has favored increased boldness in urbanized and/or invasive Guttural Toad populations, then the expression of this trait may well be evolutionarily decoupled between the sedentary tadpole life stage and the active, dispersing, and mate-searching adult life stage. Beyond an evolutionary mechanism (e.g., local adaptation), the presence of bolder urbanized adult Guttural Toads ), but the absence of such a trait in tadpoles, may be a result of circumstance, due to behavioral flexibility, learning, and/or acclimatization (e.g., Levey et al. 2009; Wright et al. 2010; Couchoux and Cresswell 2012; Vincze et al. 2016; Evans et al. 2019) . For example, several phenotypic traits of adult invasive Cane Toads are strongly influenced by prior experience and developmental environment, such as thermal acclimation (McCann et al. 2014) or antipredator response (Hudson et al. 2017) . Novel environments, such as cities or the extra-limital ranges of introduced populations, often require animals to phenotypically respond rapidly to novel or stressful ecological conditions (Grimm et al. 2008) . Behavioral flexibility and learning represent some of the quickest mechanisms to adaptively match novel selective regimes (Ghalambor et al. 2007; Lande 2009; Fox et al. 2019) . Such scenarios can be seen with Northern Mockingbirds (Mimus polyglottos), which are capable of rapidly learning and identifying individual humans and adjusting their flight distance based on their prior experience with a specific person (Levey et al. 2009 ), or experience-driven decreases in risk aversion by Burrowing Owls (Athene cunicularia) towards humans (Cavalli et al. 2018) . Therefore, it follows that Guttural Toads in urban habitats could have their behavioral traits incrementally shaped through differential experiences, rather than evolutionary mechanisms, resulting in the expression of a bolder adult phenotype in these populations. Our finding that tadpoles became less bold with increasing age provides clear evidence that boldness levels can change across ontogeny, even within a given life stage, which provides important insights for future research looking to disentangle where the divergence in adult boldness between urban and rural populations occurs . This finding was contrary to our predictions because, as we expected that older and larger tadpoles would spend more time foraging, and as there often is a correlation between foraging activity and boldness, we expected our tadpoles to also be bolder as they aged (Wilson and Krause 2012b; Urszán et al. 2015a, b; Kelleher et al. 2018 ). However, our results aligned with previous research on Common Toad tadpoles, similarly reared in a common garden environment, which were shown to decrease their activity levels (which can correlate with boldness; Kelleher et al. 2018) across larval development (Kurali et al. 2018) . Thus, it is likely that tadpoles early in development (after the consumption of the yolk sac; Gosner 1960) show increased foraging activity and boldness to enhance resource acquisition, which then lessen over time due to consequent growth and energy storage (Ramamonjisoa et al. 2019; Davenport et al. 2020 ). This change in behavioral traits across ontogeny, potentially related to energetics, will also play an important role in interpreting future research on how traits, like boldness, change between life stages and how this may or may not relate to adaptive behavioral shifts associated with urbanization or biological invasions. Overall, we showed that there is no innate difference in boldness among Guttural Toad tadpoles from rural-native, urban-native, and urban-invasive ranges, despite the welldocumented divergence of boldness in adults . Our study highlights the importance of studying the ontogenetic and environmental context of behavioral shifts in novel environments. We encourage researchers, interested in rapid evolution in urban and novel environments, to carefully choose experimental approaches that identify the underlying mechanisms driving population divergence (see Lambert et al. 2020 ). Although this study was able to provide evidence against one potential driver of human-induced behavioral change (i.e., fixed, heritable traits that persist across life stages) within this study system, there remains a number of alternative mechanisms and means from which phenotypes can arise that will require future research. Boldness behavior and stress physiology in a novel urban environment suggest rapid correlated evolutionary adaptation Challenges of a novel range: water balance, stress, and immunity in an invasive toad A multivariate approach to understanding shifts in escape strategies of urban lizards Fitting linear mixedeffects models using lme4 Fortune favors the bold toad: urban-derived behavioral traits may provide advantages for invasive amphibian populations Shrinking before our isles: the rapid expression of insular dwarfism in two invasive populations of guttural toad (Sclerophrys gutturalis) Bold New World: urbanization promotes an innate behavioral trait in a lizard Stressed tadpoles mount more efficient glucocorticoid negative feedback in anthropogenic habitats due to phenotypic plasticity Urban evolution of invasive species Personality trait differences between mainland and island populations in the common frog (Rana temporaria) Living at the edge of the front; reduced predation risk to invasive round goby in a Great Lakes tributary Behavioural responses to a changing world: mechanisms and consequences The role of habituation in the adjustment to urban life: an experimental approach with burrowing owls Tadpoles of Africa: the biology and identification of all known tadpoles in Sub-Saharan Africa. Chimaira, Frankfurt am Main Collet J, Fellous S (2019) Do traits separated by metamorphosis evolve independently? Concepts and methods Personality constraints versus flexible antipredation behaviors: how important is boldness in risk management of redshanks (Tringa totanus) foraging in a natural system Getting ahead: exploitative competition by an invasive lizard Can behaviour explain invasion success? A comparison between sympatric invasive and native lizards The non-consumptive effects of predators and personality on prey growth and mortality Amphibia: Anura: Bufonidae Bufo gutturalis Power, 1927 guttural toad introduced population The biology of metamorphosis Bufo gutturalis Power A meta-analysis indicates reduced predation pressure with increasing urbanization Behavioural syndromes in urban and rural populations of song sparrows Hourly movement decisions indicate how a large carnivore inhabits developed landscapes Beyond buying time: the role of plasticity in phenotypic adaptation to rapid environmental change Non-native species in urban environments: patterns, processes, impacts and challenges Adaptive versus non-adaptive phenotypic plasticity and the potential for contemporary adaptation in new environments Cane toads on cowpats: commercial livestock production facilitates toad invasion in tropical Australia Toads in the backyard: why do invasive cane toads (Rhinella marina) prefer buildings to bushland? A simplified table for staging anuran embryos and larvae with notes on identification Global change and the ecology of cities Is the behavioural divergence between range-core and range-edge populations of cane toads (Rhinella marina) due to evolutionary change or developmental plasticity? Captive care and breeding of marine toads, Bufo marinus Boldness related to size in the hermit crab Coenobita compressus at undisturbed, but not disturbed beach Human-modified habitats facilitate forest-dwelling populations of an invasive predator Evolutionary shifts in antipredator responses of invasive cane toads Anthropogenically induced adaptation to invade (AIAI): contemporary adaptation to human-altered habitats within the native range can promote invasions Trade, transport and trouble: managing invasive species pathways in an era of globalization Mediterranean hatchling lizards grow faster at higher altitude: a reciprocal transplant experiment Behavioural repeatability is affected by early developmental conditions in a butterfly Animal personality and behavioral syndromes in amphibians: a review of the evidence, experimental approaches, and implications for conservation Individual behavioural consistency and plasticity in an urban spider Does the spatial sorting of dispersal traits affect the phenotype of the non-dispersing stages of the invasive frog Xenopus laevis through coupling? Resource-dependent temporal changes in antipredator behavior of common toad (Bufo bufo) tadpoles Adaptive evolution in cities: progress and misconceptions Global land use change, economic globalization, and the looming land scarcity Adaptation to an extraordinary environment by evolution of phenotypic plasticity and genetic assimilation Emmeans: estimated marginal means, aka least-squares means Urban mockingbirds quickly learn to identify individual humans Differences in standard metabolic rate at the range edge versus the center of an expanding invasive population of Xenopus laevis in the West of France Behavioural responses of wildlife to urban environments Challenges of dehydration result in a behavioral shift in invasive toads The effects of anthropogenic global changes on immune functions and disease resistance Small mammals in the big city: behavioural adjustments of non-commensal rodents to urban environments Rapid acclimation to cold allows the cane toad to invade montane areas within its Australian range Invasive amphibians in southern Africa: a review of invasion pathways Growing up in a new world: trait divergence between rural, urban, and invasive populations of an amphibian urban invader Metabolic and behavioral adaptations of greater white-toothed shrews to urban conditions How do invasive species travel to and through urban environments? Allocation trade-offs impact organ size and muscle architecture in an invasive population of Xenopus laevis in Western France Local and landscape drivers of carabid activity, species richness, and traits in urban gardens in coastal California Differences in growth and foraging behavior of native and introduced populations of an invasive crayfish Differences in aggression, activity and boldness between native and introduced populations of an invasive crayfish Big city living: a global meta-analysis reveals positive impact of urbanization on body size in lizards Scientists' warning on invasive alien species Predation decreases cohort foraging activity and growth, yet increases individual size variation in prey Integrating animal temperament within ecology and evolution Competitor-induced plasticity in tadpoles: consequences, cues, and connections to predator-induced plasticity Fine-tuned phenotypes: tadpole plasticity under 16 combinations of predators and competitors A comparative study of activity levels in larval anurans and response to the presence of different predators R: a language and environment for statistical computing. R Foundation for Statistical Computing Early social environment influences the behaviour of a family-living lizard Evolution of invasive traits in nonindigenous species: increased survival and faster growth in invasive populations of rusty crayfish Aquatic mesocosms. In: Dodd CK (ed) Amphibian ecology and conservation: a handbook of techniques Behavioural adjustments for a life in the city Origin of invasive populations of the Guttural Toad (Sclerophrys gutturalis) on Reunion and Mauritius Islands and in Constantia Artificial light at night increases growth and reproductive output in Anolis lizards Urbanization and individual differences in exploration and plasticity Urban and rural habitats differ in number and type of bird feeders and in bird species consuming supplementary food No personality without experience? A test on Rana dalmatina tadpoles Behavioural consistency and life history of Rana dalmatina tadpoles Plasticity in a changing world: behavioural responses to human perturbations A comparative test of the adaptive plasticity hypothesis: relationships between habitat and phenotype in anuran larvae Rapid adaptive response to a Mediterranean environment reduces phenotypic mismatch in a recent amphibian invader Invasive toads adopt marked capital breeding when introduced to a cooler, more seasonal environment Integrating age structured and landscape resistance models to disentangle invasion dynamics of a pond-breeding anuran Habituation to human disturbance is faster in urban than rural house sparrows Introduced species: a significant component of humancaused global change Size, scaling, and the evolution of complex life cycles Rapid evolution and behavioral plasticity following introduction to an environment with reduced predation risk Rapid evolution in introduced species,'invasive traits' and recipient communities: challenges for predicting invasive potential Complex life cycles Metamorphosis and animal personality: a neglected opportunity Personality and metamorphosis: is behavioral variation consistent across ontogenetic niche shifts? Behavioral flexibility and species invasions: the adaptive flexibility hypothesis A protocol for data exploration to avoid common statistical problems We would like to thank NCC Environmental Services (Western Cape), Durban Botanic Gardens, University of KwaZulu-Natal Westville Campus, and Shongweni Dam and Nature Reserve for providing us with toads or access to their facilities for toad collection. Thanks to both N.S. Dludla for assistance during the experiment and P. Kowalski for assistance during fieldwork. We would like to thank J. Jeschke and the members of the MeaseyLab at Stellenbosch University for providing helpful comments and discussions that helped shape this research and manuscript. An anonymous reviewer, B. Carlson, and T. C. M. Bakker helped improve the quality of this manuscript.Funding This work was funded by the DSI-NRF Centre of Excellence for Invasion Biology (JBG & JM) and the Centre for Collaborations in Africa at Stellenbosch University (JBG). Also, during this work, JLR's postdoctoral salary was supplied by the Claude Leon Foundation and Natural Sciences and Engineering Research Council of Canada. The datasets and R code for this study are available from Open Source Framework (OSF) and can be found at https:// doi. org/ 10. 17605/ osf. io/ 9ptcn. Ethical approval Ethical approvals for toad captures and the experimental protocols followed animal ethics guidelines set out and approved by the Stellenbosch University Research Ethics Committee (ACU-2019-13154) and the University of KwaZulu-Natal Animal Research Ethics Committee (AREC/009/020). Field collections in the native range were undertaken under an Ezemvelo KZN Wildlife Ordinary Permit (OP 4072/2019) and a Cape Nature Permit to Export and Transport Protected Wild Animals (CN13-59-12357) for toads from the invasive range. The authors declare no competing interests.