- Research article
- Open Access
Vertical clinging and leaping induced evolutionary rate shifts in postcranial evolution of tamarins and marmosets (Primates, Callitrichidae)
BMC Ecology and Evolution volume 21, Article number: 132 (2021)
Callitrichids comprise a diverse group of platyrrhine monkeys that are present across South and Central America. Their secondarily evolved small size and pointed claws allow them to cling to vertical trunks of a large diameter. Within callitrichids, lineages with a high affinity for vertical supports often engage in trunk-to-trunk leaping. This vertical clinging and leaping (VCL) differs from horizontal leaping (HL) in terms of the functional demands imposed on the musculoskeletal system, all the more so as HL often occurs on small compliant terminal branches. We used quantified shape descriptors (3D geometric morphometrics) and phylogenetically-informed analyses to investigate the evolution of the shape and size of the humerus and femur, and how this variation reflects locomotor behavior within Callitrichidae.
The humerus of VCL-associated species has a narrower trochlea compared with HL species. It is hypothesized that this contributes to greater elbow mobility. The wider trochlea in HL species appears to correspondingly provide greater stability to the elbow joint. The femur in VCL species has a smaller head and laterally-oriented distal condyles, possibly to reduce stresses during clinging. Similarly, the expanded lesser trochanters visible in VCL species provide a greater lever for the leg retractors and are thus also interpreted as an adaptation to clinging. Evolutionary rate shifts to faster shape and size changes of humerus and femur occurred in the Leontocebus clade when a shift to slower rates occurred in the Saguinus clade.
Based on the study of evolutionary rate shifts, the transition to VCL behavior within callitrichids (specifically the Leontocebus clade) appears to have been an opportunity for radiation, rather than a specialization that imposed constraints on morphological diversity. The study of the evolution of callitrichids suffers from a lack of comparative analyses of limb mechanics during trunk-to-trunk leaping, and future work in this direction would be of great interest.
Facilitated by the collections in natural history museums from around the globe, the comparative analysis of mammalian limb long bones is an important means of gaining insight into phenotypic evolutionary disparity and diversification (e.g., [1,2,3,4,5,6,7,8,9,10,11]). Bone shape reflects adaptations to specific functional demands related to taxon ecology, as well as being constrained by body size (allometric changes) and phylogenetic history [7, 12,13,14]. The study of functional shape in respect to organismal evolution has been facilitated by the advent of three-dimensional (3D) imaging techniques and improvements and availability of software for multivariate statistics, including being able to account for the phylogenetic relatedness of analyzed specimens [13, 15,16,17,18].
Although the primate postcranial skeleton is relatively conservative in terms of the number, and the organization, of bones (e.g. in comparison to bats, cetaceans, or perissodactyls), primates have nevertheless evolved a wide variety of locomotor behaviors [19,20,21,22]. The locomotor repertoire can also be relatively diverse even within individual primate species . Tamarins and marmosets (Callitrichidae; Primates), a species-rich clade of New World primates, display striking differences in locomotor and feeding ecology despite a lack of notable morphological differences. Comparative analysis of 3D limb long bone shape has the potential to detect subtle phenotypic differences between species and is employed here to gain insight into the ecologically significant disparity and evolution within Callitrichidae.
The Callitrichidae, (also sometimes classified as sub-family Callitrichinae [24,25,26]) constitute a group of over 60 species and subspecies  living across South and Central America . All species are of a relatively small size, with head and body lengths ranging from just 130 to 370 mm  and body masses of 105 to 700 g . Callitrichids also possess secondarily pointed claws (tegulae) instead of flat nails on all digits except the hallux, and exhibit lateral sequence gaits instead of the diagonal sequences usually used by quadrupedal primates. Due to these combined characteristics, members of the Callitrichidae have been proposed as viable models for certain hypothetical stages of early primate evolution (see , for a review). Diurnality and arboreality are also collective traits across all callitrichid species .
Locomotor behavior of Callitrichidae
In contrast to the shared ecological traits described above, locomotor behavior differs greatly between callitrichid taxa [23, 30, 31]. Some species use vertical clinging and leaping (VCL) and make extensive use of vertical supports, employing trunk-to-trunk leaping to travel, forage or escape predators. VCL species, such as the spring tamarin Callimico goeldii, are often found in forest layers close to the ground . Other callitrichids preferentially use horizontal leaping (HL), and leap from and usually land on thin, flexible, terminal branches in the canopy, such as lion tamarins of the genus Leontopithecus . Across Callitrichidae a continuous spectrum exists, ranging from either VCL or HL specialists to various intermediate degrees of tropism (Table 1 and references therein). To a certain extent this locomotor diversity is linked to diversity in feeding behavior [24, 34, 35]. VCL is often used to facilitate various types of foraging on large tree trunks, including gummivory, and foraging for fungi, and insects in bark, whereas terminal branch locomotion and leaping are associated with fruit and nectar exploitation, as well as feeding on insects in the canopy, or small vertebrates [24, 36,37,38]. VCL is common in the subclade formed by Callimico, Cebuella, Mico, and Callithrix, and is also exhibited by some tamarin species [24, 39,40,41]. A predominant use of VCL is therefore not restricted to a single lineage, so it therefore is possible to study two independent adaptive specializations toward VCL (Fig. 1). Moreover, preference in the orientation and diameter of the support used while leaping appears to be decoupled from the environment, as illustrated by the behavior of mixed-species groups comprised of species with different specializations. The use of the substrate is therefore apparently a matter of species’ preference and differential advantage, rather than substrate availability [39, 42,43,44,45].
Morphological correlates to locomotor behaviour
Specific morphological adaptation to VCL has been extensively studied in Strepsirrhini [74,75,76], and callitrichid VCL kinematics, mechanics, and ecology have also been studied [24, 33, 75, 77,78,79]. Still, morphological correlates of adaptations to this type of locomotion in Callitrichidae have received limited attention, and previous work was focused on functional interpretations of fossils or poorly known species [80, 81], or on hand and feet morphology [42, 68]. The intermembral index (IMI) obtained from the ratio of forelimb length (humerus + radius): hindlimb length (femur + tibia)  has been suggested as a marker of specialization to VCL [39, 76, 82], and it has been demonstrated that a proportionally longer hind-limb (i.e. a low IMI) is linked to greater take-off speed . Garber and Leigh  showed that Callimico goeldii has extremely elongated hind-limbs, a characteristic also observed in Strepsirrhini. The authors also described a long patellar groove in Leontocebus fuscicollis and interpreted it as an adaptation to greater knee excursion during vertical clinging and leaping. Nevertheless, Garber and Leigh  measured a high IMI in Leontocebus fuscicollis, mostly because of the elongated forearm, which increases when the hand is included. They discussed it also as an adaptation to trunk-to-trunk leaping, the relative long forelimb providing a greater distance when flexing the limb to dissipate the landing impact, and the mass distribution farther from the bodyprovides an advantage during in-air rotation. Similarly, Falsetti and Cole  suggested that the greater IMI in Leontocebus fuscicollis compared to Callithrix jacchus and Saguinus oedipus might increase agility on trunks, with the greater forelimb span facilitating clinging and deceleration during landing. But, as detailed above, a lower IMI has also been suggested to be beneficial for leaping (a relatively longer hindlimb allowing for a longer distance during propulsion) indicating that the intermembral index is limited as an explanatory variable for differences in locomotor evolution among Callitrichidae. There may be different evolutionary pathways behind the independently acquired VCL behavior in Leontocebus and in the Callithrix/Callimico/Mico/Cebuella clade, potentially reflected in subtle differences of limb long bone shape between both taxa.
Bicca-Marques  showed that Leontopithecus have long and slender hands, and described a similar pattern in Leontocebus fuscicollis, Leontocebus tripartitus, and Leontocebus nigricollis. However these characteristics are discussed in respect to feeding strategies, not locomotion. Smith and Smith  showed a link between foot length, tegulae curvature, and the use of vertical supports. They also demonstrated the increasing number of plantar and palmar ridges in relation to the degree of exudativory.
Differences in functional demands posed by VCL and HL
VCL differs from HL in the mechanical demands faced by the animal, especially in respect to the limb long bones. Generation of leaping forces during trunk-to-trunk leaping mostly involves the hindlimbs [77, 78]. The flight phase during VCL also involves a change in body orientation to face the landing support, which is not necessary during HL. Landing on non-compliant supports implies high compressive forces as all the landing energy is actively dissipated during leg flexion. During HL, the majority of the landing supports are highly compliant, making it necessary to grasp and balance rather than to dissipate impact forces. The evolutionary role of primate proximal limb long bones (i.e., the humerus and femur) in locomotor behavior is considered to be very important [83,84,85,86,87]. All mechanical loads transmitted between the limb and the tree trunk support pass via these bones. The bone morphology can therefore be expected to reflect different functional demands related to locomotor behavior, and thus provide insight into how morphology relates to ecology and may reflect the evolution of the entire family.
The study of evolutionary rates (both molecular and phenotypic) is key for the understanding of macroevolution, and is the subject of much debate and methodological development [17, 88,89,90,91,92,93]. The rhythm of evolutionary change is a useful indicator of either the acquisition of a new fitness optimum, or of a strong stabilizing selection acting on the system. The adaptive radiation model  predicts a higher evolutionary rate in response to the emergence of a new function, or during the selection process in response to a new or changing environment. In the case of Callitrichidae, the development of VCL locomotion is hypothesized to coincide with an increase in the evolutionary rate of the locomotor structures involved.
Morphological differences for femur and humerus between VCL and HL species are investigated. Significant differences in shape are expected (e.g. reinforced hindlimb extensors insertion in VCL), and the size difference between the groups is also considered in order to identify any deviations from isometric shape change with increasing body size (allometry). The differences themselves are described and their functional significance discussed. Finally, as the exploitation of a new ecological niche is expected to provide the opportunity for rapid diversification, changes in the evolutionary rates across Callitrichidae are investigated in order to determine if VCL has been a driver of increased morphological diversification. All of these questions are analyzed while accounting for body size and accounting for the statistical non-independence of taxa due to their phylogenetic relatedness using multivariate statistical and phylogenetic comparative methods.
We observed a significant phylogenetic signal for centroid size for both the humerus (K = 0.832, p-value = 1e−04) and the femur (K = 0.680, p-value = 0.0022).
Differences between ecological groups
For both the humerus and femur, VCL species have a smaller centroid size compared to HL species (t-test on log-transformed centroid size: p-value = 0.0075 and 0.0128 for the humerus and femur, respectively). However there is an overlap in the size distribution of the two groups (Fig. 2) and when phylogeny is taken into account the differences are no longer significant (phylogenetic ANOVA on 1000 simulations for humerus: p-value = 0.251, for femur: p-value = 0.298).
Per taxa centroid size difference
The tamarins (Saguinus and Leontopithecus) are distinct from the other species in our sample in that they have a relatively homogenous centroid size for both humerus and femur (Fig. 3A, B). Leontopithecus specimens present large humerus and femur sizes, whereas Cebuella present small values, and Callithrix species lie between these two groups. Callimico display a medium value for the humerus, but a high value for the femur. Saguinus midas and S. niger have high femoral centroid sizes. In terms of the ratio of the humeral centroid size over the femoral one (Fig. 3C), Callimico exhibits a larger hind limb proximal segment compared to the forelimb (ratio = 0.929). To a lesser extent, this is also the case in Callithrix species (ratios between 0.933 and 0.940) and Saguinus leucopus and S. labiatus (both 0.935). Leontopithecus specimens present the highest ratios (0.956–0.957), followed by Cebuella (0.951). Leontocebus (0.940–0.946), Mico (0.939–0.944) and all other Saguinus (0.935–0.942) display medium values.
Evolutionary rates of centroid size
Evolutionary rates for humeral and femoral size are significantly higher for the inner nodes of the Leontocebus clade, and all nodes of the clade show significantly higher rates for the size ratio (Fig. 3A–C). The evolutionary rates for the centroid size of the femur are significantly higher for the node that separates Leontopithecus from marmosets and Callimico, and for the node that separates Cebuella pygmaea from Mico species. There is a significantly lower evolutionary rate for humeral size and femoral size for the basal node of Saguinus, and for the node that separates the Saguinus midas and the Saguinus goeffroyi clades. Femoral size is also significantly lower at the node separating Saguinus leucopus from the S.geoffroyi/S.oedipus clade.
We did not observe a significant phylogenetic signal for shape for either the humerus or the femur (Kmult of 0.080 and 0.087, p-value of 0.67 and 0.64 for the humerus and femur respectively). No significant deviation from isometric shape in respect to centroid size was found in the humerus, however, the p-value is close to significance. There is significant allometry in femoral shape (Table 2). The locomotor categories have a significant impact on both the humeral and femoral shape. The interaction term that accounts for both size and locomotor mode is significant for the humerus, indicating that the allometric slopes in the humerus differ significantly between the two locomotor groups.
The first axis (PC1) of the phylomorphospace computed for humeral shape (Fig. 4A) reflects 41% of total variance. Leontopithecus plots on the negative end of PC1, whereas Callithrix kuhllii and S. geoffroyi fall on the extreme positive end. Species of the genera Leontocebus and Mico plot close to the negative end, as well as Saguinus mystax. PC2 (16% of total variance) separates Saguinus bicolor on the extreme positive side and Leontocebus melanoleucus on the extreme negative. There is no obvious separation between the two locomotor categories along either axis. PC1 on humeral shape separates species with robust bones, a caudally bent diaphysis, and proportionally large epiphyses on the positive side from species with relatively more gracile bones with smaller epiphyses. PC2 on humeral shape separates the species with a diaphysis bent posteriorly around the deltoid tuberosity on the negative side from species with a more straight shaft on the positive side. For the humerus, the between-group PCA separates the locomotor groups with an accuracy of 84%. The p-value of the group differences obtained from permutation testing is 0.3812. The observed differences in shape are a straighter humeral neck in VCL species whereas HL species have a humeral head more caudally curled (Fig. 5A). The medial lip of the trochlea is expanding more medially in HL resulting in a broader cranial area of the trochlea, and the medial supratrochlear area is more expanded. The lateral supracondylar crest has a broader proximal portion in HL species compared to those that use VCL.
The phylomorphospace computed for femoral shape (Fig. 4B) displays along PC1 (34% of total variance) Callithrix kuhllii on the extreme negative end from most of the other species. Still, species of the genus Leontopithecus, Cebuella and Saguinus bicolor have lower values on PC1 than all other species, except for Callithrix kuhllii. PC2 (11% of total variation) depicts Leontopithecus rosalia on the extreme positive end, and Leontocebus tripartitus and Leontocebus weddelli on the extreme negative end. PC2 appears to tentatively separate the two locomotor categories, with VCL species trending towards lower values, whereas HL species trend towards higher values along this PC axis.
PC1 of femoral shape separates species with proportionally larger epiphyses on the negative end from species with relatively smaller ones on the positive end. PC2 of femoral shape separates species with a proportionally larger femoral head and a smaller lesser trochanter on the positive end of the axis. Please refer to Additional file 1 (Fig. S3) for graphic representation and a more detailed description.
The between-group PCA on the femur (Fig. 5B) separates the groups with an accuracy of 88%, but the p-value of the between-group difference remains non-significant (0.0977). The VCL group differs in shape from the HL by its proportionally smaller femoral head, larger lesser trochanter, and narrower distal part of the medial condyle, with the plane tangent to the most distal point of two condyles tilted toward the medial side.
Evolutionary rates of shape changes
When testing for shifts in evolutionary rates (Fig. 6) we observe significantly higher rates for the evolution of the shape of both femur and humerus in the Leontocebus clade. The node separating Callithrix kuhllii from Callithrix penicillata and Callithrix jacchus has a significantly higher evolutionary rate for femoral shape change. Both the crownward and stemward nodes relative to the aforementioned node present high evolutionary rates but are not significant (Additional file 1: Table S4). There is a significantly lower evolutionary rate for femoral shape change for the basal node of Saguinus. Values for the other nodes are reported in Additional file 1 (Additional file 1: Table S4 and Fig. S4).
Evolution of size in callitrichid proximal limb long bones
To gain insight into the evolution of callitrichid primates and specifically to investigate whether differences in locomotor ecology are reflected in proximal limb bone morphology, we compared size and shape within a comparative sample of humeri and femora. Within this clade of primates, VCL has independently evolved twice, under the hypothesis of HL being the ancestral state (see Garber  for a synthesis).
We observed that VCL species present on average a smaller centroid size for both the humerus and femur compared to HL species. It was noticed in previous work  that among callitrichids, VCL species are generally smaller in size than their HL counterparts. Still, as noticed by Garber , there is an overlap in the size range of the two ecological categories. In our dataset, the highly specialized VCL Callimico goeldii has the largest centroid size for the femur and the fourth largest centroid size for the humerus, behind three HL species. The relatively large body mass of this species was previously noted (male mean 499 gr, female mean 468 gr; ) compared to other VCL species like Callithrix jacchus (male mean 317 gr, female mean 324 gr; ).
The evolution of centroid size is highly influenced by phylogenetic inertia. When mapped across the phylogeny (Fig. 3A), humeral centroid size shows limited variation, except among the Leontopithecus, Callimico, and Cebuella clades. The variations of femoral size are greater across the phylogeny, and accordingly the ratio of humeral to femoral centroid size also varies over the phylogeny.
Intermembral Index (IMI) is often seen as a better predictor for locomotion in primates than overall size [74, 95, 96]. As a general trend among primates, highly specialized VCL primates tend to present a low IMI compared to species with quadrupedal habits, the brachiating species scoring higher [23, 97]. Additionally, IMI tends to scale positively with body size, thus the combined effect of body size and locomotor habits should be considered . Despite this general trend among primates, Jungers’ results  suggest that IMI is independent of body mass within callitrichids. This variable IMI scaling pattern in callitrichids may cause comparisons of IMI with other primate groups to be misleading. Additionally, as element proportions within each limb may vary, comparing a humerus to femur size ratio to the IMI would require extra caution.
Much like IMI, the humerus to femur size ratio is low for Callimico goeldii . We also calculate a low size ratio within the Callithrix clade. A high IMI value in Cebuella was noted by Garbino , in agreement with the size ratio that we observed. It has been suggested that both high and low IMI might be associated with VCL within Callitrichidae, indicating a divergent adaptive path . The elongated hindlimb observed in Callimico may provide a longer propulsive phase, allowing for greater speed at take off and thus increasing the potential for larger distances to be crossed while leaping, similar to the adaptation described for strepsirrhines specialized in VCL [39, 75, 98]. Alternatively, it is hypothesized that the elongated forelimb observed in Leontocebus allows for a longer deceleration phase, thus dampening the peak force during landing . We have observed a relatively short humerus compared to the femur in the genus Callithrix, suggesting an adaptive response to VCL similar to the one observed in Callimico, and, to a lesser extent, also similar to what was deduced from quantifying IMI [65, 78]. Nevertheless, the analysis of trunk-to-trunk leaping kinematics  determined that Callithrix jacchus uses a leaping pattern similar to the one used by Cebuella: both present low angle and low velocity at take-off, thus losing more height during flight when compared to a jump of a similar distance by Callimico . Thus Callithrix jacchus and Cebuella are expected to rely on multiple shorter jumps, instead of covering a long distance in a single trunk-to-trunk leap [33, 78, 98].
We found that size ratio in the whole Leontocebus genus is greater than in Callimico and Callithrix, and lower than in Cebuella. Garber  measured VCL habits (preferred jump distance) in natura that place Leontocebus weddelli in between Callimico and Cebuella.
The species of Mico included in this study have a similar ratio to Leontocebus. Moreover, the humeral shape tends to be similar suggesting an ecological similarity between Mico and Leontocebus (Fig. 4A). Correspondingly, the distribution area of the two genera is generally separate, as would happen if the two were competing for the same resources. It has been previously suggested that a competitive interaction might exist between the small-sized tamarin (Leontocebus) and Mico . Nevertheless, overlaps in distribution exist in some areas [26, 99]. In these areas, they form mixed-species troops (Mico emiliae and Leontocebus weddelli, ; Leontocebus weddelli, Mico melanurus, and Mico rondoni; ). Thus the differences in distribution may be a result of evolutionary history rather than mutual exclusion.
Evolution of shape in callitrichid proximal limb long bones
The first two PCA axes for both the humerus and the femur are similarly structured: the Leontopithecus species are at an extreme, the HL tamarins occupy a central area, and all VCL species are spread along both axes to cover a large area. The case of Callithrix kuhllii must be considered with caution as it differs strongly from the other species of the genus Callithrix. This species is represented in our sample by a single specimen of unknown origin. The morphological separation of the genus Leontopithecus is well documented: greater body weight compared to other Callitrichidae , relative long hands  and derived incisors  that are linked with its derived food-seeking strategies and peculiar manipulative abilities [35, 47, 58, 59]. The slender humerus and femur, straight humeral shaft, and the small lesser trochanter of the femur are congruent with a reduced demand on shoulder flexors and hip adductors when compared to VCL species. The narrow and shallow trochlea reduces passive elbow stability in favor of mobility, requiring a more active control of stability [102, 103].
The shape of the humeral head differs between the VCL and HL groups (Fig. 5A). The articular surface of the humeral head expands more in the disto-caudal direction in HL species, which may allow greater retraction of the humerus. In contrast, the more proximally oriented humeral head of VCL species (i.e., the orientation of the humeral head is more in line with the humeral shaft) places the greater tubercle in a relatively lower position, apparently allowing a greater degree of freedom of the humero-scapular articulation during shoulder extension. This has also been described in VCL strepsirrhine species [87, 104], although a study of joint kinematics during these two locomotor behaviors is currently lacking for Callitrichidae. The distal articular surfaces of the humerus also differ between the groups. The more expanded trochlea, particularly in its medial margin, observed in HL might increase elbow joint stability during flexion–extension while the smaller medial margin of the trochlea and the broader capitulum observed in VCL species may be assumed to result in greater pronation/supination abilities. Likewise, Szalay and Dagosto  highlighted the reduction in trochlea breadth observed in VCL strepsirrhines compared to other locomotor specialists. These authors also described the trochlea as expanding further towards the humeral shaft in VCL species. This is not the case in Callitrichidae, likely due to the difference in support substrate diameter preferences, and posture during clinging [24, 39]. Strepsirrhine species prefer small diameter vertical supports for clinging (although clawed strepsirrhines can cling onto large-sized trunks [75, 105]). Most Callitrichidae specialize in the use of large vertical supports (tree trunks), thanks to their claws. The larger the support, the greater the elbow angle . Thus, it is unlikely that the vertical clinging locomotor behavior of callitrichids would induce the same need for an expanded trochlea in the antero-cranial direction as observed in Strepsirrhini. Again, the lack of knowledge of joint kinematics during locomotion in Callitrichidae limits conclusions about this characteristic.
The lesser trochanter is expanded in the VCL group (Fig. 5B). It provides the insertion for the iliopsoas, and this leg retractor plays a central role during clinging and has been shown to be enlarged in VCL specialists like strepsirrhine primates [74, 80]. Still, it may be suggested that an expanded lesser trochanter might provide a wider area of insertion for muscles that insert in direct proximity, such as the vastus medialis and vastus intermedius knee extensors that could be of adaptive value during leaping . The actual predominance of one of these three muscles in shaping the lesser trochanter remains to be determined.
The femoral head is smaller in VCL species when compared to HL species. This reduction of the femoral head was also noted for Cebuella compared to other marmosets by Ford and Davis . These authors suggest that it may be an adaptation to the constantly flexed hindlimb posture used during clinging.
The shape of the distal part of the medial condyle and the overall orientation of the tibiofemoral articulation during extension differ between the two groups. The orientation of the femur condyle observed in the VCL group induces a more everted position of the knee. Thus the knee articulation stays horizontal and the hip is less abducted, while the feet remain at a greater distance from each other. This is congruent with greater clinging abilities on large-diameter substrates like vertical trunks, but actual articular angles and associated forces remain to be investigated.
While VCL Strepsirrhine primates are often used to discuss the morphological changes observed in Callitrichidae, this article included, caution should be used when doing so. The fact that IMI in Callitrichidae does not follow the scaling pattern observed among primates, in addition to the fact that, unlike Strepsirrhini, Callitrichidae land hand first, may result in marked functional differences. In this context, Gebo  suggested that vertical clinging and leaping behavior should actually be seen as three distinct patterns: first that of tarsier and Galago, second the callitrichids and other keeled-nailed primates, and third the indriids and Lepilemur.
Rate shifts in callitrichid limb evolution
The rates of evolution vary significantly across the phylogeny and especially among tamarins. The rates observed in the Saguinus clade are mostly low and are significantly lower than the tree average at the root of the Saguinus clade for humerus size, and for the size and shape of the femur. The clade that separates the midas group from the oedipus group is also significantly lower for both humeral and femoral size. These low evolutionary rates appear in a group that is almost exclusively composed of HL species. However, feeding diversity is dispersed throughout the phylogeny of the Saguinus clade with some species preferentially eating the reproductive parts of plants (> 70% of fruits and flowers; S. oedipus ; S. niger ; S. labiatus ; S. leucopus [70, 108]), and others consuming a greater amount of animal prey (> 30%; S. bicolor ; S. midas ; S. geoffroyi ; S. mystax ; S. imperator ). Thus, despite the relative shape homogeneity across the genus and the low evolutionary rates observed, the group exhibits considerable ecological diversity. The large tamarin morphotype (Saguinus) appears therefore to be a generalist morphology, allowing diversification through other means than locomotor niches, for example through food resources, cranial or dental morphology, or behavior. The role of this ecological segregation between closely related species in the diversification of the clade requires further investigation.
The Leontocebus clade is marked by increased evolutionary rates in all variables, and is significant for humeral size and shape, femoral shape, and the ratio between humerus and femur size (see Additional file 1: Table S4 and Fig. S4). The apparition of VCL behavior in this genus can thus be seen as a key innovation that resulted in an opportunity to exploit new resources. The genus Leontocebus shows differences in food sources between species. Leontocebus weddelli can be highly frugivore (up to 84% in the dry season ) whereas other species can be highly carnivore, such as Leontocebus fuscicollis (up to 48%; ). Moreover, the diet can vary between populations of the same species, or even within the same group, depending on the season [34, 60, 64, 111]. In opposition to what has been previously observed during specialization , we do not observe low evolutionary rates after the apparition of VCL abilities. The addition of a new food resource to the species repertoire (here vertical trunk insects and gum) might be the reason for the greater evolutionary rates observed in this group by broadening its ecological niche. The fact that this genus can form mixed groups with other Callitrichidae species tends to confirm this access to a new feeding niche [39, 44, 45, 99, 113].
We found high evolutionary rates for humeral and femoral shape changes within Callithrix, reaching significance only for the node separating Callithrix kuhllii from the C. penicillata/C. jacchus clade. The C. penicillata/C. jacchus clade is the group where one of the highest degrees of tree gouging specialization has been described and is considered to be an adaptation to the highly seasonal food resources available in the dry forests located in north-eastearn and central Brazil . In this environment, fruit is not available year-round, thus the development of a high degree of gummivory was a key adaption . The strong adaptive constraints induced by this new environment may have resulted in rapid morphological changes within the C. penicillata/C. jacchus clade. Still, as mentioned previously, only one specimen was obtained for C. kuhllii and it appears to diverge from the rest of the sample. Thus the high evolutionary rates at the node separating Callithrix kuhllii from the C. penicillata/C. jacchus clade may be due to an artifact driven by a single specimen.
Quantified shape descriptors were used to investigate the shape evolution of the humerus and the femur associated with VCL behavior, which has appeared on two separate occasions within the evolution of the Callitrichidae.
Besides the phylogenetic inertia, VCL callitrichids share anatomical similarities that are functionally relevant. The distal articular surfaces of the humerus differ between HL and VCL species, similar to what was noted for other groups previously. The narrower trochlea may contribute to greater mobility in VCL species while the wider one visible in HL species appears to provide greater elbow stability. The femur differs in VCL species by exhibiting a relatively small femoral head, an expanded lesser trochanter, and the distal portion of the condyles are more laterally oriented. Functional interpretation suggests that these differences reflect the increased role of clinging in these species. Leaping probably involves the highest peak forces the musculoskeletal system encounters. However, these species also spend a lot of time clinging. Thus adaptations to clinging potentially yield large energetic savings and unsurprisingly seem to contribute to the evolution of limb long bone shape.
The lack of available comparative analyses of limb mechanics during trunk-to-trunk leaping in Callitrichidae limits our interpretation. Moreover, we are not aware of any published analyses quantifying bone strains during clinging behavior. The understanding of the evolution of the Callitrichidae would greatly benefit from such comparative experimental studies. More work of this kind has been conducted in Strepsirrhine primates but the marked differences in clinging and landing behaviors limit analogies between the two groups. Future work would greatly benefit from more quantitative data on locomotion to compliment research on locomotor group partition, which is arbitrary by definition. Only homogeneously quantified information would allow detailed morphological correlates to be identified. Moreover, information on the reasons for locomotion (moving toward a fruit feeding site, pursuing prey, escaping a predator) should be considered. The locomotor choices might differ during different circumstances. Moving between feeding sites might likely involve lower speed movements with a low falling risk to minimize the energetic cost, in comparison to escaping a predator, where an “all-out” strategy might be preferred.
This study has demonstrated significant differences in evolutionary rates within the Callitrichidae. The transition to VCL behavior within callitrichids (specifically the Leontocebus clade) appears to have been an opportunity for radiation, as opposed to a specialization that limits morphological diversity. In contrast, the genus Saguinus has primarily maintained a preference for HL and shows low evolutionary rates of change in terms of size and morphology of the humerus and femur, respectively, despite the diversity in feeding behavior. This suggests that this genus has reached a local fitness optimum of the locomotor apparatus.
The 58 specimens analysed in this study were obtained from the Field Museum of Chicago, the American Museum of Natural History of New York, and the Museum d’Histoire Naturelle of Paris. To get a minimal estimate of intra-specific variability, up to three specimens per species were sampled, when available, from the collections. Microfocus computed tomography (µCT) scans were realized at the Prof. Luo Lab at the University of Chicago, at the Shared Materials Instrumentation Facility at Duke University, and the μCT facilities of the Montpellier Resources Imagerie at Université de Montpellier, respectively. The resolutions used vary between 12.6 and 18 µm depending on specimen size and scanning facility. Details on specimens, scanning resolution, scanning apparatus and operators are available in Additional file 2.
Scan processing and landmark digitization
The μCT scans were segmented and surface models were extracted using Amira software (Thermo Fisher Scientific, version 6.0.0 ). The inner structure of the surface models was removed using Meshlab . Final tidying and decimation of the surface models to 500,000 triangles were made under Geomagic wrap (3D Systems 2017).
Anatomical landmarks and curves were digitized using Morphodig 1.5.4 . For a detailed map and designation of digitized anatomical landmarks and curves, please refer to Additional file 1: Tables S1–2 and Figs. S1–2. All specimens were measured twice to reduce potential digitization errors. One scan for each bone was arbitrarily chosen to serve as a template to allow semi-automated placement of surface sliding landmarks onto all surface models. The specimen Callithrix argentata AMNH 184689 was chosen as a template for the humerus and the specimen Callimico goeldii FMNH 153714 for the femur. The same points and curves, as in all specimens, were placed on template and surface sliding landmarks were digitized using IDAV Landmark .
All subsequent morphometric procedures and analyses were done using R . Curves were resampled to the number detailed in Additional file 1: Table S3 using the ‘subsampl.inter’ function published in Botton-Divet et al. . All functions used for patching, sliding, and superimposition are included in the ‘Morpho’ package for R . Surface sliding landmarks were patched on all specimens using the ‘placePatch’ function, then specimens were relaxed against the template model under the minimum bending energy criterion using the ‘relaxLM’ function. Finally, the specimens were relaxed against the Procrustes consensus of the dataset for three iterations using the ‘slider3d’ function also minimizing bending energy. Slid specimens were then superimposed (generalized Procrustes analysis [120, 121]) using the ‘procSym’ function. To allow comparative analyses, a mean per species was computed by superimposing all specimens of the same species. The Procrustes consensus multiplied by the mean of centroid sizes was then used as the mean form of a species for subsequent analyses.
Ecological categories and backbone phylogeny
Based on data available in the literature, we divided species into two categories, one for species with significant use of vertical clinging and leaping (VCL), one for preferential HL species based on the number of leaps to and from vertically oriented supports. Classification sources for each species can be found in Table 1. Considering the absence of quantitative data for the use of vertical support during locomotion for many species, and the differences in the protocols used for the available quantification, the classification of these species relies on qualitative information available in the literature. We are grateful to Prof. Eckhard Heymann (German Primate Center, Göttingen, Germany) for his help and comments on this classification attempt.
To use phylogenetic comparative methods, we built a composite phylogeny using the one published by Aristide et al.  as a backbone. We used Matauschek et al.  for the species within the Leontocebus genus and the time tree of life  was used as a source for Mico melanurus, Saguinus leucopus, and Saguinus niger. Branch lengths for added taxa were computed by scaling the branches from the above sources to match the backbone phylogeny . In order to add a taxon, the branch lengths of the closest related species present in both the backbone and patch phylogenies were used to compute a scaling factor by cross multiplication. The added species may then be integrated into the backbone phylogeny with scaled branch lengths. Leontocebus nigricollis was also replaced by Leontocebus fuscus, as it is the closest relative species according to Cropp et al.  and Matauschek et al. .
Reconstruction of ancestral states is problematic given the repartition of the VCL versus HL character traits across the phylogeny (Fig. 1). An ancestral state reconstruction using parsimony criterion (function ‘asr_max_parsimony’ from the ‘castor’ R package ) leads to an unresolved state at the root of the phylogeny (a probability of 0.5 for both VCL and HL). The use of the likelihood method with stochastic character mapping using the ‘make.simmap’ function from the ‘phytools’ R package , with equal rate parameter and 1000 iterations results in posterior probabilities at the initial node of 0.45 and 0.55 for HL and VCL respectively. This result is highly influenced by the long branch that supports the Leontopithecus genus. Thus the result of such an analysis remains unreliable in this context. To overcome this we base the polarisation of the character on a functional argument. We hypothesize that clinging is a prerequisite to VCL. All Callitrichidae can cling to a certain extent, but trunk-to-trunk leaping is only adopted by some species . VCL is therefore regarded as a derived state in the group.
Statistics and assessment of evolutionary rate shifts
Differences in size between the two ecological categories were tested using two-sided t-test (‘t.test’ function in R) on log-transformed centroid size for both the humerus and femur. These differences were tested while taking phylogeny into account by using a phylogenetic ANOVA  using the ‘phylANOVA’ function from the ‘phytools’ package  with default settings.
In order to assess how strongly our data are driven by evolutionary history, and thus how strong the functional signal is, phylogenetic signal was measured. On centroid size (log), phylogenetic signal was quantified using Blomberg’s K statistic  and tested using the method from Ives et al.  with 10,000 iterations, using the ‘phylosig’ function implemented in the ‘phytools’ package . Phylogenetic signals of shapes were measured using Kmult  and tested using the ‘physig’ function from the ‘geomorph’ package (version 3.3.1 ) with 10,000 iterations.
The effect of size (allometry), locomotor categories and the interaction of these on humeral and femoral shape was tested using MANCOVA. For this analysis, data dimensionality was reduced by selecting the first PCA axes gathering 95% of the total variance (14 and 15 axes for humerus and femur, respectively).
We computed the morphological difference between the two ecological groups by using between-group PCA. We evaluated the group separation by reclassification using the ‘typprobClass’ function from the ‘Morpho’ package  with cross-validation using Wilson’s  method for small samples. The significance of the between-group differences obtained from the between-group PCA was estimated by permutation test of the Euclidean distance between group means with 10,000 iterations. We visualized associated shape deformations by computing shape coordinates on the extreme of the axes and warping a dummy bone to the computed landmark values.
To test for changes in the evolutionary rates across the whole group we used the method developed by Castiglione et al. . This method is based on Phylogenetic ridge regression  and measures the rate of evolutionary change for each branch, allowing the mean for each clade to be computed and then compared to the overall rate observed in the rest of the tree. The significance of rate difference between the considered clade and the rest of the tree is tested by randomization. We ran the ‘search.shift’ function  from the package ‘RRphylo’  with the “clade” option on all clades smaller than one half of the tree.
Availability of data and materials
All surface models are available through the Morphosource portal (https://www.morphosource.org). The identification numbers (ID) of the surface models are provided in Additional file 2.
Vertical clinging and leaping
Principal component analysis
Microfocus computed tomography
Multivariate analysis of covariance
- Kmult :
Multivariate extension of Blomberg’s Kappa as defined in 
Bertram JE, Biewener AA. Allometry and curvature in the long bones of quadrupedal mammals. J Zool. 1992;226:455–67.
Botton-Divet L, Cornette R, Houssaye A, Fabre A-C, Herrel A. Swimming and running, a study of the convergence in long bone morphology among semi-aquatic mustelids (Carnivora: Mustelidae). Biol J Linn Soc. 2017;121:38–49.
Fabre A-C, Cornette R, Peigné S, Goswami A. Influence of body mass on the shape of forelimb in musteloid carnivorans. Biol J Linn Soc. 2013;110:91–103.
Hedrick BP, Mutumi GL, Munteanu VD, Sadier A, Davies KTJ, Rossiter SJ, et al. Morphological diversification under high integration in a hyper diverse mammal Clade. J Mamm Evol. 2019;27:1–13.
Kimura T. Long bone characteristics of primates. Z Morphol Anthropol. 1995;80:265–80.
Mancini IAD, Rieppo L, Pouran B, Afara IO, Braganca FMS, van Rijen MHP, et al. Effects of body mass on microstructural features of the osteochondral unit: a comparative analysis of 37 mammalian species. Bone. 2019;127:664–73.
Milne N, Vizcaíno SF, Fernicola JC. A 3D geometric morphometric analysis of digging ability in the extant and fossil cingulate humerus. J Zool. 2009;278:48–56.
Polly PD. Limbs in mammalian evolution. In: Hall BK, editor. Fins into limbs: evolution, development and transformation. Chicago: University of Chicago Press; 2007. p. 245–68.
Polly PD, Stayton CT, Dumont ER, Pierce SE, Rayfield EJ, Angielczyk KD. Combining geometric morphometrics and finite element analysis with evolutionary modeling: towards a synthesis. J Vertebr Paleontol. 2016;36:e1111225.
Ross D, Marcot JD, Betteridge KJ, Nascone-Yoder N, Bailey CS, Sears KE. Constraint on mammalian forelimb development: insight from developmental disparity. Evolution. 2013;67:3645–52.
Schutz H, Guralnick RP. Postcranial element shape and function: assessing locomotor mode in extant and extinct mustelid carnivorans. Zool J Linn Soc. 2007;150:895–914.
Chazot N, Panara S, Zilbermann N, Blandin P, Le Poul Y, Cornette R, et al. Morpho morphometrics: Shared ancestry and selection drive the evolution of wing size and shape in Morpho butterflies: MORPHO MORPHOMETRICS. Evolution. 2016;70:181–94.
Polly PD. Adaptive zones and the pinniped ankle: a 3D quantitative analysis of carnivoran tarsal evolution. In: Sargis E, Dagosto M, editors. Mammalian evolutionary morphology: a tribute to Frederick S. Szalay. Dordrecht: Springer; 2008. p. 65–194.
Price MV. A functional-morphometric analysis of forelimbs in bipedal and quadrupedal heteromyid rodents. Biol J Linn Soc. 1993;50:339–60.
Adams DC, Rohlf FJ, Slice DE. Geometric morphometrics: ten years of progress following the ‘revolution.’ Ital J Zool. 2004;71:5–16.
Botton-Divet L, Cornette R, Fabre A-C, Herrel A, Houssaye A. Morphological analysis of long bones in semi-aquatic mustelids and their terrestrial relatives. Integr Comput Biol. 2016;56:1298–309.
Denton JSS, Adams DC. A new phylogenetic test for comparing multiple high-dimensional evolutionary rates suggests interplay of evolutionary rates and modularity in lanternfishes (Myctophiformes; Myctophidae): rates and modularity in Lanternfishes. Evolution. 2015;69:2425–40.
Mitteroecker P, Gunz P, Windhager S, Schaefer K. A brief review of shape, form, and allometry in geometric morphometrics, with applications to human facial morphology. Hystrix Italian J Mammal. 2013;24(1):59–66.
Ashton EH, Oxnard CE. Locomotor patterns in primates. Proc Zool Soc Lond. 2009;142:1–28.
D’Août K, Vereecke EE. Primate locomotion. New York: Springer; 2011.
Jenkins FA. The movement of the shoulder in claviculate and aclaviculate mammals. J Morphol. 1974;144:71–83.
Strasser E, Fleagle JG, Rosenberger AL, McHenry HM. Primate locomotion. Boston: Springer; 1998.
Granatosky MC. A review of locomotor diversity in mammals with analyses exploring the influence of substrate use, body mass and intermembral index in primates. J Zool. 2018;306:207–16.
Garber PA. Vertical clinging, small body size, and the evolution of feeding adaptations in the Callitrichinae. Am J Phys Anthropol. 1992;88:469–82.
Heymann EW, Smith AC. When to feed on gums: Temporal patterns of gummivory in wild tamarins. Saguinus mystax and Saguinus fuscicollis (Callitrichinae). 1999;18:459–71.
Garbino GST, Martins-Junior AMG. Phenotypic evolution in marmoset and tamarin monkeys (Cebidae, Callitrichinae) and a revised genus-level classification. Mol Phylogenet Evol. 2018;118:156–71.
Rylands AB, Coimbra-Filho AF, Mittermeier RA. The systematics and distributions of the marmosets (Callithrix, Callibella, Cebuella, and Mico) and Callimico (Callimico) (Callitrichidae, Primates). In: Ford SM, Porter LM, Davis LC, editors. In the smallest anthropoids. Boston: Springer; 2009. p. 25–61.
Nowak RM, Walker EP. Walker’s primates of the World. Baltimore: JHU Press; 1999.
Nyakatura JA. Early primate evolution: insights into the functional significance of grasping from motion analyses of extant mammals. Biol J Linn Soc. 2019;127:611–31.
Nyakatura JA, Heymann EW. Effects of support size and orientation on symmetric gaits in free-ranging tamarins of Amazonian Peru: implications for the functional significance of primate gait sequence patterns. J Hum Evol. 2010;58:242–51.
Youlatos D. Comparative locomotion of six sympatric primates in Ecuador. Ann Sci Nat. 1999;20:161–8.
Yoneda M. Comparative studies on vertical separation, foraging behavior and traveling mode of saddle-backed tamarins (Saguinus fuscicollis) and red-chested moustached tamarins (Saguinus labiatus) in Northern Bolivia. Primates. 1984;25:414–22.
Garber PA, McKenney AC, Mallott EK. The ecology of trunk-to-trunk leaping in Saguinus fuscicollis: implications for understanding locomotor diversity in Callitrichines. Neotropic Primates. 2012;19:1–7.
Porter LM. Dietary differences among sympatric Callitrichinae in Northern Bolivia: Callimico goeldii, Saguinus fuscicollis and S. labiatus. Int J Primatol. 2001;22:3961–92.
Rosenberger AL. Evolution of feeding niches in new world monkeys. Am J Phys Anthropol. 1992;88:525–62.
Garber PA. Locomotor behavior and feeding ecology of the Panamanian tamarin (Saguinus oedipus geoffroyi, Callitrichidae, Primates). Int J Primatol. 1980;1:185–201.
Nadjafzadeh MN, Heymann EW. Prey foraging of red titi monkeys, Callicebus cupreus, in comparison to sympatric tamarins, Saguinus mystax and Saguinus fuscicollis. Am J Phys Anthropol. 2008;135:56–63.
Rylands AB. Marmosets and tamarins: systematics, behaviour, and ecology. Oxford: Oxford University Press; 1993.
Garber PA, Leigh SR. Patterns of positional behavior in mixed-species troops of Callimico goeldii, Saguinus labiatus, and Saguinus fuscicollis in northwestern Brazil. Am J Primatol. 2001;54:17–31.
Garber PA. A comparative study of positional behavior in three species of tamarin monkeys. Primates. 1991;32:219–30.
Youlatos D, Gasc J-P. Comparative positional behaviour of five primates. In: Bongers F, Charles-Dominique P, Forget M, Théry M, editors. Nouragues. Dordrecht: Kluwer Academic; 2001. p. 103–14.
Bicca-Marques JC. Hand specialization, sympatry, and mixed-species associations in Callitrichines. J Hum Evol. 1999;36:349–78.
Buchanan-Smith H. Polyspecific association of two tamarin species, Saguinus labiatus and Saguinus fuscicollis, in Bolivia. Am J Primatol. 1990;22:205–14.
Heymann EW. The relationship between body size and mixed-species troops of tamarins (Saguinus spp.). Folia Primatol (Basel). 1997;68:287–95.
Heymann EW, Buchanan-Smith HM. The behavioural ecology of mixed-species troops of callitrichine primates. Biol Rev. 2007;75:169–90.
Porter LM. Forest use and activity patterns of Callimico goeldii in comparison to two sympatric tamarins, Saguinus fuscicollis and Saguinus labiatus. Am J Phys Anthropol. 2004;124:139–53.
Fleagle JG. Primate adaptation and evolution. New York: Wiley; 1999.
Smith RJ, Jungers WL. Body mass in comparative primatology. J Hum Evol. 1997;32:523–59.
Buchanan-Smith HM, Hardie SM, Caceres C, Prescott MJ. Distribution and forest utilization of Saguinus and other primates of the Pando Department, Northern Bolivia. Int J Primatol. 2000;21:353–79.
Porter LM, Garber PA, Nacimento E. Exudates as a fallback food for Callimico goeldii. Am J Primatol. 2009;71:120–9.
Passamani M, Rylands AB. Feeding behavior of Geoffroy’s marmoset (Callithrix geoffroyi) in an Atlantic forest fragment of South-eastern Brazil. Primates. 2000;41:27–38.
Cunha AA, Vieira MV, Grelle CEV. Preliminary observations on habitat, support use and diet in two non-native primates in an urban Atlantic forest fragment: The capuchin monkey (Cebus sp.) and the common marmoset (Callithrix jacchus) in the Tijuca forest, Rio de Janeiro. Urban Ecosyst. 2006;9:351–9.
Correa HKM, Coutinho PEG, Ferrari SF. Between-year differences in the feeding ecology of highland marmosets (Callithrix aurita and Callithrix flaviceps) in south-eastern Brazil. J Zoology. 2000;252:421–7.
Raboy BE, Canale GR, Dietz JM. Ecology of Callithrix kuhlii and a review of Eastern Brazilian Marmosets. Int J Primatol. 2008;29:449–67.
Ferrari SF, Rylands AB. Activity budgets and differential visibility in field studies of three marmosets (Callithrix spp.). Folia Primatol. 1994;63:78–83.
Fonseca GB, Lacher T. Exudate-feeding by Callithrix jacchus penicillata in Semideciduous Woodland (Cerradao) in Central Brazil. Primates. 1984;25:441–50.
Youlatos D. Locomotion, postures, and habitat use by Pygmy Marmosets (Cebuella pygmaea). In: Ford SM, Porter LM, Davis LC, editors. The smallest Anthropoids. Boston: Springer; 2009. p. 279–97.
Catenacci LS, Pessoa MS, Nogueira-Filho SLG, De Vleeschouwer KM. Diet and feeding behavior of Leontopithecus chrysomelas (Callitrichidae) in degraded areas of the Atlantic Forest of South-Bahia, Brazil. Int J Primatol. 2016;37:136–57.
Dietz JM, Peres CA, Pinder L. Foraging ecology and use of space in wild golden lion tamarins (Leontopithecus rosalia). Am J Primatol. 1997;41:289–305.
Garber PA, Porter LM. The ecology of exudate production and exudate feeding in Saguinus and Callimico. In: Burrows AM, Nash LT, editors. The evolution of exudativory in primates. New York, NY: Springer; 2010. p. 89–108.
Norconk MA. Mechanisms promoting stability in mixed Saguinus mystax and S. fuscicollis troops. Am J Primatol. 1990;21:159–70.
Kupsch D, Waltert M, Heymann EW. Forest type affects prey foraging of saddleback tamarins Saguinus nigrifrons. Primates. 2014;55:403–13.
Veracini C. Habitat use and ranging behavior of the Silvery Marmoset (Mico argentatus) at Caxiuanã National Forest (Eastern Brazilian Amazonia). In: Porter LM, Davis LC, editors. The smallest Anthropoids. Boston: Springer; 2009. p. 221–40.
Lopes MA, Ferrari SF. Foraging behavior of a tamarin group (Saguinus fuscicollis weddelli) and interactions with marmosets (Callithrix emiliae). Int J Primatol. 1994;15:373–87.
Garbino GST. How many marmoset (Primates: Cebidae: Callitrichinae) genera are there? A phylogenetic analysis based on multiple morphological systems. Cladistics. 2015;31:652–78.
Veracini C. O comportamento alimentar de Callithrix argentata (Linnaeus 1771)(Primata, Callitrichinae). In Caxiuanã. Edited by Lisboa PLB. Belém: Mus Para Emílio Goeldi. 1997:437–446.
Soini P. Ecology of the saddle-back tamarin Saguinus fuscicollis illigeri on the Rio Pacaya, northeastern Peru. Folia Primatol. 1987;49:11–32.
Smith JM, Smith AC. An investigation of ecological correlates with hand and foot morphology in callitrichid primates: hand morphology and Exudativory in Callitrichids. Am J Phys Anthropol. 2013;152:447–58.
Garber PA, Teaford MF. Body weights in mixed species troops of Saguinus mystax mystax and Saguinus fuscicollis nigrifrons in Amazonian Peru. Am J Phys Anthropol. 1986;71:331–6.
Poveda K, Sánchez-Palomino P. Habitat use by the white-footed tamarin, Saguinus leucopus: a comparison between a forest-dwelling group and an urban group in Mariquita, Colombia. Neotrop Primates. 2004;12:6–9.
Oliveira ACM, Ferrari SF. Seed dispersal by black-handed tamarins, Saguinus midas niger (Callitrichinae, Primates): implications for the regeneration of degraded forest habitats in eastern Amazonia. J Trop Ecol. 2000;16:709–16.
García-Castillo F, Defler TR. The diet of Saguinus oedipus in a dry tropical forest and the importance of Spondias mombin gum as a “Fallback Food.” Primate Conserv. 2018;67:79.
Aristide L, Rosenberger AL, Tejedor MF, Perez SI. Modeling lineage and phenotypic diversification in the New World monkey (Platyrrhini, Primates) radiation. Mol Phylogenet Evol. 2015;82:375–85.
Anemone RL. The VCL hypothesis revisited: Patterns of femoral morphology among quadrupedal and saltatorial prosimian primates. Am J Phys Anthropol. 1990;83:373–93.
Gebo DL. Vertical clinging and leaping revisited: vertical support use as the ancestral condition of strepsirrhine primates. Am J Phys Anthropol. 2011;146:323–35.
Granatosky MC, Tripp CH, Fabre A-C, Schmitt D. Patterns of quadrupedal locomotion in a vertical clinging and leaping primate (Propithecus coquereli) with implications for understanding the functional demands of primate quadrupedal locomotion: quadrupedal locomotion in Propithecus coquereli. Am J Phys Anthropol. 2016;160:644–52.
Garber PA, Blomquist GE, Anzenberger G. Kinematic analysis of trunk-to-trunk leaping in Callimico goeldii. Int J Primatol. 2005;26:223–40.
Garber PA, Sallenave A, Blomquist GE, Anzenberger G. A comparative study of the kinematics of trunk-to-trunk leaping in Callimico goeldii, Callithrix jacchus, and Cebuella pygmaea. In: Ford SM, Porter LM, Davis LC, editors. The smallest Anthropoids. Boston: Springer; 2009. p. 259–77.
Kinzey WG, Rosenberger AL, Ramirez M. Vertical clinging and leaping in a neotropical anthropoid. Nature. 1975;255:327–8.
Ford SM. Locomotor adaptations of fossil platyrrhines. In: Fleagle JG, Rosenberger AL, editors. The Platyrrhine fossil record. Cambridge: Academic Press; 1990. p. 141–73.
Ford SM, Davis LC. Marmoset postcrania and the skeleton of the dwarf marmoset, Callibella humilis. In: Ford S, Porter L, Davis L, editors. The smallest anthropoids. Boston: Springer; 2009. p. 411–48.
Falsetti AB, Cole TM. Relative growth of the postcranial skeleton in callitrichines. J Hum Evol. 1992;23:79–92.
Burr DB, Ruff CB, Johnson C. Structural adaptations of the femur and humerus to arboreal and terrestrial environments in three species of macaque. Am J Phys Anthropol. 1989;79:357–67.
MacLatchy L. Reconstruction of hip joint function in extant and fossil primates. In Primate locomotion. New York: Plenum Press; 1998. p. 111–30.
Meldrum DJ, Fleagle JG, Kay RF. Partial humeri of two Miocene Colombian primates. Am J Phys Anthropol. 1990;81:413–22.
Ruff CB, Burgess ML, Ketcham RA, Kappelman J. Limb bone structural proportions and locomotor behavior in AL 288–1 (“Lucy”). PLoS ONE. 2016;11:e0166095.
Szalay FS, Dagosto M. Locomotor adaptations as reflected on the humerus of Paleogene primates. Folia Primatol (Basel). 1980;34:1–45.
Gould SJ, Eldredge N. Punctuated equilibria: the tempo and mode of evolution reconsidered. Paleobiology. 1977;3:115–51.
Gould SJ, Gingerich PD. Smooth curve of evolutionary rate: a psychological and mathematical artifact. Science. 1984;226:994–6.
Haldane JBS. Suggestions as to quantitative measurement of the rate of evolution. Evolution. 1949;3:51–6.
Jones KE, Smaers JB, Goswami A. Impact of the terrestrial-aquatic transition on disparity and rates of evolution in the carnivoran skull. BMC Evol Biol. 2015. https://doi.org/10.1186/s12862-015-0285-5.
Kimura M. Evolutionary rate at the molecular level. Nature. 1968;217:624–6.
Simpson GG. Tempo and mode in evolution. New York: Columbia University Press; 1994.
Castiglione S, Tesone G, Piccolo M, Melchionna M, Mondanaro A, Serio C, et al. A new method for testing evolutionary rate variation and shifts in phenotypic evolution. Methods Ecol Evol. 2018;9:974–83.
Fleagle JG, Meldrum DJ. Locomotor behavior and skeletal morphology of two sympatric pitheciine monkeys, Pithecia pithecia and Chiropotes satanas. Am J Primatol. 1988;16:227–49.
Hill WCO. The anatomy of Callimico goeldii (Thomas): a primitive american primate. Trans Am Philos Soc. 1959;49:1.
Jungers WL. Body size and scaling of limb proportions in primates. In: Jungers W, editor. Size and scaling in primate biology. New York: Plenum; 1985. p. 345–81.
Demes B, Fleagle JG, Jungers WL. Takeoff and landing forces of leaping strepsirhine primates. J Hum Evol. 1999;37:279–92.
Ferrari SF. Ecological differentiation in the Callitrichidae. In: Rylands AB, editor. Marmosets and tamarins: behavior, ecology and systematics. Oxford: Oxford University Press; 1992. p. 314–28.
Ferrari SF, Martins ES. Gummivory and gut morphology in two sympatric callitrichids (Callithrix emiliae and Saguinus fuscicollis weddelli) from western Brazilian Amazonia. Am J Phys Anthropol. 1992;88:97–103.
Gusmão AC, de Souza LS, da Costa TM, Velten T, dos Santos-Filho M, Ferronato Formação de grupo misto entre Mico rondoni Ferrari, et al. Saguinus weddelli (Deville, 1849) e Mico melanurus (É Geoffroy, 1812) na Resex do Rio Cautário, Rondônia. Brasil Neotrop Primates. 2010;2015(22):50–1.
Rose MD. Another look at the anthropoid elbow. J Hum Evol. 1988;17:193–224.
Andersson KI. Elbow-joint morphology as a guide to forearm function and foraging behaviour in mammalian carnivores. Zool J Linn Soc. 2004;142:91–104.
Bardua C, Felice RN, Watanabe A, Fabre A-C, Goswami A. A practical guide to sliding and surface semi landmarks in morphometric analyses. Integr Organ Biol. 2019;1(1):16.
Cartmill M. Pads and claws in arboreal locomotion. In: Jenkins FA, editor. Primate locomotion. New York: Academic Press; 1974. p. 45–83.
Johnson LE. Biomechanics of vertical clinging and grasping in primates. PhD thesis. Duke University; 2012.
Rosin S, Nyakatura JA. Hind limb extensor muscle architecture reflects locomotor specialisations of a jumping and a striding quadrupedal caviomorph rodent. Zoomorphology. 2017;136:267–77.
Veracini C. Dados preliminares sobre a ecología de Saguinus niger na Estação Científica Ferreira Penna, Caxiuanã, Pará, Brasil. Neotrop Primates. 2000;8:108–13.
Egler SG. Feeding ecology of Saguinus bicolor bicolor (Callitrichidae: Primates) in a relict forest in Manaus, Brazilian Amazonia. Folia Primatol (Basel). 1992;59:61–76.
Terborgh J. Five New World primates: a study in comparative ecology. Princeton: Princeton University Press; 1983.
Heymann EW. Field observations of the golden-mantled tamarin, Saguinus tripartitus, on the Río Curaray, Peruvian Amazonia. Folia Primatol. 2000;71:392–8.
Collar DC, O’Meara BC, Wainwright PC, Near TJ. piscivory limits diversification of feeding morphology in Centrarchid fishes. Evolution. 2009;63:1557–73.
Garber PA. Diet, foraging patterns, and resource defense in a mixed species troop of Saguinus mystax and Saguinus fuscicollis in Amazonian Peru. Behaviour. 1987;105:18–34.
Amira SD. A highly interactive system for visual data analysis. In: Hansen CD, Johnson CR, editors. The visualization handbook. Amsterdam: Elsevier; 2005. p. 749–67.
Cignoni P, Callieri M, Corsini M, Dellepiane M, Ganovelli F, Ranzuglia G. Meshlab: An open-source mesh processing tool. In: Eurographics Italian Chapter Conference. vol. 2008. Salerno: The Eurographics Association: 2008. pp. 129–36.
Lebrun, R. MorphoDig, an open-source 3D freeware dedicated to biology. In IPC5, August 2018. Paris
Wiley DF, Amenta N, Alcantara DA, Ghosh D, Kil YJ, Delson E, Harcourt-Smith W, Rohlf FJ, St. John K, Hamann B. Evolutionary Morphing. Proceedings of IEEE Visualization 2005 (VIS’05); Minneapolis. 2005, pp. 431–438.
R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria, 2019.
Schlager S. Morpho and Rvcg - shape analysis in R. In: Zheng G, Li S, Szekely G, editors. Statistical shape and deformation analysis. London: Academic Press; 2017. p. 217–56.
Gower JC. Generalized procrustes analysis. Psychometrika. 1975;40:33–51.
Rohlf FJ, Slice D. Extensions of the Procrustes method for the optimal superimpositon of landmarks. Syst Biol. 1990;39:40–59.
Matauschek C, Roos C, Heymann EW. Mitochondrial phylogeny of tamarins (Saguinus, Hoffmannsegg 1807) with taxonomic and biogeographic implications for the S. nigricollis species group. Am J Phys Anthropol. 2011;144:564–74.
Kumar S, Stecher G, Suleski M, Hedges SB. TimeTree: a resource for timelines, timetrees, and divergence times. Mol Biol Evol. 2017;34:1812–9.
Cropp SJ, Larson A, Cheverud JM. Historical biogeography of tamarins, genus Saguinus: the molecular phylogenetic evidence. Am J Phys Anthropol. 1999;108:265–89.
Louca S, Doebeli M. Efficient comparative phylogenetics on large trees. Bioinformatics. 2018;34:1053–5.
Revell LJ. phytools: an R package for phylogenetic comparative biology (and other things): phytools: R package. Methods Ecol Evol. 2012;3:217–23.
Garland T, Dickerman AW, Janis CM, Jones JA. Phylogenetic analysis of covariance by computer simulation. Syst Biol. 1993;42:265–92.
Blomberg SP, Garland T, Ives AR. Testing for phylogenetic signal in comparative data: behavioral traits are more labile. Evolution. 2003;57:717–45.
Ives AR, Midford PE, Garland T. Within-species variation and measurement error in phylogenetic comparative methods. Syst Biol. 2007;56:252–70.
Adams DC. A generalized K statistic for estimating phylogenetic signal from shape and other high-dimensional multivariate data. Syst Biol. 2014;63:685–97.
Adams DC, Collyer ML, Kaliontzopoulou A, Sherratt E. Geomorph. Software for geometric morphometric analyses. R package version 3.0.5. https://cran.r-project.org/package=geomorph. 2017.
Wilson SR. On comparing fossil specimens with population samples. J Hum Evol. 1981;10:207–14.
Kratsch C, McHardy AC. RidgeRace: ridge regression for continuous ancestral character estimation on phylogenetic trees. Bioinformatics. 2014;30:i527–33.
Raia P, Castiglione S, Serio C, Mondanaro A, Melchionna M, Di Febbraro M, et al. RRphylo: Phylogenetic Ridge Regression Methods for Comparative Studies. R Package Version. 2019;2.
Garbino GST, Nascimento FO. Mico humeralifer (Primates: Callitrichidae). Mamm Species. 2014;911:40–7.
Karantanis N-E. Comparative positional behaviour in three captive callitrichid species: Leontopithecus chrysomelas, Saguinus imperator and Cebuella pygmaea. MSc Thesis, University College of London; 2010.
We thank E.W. Heymann and C. Veracini for their help and comments. We thank P. Berles for our bibliographic exchanges. We are grateful to the collection managers and staff for allowing access and loan of the specimens: R.S. Voss, N.Duncan, R.D. Macphee, N.B. Simmons, S. Ketelsen, E. Hoeger and M. Surovy from the AMNH; A.W. Ferguson, L. Smith, J. Phelps and L. Heaney from the FMNH; J. Cuisin, A. Verguin, C. Bens from the MNHN. We thank all those who helped during the scanning process: D. Boyer and J. Gladman from Duke University; Z.-X. Luo and A.I. Neander from the University of Chicago; R. Lebrun from the University of Montpellier. We thank the two anonymous reviewer whose comments helped greatly to improve the early version of this manuscript. We are grateful to E. Kerr for her English editing services. LBD also thanks M. Depraetere for her comments on the initial manuscript.
Open Access funding enabled and organized by Projekt DEAL. We acknowledge support by the Open Access Publication Fund of Humboldt-Universität zu Berlin. LBD received funding from the Einstein Foundation (Grant no.: IPF 2017-389). JAN received funding from the German Research Foundation (Grant no.: DFG NY 63/2-1). The funders have no role in the design of the experiment, analyses or writing of the present manuscript.
Ethics approval and consent to participate
All specimens used are museum specimens, no new specimens were obtained for the present study. Visits of the collections and loans of the specimens were submitted to collection managers authorization through forms available on their respective institution websites.
Consent for publication
Authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
Designation of the anatomical landmarks used for the humerus. Fig. S1. Location of the anatomical landmarks used for the humerus. Table S2. Designation of the anatomical landmarks used for the femur. Fig. S2. Location of the anatomical landmarks used for the femur. Table S3. Number of curve sliding landmarks per curve for humerus and femur. Fig. S3. Deformations associated with PCA axes. Table S4. Test for evolutive rate shift in shape across the phylogeny. Fig. S4. Node numbering used in the Table S4.
Additional file 2.
Details on specimens, scanning resolution, scanning apparatus and operators.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Botton-Divet, L., Nyakatura, J.A. Vertical clinging and leaping induced evolutionary rate shifts in postcranial evolution of tamarins and marmosets (Primates, Callitrichidae). BMC Ecol Evo 21, 132 (2021). https://doi.org/10.1186/s12862-021-01848-z
- Geometric morphometrics