Introduction

Arboreal environments impose strong physical constraints on vertebrate morphology and locomotor adaptations, due to the diversity of supports varying in orientation, diameter, and compliance. For arboreal mammals, the ability to navigate sloping and vertical trunks, branches, and lianas is essential, as these constitute the most available supports in trees. Climbing on vertical supports is a key aspect of arboreal locomotion in both fully arboreal and scansorial species, enabling access to the canopy from the ground, diverse food sources, enhanced vigilance, predator avoidance, and secure resting and nesting sites 1. Climbing abilities in mammals date back to the Jurassic 24 and have played a central role in the evolution of euarchontoglires, especially primates 513. Adaptations for climbing vertical supports have been identified in early Paleogene primates through postcranial fossil analyses 1417. Frequent vertical climbing may have been a prerequisite for the development of hallucal and pollical grasping in early primates 18,19, and has been proposed as an alternative scenario to the hypothesis of cautious displacements on narrow branches of angiosperm tree peripheries 2022. Several extant mammal species have been proposed as models for early primate arboreal adaptations, including small strepsirrhines 23,24, platyrrhines 13,19,25, treeshrews 26, rodents 2729, carnivorans 30, and marsupials 19,3135. While not all arboreal mammals travel on narrow terminal branches, all rely on vertical supports to navigate canopy strata. However, few comparative studies have examined vertical climbing in non-primate arboreal mammals, and most have focused on ascents 11,3638. Although behaviors such as dropping, jumping, or gliding allow arboreal animals to move to lower tree layers, the ability to safely descend sloping and vertical supports remains critical, yet largely understudied 39.

Non-primate arboreal mammals often adopt a head-first posture when descending vertical supports, as reported in some rodents 12,40,41 and procyonid carnivorans 30,42. This strategy relies on reversed feet and functional claws that enable secure gripping of tree bark during descent. Head-first descent is also common in prehensile-tailed mammals, such as certain platyrrhine primates 4345 and the kinkajou 30. In species lacking a prehensile tail, and in primates that lack claws and use their digits to wrap around supports, alternative descent strategies are required. In a study of nine small to medium-sized strepsirrhines, Perchalski et al. 46 found that lemurs under 1 kg predominantly used head-first descents regardless of support properties, whereas larger species (>1 kg) increasingly adopted tail-first postures on supports angled at 45° or more, depending on species. Both head-first and tail-first strategies involve orthograde postures with the body aligned parallel to the support 7. Perchalski and colleagues also recorded other descent behaviors, including leaping, dropping, and asymmetrical descents that did not fit into strict head-first or tail-first categories. These findings suggest that in strepsirrhines, body mass and locomotor adaptations influence the choice of descent posture. However, such behavioral variation also highlights the need for broader comparative data across arboreal mammals to fully understand the functional and ecological significance of vertical descent strategies in an evolutionary context.

From a locomotor perspective, speed and duty factor (DF, the percentage of stride duration a limb remains in contact with the support) are key indicators of arboreal locomotor performance 47. Studies across tetrapods have shown that animals generally move more slowly on narrow supports than on wider ones 48,49. Small primates 24,50,51, rodents 12,52,53, and marsupials 34,54 typically decrease speed and increase duty factor, especially for the forelimbs, during descents compared to ascents on inclined supports and compared to locomotion on horizontal supports. These mechanisms supposedly help enhance stability by allowing more cautious displacements 39,48,53. Gait patterns also reflect arboreal adaptations, particularly in primates 36,47,55,56. Vertical climbing employs both symmetrical (walk, trot, run) and asymmetrical (bound, half-bound, gallop) gaits, similar to above-branch locomotion. Asymmetrical gaits, often used by small mammals at high speeds 48,57, allow to increase the distance traveled, but often include an aerial phase 58, raising support reaction forces and potentially causing oscillations or breakage on narrow horizontal or sloping branches, complicating locomotor control and safety. Benefits of asymmetrical gaits have thus been proposed to be related to both body mass and support diameter 59. On vertical supports, such as trunks or lianas, these risks are reduced due to different force distributions and gravity orientation. Symmetrical gaits are widespread in primates and have been linked to the use of the fine branch niche, as enhanced grasping and hindlimb dominance supposedly increase stability during cautious locomotion on narrow horizontal or sloping supports 56,6062. Notably, primates adjust their symmetrical gaits with support inclination, favoring diagonal-sequence gaits during ascents and lateral-sequence gaits during descents 51,6367. Symmetrical gaits have also evolved convergently in rodents 12,52, scandentians 68, carnivorans 30, and marsupials 6972. However, most kinematic studies have focused on horizontal or moderately sloped supports (∼45°), rarely addressing vertical (∼90°) supports or analyzing both symmetrical and asymmetrical gaits comprehensively. As a result, comparative data on how vertical supports of varying diameters affect descent and ascent kinematics remain limited, hindering our understanding of this common arboreal behavior.

In addition to body mass, limb proportions influence locomotor and postural strategies 46,63. The relative length of autopodials, particularly manual and pedal digits, is closely linked to grasping ability and thus affects locomotor performance in primates 13,19,21,25,73 and small marsupials 19,32. Species with relatively short autopods may struggle to maintain a secure grip during vertical descent. Fore- and hindlimb length ratios have also been proposed as predictors of postural adaptation in vertical locomotion, with shorter forelimbs generating greater hindlimb traction during head-first descents 36,46. Consequently, a lower intermembral index (forelimb-to-hindlimb length ratio) increases the risk of forward pitching when descending vertically head-first compared to tail-first 46.

Variation of weight distribution along the body, notably caused by head mass differences, may also influence postural strategies during vertical descent (i.e., head-first vs. tail-first). Primates exhibit the highest brain-to-body mass ratio among mammals 74. This relative brain expansion, traced back to early primate evolution, is a key synapomorphy alongside a divergent hallux and pollex, nails replacing claws, elongated hindlimbs and autopods, and large convergent eyes 21,75,76. However, the evolutionary sequence and interplay between locomotor and cognitive specializations remain unresolved, despite consensus on the small body size 13,77 and arboreality 7880 of primate ancestors. While brain enlargement likely increases relative head mass, its effects on arboreal locomotion have yet to be investigated in primates or other mammals. Given that head mass influences both the center of mass and moment of inertia, greater head mass may increase toppling risk during vertical head-first descent.

In this study, we aim to test (a) whether postural and kinematic adjustments during vertical descents versus ascents vary consistently across mammalian groups; (b) whether these strategies correlate with specific morphological parameters; and (c) whether such correlations can be used to infer descent postures in fossils, shedding light on the evolutionary history of vertical arboreal locomotion, particularly in primates. To address these questions, we analyzed locomotor sequences of descents and ascents on vertical supports in 21 small- to medium-sized arboreal mammal species, including taxa proposed as modern analogues of early primates. We first examined the effect of support diameter (small, medium, large) on postural strategies and kinematics (speed, duty factor, gait) during descents relative to ascents. We then assessed the relationships between descent strategies and morphological parameters (body mass, limb proportions, and relative head mass) for each species. These data allowed us to build a model incorporating seven significant morphological predictors, which we applied to infer vertical descent behavior in 13 extinct euarchontoglires, including adapiforms (early strepsirrhines), omomyiforms (early haplorhines), a sciuromorph (early rodent), and plesiadapiforms (stem primates).

We predict that:

H1) Postural strategies: species under 1 kg and/or possessing a prehensile tail will exhibit head-first descent postures more frequently than larger species, which will favor tail-first postures 46.

H2) Kinematic adjustments: animals will adopt more cautious locomotor patterns during vertical descents compared to ascents, characterized by reduced speed, increased duty factors, and a higher proportion of symmetrical gaits, particularly on narrow supports 48,52,53,63,68.

H3) Morphological correlates: the intermembral index will positively correlate with the proportion of head-first descents, whereas relative autopod length and head size will correlate negatively with head-first descent frequency 21,36,46.

H4) Evolutionary implications: early euarchontoglires, due to their small size and generalized locomotor abilities, likely employed a high proportion of head-first descents 12,53, whereas early euprimates may have increasingly favored tail-first descents associated with larger body and brain, and elongated hindlimbs and autopods 21,75,76.

Results

Descent strategies on vertical supports among arboreal mammals

We recorded the locomotion of 57 adult individuals from 21 species: 6 strepsirrhine primates, 5 platyrrhine primates, 1 scandentian, 3 rodents, 3 carnivorans, and 3 marsupials (Fig. 1, Table S1). We categorized supports in three diameters: small (i.e., twigs), medium (i.e., branches), and large (i.e., boughs and very large trunks) relative to the size of the animals’ hands and feet (see section Materials and Methods for more details). Whenever possible, we analyzed up to 10 ascents and 10 descents per individual on each support type. In total, we analyzed 1,400 ascents and 1,390 descents.

Strategies of descent on vertical supports of various diameters by species.

A. Photographs illustrating the three different strategies of descent identified in this study. Black arrows represent the axis of the body regarding the support and the direction of the movement, and grey arrows represent the direction of the gaze. B. Proportions (in %) of each descent type (yellow = head-first descent, light green = side descent, and dark green = tail-first descent) occurring on each support category (vertical branches of large, medium, and small diameters represented by brown icons at the top of the graph), by species. Proportions of all individuals of the same species (n) were averaged. See Table S1 for details on the individuals studied, and Table S2 for the mean body masses (BM) references. NA = no data for the given condition. Results of statistical tests between support diameters (ANOVAs or MANOVAs on bootstrapped samples, confirmed by post hoc tests) are represented with annotations defined as NS (non-significant): p > 0.05, *: p < 0.05, **: p < 0.01, and ***: p < 0.001. See SuppTables A for associated p-values. C. Associated phylogeny of the species studied with branch length obtained from http://timetree.org 129. D. Photographs of individuals of various species descending vertical supports of different diameters.

We identified three distinct strategies of descent: head-first descent, side descent, and tail-first descent (Fig. 1A) which partly overlap with Perchalski’s ethogram 46. During a head-first descent the body is positioned head-down, parallel to the support, hands and feet are anchored to the support, and the face is directed downward (see Trichosurus vulpecula in Fig. 1A). During a side descent the body is rotated perpendicularly to the support and the face is also directed downward (see Plecturocebus cupreus in Fig. 1A). During a tail-first descent the body is positioned head-up and parallel to the support similarly to ascents, and the face is either directed horizontally or downward (see Eulemur mongoz in Fig. 1A). We found a clear difference between primates and non-primates in the occurrence of each descent strategy on each support type (Fig. 1B). Most non-primate species, although very diverse in body mass and morphology, showed only head-first descents independently of the orientation and diameter of the supports. The racoons, the largest species of this study (∼ 6 kg), were an exception and mainly exhibited tail-first descents. The other studied carnivorans, i.e., the coatis and the kinkajous (the latter being the only prehensile-tailed of this study along with the brushtail possum Trichosurus vulpecula), as well as all marsupials, rodents, and treeshrews exhibited only head-first descents independently of the support diameter (Fig. 1B, D). In sharp contrast, primates rarely used head-first descents on vertical supports except for lorises (Xanthonycticebus pygmaeus, strepsirrhine) and tamarins (Saguinus imperator and S. Oedipus, platyrrhine), the smallest primate species of this study (∼ 360 g to ∼ 470 g respectively). The other strepsirrhines (i.e., lemurs) descended vertical supports mainly using tail-first descents with higher proportions on small supports, and the other platyrrhines generally used side descents with also increased proportions on small supports. The difference between phylogenetic groups is substantiated by a significant phylogenetic signal contained in the overall proportion of head-first descent strategy on vertical supports (k=0.58, p = 6.6.10−4, effect size= 2.57).

Kinematic adaptations during ascents and descents on vertical supports

For each complete ascent and descent sequence of locomotion, characterized by successive lift-off, touchdown, and following lift-off of all four limbs, we quantified: the absolute speed, the relative speed (in relation to the individual’ body length), the forelimbs and hindlimbs’ duty factors (Fig. 2), and the gait type (categorized as symmetrical or asymmetrical, Fig. 3).

Kinematics of locomotion during ascents and descents on vertical supports of varying diameters.

A. Mean absolute speeds by species during ascents and each descent strategy on the vertical supports of large, medium, and small diameters (represented by brown icons on the top left of the plots), with box plot representations of the medians for all species combined in ascents and each descent strategy. A. Mean relative speeds (based on individuals’ body length) by species during ascents and each descent strategy on each support type, with box plot representations of the medians for all species combined in ascents and each descent strategy. C. Mean duty factors by species, differentiating forelimbs (left box plots) and hindlimbs (right box plots), during ascents and each descent strategy on each support type, with box plot representations of the medians for all species combined in ascents and each descent strategy. Results of statistical tests between ascents and each descent strategy (two-sided Mann Whitney U tests) are represented with annotations defined as NS: p > 0.05, *: p < 0.05, **: p < 0.01, and ***: p < 0.001. See SuppTables B, C and D for associated p-values. Percentages correspond to the relative variation of each kinematic variable mean in each descent condition compared to ascents.

Proportions (in %) of gait types by support diameters for ascents and each descent strategies, for all phylogenetic groups.

Results of statistical tests on the fraction of symmetrical versus asymmetrical gaits by group between ascents and descents on each support (two-sided Wilcoxon signed rank tests) are represented with annotations defined as NS: p > 0.05, *: p < 0.05, **: p < 0.01, and ***: p < 0.001. NA = not applicable. DSDC: Diagonal sequence - diagonal couplet, LSDC: Lateral sequence – diagonal couplet, LSLC: Lateral sequence – lateral couplet. See SuppTables F for associated post hoc corrected p-values. See Fig. S1 for detailed proportions of gait types by species and support diameters.

We found a significant general tendency for all species to reduce their absolute speed by approximately 30 to 43 % when descending head-first on vertical supports compared to ascents, regardless of support diameter (Fig. 2A, SuppTables B). However, when descending sideways or tail-first, primates did not significantly reduce their speed compared to ascents, except for side descents on medium supports (−16,36 %). On the contrary, in primates the relative speeds decreased during side and tail-first descent on all supports by between 17 and 46 %, indicating that animals tended to decrease the distance travelled compared to ascents (Fig 2B, SuppTables C). But during head-first descents the relative speeds tended to increase by ∼ 10 % on medium supports, while it decreased on small and large supports. It is interesting to note that the smallest species of our dataset, i.e. rodents and treeshrews (< 200 g), exhibited the highest relative speeds compared to other taxa, thus travelling larger distances in proportion to their size (Fig. 2B). On the other hand, all studied animals tended to increase both forelimbs and hindlimbs duty factors during descents compared to ascents, particularly during head-first descents on all supports by ∼ 4%, and during side descents on medium and small vertical supports, but not when descending tail-first (Fig. 2C, SuppTables D).

We found that gait contains a low, yet significant phylogenetic signal (Kmult = 0.238; p = 4.31.10−2; Effect Size: 1.72), suggesting the impact of shared evolutionary history of related species on locomotor characteristics. We additionally found that there is a significant similarity between phylogenetically related species and dissimilarity across phylogenetically distant groups, allowing us to merge the gait data at the phylogenetic group level (as defined in Fig. 1C) to obtain a more general picture of the large dataset (MANOVA: p = 2.09.10−2, Wilks’ λ = 3.36.10−3, df = 5, SuppTables E). Considering the difference in the proportions of descent strategies obtained between strepsirrhines and platyrrhines, we decided to keep these two subgroups separate rather than grouping all primates together. When analyzing the effect of vertical descent strategies on gait types by phylogenetic groups, we found that animals exhibited a mosaic of gait types depending on the species, support diameter and direction of movement. Interestingly, the proportion of asymmetrical gaits significantly increased during descents on all supports for all groups studied (Fig. 3, SuppTables F). In effect, during vertical descents, the animals exhibited more bounds, gallop, and related asymmetrical gaits with relatively high duty factors such as ambling 81, not only during side and tail-first descents, but also during head-first descents (Fig. S1). Primates and marsupials overall exhibited a higher proportion of symmetrical gaits compared to the other studied mammals and particularly diagonal and lateral sequence - diagonal couplets (DSDC and LSDC respectively), especially during vertical ascents (Fig. 3, see Fig. S1 for detailed gait types by species and support diameters). We also documented symmetrical gaits in rodents and carnivorans. Primates, marsupials, and rodents exhibited more symmetrical gaits during ascents on vertical supports of small and medium diameter compared to large vertical supports. Carnivorans, especially raccoons, exhibited a relatively small proportion of symmetrical gaits on large and medium supports during ascents, and treeshrews used only asymmetrical gaits during both ascents and descents on vertical supports of medium and large diameter. Furthermore, strepsirrhines exhibited some lateral sequence – lateral couplets (LSLC) only during vertical descents, and rodents and carnivorans exhibited more DSDC during descents and more LSDC during ascents.

Correlations between vertical descent strategies and morphological features

To investigate links between vertical descent strategies and morphological parameters, we collected ten morphological variables including body mass (Table S1), eight intrinsic body proportions calculated from measurements taken on each studied animal (Table 1), and relative head mass of each species, approximated by its encephalization quotient (EQ) (Table S2). We used the EQ 82 as it quantifies whether the brain mass of a given species is lower (<1) or higher (>1) than expected from a general brain-to-body allometric relationship. We calculated the mean EQ for each species from their mean endocranial volumes, collected either from the literature or from our own measurements for previously documented species (Tables S2, S3). To ensure that the EQs of the different species studied are comparable meaningful, we tested the correlation between the brain and body masses in our dataset and found a significant and consistent positive allometry for all species (Pearson correlation: R = 0.96, p = 5.0.10−12), and similar allometric coefficients when restricting the analysis to phylogenetic groups (Fig. S2).

Definitions of body measurements performed on extant and extinct individuals and calculation of limbs proportions.

Among the ten morphological variables analyzed, we found that the arm and leg relative lengths correlate with each other and with the mean EQ, and the leg and tail relative lengths correlate with each other (SuppTables G). Moreover, we found that body mass, EQ, relative arm and leg lengths, and relative tail length significantly and positively correlate with the proportion of head-first descent strategy on vertical supports, while the intermembral index significantly and negatively correlates with the proportion of head-first descent strategy (all diameters combined) (Fig. 4 A, B, C, Fig. S3 A, B, C, SuppTables H). Also, in primates, only the relative length of the hallux correlates significantly and negatively with the proportion of side descents (Fig. 4 D, SuppTables I). Primates stand out from the other studied groups in having relatively longer legs (> 70% of body length, Fig. 4C), arms (> 70% of body length, Fig S3A), tails (> 90% of body length, except for the lorises who have a very reduced tail, Fig. S3C), and high EQs (> 1.5, Fig. 4B). In comparison, carnivorans, albeit being the largest studied species (> 3kg, Fig. 4A), possess relatively shorter limbs (∼ 60% of body length, Fig. 4C, Fig. S3A), and medium EQs (> 1 and < 1.2, Fig. 4B). The studied rodents, scandentians, and marsupials also all possess relatively shorter limbs (< 70% of body length, Fig. 4C, Fig. S3A), and tails (< 90% of body length, Fig. S3C), and show small to medium body masses (< 3 kg), and low EQs (≤ 1).

Morphological correlates with vertical descents and prediction of descent strategy in extinct species.

Spearman correlations and associated corrected p-values between the proportion of head descents by extant species and A. the logarithm of mean body mass, B. the mean encephalization quotient (EQ), C. the leg to body length ratio (in %), and D. between the proportion of side descent in extant primate species only and the hallux to foot length ratio (in %). Each plot includes also the head descent proportion (A., B., C.) and side descent proportion (D.) of fossils predicted with the Schafer’s multiple imputation procedure. See Table 1 for body measurements definitions and Tables S2 and S3 for mean body masses, ECV and calculated EQ by species. See Fig. S3 for the plots of the percentage of head descent against the six other morphological variables analyzed by species, and SuppTables H and I for Spearman rhos and associated corrected p-values between each morphological variable and the proportion of respectively head (SuppTables H) and side descents (SuppTables I) on vertical supports. See Table S4 for calculated morphological variables of extinct species and SuppTable J for the resulting predictions of their descent behavior.

Integrating extant behavioral and morphological data to hypothesize vertical descent strategy in early euarchontoglires

Based on these results, we attempted to infer the likely descent strategies of extinct species, using their known body proportions. We selected 13 fossils relevant to the understanding of early primate evolution including 7 euprimates, 5 plesiadapiforms and 1 rodent with relatively complete postcranial and/or cranial remains (Table S2, S4). We used a Multiple Imputation (MI) method to predict descent behaviors given partial fossil morphological data. Following the benchmark of Clavel and colleagues 83 which compared several MI methods for the precise problem of morphometric fossil estimation, we used Schafer’s MI procedure from the R package norm, with 5000 runs of 50 steps Markov-Chains. We calculated the prediction of head-first descent proportion on vertical supports (all diameters combined) for all fossils, using all extant and extinct species, and retaining all six morphological variables that significantly correlate with head-first descent (Fig. 4A, B, C, Fig. S3A, B, C). We also computed the prediction of side descents only for primate and plesiadapiform fossils, retaining all extant primate species and the respective correlated variable, i.e., the relative length of the hallux (Fig. 4D).

Fossil species were mostly reconstructed as employing a relatively high proportion of vertical head-first descent strategy (> 75%), except for Darwinius masillae and Europolemur kelleri, two medium-sized (660g and ∼1.5kg respectively) adapiformes (early strepsirrhines primates), that were estimated at around 30% and 60% of vertical head-first descent strategy respectively (Fig. 4A, C, Fig. S3A, B, C, SuppTables J). Darwinius remains are very well preserved and we were able to calculate all the morphological variables for this species. For Europolemur, however, we could only include the estimated body mass and foot proportions in the analysis. The other adapiforms of this study, Adapis parisiensis and Notharctus tenebrosus, were reconstructed with high head-first descent proportions (∼ 83% and ∼ 86% respectively) despite being larger (∼ 1.5kg and ∼ 9kg respectively). However, due to the lack of their postcranial remains described in the literature, we could only include their reconstructed body mass and EQ in this study, the latter being relatively low (0.56 and 0.26 respectively) compared with extant strepsirrhines. Interestingly, the omomyiforms (early haplorrhines) studied here, Archicebus achilles, Microchoerus erinaceus, and Tetonius homunculus were reconstructed with very high head-first vertical descent proportions (> 90%, Fig 4A, B, SuppTables J), but they are also known mainly from their cranial remains and only body mass and EQ were used in the head-first descent model. Yet, their EQs are also very low (< 1) compared with those of extant haplorrhines (Fig 4B, Table S2, SuppTables J). Concerning the plesiadapiformes and the fossil rodent studied, which display variable body masses ranging from the small Carpolestes simpsoni (∼ 100 g) to the medium-sized plesiadapids (∼ 2 kg, Fig. 4A, Table S2), they were all reconstructed as exhibiting a relatively high proportion of head-first descent strategy (> 75%, SuppTables J). Those species possess relatively elongated limbs whose proportions fall between extant non-primates and primates (Fig. 4C, Fig S3A, B, Table S4), but retain a low EQ (Fig. 4B, Table S2). Concerning the reconstruction of side descent strategy based on hallucal proportions, our model reconstructed the omomyiform Archicebus achilles with a low side descent proportion (23%, Fig 4D, Table S4, SuppTables J), despite its small size (∼ 25g) but relatively very long hallux (50% of foot length), similar to extant strepsirrhines. The same holds for the adapiforms Darwinius masillae and Europolemur kelleri that were estimated with moderate side descent proportions (33% and 55% respectively, Fig. 4D, SuppTable J). Finally, the plesiadapiforms Carpolestes simpsoni and Plesiadapis insignis, with hallucal proportions similar to the studied platyrrhines, were reconstructed with relatively high side descent proportions (respectively 58% and 75%, Fig. 4D, SuppTable J).

Finally, we investigated if the proportion of head-first descent on vertical supports follows a particular evolutionary pattern. Based on the literature, we included extinct taxa in the phylogeny of extant species (Fig. 5) and compared several quantitative models of evolution using the Akaike Information Criterion (AIC). We considered simple models commonly included in studies of morphological evolutionary dynamics: a single-rate Brownian Motion model (BM), the Early Burst model (EB), and a single-peak Ornstein-Uhlenbeck model (OU). We also considered models with two optima: a multiple-rate BM model (BMM) and a multiple-peak OU model (OUM). We considered various relevant phylogenetic nodes where peak shifts or rate changes could have occurred. We found that OUM models with a peak-shift for ‘euprimates + plesiadapiforms’ or ‘euprimates only’ yielded the best fits (AIC= 9.46 and 9.48 respectively) both clearly favored from the other evolutionary models tested (SuppTables K). These models have a head-first descent proportion theta peak of ∼44% and ∼45% respectively, for a theta peak of ∼95% and ∼93% for other mammals groups respectively. After rescaling the phylogenetic tree height to 1, the strength of return to the optimum is quite high (alpha = 7.49 and 7.51 for the ‘euprimates + plesiadapiforms’ and ‘euprimates only’ models, respectively) and the phylogenetic half-life is small (0.093 and 0.092, respectively) confirming that these OUM models fit our dataset best 84.

Phylogeny of studied species, including extant and extinct taxa, and representation of the best evolutionary model of head-first descent proportion.

We calculated a consensus tree for extant species based on 1000 trees from vertlife.org, to which we added extinct taxa using age and phylogenetic topology from literature. See Table S2. for estimated age of extinct species and associated references. The best evolutionary model of head-first descent proportion is an Ornstein-Uhlenbeck evolutionary model with two optimums (OUM) with a peak shift occurring at the node ‘plesiadapiformes + euprimates’. The associated theta peaks of head-first descent proportions of the two subgroups ‘plesidapiforms + euprimates’ vs ‘others’ are represented in red and black respectively.

Discussion

The variability of vertical descent strategies in primates compared to other mammals

To date, this study is the first comparative quantification of ascent and descent strategies on vertical supports across a broad sample of small- to medium-sized arboreal mammals. Overall, our results reveal that vertical descent strategies are strongly influenced by phylogenetic relationships and not solely by body mass, contrary to our initial hypothesis (H1). Specifically, not all species >1 kg systematically descended tail-first. Primates exhibited greater variability in descent strategies—including head-first, side, and tail-first postures—whereas scandentians, rodents, carnivorans, and marsupials primarily employed head-first descents, with the exception of raccoons. Both head-first and tail-first postures involve the body remaining parallel to the support, as previously documented in rodents 12,40,41, procyonids 30,42, prehensile-tailed platyrrhines 4345, and non-prehensile tailed strepsirrhines and marsupials 46,71. However, we also identified a distinct “side descent” posture, observed exclusively in primates, in which the body is rotated perpendicular to the support while the face remains directed downward. This posture likely partly corresponds to the “other descent strategies” described by Perchalski et al. (2021). We quantified side descents in both strepsirrhines and platyrrhines, with non-callitrichid platyrrhines using it most frequently (>70%) across all support diameters. In contrast, the smallest primates (i.e., pygmy slow loris N. pygmaeus, ∼360 g; and tamarins S. imperator ∼475 g, S. oedipus ∼410 g) primarily employed head-first descents, while larger strepsirrhines, particularly lemurs, used tail-first descents in >80% of trials across all support diameters. Primates generally increased the use of tail-first and side descents on small vertical supports, suggesting that narrow supports impose additional constraints favoring more upright postures. This marked behavioral divergence among phylogenetic groups underscores the role of evolutionary history and supports recent findings emphasizing the influence of shared ecological niches and postural adaptations in shaping primate locomotor performance 63,67,8587.

The kinematic adjustments of arboreal mammals during vertical descents

In terms of locomotion, we found that all species reduced their speed and increased both forelimb and hindlimb duty factors during head-first descents compared to ascents on vertical supports, partially supporting our second hypothesis (H2). This pattern, previously observed in small primates 24,50,51, rodents 12,52,53, and marsupials 34,54,71, suggests that head-first descent imposes greater demands on maintaining secure grips, leading to more cautious movements regardless of support diameter, characterized by increased limb contact times and reduced velocity 48,53. However, when employing side or tail-first descent strategies, primates did not significantly reduce absolute speed nor markedly increase duty factors, but instead decreased relative distance traveled per cycle. This unexpected result suggests that tail-first and side descent strategies are highly effective, allowing primates to maintain comparable speeds while adopting flexible manual and pedal postures 19,88. During vertical ascent, locomotion primarily relies on hindlimb propulsion 8,89, whereas head-first descent requires increased forelimb braking and frictional control at the extremities 8,50,51,64,90,91. Our results suggest that the upright descent strategies employed by primates mitigate these braking demands on the forelimbs while preserving static stability and high descent velocity.

Moreover, we found that gait patterns are strongly constrained by phylogenetic history, consistent with previous findings in primates 63,67,8587. Primates, marsupials, rodents, and carnivorans exhibited both symmetrical and asymmetrical gaits, though primates and marsupials employed the highest proportions of symmetrical gaits, particularly DSDC and LSDC gaits. This aligns with earlier reports of habitual DSDC gait use in primates, some marsupials, and carnivorans, specifically the kinkajou Potos flavus 49,56,60,61,71,72,92. Interestingly, only strepsirrhines increased their use of LSLC gaits during descents, as previously observed in platyrrhines on sloped supports 51,6367, whereas rodents and carnivorans favored more DSDC gaits during descents compared to more LSDC gaits in ascents. Both DSDC and LSLC gaits have been proposed to enhance stability by facilitating precise fore- and hindfoot placement, especially during horizontal walking on narrow branches 49,60,93. On the other hand, LSDC gaits promote stability by keeping the location of the center of mass relative to the support polygon of the limbs 60,61,66,9395. Most previous studies have focused on symmetrical gaits, particularly DS, due to their proposed adaptive significance for fine-branch locomotion in primates 56,60,63. However, several studies have also reported increased use of asymmetrical and LS gaits on sloped and vertical supports in primates and other arboreal mammals, emphasizing the importance of behavioral flexibility rather than strict kinematic specialization in arboreal environments to cope with complex 3D environments 12,53,63,64,92,96.

We found that the proportion of symmetrical gaits increased on small vertical supports, partially validating our second hypothesis (H2), but only during ascents. During descents, the proportion of asymmetrical gaits significantly increased across all supports and species, regardless of descent strategy. This key finding suggests that vertical locomotion requires distinct kinematic adjustments: symmetrical gaits enhance climbing efficiency, whereas asymmetrical gaits facilitate stability during descents. Symmetrical gaits, characterized by regular stance and swing phases, allow for controlled force application, lower reaction forces, control and transfer of the moments and torques imposed on the body axes, and minimized mediolateral deviations, keeping the center of mass within the support polygon throughout the stride 56,6062. Conversely, asymmetrical gaits permit higher speeds, reduce metabolic cost, provide an energetic input to locomotion through the sagittal movement of the spine, and stabilize locomotion by reducing peak forces and center of mass oscillations through shorter, more frequent strides 24,52,53,5759. Overall, while symmetrical gaits enhance stability on narrow supports during horizontal or oblique locomotion and during vertical climbing 11,52,62,64,69,72, asymmetrical gaits may offer superior control and security during vertical descents by enabling simultaneous grasping with both fore- and hindlimbs 5153,5759,68.

The correlation between morphological adaptations and vertical descent strategies

We identified a clear influence of body proportions on head-first descent behavior in our dataset, with body mass, relative forelimb, hindlimb, and tail lengths, and relative head size all significantly affecting descent strategies. Heavier animals with relatively larger heads as well as elongated limbs and tails were more likely to avoid head-first descents, partially supporting our third hypothesis (H3). As predicted, the intermembral index (IMI) was negatively correlated with the proportion of head-first descents 36,46. This relationship was largely driven by strepsirrhines, which exhibited particularly low IMI values compared to other taxa. Strepsirrhines also performed the fewest head-first descents and the highest proportion of tail-first descents, highlighting the impact of hindlimb elongation on vertical locomotor performance by shifting the center of mass posteriorly. In contrast, platyrrhines displayed higher IMI values (> 80), and all non-primate species exhibited relatively shorter limbs (forelimbs < 55% and hindlimbs < 70% of body length) compared to primates. These findings support previous hypotheses that elongated hindlimbs may increase the risk of imbalance during head-first descent, by simultaneously increasing the forelimbs’ braking demand and the hindlimbs’ traction requirements 36,46.

We also detected a significant relationship between tail length and descent strategy: primates generally possessed proportionally longer tails (> 90% of body length) than the non-primate species. Although tails primarily function in balance during horizontal arboreal locomotion 9799 and serve to stabilize the body during leaping or aerial maneuvers 100, increased tail length may also contribute to a posterior shift of the center of mass that discourages head-first descent. Interestingly, we also found a negative correlation between relative head size (approximated by the EQ) and the proportion of head-first descents, suggesting that larger head mass may further bias animals toward more upright postures during vertical descent. To our knowledge, this study is the first to explicitly assess the influence of body mass distribution—and particularly relative head mass—in the context of arboreal locomotion. Within our dataset, primates, especially platyrrhines, exhibited higher encephalization quotients (EQ > 1.5) and engaged in head-first descent postures less frequently than other mammals, which possessed lower EQ values (< 1.2) and relied more on head-first strategies. A larger brain not only demands increased cerebral blood flow, which may impose physiological limitations during acrobatic arboreal behaviors, but also increases head mass, thereby raising the moment of inertia and shifting the center of mass anteriorly. These factors support our hypothesis that an enlarged head may mechanically constrain animals to avoid head-first vertical descent by increasing the risk of forward toppling during head-first postures, particularly in primates. However, it is noteworthy that certain larger arboreal primates (> 5 kg), such as howler monkeys 43 and orangutans 101, both characterized by higher encephalization and, in the case of orangutans, markedly elongated forelimbs and elevated intermembral indices 102, have been observed performing head-first descents in their natural environments. Although these behaviors appear relatively infrequent, they underscore the need for further investigation into how brain expansion and associated morphological changes interact to influence arboreal locomotor strategies across taxa.

Contrary to our expectations (H3), we did not find a significant correlation between the relative size of the autopods and the proportion of head-first descents. However, within primates, we found a negative correlation between the hallux length and the proportion of side descents, with the longer the hallux the less side descents employed. Digital elongation of the hand and foot in primates has long been interpreted as enhancing grasping efficiency on arboreal supports 13,19,21,25,73,103. During vertical ascent the superior point of attachment (i.e., the hands) ensures a critical hold, whereas in head-first descents both hands and/or feet may ensure holding. In contrast, tail-first or side descents may allow primates to rely more extensively on the foot for secure hold, supporting most of the weight. Strepsirrhines possess a highly specialized foot with a very divergent elongated hallux, while some platyrrhines possess a less elongated hallux but bear claw-like nails on their distal phalanges which they use to cling to tree bark similarly to other mammals. These differences in hand and foot morphology among primates are reflected in their vertical descent behavior: strepsirrhines more frequently adopt tail-first descents, relying more heavily on the foot as a secure point of support, whereas platyrrhines more often employ side descents, using both hands and feet to grasp the support efficiently. Taken together, our results highlight that efficient vertical descent is a complex phenomenon governed by the interplay of kinematic adjustments, locomotor flexibility, body mass distribution, limb proportions, grasping ability, and shared evolutionary history among species. Primates exhibit particular adaptations for vertical descent, including upright postures and a remarkable capacity for locomotor flexibility, which may derive from their ancestral adaptation to navigating supports of varying size and orientation 63. In contrast, other arboreal mammals appear to have been more restricted to the use of less diversified supports during their evolutionary histories.

Implications for the reconstruction of early primates locomotion

Our inference model reconstructed fossil species as employing relatively high proportions of head-first descent (> 75%), except for the two medium-sized adapiforms Darwinius masillae and Europolemur kelleri (estimated at ∼30% and 60%, respectively), which already displayed relatively long legs and tail. This validates our fourth hypothesis (H4) and suggests that postural adaptations for vertical descent in primates evolved progressively. Although some species already reached relatively high body mass, such as the adapiform Notharctus tenebrosus (∼9 kg, 104, the plesiadapiform Plesiadapis sp. (∼2 kg, 105), or the rodent Ischyromys typus (∼1.3 kg, 106), they retained low EQs and moderate limb elongation. Interestingly, the plesiadapoids Carpolestes simpsoni and Plesiadapis insignis, whose phylogenetic positions within Euarchonta remain debated 78,107109, have been reconstructed with a high proportion of side descent, associated with their slightly elongated hallux similar to extant platyrrhines. Plesiadapiforms are central to the debate on primate origins as the family Plesiadapidae is often regarded as a stem or sister group to euprimates 108,110,111. However, plesiadapiforms are mainly reconstructed as generalized quadrupedal walkers and climbers 20,110,112,113. In contrast, early primates (including both omomyiforms and adapiforms) are all reconstructed as vertical climbers and leapers at variable degrees 114116. Nonetheless, our evolutionary analysis corroborates current phylogenetic hypotheses suggesting that euprimates and plesiadapiforms share an adaptive regime distinct from the other groups studied [best OUM model: euprimates + plesiadapiforms vs. other], highlighting that descent strategy on vertical supports was tightly constrained throughout evolution.

These results lead us to suggest a refinement of the existing scenario on the order of acquisition of early primate specializations. Considering that early euarchontoglires were probably small to very small (≤ 100g 109) with shorter hindlimbs, autopods, and reduced brain size, it is plausible that they primarily used mostly head-first descents and asymmetrical gaits on vertical supports 12,53. As euprimates evolved enhanced grasping abilities, elongated hindlimbs and tail, and eventually larger brains, they likely began to adopt side and tail-first vertical descent postures. The evolutionary success of early primates likely resided in their ability to colonize environments which are mechanically more challenging and difficult to access by other animals, such as horizontal terminal branches and small vertical supports. The frequent use of vertical supports, such as trunks, lianas, climbing plants, or even tall grass thickets, favored the early acquisition of a divergent hallux and pollex, improving grasping ability 15,18,19. These supports were abundant in the tropical forests of the northern hemisphere during the Paleogene, the likely origin place and time of primates 78,117,118. The use of these supports may have facilitated rapid access to the canopy, provided platforms for scanning for food or predators, efficient escape routes and rapid changes of height, and possibly privileged access to a greater diversity of food sources 119. Autopodial innovations were followed hindlimb elongation, likely enhancing leaping abilities as inferred in various early euprimate fossils 16,120. Altogether, these hindlimb specializations, followed by a progressive increase of the body and tail lengths, and brain size, probably altered body mass distribution. This may have further promoted side- and tail-first vertical descent strategies, allowing efficient vertical navigation while maintaining environmental awareness. Such adaptations may have thus favored primates to frequently adopt upright postures when negotiating vertical supports in an early stage of their evolutionary history, while increasing the use of symmetrical gaits to negotiate other types of supports such as narrow horizontal branches. Building on previous hypotheses on autopodial grasping evolution of primates 21,23,121,122, upright orthograde postures would have also freed the hand from locomotion, while the foot remained specialized for efficient support holding. Functional differences between hand and foot have emerged early in primates 5,90,123,124, with the anchoring foot playing a fundamental role in body support, and the hand redirected toward other activities, such as food grasping and manipulation, substantiated by the great diversity of manual postures in extant species 19,125,126.

Limitations of experimental conditions and future directions

Every study that includes experiments with animals housed in zoos displays some limitations. Some enclosures required enrichment with new supports and in some cases the animals did not readily use the support provided (see the Materials and Methods section). It is therefore possible that this experimental bias may have had an impact on the proportion of descent behaviors and the kinematics quantified for some species. For instance, although raccoons and coatis have been observed to descend head-first on small and medium supports in the wild (personal observations) they did not use them in their enclosures. Thus, quantifying vertical descent strategies as well as the relative frequency of vertical support use in the wild would provide more conclusive evidence regarding the vertical descent repertoire in those species. Moreover, comparing kinematics of vertical ascents and descents in more species of larger body mass (> 5 kg) with varying morphological proportions, such as apes and large carnivorans, would help to better understand the complex parameters at play regarding mechanical adaptation to arboreal locomotion in larger mammals. Such studies will also allow to develop reconstruction models better adapted to large fossils, particularly those primates with relatively lower EQ than extant species, thereby refining our understanding of the evolution of those adaptations within this group.

Materials and Methods

Experimental design

Animals studied

We collected data on animals housed in zoological parks except for the forest dormice Dryomys nitedula which were wild caught (Table S1). In zoos, individuals were video recorded directly in their enclosures. All the enclosures were large enough to allow unconstrained locomotion of the animals. We placed non-treated wooden branches of variable diameters and orientations in each enclosure before filming to ensure equal access to a variety of support types. All studied individuals were adults, in good shape, and did not display any stereotypical behavior before or during the experiments. The two individuals of Dryomys nitedula were caught in their natural habitat in Vasilika, Greece, with the use of specialized wooden nest boxes positioned in trees (permit by the Hellenic Ministry of Environment and Energy 8817/242). They were immediately transferred and housed to glass terrariums (L= 100; W= 40; H= 50 cm) enriched with tree branches at room temperature and with normal photoperiod (complying to the regulations of the Ethics committee of the School of Biology, Aristotle University of Thessaloniki, Greece). After observations and filming directly in these terrariums, the animals were released at the same location they were originally captured. Animal observation was performed in compliance with the International Primatological Society (IPS) Ethical Guidelines for the Use of Nonhuman Primates in Research and the Association for the Study of Animal Behavior (ASAB) and the Animal Behavior Society (ABS) ethical guidelines for the use of animals in research 127.

Behavioral data collection

We focused on vertical supports (90°± 5°). We categorized support diameters into three categories: small (support circumference < hand and foot lengths), medium (support circumference ≈ hand and foot lengths), and large (support circumference > hand and foot lengths) following previously described methods 19, resulting in three support types tested for each individual. Sequences were recorded using a portable video camera (Panasonic HC-V770 camcorder 120fps, full HD 1080p). To ensure both a complete view and close ups of the animals during locomotion, we also used three small action video cameras (Mobius ActionCam 60fps, HD 720p) installed in three different view angles (frontal, lateral and ventral) and in relatively close distance (between 20 and 50cm depending on the size of the animals) from each support type. After a habituation period, the individuals were recorded using an alternation of scan sampling and focal-animal sampling methods 128 in continuous recording sessions of 10 to 30 min. Animals were filmed as they moved freely over supports, and we stimulated them by providing small pieces of food placed randomly along the supports when necessary. For nocturnal species (Table S1), observations were conducted either in artificial or real nocturnal conditions with addition of red lighting to enable recording without disturbing the animals.

Video analysis

We analyzed videos using Adobe Premiere Elements 12. All video analysis was performed by the same person (SLDT) to avoid inter-observer biases. We focused on locomotion only and excluded other behaviors. We analyzed a minimum of three and up to ten passages in ascent and ten passages in descent for each individual on each support type (i.e. small, medium, and large), resulting in up to thirty ascents and thirty descents analyzed for each individual, with the following exceptions: we could not obtain data for Tupaia belangeri on small supports, for Graphiurus murinus on small supports; for Typhlomys chapensis on medium supports, for Procyon lotor on small supports, and for Nasua nasua on the medium and small supports. For Tupaia belangeri, Graphiurus murinus, and Typhlomys chapensis, these supports were not available during the data acquisition process. However, for Procyon lotor and Nasua nasua, these supports had been added and were available in the enclosures, but were not used by the animals.

Kinematic analysis of locomotion

To quantify the absolute speed of the individuals for each sequence of locomotion we extracted two snapshots and their associated frame number from the beginning and the end of the movement, using videos with a lateral view relative to the support. We merged both images using the product layer functionality in Adobe Photoshop CS6. We then measured in ImageJ the distance travelled by the individuals during each stride using markers of known length positioned in each support as references. We calculated the associated speed as the ratio of this distance to the duration separating the snapshots.

We also calculated the relative speed for each sequence of locomotion as the ratio of the absolute speed to the body length of each animal studied, defined as the total length between the top of the nose to the base of the tail. We used snapshots of each individual to measure body lengths in ImageJ, using a scale marker as a metric reference positioned on top of the branches.

We analyzed the strides for each sequence of locomotion by recording the frame number of the first touchdown, the following lift-off and the subsequent touchdown of each limb on the support. We calculated the associated stride duration of each limb, using the known frame rate of the associated videos. We then calculated the duty factor (DF) of each limb as the percentage of time spent in contact with the support relatively to the total stride duration. Forelimbs and hindlimbs duty factors were averaged to obtain the mean DF.

We classified the gait type following previously described methodology 57. We used the midpoints of each limb contact intervals (i.e., the temporal midpoint between first touchdown and following lift-off). We categorized gaits as symmetrical when the temporal lag between left and right forelimbs and hindlimbs midpoints were both ≥ 40% and ≤ 60% of the total stride duration, and asymmetrical otherwise. Further categorization among symmetrical gaits were made using the limb phase, which corresponds to the percentage of total stride duration separating a hindlimb touchdown from its ipsilateral forelimb touchdown. We thus calculated the left and right limb phases using the midpoints of forelimbs and hindlimbs contact intervals. Left and right limb phases were then averaged to obtain the mean limb phase for each stride. We categorized as lateral sequence – lateral couplet (LSLC) when the mean limb phase was ≥ 0% and < 25%, as lateral sequence – diagonal couplet (LSDC) when ≥ 25% and < 50%, as diagonal sequence – diagonal couplet (DSDC) when ≥ 50 and < 75%, and as diagonal sequence – lateral couplet (DSLC) when ≥ 75% and ≤ 100%. Classification of asymmetrical gaits used the stance period, which corresponds to the percentage of limb pair total contact duration separating left and right limbs of the same pair. We classified as bounds the gaits when both forelimb and hindlimb stance periods were ≤ 10 %, as half-bounds those when either (i) the forelimbs stance period was ≥ 10% and the hindlimbs stance period was ≤ 10%, or (ii) when the forelimbs stance period was ≤ 10% and the hindlimbs stance period was ≥ 10%, and as gallop when both forelimb and hindlimb stance periods values were ≥ 10 %.

Morphological data collection

For each studied individual, we extracted three to five lateral view screenshots from the videos. On these images we measured lengths in ImageJ allowing to calculate eight intrinsic body proportions (Table 1). We calculated the average proportions over screenshots for each individual, and ultimately the average proportions by species.

The mean body mass of each species studied was gathered from the literature (Table S1). We could have used the actual body masses of certain individuals studied in zoos, but we could not have access to this information for most species, so we decided to use the standard reference values from literature. Additionally, we calculated the mean encephalization quotient (EQ) of these species as a proxy for their relative head to body size ratio and to allow comparison with fossil specimens for which the EQ is commonly studied and provided in the literature (Table S2). Prior to the calculation of the EQ in extant species, we collected the mean endocranial volumes (ECV) either from the literature or from own measurements in the case of species that were not previously documented (Table S3). For these species we quantified the ECV from preserved skulls of specimens housed in the anatomical collections of the Museum für Naturkunde of Berlin, by filling their endocranial cavity with chia seeds of known volumic mass to infer their volume. We thus provide hereby the mean endocranial volumes for the treeshrew Tupaia belangeri, the rodents Dryomys nitedula, Graphiurus murinus and Typhlomys cinereus, and the marsupial Petaurus breviceps (Table S3).

For comparative purposes, we collected from the literature the reconstructed body mass and EQ of key fossil specimens of euarchontoglires with relatively complete postcranial and/or cranial remains (Table S2). Whenever possible, we also measured the same postcranial lengths as for extant species, using published photographs (see Table S2 for associated references, and Table S4 for body proportions calculated in fossils).

Construction of phylogenetic tree including fossil taxa

For illustration purposes in Fig. 1C, we used a phylogeny of extant species from timetree.org 129 which derives from an aggregation of molecular phylogenetic data.

For phylogenetic comparative analysis (phylogenetic signals and evolutionary model comparisons) we generated 1000 trees from the tip-date DNA-only version of mammalian trees of the studied extant species from vertlife.org 78. We computed a consensus tree using the consensus function of the R package ape and made it strictly ultrametric by extending the tips with the force.ultrametric function from the R package phytools. To build and compare evolutionary models of the fossils’ descent behaviors, we supplemented this tree of extant species with the extinct taxa studied (Fig. 5) using estimated age and phylogenetic topology from literature (Table S2).

Statistical analysis

Statistical analysis was performed with RStudio version 1.3.1093 (R version 3.6.3).

Proportions of descent strategies

To assess the effect of support diameter (e.g., small vs. medium vs. large) on the proportion of descent types between individuals of the same species, we performed either analysis of variance (ANOVA) or multivariate analysis of variance (MANOVA) tests (respectively aov and manova functions in R) depending on the number of descent strategies exhibited (Fig. 1B, SuppTables A). To even out the testing power among species for which we had a different number of individuals, we bootstrapped 10 individuals and 10 descent events on each support category to perform the tests. P-values were then Benjamini-Hochberg corrected.

To assess the extent of phylogenetic information contained in the proportion of head-first descents on vertical supports, we computed Blomberg’s K phylogenetic signal 130 using the physignal function from the geomorph package in R with 50000 iterations.

Locomotion

To test for the global differences of speed, relative speed, and duty factor between ascents and each descent type, we performed two-sided Mann Whitney U tests on 40 bootstrapped events (20 ascents and 20 descents) for all species on each support diameter. P-values were then Benjamini-Hochberg corrected. We also calculated the mean percentage of variation between ascents and each descent category (Fig. 2, SuppTables B, C, D). For tests over duty factor, we considered both forelimbs and hindlimbs combined.

We calculated the frequency of all possible gait type values occurring for each individual on each support category differentiating ascents and the different descent strategies when relevant. These frequencies were PCA-transformed to obtain linearly independent values which serve as low-dimensional representation of the overall gait data. Using the same physignal function as described above for phylogenetic signal computation of head-first proportion on vertical supports, we computed the Kmult statistic 131 with 50000 iterations on the PCA-transformed data of gaits, keeping only the first principal components that retain 95% of variance, to assess the extent of phylogenetic information contained in gaits.

To check for enough similarities between species of the same phylogenetic group as defined in Fig. 1C, and dissimilarities between groups to justify pooling data at the group level, we performed MANOVAs with Benjamini-Hochberg corrected post hoc pairwise comparisons on PCA-transformed frequencies of gaits, keeping principal components explaining at least 95% of variance, using the manova function from the stats package in R and Wilks’s test (SuppTables E).

To test for changes of the fraction of symmetrical versus asymmetrical gaits between ascents and descents for each group of species and on each support diameter separately, we performed two-sided Wilcoxon ranked signed tests on bootstrapped samples. For each group and each support, we took 100 bootstrapped species and sampled 20 events (10 ascents and 10 descents) at random. The resulting p-values were Benjamini-Hochberg corrected (Fig. 3, SuppTables F).

Morphological correlates

To control that brain allometry is homogeneous among all phylogenetic groups, to be able to compare EQ between species, we computed Pearson correlations between the Log transformed brain and body masses, over all species and by phylogenetic group (Fig. S2).

To assess the level of correlation among the 10 morphological variables, we computed Spearman correlations on all pairs. P-values were then Benjamini-Hochberg corrected (SuppTables G).

We computed Spearman correlations to quantify the relation between morphological variables and the proportions of head-first descent (for all species) or side descent (for primates) on vertical supports (all diameters combined) and the resulting p-values were Benjamini-Hochberg corrected (Fig. 4, Fig. S3, SuppTables H, I).

Estimation of extinct species vertical descent strategy

We isolated the six morphological variables for which we obtained a significant correlation with head-first descent proportions on vertical supports (body mass, EQ, relative arm length, relative leg length, intermembral index and relative tail length). We used the Schafer’s Multiple Imputation procedure from the R package norm to compute estimations of the fossils’ proportions of head-first descent given their respective morphological features (Fig. 4 A, B, C, Fig. S3 A, B, C, SuppTables J). We simulated 5000 runs of Markov Chains with 50 steps to infer the head-descent proportions of these 13 extinct species. We conducted the same procedure to then infer the proportion of side descent in euprimates and plesiadapiforms fossils by restricting the set of extant species to primates, and using the only morphological variable which correlated with side descents (relative hallux to foot length, Fig. 4 D, SuppTables J).

Evolutionary model comparison for head-first descent strategy

To evaluate if the proportion of head-first descent strategy on vertical supports follows a specific evolutionary model, we compared different quantitative models of trait evolution in all extant and extinct taxa studied. We first fitted three simple models: Early Burst (EB), Brownian Motion (BM) and Ornstein-Ulhenbeck (OU) using the mvEB, mvBM and mvOU functions from the R package mvMorph respectively 132. Then, we considered models with two distinct optima within the phylogeny. We fitted the Brownian Motion and Ornstein-Ulhenbeck models with multiple optima (respectively BMM and OUM) using the mvBM and mvOU functions. We tested a change of the evolutionary optimum at different phylogenetic nodes: I) euprimates, II) euprimates + plesiadapiforms, III) euarchontans (euprimates + plesiadapiforms + scandentians), and IV) euarchontoglires (euarchontans + rodents). We compared models using the Akaike Information Criterion (AIC) and estimated the phylogenetic half-life in the OUM models that best fitted head-first descent strategies in our sample, using the halflife function from mvMorph, after rescaling the phylogenetic tree height to 1 (SuppTables K).

Data and materials availability

All aggregated data needed to evaluate the conclusions of this study are included in this paper and its supplementary materials. Raw video data are available upon reasonable request which should be addressed to SLDT and DY. The R scripts used in this study are also available upon request and should be addressed to SLDT.

Acknowledgements

We greatly thank the staff of the Zoological and Botanical Parc of Mulhouse, the Zoological Park of Paris, and the Spaycific’Zoo (France) for access to their animals and for their assistance during material installation and data collection. We thank Dr. A. Panyutina and Dr. A. Kuznetsov (Lomonossov Moscow State University and Russian Academy of Sciences) and the Moscow Zoo staff (Russian Federation) for video recordings of Graphiurus murinus, Typhlomys chapensis, and Tupaia belangeri. We thank Dr. N.E. Karantanis and the staff of the Nowe Zoo Poznan (Poland) for the video recordings of Xanthonycticebus pygmaeus. We thank the staff of the Laboratoire d’Ecologie Generale in Brunoy (France) for the video recording of Caluromys philander. We thank C. Funk from the Museum für Naturkunde of Berlin for access to the collection of mammal specimens. We thank Dr. A. Llamosi for his help with statistical analysis.

Additional information

Fundings

Research funding was from the Fyssen Foundation and the Alexander von Humboldt Foundation postdoctoral grants (to SLDT) and from the German Research Foundation grant DFG NY63 2/1 (to JAN). Travel funding for behavioral data collection in Poland and The Russian Federation was from Erasmus International Mobility funds (to DY).

Author contributions

Conceptualization: SLDT

Methodology: SLDT, DY, JAN

Investigation: SLDT, DY, JAN

Visualization: SLDT

Funding acquisition: SLDT, DY, JAN

Writing—original draft: SLDT

Writing—review & editing: SLDT, DY, JAN

Additional files

Supplementary file 1

Supplementary file 2