Introduction

Laterality refers to the tendency of animals to use the structures on one side of the body more frequently or efficiently than the other side (Bianki 1988; Bisazza et al. 1998). It relates to the ability of individuals to perform different functions simultaneously, which can influence their fitness (Rogers 2000; Rogers and Andrew 2002; Ferrari et al. 2015), for example, feeding while being attentive to the presence of predators (Rogers et al. 2004), handling food while holding a tool (McGrew and Marchant 1992), and swimming or flying while sleeping (Mascetti 2016). Laterality is associated with a functional difference between the brain’s hemispheres (Estévez-González 1992; Vallortigara et al. 1999). The right hemisphere of the brain is associated with rapid behavioural responses that individuals should have when facing stimuli such as predators or competitors’ attacks (Davidson 1992; Forrester et al. 2018). On the contrary, the left hemisphere of the brain is associated with analytical functions such as tool use, remembering cues in the habitat, and attacking strategies (Robins 2006; Forrester et al. 2018). In general, laterality reflects the evolutionary history of the species and multiple ecological factors that individuals experience during ontogeny (Rogers et al. 2013; Ferrari et al. 2015; Güntürkün et al. 2020), and its study is one of the diverse approaches to understanding the function of the animal mind and its cognitive capacity (Brown and Magat 2011; Rogers 2021).

Due to selective pressure, the predominant use of one side of the body by individuals can be significantly higher or lower than expected by chance (Bisazza et al. 1998; Frasnelli and Vallortigara 2018). The strength of laterality in animals (e.g., the percent of right-handed individuals in a population) could be influenced by the characteristics of their microhabitat and the predominant locomotion movements (Malashichev 2006, 2016). Specifically, the size, slope, and type of substrate that individuals use can affect the kinematics of their movements, thus selecting for specific morphologies, behaviors, and adaptive locomotion patterns (Losos 2009; Herrel et al. 2013; Moen et al. 2016; Fratani et al. 2020). Concerning locomotion type, asynchronous (i.e., non-simultaneous) movements of limbs such as walking and climbing, involve alternating muscle contractions and are associated with higher neural asymmetry compared to synchronous (i.e., simultaneous) movements of limbs, such as swimming and hopping (Syme and Josephson 2002; Malashichev 2006). Therefore, it is expected that species with asynchronous movement exhibit a higher strength of laterality than species with synchronous movements of their limbs (Giljov et al. 2015; Malashichev 2016).

On account of the high diversity of ecological characteristics and life history traits in anurans, Malashichev and Robins (2018) used anuran species as a model to test the expected relationships between the strength of laterality, microhabitat preferences, and locomotion type in animals (Fig. 1). These authors hypothesized that if anurans are aquatic, they tend to exhibit synchronous movements of their limbs during swimming (but see Abourachid and Green 1999 for some exceptions); if anurans are terrestrial, they show both synchronous (hopping) and asynchronous (walking) movements of their limbs; and if anurans are arboreal, they tend to have asynchronous movements of their limbs (walking, climbing) (Abourachid and Green 1999; Malashichev 2006; Sullivan 2018; Fig. 1). Therefore, anurans living in aquatic microhabitat may show less strength of laterality than anurans living in terrestrial and arboreal microhabitats (Malashichev and Robins 2018). Fossorial species exhibit a different pattern, displaying synchronous movements in the forelimbs, but asynchronous movements in the hindlimbs during dragging (Malashichev 2006; Malashichev and Robins 2018). In fact, in the fossorial species Ceratophrys ornata and Pelobates fuscus there is not laterality in the use of forelimbs, but in hindlimbs (Malashichev 2006). Regarding semi-arboreal species, some of them may be predominantly hoppers, while others are walkers or climbers (Duellman 2001, 2005; Vargas-Salinas et al. 2019).

Fig. 1
figure 1

Graphical representation of the hypothesized relationship between the strength of laterality (percentage of right-handed individuals) in the snout-wiping test for anuran species, according to microhabitat (aquatic, semi-aquatic, fossorial, terrestrial, semi-arboreal, arboreal) and predominant locomotion type (synchronous, asynchronous). The figure is adapted from Malashichev and Robins (2018). The grey areas enclose species belonging to each microhabitat type. The colored lines on the right side of the panel indicate the locomotion types that can be present in species using a given microhabitat type. For example, the microhabitat “Terrestrial” can include species that are predominantly walkers or hoppers. The triangle on the left side of the panel represents the hypothesized highest strength of laterality at the base (toward species with asynchronous movement of limbs) and a reduction in the strength of laterality towards the top (species with synchronous movement of limbs). Black triangles indicate the strength of laterality for 14 species according to published literature (Bisazza et al. 1996; Rogers 2002; Malashichev 2006; Robins and Rogers 2006; Sovrano 2007; Sullivan 2018; Robins and Rogers 2021). Open red circles indicate the strength of laterality for 19 Neotropical species (this study)

The model proposed by Malashichev and Robins (2018) offers testable predictions about the strength of laterality in anurans; however, there are limitations to it. First, it is supported by only 12 anuran species, most of them from temperate regions (Fig. 1). Second it does not include body size, which is a pivotal phenotypic trait in many ecological and evolutionary aspects of anurans (Wells 2007; Amado et al. 2019). In the present study, we evaluated the strength of laterality in forelimbs of 20 neotropical anuran species. Consistent with the model by Malashichev and Robins (2018), we expected that terrestrial and arboreal species with asynchronous movements would exhibit greater strength of laterality in forelimbs compared to aquatic species. In addition, we pooled our data with those from previous studies on anurans (e.g., Bisazza et al. 1996; Malashichev 2006; Malashichev and Deckel 2006; Robins and Rogers 2021) to test a possible relationship between strength of laterality in forelimbs, with microhabitat differences, locomotion type and body size. We are aware that the absence of laterality in forelimbs does not preclude laterality in other body parts and behaviours (Rogers 2014; Güntürkün et al. 2020). However, investigating relationships regarding forelimb laterality is an important step for understanding how behaviour and habitat use relates to the cognitive capacity of anurans and other vertebrates.

Materials and methods

Data collection

Between August 2019 and April 2023, we conducted 35 field trips to 10 localities in Colombia (seven in the central Andes, one in the Pacific region, one in the Sierra Nevada of Santa Marta in the Caribbean region, and one in the Valle of Magdalena Medio). These localities are characterized by their humid habitats and high diversity of anurans, including aquatic, semi-aquatic, fossorials, terrestrial, semi-arboreal and arboreal species. We collected data from all available microhabitats. In each locality, we searched randomly for anurans during diurnal (06:30 to 10:30 and 16:00 to 17:30) and nocturnal (18:30 to 02:00) visual encounter surveys without restrictions (Heyer et al. 1994; Angulo et al. 2006). Species identification was based on morphology using local/regional field guides and taxonomic descriptions (e.g., Vargas-Salinas and Bolaños-Lizalda 1999; Acosta-Galvis et al. 2006; Bolívar et al. 2010; Domínguez et al. 2010; Pérez-González et al. 2016; Vanegas-Guerrero et al. 2016; Vargas-Salinas et al. 2019). We followed the taxonomic nomenclature of Frost (2023). We only captured mature individuals for our experiment. They were recognized by their advertisement calls in males, the presence of eggs in the belly of females, or determined by body size. To avoid pseudoreplication, within the same species we distinguish among individuals based on their unique coloration patterns (Donnelly et al. 1994; Ferner 2010). So, during the experiment, each individual was tested only once. For species whose individuals show no difference in coloration pattern, individuals were collected and placed in wet bags until additional individuals were found. After capturing an individual, we recorded its body size (snout-vent length, SVL) using a digital caliper and collected additional information about its microhabitat (aquatic, semi-aquatic, fossorial, terrestrial, semi-arboreal, arboreal). The locomotion type (synchronous, asynchronous) was determined based on literature and personal observations in the field.

The methodology described below was carried out in the field at a pre-selected site prior to capturing the individuals. To assess the strength of laterality in forelimbs of anurans (Fig. 2), we followed the methodology proposed by Bisazza et al. (1996). Prior to the experiment, we placed each individual on a horizontal surface (Fig. 2a), and used a small opaque plastic container to enclose it for five minutes (Fig. 2b). Then, we handled one individual at a time, while being handled, the subject’s limbs were restricted from movement to control for possible biases in the use of their limbs for snout-wiping during the experiments (Malashichev and Nikitina 2002). We placed a wet microfiber cloth covering its head and eyes (Fig. 2c). The size of the cloth was adjusted to the head size of each individual. The experiments began by placing the individual on a horizontal surface and recording the first forelimb it used to remove the cloth from its head (Fig. 2d). The percentage of individuals using the right forelimb to remove the cloth served as a proxy for the strength of laterality (Robins et al. 1998; Malashichev and Nikitina 2002). After the experiments, all individuals were released back into the site of capture. For each species, the experiment was conducted in a single locality to control for possible interpopulation variation in the strength of laterality (Frasnelli and Vallortigara 2018; Vallortigara, personal communication).

Fig. 2
figure 2

Methodology used to quantify the percentage of forelimb right-handed anurans (Bisazza et al. 1996). After placing an individual on a horizontal and smooth platform a, he was covered with an opaque plastic container to favor coupling to the experimental arena b. After five minutes, the plastic container was removed, the individual was taken, and a wet microfiber cloth was placed over his head, covering his eyes c. We recorded the first forelimb used by the individual to remove the microfiber cloth d. The red arrow in the last image indicates the movement that would be recorded for the focal individual (in this example: left forelimb). Two videos (Appendix I) shown the use of forelimbs by individuals of Leucostethus brachistriathus (Dendrobatidae) and Atelopus laetissimus (Bufonidae) under snout-wiping testing conditions, as well as Andinobates bombetes (Dendrobatidae) under natural conditions

For expanding our dataset, we performed a systematic review of the literature on the strength of laterality of forelimbs in anurans. The literature search was conducted using Google Scholar, Web of Science, and Scopus platforms using the keywords “lateralization”, “asymmetry”, “locomotion”, “frogs”, and “behavior”. We searched for papers written in English and Spanish languages.

Statistical analysis

To determine whether populations of our target species tend to have a higher or lower percentage of right or left-handed individuals than expected by chance (50%), we performed a two-tailed binomial test (Wagner-Menghin 2005; Whitlock and Schluter 2015). To test for a relationship between the strength of laterality in forelimbs and the microhabitat type, locomotion type, and body size, we transformed the percent of forelimb right-handed individuals per species (strength of laterality) into a “delta value of laterality”. This delta value represents the absolute difference between the recorded percentage of right-handed individuals and the expected percentage of right-handed individuals under the null hypothesis (i.e., 50%). This transformation was necessary because laterality in forelimb can be two-tailed biased (left- or right-handed), and analyzing the statistical data using the percentage of right-handed individuals alone could mask patterns if there are right-handed and left-handed species.

Given that species cannot be considered as independent statistical units (Felsenstein 1985; Revell and Harmon 2022), all of the statistical analysis were corrected for the phylogenetic relationship between species (Jetz and Pyron 2018). These analyses were conducted by pooling delta values of laterality obtained from our experiments and those found in literature. Genetic information was not available for the common toad Rhinella horribilis; therefore, we excluded this species from the following analyses. Due to the small sample size (number of species), which limited our ability to include all potential explanatory variables (body size, microhabitat type, locomotion type, and interactions), we initially tested the relationship between microhabitat type and delta values of laterality using a phylogenetic ANOVA (Revell 2012). In this analysis, the number of levels in the categorical variable “Microhabitat type” was reduced from six (Fig. 1) to three: (i) aquatic, comprising aquatic and semiaquatic species; (ii) terrestrial, referring to both terrestrial and fossorial species, and (iii) arboreal, putting together arboreal and semi-arboreal. The relationship between body size and locomotion type with delta values of laterality in forelimb was tested using a phylogenetic ANCOVA (Fuentes-G et al. 2016). These two comparative phylogenetic analyses were performed using the functions ape (Paradis et al. 2004), geiger (Harmon et al. 2008), phangorn (Schliep 2011), nlme (Pinheiro et al. 2023), and phytools (Revell 2012) with the packages “geiger”, “ape” and “phytools” (Harmon 2018; Revell and Harmon 2022) from R software (v. 3.6.1, R-Development-Core-Team, 2019).

Results

We conducted our laterality experiments on 338 individuals belonging to 19 anuran species: one semi-aquatic, nine terrestrial, one semi-arboreal, and eight arboreal (Table 1). The results of the binomial test indicate that populations of four species (Fig. 3) exhibit a higher proportion of left- or right-handed individuals than expected by chance (Table 1): most individuals of Andinobates bombetes, and Dendropsophus microcephalus used their right hands, while in Leucostethus brachistriatus and Atelopus laetissimus most individuals used their left hands. The other 16 tested species did not show a laterality tendency that significantly differs from that expected by chance (Table 1).

Table 1 Summary of information about the strength of laterality (percentage of forelimb right-handed individuals), ecological and phenotypic characteristics for 33 anuran species. Body size (snout-vent length, SVL) is expressed in mean and range values. A methodology similar to this study was used for all species (see methodology text). The results of the Binomial test (i.e., p-value) are shown for the 20-target species of this study. The absence of values in some cells is due to the unavailability of detailed information in the papers used as source of information. Delta values of laterality refer to the difference (absolute value) between the recorded percentage of forelimbs right-handed individuals and the expected percentage of right-handed individuals under a null hypothesis (i.e., 50%).
Fig. 3
figure 3

Images of four anuran species that exhibit a statistically significant difference or tendency (at alpha = 0.05) to use the right or left forelimb more frequently in snout-wiping experiments than expected by chance: a Andinobates bombetes, b Leucostethus brachistriatus, c Atelopus laetissimus, and d Dendropsophus mirocephalus

In the systematic review, we found 26 papers (published between 1992 and 2023) in which authors employed a similar methodology to study anurans as in our study (Table 1). In these papers, the strength of laterality was recorded for 14 anuran species, and included information on the type of microhabitat in which they were reported (one aquatic, four semi-aquatic, two fossorial, five terrestrial, and two arboreal, Table 1). The phylogenetic ANOVA revealed that delta values of laterality did not differ significantly between aquatic, terrestrial, and arboreal species (phylANOVA F = 3.09, p = 0.37, N = 32 species). Results of the phylogenetic ANCOVA (Table 2) indicates that species with asynchronous limb movements exhibit higher strength of laterality compared to species with synchronous limb movements, while body size is not related to delta values of laterality. The interaction between body size and locomotion types did not have a significant effect on delta values of laterality (Table 2).

Table 2 Results of phylogenetic ANCOVA indicating the relationship between the strength of laterality (delta value of laterality), body size (snout-vent length, SVL), and locomotion type (synchronous, asynchronous) in 33 anuran species. The symbol * indicates an interaction between explanatory variables. This analysis includes the phylogenetic relationship among species following Jetz and Pyron (2018). The toad Rhinella horribilis was not include in this analysis because unavailability of genetic sequences

Discussion

We found that four out of the 19 anuran species tested in this study exhibited a higher proportion of left- or right-handed individuals than expected by chance. These four species were Andinobates bombetes, Leucostethus brachistriatus (Dendrobatidae), Atelopus laetissimus (Bufonidae), and Dendropsophus microcephalus (Hylidae). When we pooled our results with those obtained in similar studies (total = 32 species excluding R. horribilis), we found that anurans with asynchronous movements showed a higher strength of laterality in the use of their forelimbs compared to anurans with synchronous movements. The strength of laterality in forelimbs was not related to microhabitat or body size of the species.

In quadrupeds, asynchronous limb movements are associated with greater functional asymmetry in the neural systems compared to organisms with synchronous limb movements (Mikliaeva et al. 1987; Forward et al. 1962; Malashichev 2006). In turn, this would mean that a species with a tendency to use limbs unilaterally will exhibit a higher strength of laterality (Mikliaeva et al. 1987; Malashichev 2006). The results of the phylogenetic ANCOVA performed with the pooled data of 32 anuran species, support this statement. Moreover, the four species in which we detected a significant strength of laterality (D. microcephalus, A. bombetes, L. brachistriatus, A. laetissimus) exhibit asynchronous movements in the use of their forelimbs.

In the tree frog D. microcephalus (right-handed), the individuals inhabit areas with temporary or semi-permanent bodies of water. While individuals jump or climb grasses and thin branches of shrubs, it is possible to observe the asynchronous movement of their limbs. In the case of the two dendrobatid species, A. bombetes (right-handed) and L. brachistriatus (left-handed), they inhabit leaf litter mounds and cavities under trunks or rocks on the forest floor (Lötters et al. 2007; Vargas-Salinas et al. 2014). Due to the numerous physical obstacles in and under such substrates (e.g., small stones, fallen branches, irregularities of the ground, leaf litter) these species predominantly engage in small walks (personal observations). Additionally, as is typical in dendrobatids, these species exhibit active foraging, conducting short walks in search of food (Bourne 2001; Staudt et al. 2010; Duarte-Marín et al. 2020) and engage in physical combat where asynchronous limb movements are evident (Pröhl 2005; Lötters et al. 2007; Wells 2007; Duarte-Marín et al. 2020).

Like dendrobatids, harlequin toads (genus Atelopus, Bufonidae) walk on the ground or climb understory vegetation (Lötters 1996; personal observations). Males of Atelopus also engage in visual displays (semaphoring) to attract mates (Lindquist and Hetherington 1998; Hödl and Amézquita 2001; Rößler et al. 2019). Both walking and visual displays involve asynchronous movements of the limbs. Therefore, in line with the expected relationship between asynchronous movements and the strength of laterality in anurans (Malashichev and Robins 2018; Fig. 1), 18 out of 24 individuals (75%) of A. laetissimus used their left forelimb to remove the microfiber in our experiments. The 80% of A. spurrelli individuals also used their left forelimb for snout-wiping, although this tendency was not statistically significant in the Binomial test, possibly because the small sample size (N = 10 individuals). The natural history characteristics of species in the family Dendrobatidae and the genus Atelopus (Bufonidae) suggest that they frequently exhibit asynchronous movements of their forelimbs, and hence, it is expected that they exhibit a clear tendency to be left- or right-handed. Future studies are required to corroborate this hypothesis.

According to Malashichev and Robins (2018), anurans in aquatic microhabitats exhibit a lower interspecific variability in the strength of forelimb laterality than anurans in terrestrial and arboreal microhabitats (Fig. 1). However, the results of the phylogenetic ANOVA indicate that there is no relationship between the microhabitat type and the strength of forelimb laterality in anurans. This discrepancy with the predictions of the model by Malashichev and Robins (2018) may occur because these authors only had information for 12 species and linked a specific locomotion type to a specific microhabitat type. For instance, they mention that terrestrial species tend to be hoppers while arboreal species are climbers. Since we included a greater number of species to test such a model (i.e., 32 species), a mismatch between locomotion type and microhabitat type arises. That is, within one microhabitat type, there are species with synchronous and species with asynchronous locomotion types (Enríquez-Urzelai et al. 2015).

Body size was unrelated to the strength of laterality, which was unexpected because in vertebrates, body size is known to be related to many physiological, ecological, and life history traits of species (Schmidt-Nielsen 1984; Wells 2007). Malashichev and Robins (2018) mention that Bufotes viridis (a terrestrial anuran), does not exhibit the same strength of laterality in forelimbs as other bufonids because individuals of this species are small compared to other representatives of the family Bufonidae. Hence, B. viridis tend to jump rather than walk. With this statement, those authors suggest that body size could be related to locomotion type in anurans. However, Enríquez-Urzelai et al. (2015) and Citadini et al. (2018) reported for 17 Mediterranean and 64 tropical anuran species, respectively, that there was no relationship between body size and locomotion type. The absence of this relationship can be attributed to the allometry between body size and limb size (Collar et al. 2011; Enríquez-Urzelai et al. 2015). Arboreal and aquatic species have relatively longer fore and hindlimbs than terrestrial and fossorial species regardless of body size differences (Enríquez-Urzelai et al. 2015; Fabrezi et al. 2017; Citadini et al. 2018; Moen 2019). Given that we find a relationship between the strength of laterality in forelimb with locomotion type, but not with body size, our results support the findings by Enríquez-Urzelai et al. (2015) and Citadini et al. (2018). That is, the strength of laterality in anurans may be more related to relative size of limbs or a set of morphological traits that define body architecture of species (i.e., ecomorphs) rather than just body size (Losos 2009; Moen et al. 2016, Moen and Wiens 2017).

We are aware of potential limitations in our study. First, laterality was not assessed at the individual level (e.g., repeated use of the same forelimb for snout wiping by an individual), but rather at the population level. It is possible that each individual has a strong tendency to be right-or left-handed, but if a similar number of individuals in the sample are right- or left-handed, statistical tendences may not be detected in studies like ours. However, it has been found that laterality recorded at the individual level is correlated with laterality at the population level (Frasnelli and Vallortigara 2018). Second, we tested the potential relationship between the strength of laterality in forelimbs and body size (measured as SVL), which has not been quantitatively tested before in a set of anuran species. However, species with similar body size can exhibit divergent body architectures (ecomorphs), which in turn determine performance in various tasks (Losos 2009; Muschick et al. 2012; Moen et al. 2013). Considering that strong selection processes have driven morphological convergences in anurans that occupy similar microhabitats (Moen et al. 2016, 2017), our comparison of forelimb laterality between aquatic, terrestrial and arboreal species serves as a proxy for comparing anuran ecomorphs. Nonetheless, within the same ecomorph there can still be variability in traits; further studies on laterality that directly examine sets of morphological traits of ecomorphs and include a large number of species may reveal new relationships (Moen 2019).

This study documents some of the highest strengths of laterality recorded to date in forelimb use for snout-wiping in anurans (Robins et al. 1998, Malashichev 2002, 2006), and to the best of our knowledge, is the first one for neotropical anuran species (but see Anderson et al. 2021). Our results mostly align with the predictions of the model hypothesized by Malashichev and Robins (2018) and demonstrate that despite advances in our understanding of the interplay between ecology and laterality in animals (Rogers and Andrew 2002; Csermely and Regolin 2013; Rogers et al. 2013), there are still many gaps that need to be addressed. For instance, the right hemisphere and the left hemisphere of the brain in vertebrates control different types of activities and abilities (Robins 2006; Csermely and Regolin 2013). Therefore, the fact that some species, even within the same family (e.g., Dendrobatidae), tend to be right-handed while others tend to be left-handed is intriguing. Testing possible causes of these contrasting results is beyond the scope of the current study, but an interaction between aposematism, predation risk, and tendencies in forelimb laterality cannot be ruled out (Anderson et al. 2021). Exploring the relationship between the strength of laterality and the asynchronous use of limbs for capturing and manipulating prey (Gray et al. 1997; Csermely and Regolin 2013) is also an aspect that has been overlooked in anurans and requires further studies (but see Bizzasa et al. 1996, 1997). Overall, the strength of laterality in animals can be influenced by numerous factors (Rogers et al. 2013; Malashichev and Robins 2018), and with this study, we hope to encourage future research on laterality and its relationship with the evolutionary ecology of neotropical anurans.