Post by Admin on Feb 27, 2021 2:56:36 GMT
Phalangeal curvature
Phalangeal curvature reduces diaphyseal strain associated with suspension (32, 33) and is therefore an important correlate of primate positional behavior (19, 30, 36, 37). To clarify how phalangeal curvature in the Ar. ramidus hand compares to our comparative sample, we quantified the curvature of the PP3 in a large sample of anthropoids and fossil taxa using the included angle method (Fig. 3). The overall pattern of phalangeal curvature in the extant and fossil sample is consistent with that of previous studies (19, 36, 37). The phalanges of Hispanopithecus, Danuvius, and Ar. ramidus are characterized by pronounced curvature, falling within the ranges of variation of the most suspensory anthropoids (Pan, Pongo, hylobatids, and atelines; Fig. 3). The phalanges of the remaining Miocene hominoid taxa (Pierolapithecus, Griphopithecus, Ekembo, Sivapithecus, and Oreopithecus) as well as Au. afarensis and Au. sediba fall within the overlapping ranges of multiple anthropoid taxa including Pan, hylobatids, cercopithecins, and colobines.
Fig. 3 Proximal phalangeal curvature in anthropoid primates.
The gray bar indicates the range of variation among the Miocene fossil hominoids included here. The curvature of the Ar. ramidus PP3 falls within the ranges of variation of P. troglodytes, Pongo, hylobatids, and atelines; between the highly suspensory Miocene hominoids Danuvius and Hispanopithecus; and above the ranges of variation of all other taxa.
Evolutionary modeling
To test alternative hypotheses about the relationship between hand morphology and locomotor behavior, we reconstructed the adaptive landscape of primate hand evolution using phylogenetic comparative methods [(38, 39); see also (5, 22, 40)]. These methods allow us to translate adaptive hypotheses into explicit evolutionary models tested against comparative data in a maximum likelihood framework. They also allow us to estimate evolutionary parameters such as the optimal hand morphology for a given selective regime (i.e., locomotor category).
Our evolutionary modeling analyses use the first three PCs from our PCA to reduce the dimensionality of the dataset. In our a priori adaptive hypotheses (fig. S9), we assigned Australopithecus and Homo to their own selective regime and Ar. ramidus to the same regime as Pan given the results of our morphometric analyses. The first hypothesis (H1) recognizes taxa that are mostly terrestrial (Gorilla beringei, Papio, Theropithecus, and Mandrillus), those that are semiterrestrial (Pan, G. gorilla, Semnopithecus entellus, Macaca mulatta, Macaca nemestrina, Macaca nigra, Cercocebus, and Chlorocebus aethiops), those that are mostly arboreal (all other non-hominin taxa), and those who frequently use their hands in nonlocomotor contexts including Homo and Australopithecus. The second hypothesis (H2) is identical to H1 but separates mostly arboreal quadrupedal taxa from those that are more suspensory (Pongo, hylobatids, Ateles, and Brachyteles). The third hypothesis (H3) includes regimes for digitigrade taxa (M. mulatta, M. nigra, Papio, Theropithecus, Mandrillus, and Cercocebus), knuckle-walking taxa (Pan and Gorilla), suspensory taxa (Pongo, hylobatids, Ateles, and Brachyteles), nonlocomotor taxa (Homo and Australopithecus), and palmigrade taxa (all remaining anthropoids). We tested an alternative hypothesis (H4) based on H3 in which we placed Pan and Ardipithecus in a regime with the suspensory taxa and Gorilla in its own terrestrial knuckle-walking regime (fig. S9). The assignment of cercopithecoid taxa into locomotor or hand posture categories was based on published literature (41). Next, we used SURFACE to identify selective regimes based on morphological similarity without identifying them a priori (39), which has been used previously to test hypotheses of human and primate evolution (5, 22, 40). Last, we tested the relative support of our results with two simpler alternative evolutionary hypotheses including evolution by genetic drift (i.e., Brownian motion; HBM) and adaptation with a single optimum [i.e., Ornstein-Uhlenbeck (OU); HOU1]. Support for HBM or HOU1 would suggest that our data do not reflect an adaptive signal related to locomotor behavior.
Our results show that hand morphology is best described by 11 selective regimes with all a priori adaptive hypotheses performing worse than the SURFACE model (Fig. 4 and tables S6 and S7). The selective regimes identified by SURFACE likely reflect a combination of locomotion and other factors (Fig. 4 and table S7). For example, the best supported a priori model fit using the “ouch” package in R included a convergent suspensory regime containing atelines and Asian hominoids, but SURFACE separated them into distinct regimes that likely reflect their distinct evolutionary histories (e.g., Ateles and Brachyteles pollical reduction) and body sizes. Likewise, we found support for terrestrial knuckle-walking (Pan and Gorilla) and digitigrady (M. mulatta, M. nigra, Papio, Theropithecus, Mandrillus, and Cercocebus) regimes that were subsequently separated by SURFACE. Note that SURFACE identified three regimes representing the spectrum of nonsuspensory quadrupedal hand postures and substrate preferences among cercopithecoids and platyrrhines ranging from arboreal palmigrady to terrestrial digitigrady, mostly reflecting variation in metacarpal and phalangeal lengths relative to the hand geometric mean. Colobus was placed in its own regime nearest the arboreal palmigrady optimum reflecting evolution in the direction of the Ateles and Brachyteles regime associated with pollical reduction. SURFACE placed Ar. ramidus into a regime with Pan paniscus and P. troglodytes, suggesting adaptation to an optimal hand shape related to shared aspects of positional behavior, body size, and evolutionary history. SURFACE detected an evolutionary shift in hand morphology between the Homo-Pan ancestral phenotypic optimum, including Ardipithecus and Pan, and all later hominins primarily related to decreased intrinsic phalangeal length and increased MC1 length.
Fig. 4 The evolution of anthropoid hand shape according to SURFACE.
(A) Phylogenetic tree with branches painted according to selective regime. (B) Species means for PC1 and PC2 (small circles) and estimated phenotypic optima (large circles). (C) Species means for PC1 and PC3 (small circles) and estimated phenotypic optima (large circles). Ar. ramidus was placed in the same selective regime as P. troglodytes and P. paniscus. In contrast, all later hominins were placed in a selective regime with modern humans, which results in an evolutionary shift in hominin hand morphology from a Pan-like ancestor between ca. 4.4 and 3.2 Ma ago.
In addition to testing adaptive hypotheses, the OU model allows for the estimation of evolutionary parameters such as the phylogenetic half-life (t1/2; table S7), which is the average time it takes to evolve half the distance from the ancestral optimum to a new optimum given a regime shift. The phylogenetic half-life (t1/2) therefore quantifies the rate of adaptation. The estimate of t1/2 for PC1, which is the axis separating the ancestral hominin optimum from the primary hominin optimum, is 722 ka (thousand years) (table S7), which is relatively fast evolution compared to the total length of the tree [46.8 Ma (million years)]. As discussed above, PC1 is primarily associated with the intrinsic lengths of the metacarpals and phalanges. The estimates of t1/2 for PC2 and PC3 are 1.26 Ma and 5 ka, respectively. PC2 and PC3 are associated with metacarpal length, joint, and midshaft dimensions. The estimates of the phylogenetic half-life for each PC suggest that anthropoid metacarpal length and midshaft dimensions (i.e., MC1 relative to MC5) have evolved more slowly than phalangeal lengths and articular dimensions (e.g., the ulnar carpometacarpal and MCP joints). We used a Monte Carlo simulation method (42) to evaluate the statistical power to distinguish between alternative evolutionary models given our dataset and phylogeny (fig. S10). The simulation results show that all models can be distinguished from each other and support information criteria metrics [Akaike information criterion (AIC), small sample size corrected Akaike information criterion (AICc), and Schwarz information criterion (SIC)]. Previous studies have criticized the use of SURFACE for multivariate evolutionary modeling (43). We include the use of SURFACE for explicit comparison to recent similar studies (5, 22, 40) and as an additional perspective to our a priori models fit using the ouch package (38) in R. Overall, our use of Monte Carlo simulations provides additional support for the model fit by SURFACE.
We conducted a SURFACE analysis using the first three PC scores derived from our PCA on a modified dataset including H. laietanus (fig. S11 and table S8). SURFACE identified nine select regimes including a convergent regime containing H. laietanus, Pongo, Lagothrix, and Ateles. The locations of the evolutionary shifts estimated by SURFACE resemble our initial analysis using a larger dataset but with some differences. Ar. ramidus is placed in the same selective regime with P. troglodytes and P. paniscus, but expectedly, the analysis on the reduced dataset also places Au. afarensis in this regime since it excludes MC1 data. The Au. afarensis hand clusters with Au. sediba, H. naledi, and H. neanderthalensis in the PCA (fig. S5) but is distinguished from them along PC3 (fig. S11), which is associated with MCP and IP joint dimensions and MC4 length (table S2). Instead, SURFACE detects an evolutionary shift at the base of the Au. sediba and Homo clade. In addition, colobines, cercopithecins, and most papionins are placed in a single regime with evolutionary shifts toward a more terrestrial digitigrady regime in Papio and Theropithecus. Overall, the results of our SURFACE analysis on a modified dataset are consistent with our initial analyses because it recovers an evolutionary shift at the base of the great ape clade associated with suspensory adaptations.
Phalangeal curvature reduces diaphyseal strain associated with suspension (32, 33) and is therefore an important correlate of primate positional behavior (19, 30, 36, 37). To clarify how phalangeal curvature in the Ar. ramidus hand compares to our comparative sample, we quantified the curvature of the PP3 in a large sample of anthropoids and fossil taxa using the included angle method (Fig. 3). The overall pattern of phalangeal curvature in the extant and fossil sample is consistent with that of previous studies (19, 36, 37). The phalanges of Hispanopithecus, Danuvius, and Ar. ramidus are characterized by pronounced curvature, falling within the ranges of variation of the most suspensory anthropoids (Pan, Pongo, hylobatids, and atelines; Fig. 3). The phalanges of the remaining Miocene hominoid taxa (Pierolapithecus, Griphopithecus, Ekembo, Sivapithecus, and Oreopithecus) as well as Au. afarensis and Au. sediba fall within the overlapping ranges of multiple anthropoid taxa including Pan, hylobatids, cercopithecins, and colobines.
Fig. 3 Proximal phalangeal curvature in anthropoid primates.
The gray bar indicates the range of variation among the Miocene fossil hominoids included here. The curvature of the Ar. ramidus PP3 falls within the ranges of variation of P. troglodytes, Pongo, hylobatids, and atelines; between the highly suspensory Miocene hominoids Danuvius and Hispanopithecus; and above the ranges of variation of all other taxa.
Evolutionary modeling
To test alternative hypotheses about the relationship between hand morphology and locomotor behavior, we reconstructed the adaptive landscape of primate hand evolution using phylogenetic comparative methods [(38, 39); see also (5, 22, 40)]. These methods allow us to translate adaptive hypotheses into explicit evolutionary models tested against comparative data in a maximum likelihood framework. They also allow us to estimate evolutionary parameters such as the optimal hand morphology for a given selective regime (i.e., locomotor category).
Our evolutionary modeling analyses use the first three PCs from our PCA to reduce the dimensionality of the dataset. In our a priori adaptive hypotheses (fig. S9), we assigned Australopithecus and Homo to their own selective regime and Ar. ramidus to the same regime as Pan given the results of our morphometric analyses. The first hypothesis (H1) recognizes taxa that are mostly terrestrial (Gorilla beringei, Papio, Theropithecus, and Mandrillus), those that are semiterrestrial (Pan, G. gorilla, Semnopithecus entellus, Macaca mulatta, Macaca nemestrina, Macaca nigra, Cercocebus, and Chlorocebus aethiops), those that are mostly arboreal (all other non-hominin taxa), and those who frequently use their hands in nonlocomotor contexts including Homo and Australopithecus. The second hypothesis (H2) is identical to H1 but separates mostly arboreal quadrupedal taxa from those that are more suspensory (Pongo, hylobatids, Ateles, and Brachyteles). The third hypothesis (H3) includes regimes for digitigrade taxa (M. mulatta, M. nigra, Papio, Theropithecus, Mandrillus, and Cercocebus), knuckle-walking taxa (Pan and Gorilla), suspensory taxa (Pongo, hylobatids, Ateles, and Brachyteles), nonlocomotor taxa (Homo and Australopithecus), and palmigrade taxa (all remaining anthropoids). We tested an alternative hypothesis (H4) based on H3 in which we placed Pan and Ardipithecus in a regime with the suspensory taxa and Gorilla in its own terrestrial knuckle-walking regime (fig. S9). The assignment of cercopithecoid taxa into locomotor or hand posture categories was based on published literature (41). Next, we used SURFACE to identify selective regimes based on morphological similarity without identifying them a priori (39), which has been used previously to test hypotheses of human and primate evolution (5, 22, 40). Last, we tested the relative support of our results with two simpler alternative evolutionary hypotheses including evolution by genetic drift (i.e., Brownian motion; HBM) and adaptation with a single optimum [i.e., Ornstein-Uhlenbeck (OU); HOU1]. Support for HBM or HOU1 would suggest that our data do not reflect an adaptive signal related to locomotor behavior.
Our results show that hand morphology is best described by 11 selective regimes with all a priori adaptive hypotheses performing worse than the SURFACE model (Fig. 4 and tables S6 and S7). The selective regimes identified by SURFACE likely reflect a combination of locomotion and other factors (Fig. 4 and table S7). For example, the best supported a priori model fit using the “ouch” package in R included a convergent suspensory regime containing atelines and Asian hominoids, but SURFACE separated them into distinct regimes that likely reflect their distinct evolutionary histories (e.g., Ateles and Brachyteles pollical reduction) and body sizes. Likewise, we found support for terrestrial knuckle-walking (Pan and Gorilla) and digitigrady (M. mulatta, M. nigra, Papio, Theropithecus, Mandrillus, and Cercocebus) regimes that were subsequently separated by SURFACE. Note that SURFACE identified three regimes representing the spectrum of nonsuspensory quadrupedal hand postures and substrate preferences among cercopithecoids and platyrrhines ranging from arboreal palmigrady to terrestrial digitigrady, mostly reflecting variation in metacarpal and phalangeal lengths relative to the hand geometric mean. Colobus was placed in its own regime nearest the arboreal palmigrady optimum reflecting evolution in the direction of the Ateles and Brachyteles regime associated with pollical reduction. SURFACE placed Ar. ramidus into a regime with Pan paniscus and P. troglodytes, suggesting adaptation to an optimal hand shape related to shared aspects of positional behavior, body size, and evolutionary history. SURFACE detected an evolutionary shift in hand morphology between the Homo-Pan ancestral phenotypic optimum, including Ardipithecus and Pan, and all later hominins primarily related to decreased intrinsic phalangeal length and increased MC1 length.
Fig. 4 The evolution of anthropoid hand shape according to SURFACE.
(A) Phylogenetic tree with branches painted according to selective regime. (B) Species means for PC1 and PC2 (small circles) and estimated phenotypic optima (large circles). (C) Species means for PC1 and PC3 (small circles) and estimated phenotypic optima (large circles). Ar. ramidus was placed in the same selective regime as P. troglodytes and P. paniscus. In contrast, all later hominins were placed in a selective regime with modern humans, which results in an evolutionary shift in hominin hand morphology from a Pan-like ancestor between ca. 4.4 and 3.2 Ma ago.
In addition to testing adaptive hypotheses, the OU model allows for the estimation of evolutionary parameters such as the phylogenetic half-life (t1/2; table S7), which is the average time it takes to evolve half the distance from the ancestral optimum to a new optimum given a regime shift. The phylogenetic half-life (t1/2) therefore quantifies the rate of adaptation. The estimate of t1/2 for PC1, which is the axis separating the ancestral hominin optimum from the primary hominin optimum, is 722 ka (thousand years) (table S7), which is relatively fast evolution compared to the total length of the tree [46.8 Ma (million years)]. As discussed above, PC1 is primarily associated with the intrinsic lengths of the metacarpals and phalanges. The estimates of t1/2 for PC2 and PC3 are 1.26 Ma and 5 ka, respectively. PC2 and PC3 are associated with metacarpal length, joint, and midshaft dimensions. The estimates of the phylogenetic half-life for each PC suggest that anthropoid metacarpal length and midshaft dimensions (i.e., MC1 relative to MC5) have evolved more slowly than phalangeal lengths and articular dimensions (e.g., the ulnar carpometacarpal and MCP joints). We used a Monte Carlo simulation method (42) to evaluate the statistical power to distinguish between alternative evolutionary models given our dataset and phylogeny (fig. S10). The simulation results show that all models can be distinguished from each other and support information criteria metrics [Akaike information criterion (AIC), small sample size corrected Akaike information criterion (AICc), and Schwarz information criterion (SIC)]. Previous studies have criticized the use of SURFACE for multivariate evolutionary modeling (43). We include the use of SURFACE for explicit comparison to recent similar studies (5, 22, 40) and as an additional perspective to our a priori models fit using the ouch package (38) in R. Overall, our use of Monte Carlo simulations provides additional support for the model fit by SURFACE.
We conducted a SURFACE analysis using the first three PC scores derived from our PCA on a modified dataset including H. laietanus (fig. S11 and table S8). SURFACE identified nine select regimes including a convergent regime containing H. laietanus, Pongo, Lagothrix, and Ateles. The locations of the evolutionary shifts estimated by SURFACE resemble our initial analysis using a larger dataset but with some differences. Ar. ramidus is placed in the same selective regime with P. troglodytes and P. paniscus, but expectedly, the analysis on the reduced dataset also places Au. afarensis in this regime since it excludes MC1 data. The Au. afarensis hand clusters with Au. sediba, H. naledi, and H. neanderthalensis in the PCA (fig. S5) but is distinguished from them along PC3 (fig. S11), which is associated with MCP and IP joint dimensions and MC4 length (table S2). Instead, SURFACE detects an evolutionary shift at the base of the Au. sediba and Homo clade. In addition, colobines, cercopithecins, and most papionins are placed in a single regime with evolutionary shifts toward a more terrestrial digitigrady regime in Papio and Theropithecus. Overall, the results of our SURFACE analysis on a modified dataset are consistent with our initial analyses because it recovers an evolutionary shift at the base of the great ape clade associated with suspensory adaptations.