(C) PLOS One This story was originally published by PLOS One and is unaltered. . . . . . . . . . . Bipedalism and brain expansion explain human handedness [1] ['Thomas A. Püschel', 'Institute Of Human Sciences', 'School Of Anthropology', 'Museum Ethnography', 'University Of Oxford', 'Oxford', 'United Kingdom', 'Rachel M. Hurwitz', 'Chris Venditti', 'School Of Biological Sciences'] Date: 2026-05 Abstract Humans exhibit a striking and near-universal population-level right-hand preference, an evolutionary singularity unmatched among primates. Despite its pervasiveness, the origins of this lateralization remain poorly understood. Here, we combine phylogenetic comparative methods with meta-analysis to investigate manual lateralization across 41 anthropoid species (n = 2,025), testing longstanding eco-evolutionary hypotheses for handedness direction (mean handedness index, MHI) and strength (mean absolute handedness index, MABSHI). Our models reveal significant phylogenetic signal for both traits and identify Homo sapiens as an evolutionary outlier, exhibiting exceptional rightward bias and strength relative to phylogenetic expectations. However, this outlier status disappears when brain size (endocranial volume) and intermembral index are included, suggesting these factors are central to the emergence of human handedness. We also show that high MABSHI evolved early in hominin evolution, while MHI increased to unparalleled levels with the appearance of the genus Homo. Our findings identify bipedalism and neuroanatomical expansion as likely key drivers of uniquely human lateralization, while also revealing broader ecological patterns shaping handedness across primates. This work provides a framework for disentangling human-specific adaptations from general primate trends in the evolution of behavioral asymmetries. Citation: Püschel TA, Hurwitz RM, Venditti C (2026) Bipedalism and brain expansion explain human handedness. PLoS Biol 24(4): e3003771. https://doi.org/10.1371/journal.pbio.3003771 Academic Editor: Catherine Hobaiter, University of St Andrews, UNITED KINGDOM OF GREAT BRITAIN AND NORTHERN IRELAND Received: October 23, 2025; Accepted: April 8, 2026; Published: April 27, 2026 Copyright: © 2026 Püschel et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Data Availability: All relevant data are within the paper and its Supporting information files. Funding: This work was funded by a Leverhulme Trust (https://www.leverhulme.ac.uk) Research Leadership Award to CV (RL-2019-012). The funder had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist. Abbreviations: ECV, endocranial volume; IMI, inembral index; MABSHI, mean absolute handedness index; MHI, mean handedness index Discussion We present a comprehensive phylogenetic comparative meta-analysis of handedness in anthropoids that accounts for shared evolutionary history, effect variation across multiple sources, and potential publication bias, thus addressing many previous limitations. By testing multiple hypotheses about the evolution of handedness simultaneously, we offer new insights into the role of evolutionary selection in shaping both the direction (MHI) and strength (MABSHI) of manual lateralization. Our results should be understood with the caveats which come from collating data from multiple sources, and that handedness as determined by standardized protocols may not fully reflect the complexity of primate hand use and manual preferences in the wild [4,34,35]. Nevertheless, the meta-analytical framework applied here mitigates to some extent data heterogeneity by modeling between-study variance and applying inverse-variance weighting, thereby giving greater influence to more precise estimates [42]. Thus, by integrating behavioral data using a meta-analytical and phylogenetic comparative approach [37], we show that humans occupy an extreme position within the primate order with respect to both the strength and direction of handedness. Our results highlight significant phylogenetic signal in both MHI and MABSHI, indicating that handedness traits have evolved with descent along the branches of the primate phylogenetic tree. However, the data do not fully support any single existing hypothesis, with most predictors only showing significance when H. sapiens is included. This suggests that many hypotheses regarding handedness may be overly anthropocentric. As such, future research should seek to clearly distinguish between human-specific and broader primate explanations for manual lateralization as it extends to hand preferences, manipulative abilities, and communicative influence [43]. A striking finding of this study is the identification of H. sapiens as an outstanding outlier in MABSHI but especially MHI. Humans display a pronounced right-handed bias (MHI = 0.76), which contrasts sharply with the phylogenetic prediction of the reduced model excluding humans (MHI = 0.0) (Fig 3c). Likewise, humans show extreme handedness strength (MABSHI = 0.94), with MABSHI values near the highest observed among anthropoids, although some arboreal species, such as spider monkeys and langurs, also exhibit high levels of lateralization strength. These results strongly support the hypothesis that human exceptionalism in handedness is likely due to to strong, human-specific selective pressures. When considering the broader primate phylogeny, several predictors show credible effects for both MHI and MABSHI include factors such as IMI, ECV, diet, and social system. However, when H. sapiens is excluded from the analysis, only a limited number of predictors retain credible effects, with most posterior intervals overlapping zero. Notably, IMI, substrate use emerged as key factors. Lower IMI (indicative of more quadrupedal or leaping locomotion) correlate with leftward bias in MHI, and arboreal species generally exhibit higher MABSHI values. Arboreal species, needing flexible stabilization and precise movements in the canopy, often use one hand for support (left or right being arbitrary) and benefit from higher MABSHI for enhanced motor control. Terrestrial species, facing comparatively simpler and more uniform substrates, exhibit less complex locomotor patterns, requiring less specialized and lateralized motor skills owing to more repetitive and straightforward movements. This highlights the importance of locomotor and ecological factors in shaping handedness across non-human primates, whereas our large brain size emerges as an important link to human’s strong directionality in handedness. These results imply that our unusual gait was the main initial driver of our exceptional handedness strength with our large brain more linked to the directionality. In humans, the evolution of bipedalism and the subsequent freeing of the hands may have intensified selective pressures for stronger hand preferences. This finding aligns with previous research linking bipedal posture to the evolution of more pronounced handedness in fossil hominins, further suggesting that the evolutionary trajectory of human handedness is rooted in our unique locomotor adaptations [44–47]. Our predictions for other hominin species using our reduced models provide a more precise temporal perspective on the evolution of handedness direction and strength in human evolution (Fig 4). Our results indicate consistently high handedness strength (MABSHI) among hominins from early on in the lineage. Although hominins became increasingly terrestrial throughout their evolution [48], a locomotor behavior that reduces MABSHI among other anthropoids, the persistence of strong handedness strength is likely explained by our unique mode of locomotion—bipedalism—which freed the upper limbs and enabled specialized manual behaviors [49], maintaining selective pressure for lateralized hand use. This consistently high MABSHI level may also offer some support for those views that consider that that adaptations for bipedalism arose in an arboreal context [50], as more arboreal species tend to show higher handedness strength values. In contrast, handedness direction (MHI) follows a different evolutionary pattern (Fig 4). Early hominins such as Ar. ramidus and Au. afarensis exhibit MHI values relatively similar to other great apes. It is with the emergence of the genus Homo, and particularly the onset of significant encephalization, that we observe a marked increase in MHI values, reaching levels that are highly unusual among anthropoids. This pattern suggests a novel link between the evolution of directional handedness and increasing brain size [51], highlighting a link between these two human evolutionary hallmarks. In this context, the emergence of pronounced right-handedness bias in humans may be viewed as part of a broader suite of neurological and behavioral specializations tied to the unique cognitive trajectory of our lineage. Importantly, this connection between encephalization and handedness direction has not been demonstrated in such a temporally resolved, phylogenetically informed framework before, offering new insights into the deep evolutionary roots of human lateralization. The intriguing result for H. floresiensis, which shows relatively low MHI despite its placement within the genus Homo, may be explained by its unusual combination of small brain size and a unique locomotor repertoire blending bipedalism with arboreality. While the pelvis and lower limbs of H. floresiensis exhibit clear adaptations for upright walking, features such as long feet, an elongated forefoot, and curved phalanges suggest a locomotor pattern more similar to Australopithecus, including climbing behaviors [52]. This unique morphology may have reduced handedness direction evolution in this lineage. However, further analyses are necessary to confirm this hypothesis. While only studies that employed the standardized “tube task” to calculate MHI and MABSHI values were included, and analyses were designed to minimize individual study bias and account for sample heterogeneity, some limitations must be acknowledged. Direct comparisons between human and non-human primate handedness are inherently constrained by differences in data collection, including potential variation in exact tasks, trial numbers, and sample sizes across species, as well as by the fact that human handedness data are often derived from questionnaires or observations of adults with strongly established, and likely culturally reinforced, manual preferences [16,43]. Importantly, handedness is also context-dependent: in humans, strong right-handedness is most consistently observed in goal-directed object manipulation, whereas its strength may weaken or even reverse in other behavioral contexts (e.g., repetitive gross motor or social actions) [43,53–56]. This further complicates cross-species comparisons, as different experimental paradigms may capture distinct aspects of lateralized behavior [36]. In non-human primates, for example, stronger right-handed biases are often reported in tasks involving complex bimanual coordination, including tool use and gestural behaviors [56,57]. These limitations reflect the nature of the available data rather than shortcomings of the present study. Nevertheless, by restricting analyses to a single, widely used behavioral paradigm in non-human primates and applying meta-analytical methods that explicitly model between-study variance [42], we aim to maximize comparability across taxa. Concerns regarding the independence of data points are warranted, as the reuse of subjects across studies cannot be entirely ruled out in the existing literature. Where possible, we minimized this risk through careful screening of study samples at the individual level, which applied to the vast majority of our dataset. In any case, any residual non-independence would be expected to inflate variance rather than systematically bias effect estimates, rendering our conclusions conservative. Furthermore, the effects of culture are difficult to fully account for in comparative analyses of handedness. Although some non-human primates exhibit forms of cultural transmission [58–60], humans are unique in displaying cumulative cultural evolution, which may amplify or stabilize behavioral asymmetries. This raises the possibility that greater variation in the strength and direction of handedness exists than is captured in the available data [61]. However, the predominance of right-handedness in humans is a robust cross-cultural phenomenon [8,62], and no left-dominant human society has been documented, making it unlikely that culture alone can account for the evolutionary trajectory of handedness in our species. Taken together, it is plausible that bipedalism and encephalization did not operate as independent drivers of handedness but acted synergistically. The initial adoption of an upright gait freed the upper limbs, creating novel opportunities for tool use, gestural communication, and other fine motor behaviors in which lateralization would have conferred performance advantages. Although terrestrial species generally exhibit lower handedness levels, humans may represent an alternative expression of the arboreal handedness pattern, despite being primarily terrestrial. The freeing of our hands from locomotor constraints allowed them to be deployed asymmetrically more frequently and extensively than in any other anthropoid lineage. This likely had precursors in an ancestor engaged in significant arboreal activities, as evidenced by the hominin fossil record, which indicates that arboreal behaviors were crucial not only at the origin of bipedalism but also in later lineages [63,64]. Concurrently, increases in brain size and associated cortical reorganization may have promoted greater hemispheric specialization, thereby enhancing the neural efficiency of such lateralized behaviors, especially after the emergence of the genus Homo. In this context, the initial locomotor change prompted by bipedalism can be seen as providing ecological and anatomical opportunities for manual specialization, while encephalization may have later reinforced and further canalized population-level patterns of lateralization. Furthermore, culture may have acted concurrently with or amplified the effects of this emerging trajectory of hominin right-handedness. While our analysis provides substantial evidence for the role of both ecological and anatomical factors in the evolution of handedness, it also raises important questions for future research. Notably, our findings underscore the need to refine hypotheses that distinguish human-specific factors, such as cumulative cultural evolution, bipedalism, or ECV, from broader primate trends in terms of handedness. Additionally, expanding this analysis to include non-primate taxa which display clear limb lateralization, such as parrots or kangaroos, would be valuable in investigating the potential for convergent evolution in handedness across species [65–67]. Finally, expanding the fossil sample size, by, e.g., using other handedness proxies such as dental wear analyses, could significantly enhance our understanding of handedness evolution in extinct hominins and contribute to more refined and robust phylogenetic models. In conclusion, this study highlights the complexity of handedness evolution and demonstrates the importance of considering both ecological and anatomical factors in understanding the selective pressures that have shaped human lateralization. Our findings support a more nuanced interpretation in which human handedness, particularly in the context of goal-directed object manipulation, reflects both continuity with broader primate trends and a uniquely pronounced degree of specialization, with humans representing an extreme at the upper end of the handedness spectrum. Future research will benefit from a more comprehensive approach that considers both the primate and broader animal phylogenies, as well as the distinct evolutionary trajectory of H. sapiens. Materials and methods We combined the two most recent and comprehensive datasets available for primate handedness to carry out our subsequent analyses (S1 Data). Soto and colleagues 2022 [36] corresponds to a meta-analysis of hand preferences in coordinated bimanual tasks in non-human primates that compiled multiple sources on non-human handedness, whilst Caspar and colleagues 2022 [15] consists of a phylogenetic comparative analysis that combined both new experimental data, as well as published sources on hand preference. To enable meaningful comparisons, we only considered studies that used the so-called “tube task” as this experimental paradigm corresponds to a simple and widely applicable test to determine primate hand preferences [13,68]. Importantly, while human handedness is often determined via questionnaire, these data use results of humans engaging in the “tube task” [43]. We only included data on MHI and MABSHI, as these two variables captured the main aspects of our phenomena of interest, namely the direction and strength of the manual asymmetries. The MHI is calculated by averaging the handedness scores of all individuals in a group, where handedness scores range from −1 (indicating complete left-handedness) to +1 (indicating complete right-handedness), with 0 representing ambidexterity. The MHI provides insight into the overall direction of hand preference within a population or sample. The MABSHI quantifies the strength of hand preference in a group, regardless of the direction of the preference (i.e., left, or right). It is calculated by taking the absolute value of the handedness scores of all individuals in the group, then averaging these absolute values. This is useful for determining the consistency or strength of hand preference within a population, irrespective of whether the preference is for the left or right hand. Therefore, we extracted all the information available for these two variables (i.e., effect sizes, variance measures and sample sizes) from [36] and [15] being careful of not duplicating sources that were referenced in both studies. The final dataset included information from 2,025 individuals across 41 species of anthropoid primates. It is important to bear in mind that in the present study, we exclusively focus on the evolutionary underpinnings of population-level hand preferences, as opposed to studies or analyses that have focused on the individual-level of variation within species. We downloaded 100 trees from http://vertlife.org/ comprising the 41 anthropoid species for which we have data to be used in throughout our analyses. These phylogenies correspond to Bayesian-inferred trees built using a “backbone-and-patch” approach based on a 31-gene super-matrix [69]. These phylogenies were used in our subsequent comparative analyses. To assess the hypotheses proposed to explain primate and/or human hand preference (S1 Table) we compiled diverse relevant covariates from the literature. These variables are BM, body mass [kg]; ECV, endocranial volume [cm3]; FRUIT, percentage of fruits in diet; DIET, percentage of fruits and animals in diet; IMI, intermembral index; DIM, body mass sexual dimorphism; TOOL, tool use [0 = absence, 1 = presence); SUBS, substrate preference [0 = arboreal, 1=both, 2 = terrestrial]; SOC-SYS, social system [0 = solitary, 1 = pair, 2 = group]; CL, intra-sexual competition levels sensu [1–4,70]; EXT, extractive foraging [0 = absence, 1 = presence]; social learning sensu [71] [0 = absence, 1 = presence] (Fig 1; S1 Table; S1 Data). In the case of IMI and body mass sexual dimorphism (DIM), there were a limited number of species for which we did not find adequate data in the literature (see S1 Data for more details). Therefore, we applied a phylogenetic imputation procedure known as “PhyloPars” to deal with this missing data. Phenotypic covariance was assumed to be equivalent among species, and we also assumed Brownian motion in our imputation procedure. We used the maximum credibility clade phylogeny from our 100 trees when carrying out this procedure. Missing observations were incorporated by maximizing the log-likelihood of the covariance parameters using relevant covariates (i.e., body mass for IMI and body mass sexual dimorphism), thus allowing us to predict means and covariances for missing values at the tips of the phylogenetic tree [72]. This phylogenetic imputation procedure was carried out using the R package “Rphylopars” v.0.3.9 [73]. All the covariates that were considered in our subsequent modeling procedures were used as shown in S1 Table to test the 10 hypotheses proposed to explain handedness. Continuous predictors were logged 10 and scaled prior to the modeling steps, while the few proportion variables were logit-transformed to ensure numerical stability. All our Bayesian meta-analytical phylogenetic comparative analyses were done using the “brms” package, which provides an R interface to fit Bayesian generalized linear and non-linear multivariate models using Stan [74]. We applied generic weakly informative priors as all our continuous predictors were standardized (i.e., scaled, and centered). Continuous and categorical predictors were modeled using a normal distribution with mean 0 and variance 1. Prior predictive checks show that using wider priors (e.g., mean 0 and variance 5 or 10) did not influence our obtained results, and therefore narrower priors were preferred for computational efficiency. All predictors shown in S1 Table were treated as fixed effects. Our random effects comprised non-phylogenetic species effects (i.e., specific effects that would be independent of the phylogenetic relationship between species such as environmental or niche effects), as well as phylogenetic species-effects assumed to be sampled from a normal distribution with a mean of 0 and covariance matrix proportional to the phylogenetic correlation matrix among taxa obtained from our phylogenetic trees. Between-study heterogeneity standard deviation (τ) was explicitly modeled via the group-level effects per study. Hence, we used (τ2), which corresponds the variance of true effects, as a measure of between-study heterogeneity. All random effects were modeled using half-Cauchy distributions with location 0 and scale 0.05, as our random-effect parameters should always be non-negative but can be close to zero [75]. As our models were also meta-analytical, each dependent variable (i.e., MHI or MABSHI) had an associated measurement error defined by the standard error of every specific study. We considered an effect to be credible (i.e., “significant”) for handedness patterns when the model’s 95% credible intervals for the intercept and corresponding slopes did not include zero. Phylogenetic signal was measured using Lynch’s h2, which is equivalent to Pagel’s λ in the context of phylogenetic generalized linear mixed models [76,77]. We ran two independent chains per model using a warmup period of 4,000 out of a total of 12,000 iterations. Convergence was assessed both visually, by looking at the obtained trace plots, as well as by using rstan’s standard convergence and efficiency diagnostic metrics for Markov chains. The Gelman-Rubin values were always one for all the estimated effects in all our models [78], thus confirming that our chains were well-mixed. Bulk and Tail Effective Sample Sizes were >1,000 for all the effects estimated by all our models, hence indicating good sampling efficiency, as well as that all our estimates were reliable. We repeated all our analyses separately for both variables, MHI and MABSHI, as different hypotheses may disparately apply to manual strength and/or direction asymmetry. In addition, we also repeated all our analyses excluding H. sapiens, to assess how our highly derived handedness patterns may influence the overall pattern observed across anthropoids. To account for topological differences, we repeated all our analyses using the 100 phylogenies previously mentioned. Values reported correspond to the average of the 100 models ran for each one of our analyses. We compared those models showing at least one credible predictor using an efficient approximate leave-one-out cross-validation approach (LOO-CV) [79]. Observations that were too influential in each model and that could not be accurately approximated were estimated using an exact cross-validation. To assess the fitness of our models, we carried out posterior predictive checks by comparing 1,000 datasets simulated from each one of our models with our original data (i.e., MHI and MABSHI) (S1 and S2 Sections). All our models showed simulated datasets that resembled our original data, thus indicating the good fit of our models. We also ran additional models including, as fixed effects, all predictors that showed a credible effect in at least one of the tested hypotheses. We then applied a model reduction procedure, iteratively removing predictors whose 95% credible intervals overlapped zero. This process resulted in four reduced models (two for MABSHI and two for MHI, each including and excluding H. sapiens) that retained only covariates with credible effects. The idea was to assess what covariates played a key role in the observed handedness patterns. To evaluate the contribution of these covariates on the variance of MHI and MABSHI, we used an R2 metric that is suitable for Bayesian multilevel models [80]. To assess H. sapiens’ apparent singularity, we used our reduced models to predict the expected levels MHI and MABSHI for our own species using a “phylogenetic outlier test” [41], while also considering a meta-analytical component. We computed a posterior distribution of predictions for the dependent variables (i.e., MHI or MABSHI) for a new observation (i.e., in our case H. sapiens) given the posterior distribution of the reduced models, the standard error for our species, the phylogenetic effects and the relevant covariates present in the reduced models. Predictive distributions that deviate strongly from the known value (i.e., outliers) provide evidence that the species has undergone a substantial amount of evolutionary change which cannot be accounted for by its phylogenetic position, branch lengths, and evolutionary change in the independent variable. The implication is that the trait has adaptive value for the species in ways not shared by its close relatives. This test was used to evaluate the idea that human handedness levels correspond to an evolutionary singularity (i.e., a derived character unique to us). This allowed us to assess the relative contribution of the covariates present in the reduced models when explaining H. sapiens handedness singularity. To predict hominin MHI and MABSHI values, we used the reduced human models to generate predictions for hominins, while marginalizing over measurement error and non-phylogenetic species effects. We assumed human-level standard errors for our meta-analytical component when computing our predictions. Data on hominin IMI was obtained from [81], whilst ECVs were obtained from [51]. Hominins were categorized as “group” in social system, whereas for substrate we classified Ardipithecus as “both,” whilst all the remaining hominin species were classified as “terrestrial.” We grafted the maximum credibility clade phylogeny from [51] to a consensus phylogeny obtained from the previously used 100 anthropoid trees. We removed the hominin species for which we did not have any data, which left us with the following species: Ar. ramidus, Au. afarensis, H. ergaster, H. erectus, H. floresiensis, and H. neanderthalensis. Acknowledgments We are grateful to Elodie Freymann for her work with Fig 1 in this article. [END] --- [1] Url: https://journals.plos.org/plosbiology/article?id=10.1371/journal.pbio.3003771 Published and (C) by PLOS One Content appears here under this condition or license: Creative Commons - Attribution BY 4.0. via Magical.Fish Gopher News Feeds: gopher://magical.fish/1/feeds/news/plosone/