The evolution of the vestibular apparatus in apes and humans
Abstract
Phylogenetic relationships among extinct hominoids (apes and humans) are controversial due to pervasive homoplasy and the incompleteness of the fossil record. The bony labyrinth might contribute to this debate, as it displays strong phylogenetic signal among other mammals. However, the potential of the vestibular apparatus for phylogenetic reconstruction among fossil apes remains understudied. Here we test and quantify the phylogenetic signal embedded in the vestibular morphology of extant anthropoids (monkeys, apes and humans) and two extinct apes (Oreopithecus and Australopithecus) as captured by a deformation-based 3D geometric morphometric analysis. We also reconstruct the ancestral morphology of various hominoid clades based on phylogenetically-informed maximum likelihood methods. Besides revealing strong phylogenetic signal in the vestibule and enabling the proposal of potential synapomorphies for various hominoid clades, our results confirm the relevance of vestibular morphology for addressing the controversial phylogenetic relationships of fossil apes.
eLife digest
Humans, gorillas, chimpanzees, orangutans and gibbons all belong to a group known as the hominoids. This ‘superfamily’ also includes the immediate ancestors and close relatives of these species, however in many instances the evolutionary relationships between these extinct ape species remain controversial.
While DNA can help evolutionary biologists to work out how living species are related to one another, fossils are typically the principle source of information for extinct species. Inferring evolutionary relationships from fossils must be done with caution, but the bony cavity that houses the inner ear – which is involved in balance and hearing and fairly common in the fossil record – has proven useful for tracing the evolution of certain groups of mammals. However, no one had previously looked to see if this structure could give insights into the evolutionary relatedness among living and extinct hominoids.
Urciuoli et al. have now used a 3D imaging technique to capture the complex shapes of the inner ear cavities of 27 species of monkeys and apes, including humans and two extinct apes (Oreopithecus and Australopithecus). The results confirmed that the shape of these structures most closely reflected the evolutionary relationships between the species and not, for example, how the animals moved.
Urciuoli et al. went on to identify features of these bony chambers that were shared within several hominoid groups, and to estimate what the inner ears of the ancestors of these groups might have looked like. The results for Australopithecus, for example, were consistent with it being most closely related to modern humans than other apes, while those for the enigmatic Oreopithecus supported the view that it was a much older species of ape that converged in some respects with other apes still alive today.
The findings highlight the potential of the inner ear for reconstructing the early branches of our family tree. They also offer the prospect of refining the controversial evolutionary relationships within the impressive diversity of extinct ape species.
Introduction
Catarrhine primates (Old World anthropoids) include two extant subclades: cercopithecoids (Old World monkeys) and hominoids (apes and humans). Based on molecular (e.g., Springer et al., 2012) and paleontological (e.g., Harrison, 2013; Stevens et al., 2013) data, both groups diverged during the late Oligocene (≥25 Ma), but experienced very different evolutionary histories. Hominoids first radiated in the early Miocene of Africa (Harrison, 2010; Begun, 2013; Begun, 2015) and subsequently diversified into Eurasia during the middle and late Miocene (Alba, 2012; Begun, 2015). Their diversity and geographic distribution (humans excluded) was much greater during the Miocene than at present, being currently restricted to a few genera in southeastern Asia and Africa. In contrast, extant cercopithecoid lineages started to diversify later and experienced a major radiation during the late Miocene (Jablonski and Frost, 2010), being currently much more diverse and widely distributed than hominoids in both Africa and Asia.
The decimated current diversity of hominoids, coupled with the fragmentary nature of their fossil record, abundant homoplasy (e.g., Larson, 1998), and the lack of known fossil hylobatids prior to the latest Miocene (Harrison, 2016) make it difficult to confidently infer the phylogenetic relationships of extinct hominoids and thus reliably infer the morphotype of the last common ancestor (LCA) of various hominoid subclades. This is required not only to adequately understand the evolutionary history of the group as a whole, but also to reconstruct the LCA of chimpanzees and humans, from which the earliest hominins evolved during the late Miocene. The contribution of Miocene apes to debates about hominoid evolution is thus diminished by the numerous controversies about the phylogenetic position of the former.
For example, putative stem hominoids from the early Miocene of Africa (Begun, 2013; Begun, 2015), such as the proconsulid Ekembo, lack most of the synapomorphies of crown hominoids, such that some authors still contend that they might represent stem catarrhines instead (Harrison, 2010; Harrison, 2013). Even more uncertain is the position of dendropithecids (e.g., Micropithecus) and other small-bodied catarrhines from the early to middle Miocene of Africa, which are generally interpreted as stem catarrhines (Harrison, 2010; Harrison, 2013) but might include stem hominoids as well (Alba et al., 2015; Begun, 2015). The same controversy applies to the European middle to late Miocene Pliobates, recovered as a stem hominoid more derived than proconsulids (Alba et al., 2015) or alternatively as a member of the stem catarrhine pliopithecoid radiation (Nengo et al., 2017). Regarding undoubted hominoids, the distinction between stem and crown taxa is by no means less controversial, in part due to the virtual lack of fossil hylobatids since their origin in the early Miocene until the latest Miocene (Harrison, 2016). This is best exemplified by the late Miocene Oreopithecus from Italy, variously considered a hominid (Moyà Solà and Köhler, 1997; Harrison and Rook, 1997) or a stem hominoid (Nengo et al., 2017). Finally, the phylogenetic placement of Miocene apes from Eurasia is also controversial. For example, most Asian forms have been classically considered pongines (Begun, 2013), but most recently the late Miocene Lufengpithecus from China has been reinterpreted as a stem hominid (Kelley and Gao, 2012; Begun, 2015). Even more controversial is the position of the middle to late Miocene European dryopithecines, interpreted as either stem hominids (Alba, 2012; Alba et al., 2015) or hominines (Begun, 2013; Begun, 2015), and further controversies apply when trying to decipher the phylogenetic relationships among various members of this extinct group.
Deciding among phylogenetic hypotheses for Miocene apes has consequences for our current understanding of hominin origins, from calibrating molecular data to estimate their divergence time to the reconstruction of the ancestral locomotor repertoire from which the earliest bipeds arose. For example, the recently described dryopithecine Danuvius from Germany has been used to argue that bipedal and suspensory adaptations characterized the last common ancestor of crown hominids (Böhme et al., 2019). However, without a phylogenetic analysis supporting a more basal branching of Danuvius compared with the older dryopithecine Pierolapithecus (Moyà-Solà et al., 2004), which was orthograde but lacked adaptations to both bipedalism and suspension (Alba et al., 2010; Alba, 2012), the implications for the ancestral condition of the group remain moot.
The recent recovery of enamel proteome sequences from the early Pleistocene (1.9 Ma) fossil pongine Gigantopithecus (Welker et al., 2019) offers some hope that molecular data will become available for Miocene apes sometime in the future. In the meantime, to better resolve the phylogeny of Miocene hominoids, the search for morphological features not very prone to homoplasy is crucial. Anatomical structures that develop early during development and are not remodeled thereafter (such as the enamel-dentine junction of teeth) represent the best candidates (Corruccini, 1987; Skinner et al., 2008). In this regard, the inner ear is also a very promising anatomical area (Spoor and Zonneveld, 1998; Spoor et al., 2003), even if thus far its phylogenetic implications have been mostly explored for fossil hominins only (Quam et al., 2016; Conde-Valverde et al., 2018; Ponce de León et al., 2018; Beaudet, 2019a; Beaudet et al., 2019b), while its application to fossil apes has been mainly devoted to locomotor inferences (David et al., 2010; Malinzak et al., 2012; Rook et al., 2004; Ryan et al., 2012).
Housed in the highly mineralized petrosal bone, which is frequently preserved in the fossil record, the inner ear is composed of a series of endolymph-filled membranes (encased in the corresponding bony labyrinths), namely the cochlea or cochlear duct (involved in hearing) and the vestibular apparatus (devoted to balance and vision stability). The vestibule consists of three (anterior, posterior, and lateral) semicircular canals (SCs) and two macular organs (utricle and saccule). The approximately orthogonal SCs sense angular accelerations and decelerations of the head, while the maculae perceive linear accelerations and thus provide gravity reference (Spoor and Zonneveld, 1998; Rabbitt et al., 2004; Johnson Chacko et al., 2018; Cheung and Ercoline, 2018).
Differences in the relative size and morphology of the SCs have been correlated with locomotor agility (Spoor et al., 1994; Spoor et al., 2007; Walker et al., 2008; Silcox et al., 2009; Ryan et al., 2012; Perier et al., 2016) and positional behavior (Le Maître et al., 2017), albeit not without criticism (Rae et al., 2016; Benson et al., 2017; Coutier et al., 2017). On the other hand, the bony labyrinth morphology has been considered of great importance for phylogenetic reconstruction in various mammals including hominins (Quam et al., 2016; Conde-Valverde et al., 2018; Ponce de León et al., 2018; Beaudet, 2019a; Beaudet et al., 2019b) and nonhuman primates (Lebrun et al., 2010; Lebrun et al., 2012; Grohé et al., 2016; Mennecart et al., 2017; Costeur et al., 2018; Schwab et al., 2019). While previous research in hominoids has yielded encouraging results for phylogenetic reconstruction (Spoor and Zonneveld, 1998; Spoor et al., 2003; Rook et al., 2004; Gunz et al., 2012), according to some authors phylogeny may not be a major component of ape vestibular morphology (Le Maître et al., 2017).
Determining if and to what extent inner ear anatomy reflects phylogeny among extant hominoids is central for assessing the potential of this anatomical area to more confidently resolve the controversial phylogenetic relationships of fossil apes. To provide insight into this question, we use phylogenetically-informed statistical analyses to test the significance and quantify the amount of phylogenetic signal of vestibular shape captured by three-dimensional geometric morphometrics (3DGM) in living hominoids and a broader sample of extant anthropoids. To capture vestibular shape, we mostly rely on a landmark-free, deformation-based 3DGM approach that takes the whole surface into account (Glaunés and Joshi, 2006; Durrleman et al., 2012b; Durrleman et al., 2012a) and enables integrating the spatial trajectory of the semicircular canals with their thickness and volumetric variations, the latter two being more difficult to assess based on mainstream landmark-based 3DGM. Since our results reveal the presence of strong phylogenetic signal in the vestibular morphology, we also employ maximum likelihood methods (Felsenstein, 1988; Schluter et al., 1997) to reconstruct the ancestral vestibular morphology for the LCA of main hominoid subclades (crown hominoids, hominids, and hominines), with the aim to identify phylogenetically informative characters that can be used in formal cladistic analysis. To test the reliability and illustrate the usefulness of our approach from a phylogenetic viewpoint, we also include two extinct hominoid taxa: the early hominin Australopithecus and the aforementioned controversial late Miocene ape Oreopithecus. The well-known phylogenetic placement of Australopithecus as the sister-taxon of humans predicts that its vestibular morphology will be somewhat derived towards the modern human condition. On the other hand, our analysis will enable testing the competing phylogenetic hypotheses for Oreopithecus (stem hominoid vs. hominid), thereby illustrating the potential of our method for clarifying the controversial affinities of extinct apes.
Results
Semicircular canal shape variation among anthropoids
In spite of a similar spatial configuration of the three SCs, their trajectory, stoutness and relative proportions are quite variable among anthropoids (Figure 1). A bgPCA performed among major anthropoid groups (platyrrhines, cercopithecoids, hylobatids, and hominids) enables their accurate distinction (Figure 2, Figure 2—figure supplement 1a, Figure 2—source data 1), as shown by classification results (more than 95% individuals correctly classified after cross-validation; Table 1). In particular, bgPC1 discriminates hylobatids from hominids and all the remaining taxa (Figure 2, Figure 2—figure supplement 1c). A landmark-based analysis applied to the same sample yields very similar results except for hylobatids, due to the reasons explained in the next section. Shape variation in the analyzed sample accounts for a strong phylogenetic signal (Kmult = 1.248, p<0.001), and this also holds for the first three bgPCs separately (see below).
When analyzed individually, we recover a similarly strong phylogenetic signal (λ = 1 and K = 1.15; Table 2) for bgPC1 (68.8% of total variance). bgPC1 captures differences in thickness and cross-section of the SCs, and is also driven by the development of the macular organs relative to the canals. Great apes fall on positive values (Figure 2, Figure 2—figure supplement 1a) due to their stout and flattened SCs combined with an extensive vertical compression of the anterior canal, a more anterosuperior insertion of the lateral canal into the vestibule, and a greater volume of the vestibular recesses relative to the canals (Figures 3a–e, 4a–e and 5a–e). Hylobatids, as well as colobine and papionin cercopithecoids, showing slender and elongated canals but maintaining well developed ampullae, largely overlap on negative values, while cercopithecins and platyrrhines display intermediate values due to their slightly inflated SCs (Figures 2, 3f–h, 4f–h and 5f–h). bgPC2 (19.6% of total variance) also bears strong phylogenetic signal (λ = 0.91; K = 1.51; Table 2), with variance accumulating among rather than within clades (as indicated by K > 1). This axis separates platyrrhines—especially Ateles (Figure 1i)—from other anthropoids due to the more reduced lateral canal in the former, which is inversely proportional to anterior canal development and vertical elongation (Figure 2a). In contrast, Gorilla (Figures 1a, 3a, 4a and 5a) occupies the opposite end of the distribution due to its large lateral canal and reduced anterior one, the latter being also vertically compressed, whereas the remaining hominoids show intermediate values along bgPC2. bgPC3 (11.6% of variance), which is driven by both trajectory and relative size of the SCs (Figure 2b), still displays a strong and significant phylogenetic signal (λ = 1 and K = 1.16; Table 2). Hylobatids (Figures 1f–h, 3f–h, 4f–h and 5f–h) display the highest scores for bgPC3, as a result of the right to acute angle formed by the apex of the common crus (CC), the latter being also shorter, a more anterosuperiorly projecting and long anterior canal, an obtuse angle between the planes identified by the anterior and lateral canals, and a more developed lateral canal relative to the posterior one, which is also posteriorly oriented. Most of the taxa fall within moderate positive and moderate negative values, with some cercopithecoids (Cercopithecus, Macaca, Papio and Nasalis, among others) and Cebus occupying the negative end of the distribution. African apes, cercopithecines, and Nasalis occupy an intermediate position, being characterized by a well-developed and anteriorly inclined posterior canal, an obtuse to right angle of the CC apex, and an obtuse angle between the vertical canals, combined with a larger lateral canal. Orangutans fall on moderately positive values and are distinguished from other great apes and humans by an anterosuperiorly projecting anterior canal (less than in hylobatids).
When fossil specimens are plotted a posteriori onto the tangent space identified by extant taxa, Oreopithecus (BAC 208; Figures 1n and 6a) falls on moderately positive scores for bgPC1 (where the distributions of hominids, cercopithecins, and platyrrhines overlap), while the two Australopithecus individuals (StW 573 and StW 578; Figures 1o,p and 6b–c) fall within the range of living great apes and humans (Figure 2, Figure 2—figure supplement 1a). This is due to the volumetric proportions of their SCs and the possession of voluminous vestibular recesses (although the latter character is less pronounced in Oreopithecus). In BAC 208, SC volume is greater on the lateral canal, while the two Australopithecus specimens possess stouter vertical canals. When the relative size of the SCs is taken into account, StW 573 shows more evenly developed canals than StW 578 and BAC 208, which both display a smaller lateral canal. Their position along the bgPC2 axis reflects these differences, with StW 573 falling on positive scores (close to the mean value for hominids and cercopithecoids; Figure 2a, Figure 2—figure supplement 1b) and the other two specimens occupying moderately negative values (within the range of extant catarrhines and approaching that of platyrrhines; Figure 2a, Figure 2—figure supplement 1b). This is caused by the comparatively smaller lateral canal and by the large vertical canal in StW 578 and in BAC 208. Due to the acute angles between the planes identified by the anterior and lateral canals and that between the planes of the posterior and anterior canals, Oreopithecus falls at the negative end of the extant anthropoid distribution for the bgPC3 (Figure 2b, Figure 2—figure supplement 1c). On the other hand, the two Australopithecus specimens occupy more intermediate values because of the possession of a right angle between the planes of the aforementioned canals.
When the bgPCs are considered at the same time (Figure 2), the australopith specimens fall well within the great ape and human range. This is further supported by their posterior probabilities of group membership based on the proximity of fossil specimen scores to groups centroids, with StW 573 and 578 being classified as hominids with p=0.678 and p=0.190, respectively (Table 3). StW 573 falls close to Pan and Homo, whereas StW 578 occupies an intermediate position between humans and orangutans due to its stouter volumetric proportions. When the posterior probabilities are computed using the centroids of the hominoid genera, StW 573 is classified as Pan as first option (p=0.368) and as Homo as second (p=0.264), while StW 578 is more clearly classified as Homo with p=0.727 (Table 4). These results suggest that both australopith specimens show vestibular similarities with extant humans, but that StW 573 display a more plesiomorphic (chimpanzee-like) morphology. In turn, Oreopithecus shows a mosaic of vestibular features (Pan-like volumetric proportions, small lateral canal, and acute angles between the anterior and both the posterior and the lateral canals) that does not match the condition of any extant taxon (Figure 2). The posterior probabilities indicate closest similarities with cercopithecoids, followed by hominids and platyrrhines, although in all instances it falls outside the variability of the extant members of these groups (p<0.05; Table 3). When comparisons are restricted to hominoid genera, Oreopithecus appears more similar to humans than to any ape genus, but again with a posterior probability that indicates significant differences with all of them (p<0.05; Table 4).
The multivariate regression between shape (deformation fields) and size (log-transformed volume of the vestibule) shows a significant correlation (i.e., allometry) at p<0.001, but nevertheless explains only a limited portion of the variance (R2 = 0.192). Bivariate regressions of the bgPCs against log-transformed cube root of vestibular volume reveal a significant correlation only for bgPC1 (R2 = 0.635, p<0.001, Table 2). A visual inspection of the scatter of points (Figure 7a) suggests that allometry for bgPC1 is more marked in hominids. This is confirmed when separate regressions are performed for hominids (R2 = 0.480, p<0.001) and the rest of the sample (R2 = 0.058, p<0.01), with the former displaying also a higher slope (Table 5). When phylogeny is considered by means of PGLS regression (Table 5), the correlation for the whole sample is still significant but explains much less variance (R2 = 0.261, p<0.01), becoming non-significant for hominids and the rest of the sample separately. bgPC3 shows a low yet significant correlation with volume for non-hominids, which becomes non-significant after PGLS correction (Table 5). Both Australopithecus and Oreopithecus overlap with the hominid scatter of points, with the two australopith specimens falling above the hominid regression line, whereas BAC 208 falls slightly below (although well above that of non-hominid anthropoids).
The bivariate regression between the log-transformed cube root of the SC volume and SC length shows in all instances a significant correlation that nevertheless only explains a limited (ca. 20–30%) amount of variance (Table 5, Figure 7b). Isometry cannot be rejected for anthropoids as a whole, but a negatively allometric relationship emerges (revealing that length increases faster than volume) when hominids and other taxa are analyzed separately (Table 5). The latter is confirmed by PGLS regressions for the whole sample and the two groups separately, which further explain a higher proportion of variance (Table 5), although the hominid regression is not significant with all probability due to small sample size. The bivariate plot (Figure 7b) shows an allometric grade shift between hominids and the remaining anthropoid taxa, which is confirmed by ANCOVA results—indicating no significant differences (F = 0.705, p=0.403) between the allometric slopes but significantly different intercepts (F = 263.26, p<0.001) between the two groups. This indicates that hominids possess more voluminous (i.e., stouter) canals than other anthropoids at equal lengths, with only minimal overlap. All the fossil specimens display hominid-like volumetric proportions (Figure 7b): StW 573 and 578 fall slightly above the hominid regression line, whereas BAC 208, due to its slenderer semicircular canals, falls below (although much closer than to the non-hominid regression line).
Exploration of a preexisting group structure in the tangent space of the vestibular shape
Recently, caution has been advised regarding the use of between-group PCA (bgPCA) applied to 3D geometric morphometric (GM) data, as it could produce spurious grouping when there are fewer groups than variables (Cardini et al., 2019). However, the same study also highlighted that the presence of a strong covariance among the variables (as found in many biological structures) largely reduces the magnitude of the problem. Interestingly, due to the properties of diffeomorphisms, the set of momenta is expected to be highly correlated, as close momenta tend to covary.
Prior to computing the bgPCA, we explored the principal components resulting from the vestibular shape GM analyses to investigate the presence of a preexisting group structure, which was found to be similar to that showed by bgPCA for both extant and fossil taxa (Figure 8). We used hierarchical clustering analysis (HCA) on the deformation fields for assessing the probability of correct classification of individuals according to the groups used in the bgPCA. The confusion matrix resulting from the HCA shows that most individuals are correctly identified in the corresponding groups (Table 6). Only in the deformation-based analyses hylobatids show a low percentage of classification (24% of the individuals), mostly due to the great similarity in the volumetric proportions and surface shape of the vestibule between this group and cercopithecoids.
Another way to ascertain that the group separation observed in our bgPCA is not the result of any bias is to compare different kind of analyses with a same data set and the same grouping factor. We thus compared our deformation-based results with those obtained from a configuration of 3D semilandmarks commonly used to investigate the vestibular shape (Gunz et al., 2012). Both analyses, based on the full primate sample, yielded a similar group separation (compare Figure 2 with Figure 9), as shown by the components resulting from the bgPCA. These results are also coherent with the biological reality, enabling the discrimination of major anthropoid clades in agreement with their phylogenetic relationships. In the landmark-based approach, hylobatids largely overlap with great apes in both bgPC1 (occupying an intermediate position) and bgPC2 (Figure 9a), but can be distinguished from them (and other anthropoids) to a large extent based on bgPC3 (Figure 9b). Shape variation in the analyzed sample, as captured by landmark-based 3DGM, accounts for a strong phylogenetic signal (Kmult = 0.973, p<0.001), and this also holds for the first three bgPCs separately (see below).
Shape differences along bgPC1 (53.3% of total variance) in the landmark-based approach embed a strong phylogenetic signal (λ = 1 and K = 1.26; Table 7). This component correlates with the insertion of the lateral canal on the vestibule, the size and shape of the posterior canal, and the roundness of the SCs. Great apes and humans fall on negative values for the first axis (Figure 9c), as they are characterized by smaller canals compared to the size of the vestibular recesses and less rounded SCs (particularly the anterior one, which is vertically compressed). Hylobatids stand on intermediate scores that largely overlap with the hominid range (Figure 9) due to a combination of long SCs, a vertically compressed anterior canal, and well separated lateral and posterior canals (since posterior canal is posteriorly displaced and the lateral canal inserts anteriorly in the vestibule). In contrast, Old world monkeys tend to be located on positive values of bgPC1 and display a protruding lateral canal that intersects the plane defined by the posterior canal.
bgPC2 (29.3% of variance) separates platyrrhines—especially Ateles found in the most negative scores (Figure 1i)—from other anthropoids due to the more reduced lateral canal in the former, which is inversely proportional to anterior canal development and vertical elongation (Figure 9d). This pattern is shared, although to a lesser extent, by humans (Figure 1d) and Theropithecus (Figure 1j), which possess more developed anterior and posterior canals relative to the lateral one. In contrast, Gorilla (Figure 1a) occupies the positive end of the distribution due to its large lateral canal and reduced anterior one, the latter being also vertically compressed, whereas the remaining hominoids show intermediate values along bgPC2. The identified phylogenetic signal displayed by bgPC2 is slightly reduced as compared to bgPC1, but still significant and very high (λ = 0.9 and K = 1.08;).
The bgPC3 (17.4% of variance) is driven by both trajectory and size of the SCs, especially of the posterior one (Figure 9e), and displays less phylogenetic signal than bgPC1 and bgPC2 (λ = 0.93 and K = 0.53, Table 7). Along this axis Homo, Ateles and Theropithecus (Figure 9b) occupy the positive end, as due to their large and anteriorly inclined posterior canal that protrudes laterally, an obtuse angle of the CC apex, large posterior canal, and an angle between the anterior and posterior canals that is close to 90°. Great apes and the majority of non-hominoid taxa fall in an intermediate position (Figure 9b), differing from the aforementioned genera by the less obtuse angle in the CC apex and a larger lateral canal. Hylobatids and, to a lesser extent, Trachypithecus show the lowest scores for bgPC3, as the result of the right to acute angle formed by the apex of the CC, more developed anterior and lateral canals relative to the posterior one, and a posteriorly tilted posterior canal.
When plotted onto the tangent space of extant taxa, the two Australopithecus specimens overlap with the range of great apes and humans for all the bgPCs (Figure 9, Figure 9—source data 1). They show a vertically compressed anterior canal together with well-separated lateral and posterior canals. The lateral canal is also moderately sinuous and its ampullar portion bends upwards, thus resulting in a negative score for bgPC1. The two specimens can be distinguished from one another by means of bgPC2 (Figure 9a), with StW 573 falling on positive values and StW 578 occupying negative ones. This is explained by the smaller lateral canal relative to the vertical SCs in StW 578, which therefore overlaps with extant human variation in both bgPC1 and bgPC2, whereas StW 573 overlaps instead with chimpanzees and bonobos (Figure 9a). Along bgPC3, both specimens display intermediate values due to the large posterior canal and for the almost right angle between the planes of the anterior and posterior SCs, overlapping with all extant anthropoids except hylobatids.
In turn, along bgPC1 Oreopithecus displays more positive values than hominoids and falls well within the range of non-hominoid anthropoids (it only slightly overlaps with the positive end of the hominoid distribution; Figure 9) due to its more coplanar lateral canal that almost intersects the plane defined by the posterior canal. Furthermore, due to its small lateral canal and fairly short CC, BAC 208 displays an intermediate value for bgPC2, within the range of extant catarrhines and slightly above the positive end of platyrrhine distribution (Figure 9a). Finally, along bgPC3 Oreopithecus falls on a positive score, differing from hylobatids, as a result of the acute angles found between anterior canal plane and those defined by both the posterior and lateral canals (Figure 9b).
Overall, the two 3DGM techniques used in this paper generally yielded similar results except for hylobatids and Oreopithecus (along bgPC1 alone). This is attributable to differences in the underlying methodological assumptions of each method when computing shape variation. In particular, our 3DGM landmark protocol measures the spatial trajectory of SCs based on their midline skeleton (Gunz et al., 2012) and hence it does not capture differences in volumetric proportions. In contrast, by comparing surfaces as a whole (Durrleman et al., 2012b; Durrleman et al., 2012a), the deformation analysis is particularly sensitive to volumetric differences. In addition, the amount of identified phylogenetic signal is very similar for both techniques, affecting the entire variance. Together with the results of the HCA, we confidently show that the separation found between the groups in the bgPCA of this study already exists in the shape data and that is not a spurious effect produced by the bgPCA method itself.
Hominoid phylomorphospace of the vestibule
The phylomorphospace approach applied to vestibular shape variation in hominoids infers different branch lengths for hominids and hylobatids from the ancestral morphology estimated for crown hominoids, which falls much closer to great apes and humans than to hylobatids for both bgPC1 and bgPC3 (Figure 10). According to our reconstructions based on the extant taxa (Figure 10), the crown hominoid LCA vestibular morphology (Figure 11a) would be characterized by equally developed and slightly inflated SCs, a fairly vertically compressed anterior canal, and by the slender portion of the lateral canal connecting more anteromedially with the vestibule. The estimated morphology for the LCA of hylobatids (Figure 11b) resembles to some extent that of the crown hominoid LCA (slightly vertically compressed anterior canal and lack of intersection between the lateral and posterior canals), combined with more monkey-like features (markedly slender canals with inflated ampullae), and others exclusive to hylobatids (an obtuse angle between a slightly anteriorly protruding anterior canal and a small posterior canal relative to the others). Hylobatid genera are generally less diverging from one other than great apes and humans. Hoolock (Figures 1f, 3f, 4f and 5f) apparently displays the most primitive morphology among hylobatids (with equally developed, rounded, almost orthogonal canals, and less anteriorly protruding anterior canal), while Hylobates (Figures 1h, 3h, 4h and 5h) and, to a lesser extent, Symphalangus (Figures 1g, 3g, 4g and 5g) show the slenderest SCs and an extremely anteriorly protruding anterior canal (as noted by Le Maître et al., 2017).
The reconstructed morphologies for LCAs of crown hominids (Figure 11c) and, to a lesser extent, crown hominines (Figure 11d) and the Pan-Homo clade (Figure 11e) are not very far from the crown hominoid ancestral condition (Figure 11a) for any of the first three bgPCs. The crown hominid LCA is characterized by a lateral insertion of the slender portion of the lateral canal on the vestibule, a moderate medioventral displacement of the posterior canal, and an increased vertical compression of the anterior canal, in combination with thick and bulgy canals and well-developed vestibular recesses. The hominine and Pan-Homo clade LCAs possess stouter canals and are very similar to one another, distinguished only by the size of the lateral and anterior canals. The LCA of the Pan-Homo clade shows a larger and less vertically compressed anterior canal, and a smaller lateral one, which also connects more anteriorly with the vestibule. The hominin (Australopithecus-Homo) LCA (Figure 11f) is closer to Homo, being characterized by the stoutest volumetric proportions, with larger vertical canals relative to the Pan-Homo clade LCA, yet smaller than those found in humans. Its anterior canal is more vertically compressed than in Homo, rather resembling the morphology of Pan, while the posterior canal is rounded, thus being intermediate between the human (laterally projecting) and the chimpanzee (laterally compressed) morphology.
Among crown hominids, Pan (Figures 1b–c, 3b–c, 4b–c and 5b–c) more closely resembles the morphology of the inferred LCAs of crown hominids and hominoids (Figure 11a,c) in the moderately inflated and equally developed SCs and in the degree of the vertical compression of the anterior canal. Pongo (Figures 1e, 3e, 4e and 5e) occupies the positive end along bgPC1 (Figure 10) due to the possession of relatively small but extremely stout canals (especially the anterior one and the common crus). It also exhibits the most vertically compressed anterior canal and a ‘triangular’ lateral canal (i.e., showing straight slender portions of the bony labyrinth close to the ampulla and to the connection with the vestibule, as previously outlined by Spoor and Zonneveld, 1998). Finally, Gorilla and Homo (Figures 1a, 3d, 4d and 5d) are derived in opposite directions along bgPC2 (Figure 10a), as the former exhibits increased lateral canal radius with a flattened cross-section, while humans retain quite rounded canals shape and cross-section, and show a reduction of the lateral canal, as opposed to more developed anterior and posterior canals.
Australopithecus is closer to humans than to any great ape, due to the possession of large vertical canals and stout volumetric proportions, being closest to the LCA of the Australopithecus-Homo clade and only slightly more derived in the same direction as Pongo due to the stouter canals (Figure 10). In contrast, Oreopithecus appears more plesiomorphic than humans and extant great apes in volumetric proportions, resembling those found in the LCA of crown hominoids. Nevertheless, Oreopithecus is clearly distinct from hylobatids in both SC stoutness and shape, being most clearly distinguished from gibbons and siamangs based on the acute angles defined the anterior canal with both the posterior and lateral canals (Figure 10b). Oreopithecus is also derived in terms of SC size (with the vertical canals much larger than the lateral one), similarly to humans albeit to a greater extent, and opposite to gorillas, due to the remarkably smaller lateral canal.
Discussion
Previous research on the morphology of the vestibular apparatus among extant mammals has focused on its relationship to positional behavior (Spoor et al., 2007; Perier et al., 2016; Le Maître et al., 2017), particularly in order to make locomotor inferences in extinct species (Spoor et al., 1994; Walker et al., 2008; Silcox et al., 2009; Ryan et al., 2012). However, the phylogenetic signal embedded in vestibular morphology has not been adequately quantified among hominoids, because previous attempts were either exploratory (Gunz et al., 2012) or based on a restricted sample (Le Maître et al., 2017). Our results indicate that main anthropoid groups can be distinguished based on vestibular shape variation, and that there are important differences not only between hylobatids, great apes, and humans, but also among extant great ape genera. A significant phylogenetic signal is found to affect the entire variance of the anthropoid sample. Thus, the shape of the SCs is overall informative from a phylogenetic viewpoint—as hypothesized for strepsirrhine primates (Lebrun et al., 2010) and carnivorans (Schwab et al., 2019), but in contrast to previous results for hominoids (Le Maître et al., 2017) and some other mammals (Grohé et al., 2016; Costeur et al., 2018).
Based on the analysis of the shape of the vestibular apparatus, we identify several potential hominoid synapomorphies (Table 8; Figure 12a–g), including among others a posteromedially displaced posterior canal and a straight segment between the lateral-most point of the lateral canal and its anteromedially situated insertion on the vestibule. These features result in an anteromedially located lateral canal (i.e., the plane defined by the posterior canal is always separated from the trajectory of the lateral canal, even when the latter is well developed, as in Gorilla and Hylobates). This would imply an increased sensitivity for angular accelerations occurring along the coronal plane, which has been correlated with orthogrady in extant hominoids (Le Maître et al., 2017).
The most evident character shared by extant apes and humans (even if somewhat variable in the latter and in Hoolock), and further displayed by the extinct genera analyzed here, is the vertical compression of the anterior canal (Figure 12a), as noted for great apes only in a previous analysis (Spoor and Zonneveld, 1998). Since the subarcuate fossa arguably constrains the shape of the anterior canal (Jeffery et al., 2008), the hominoid morphology might be related to the absence of the fossa in great apes and siamangs (Moyà-Solà and Köhler, 1993; Gannon et al., 1988; Spoor and Leakey, 1996). However, the combination of a well-developed fossa and marked vertical compression of the anterior SC found in Hylobates argues against this hypothesis. The latter is further rejected by the cercopithecoid morphology, characterized by a rounded anterior SC even in the largest terrestrial genera (Papio, Theropithecus, and Mandrillus), which unlike other cercopithecoids display a much reduced (or even absent) subarcuate fossa (Gannon et al., 1988; Spoor and Leakey, 1996).
As a result of vertical compression, all hominoids display the anterior canal projected anterosuperiorly to some extent (less accentuated in Homo and Hoolock). Such a projection of the anterior canal in Hylobates has been interpreted as a hylobatid synapomorphy (Le Maître et al., 2017). Our results raise doubts about the latter view and indicate instead that the anterior projection of this canal is variable within both hylobatids and hominids, and that the LCA of crown hylobatids might have not shown the patent elongation of the anterior canal that is found in Hylobates and Symphalangus. On the other hand, here we identify two potential synapomorphies for hylobatids (Table 8; Figure 12h,i): an obtuse angle between the planes defined by the anterior and posterior canals, and a posteriorly-inclined posterior canal, which is smaller relative to the anterior and lateral ones.
In turn, bgPC1 clearly separates hominids (great apes and humans) from cercopithecoids and hylobatids and enables the identification of some potential hominid synapomorphies (Table 8; Figure 12j,k). In particular, hominids differ from other anthropoids, including hylobatids, by derived volumetric proportions of the SCs (stouter canals relative to their length, or shorter canals relative to their volume) even when size-scaling considerations are taken into account, as well as by the possession of more extensive vestibular recess for a similar size of the SCs (Table 8; Figure 12j,k), as reflected in bgPC1. In particular, cercopithecoids and hylobatids completely overlap due to the possession of slender SCs, while hominids as a whole (even if more markedly orangutans) differ by their swollen and relatively shorter SCs (with only few cercopithecins falling within the hominid range). This might be related to the fact that hylobatids and cercopithecoids, unlike great apes, are swift moving animals that make fast and large head movements, thus requiring a limited duct sensitivity to avoid overstimulation and a quick response to angular displacement (Spoor and Zonneveld, 1998). This hypothesis is supported by biophysical models suggesting that the length of the membranous ducts is inversely proportional to their sensitivity and that a larger lumen of the ducts correlates with a reduced steadiness of the response to an external angular stimulus (i.e., the abrupt change of the position and/or posture) (Muller, 1994; Rabbitt et al., 2004). Therefore, species with shorter and thicker ducts (such as hominids) require more time to perceive and react to sudden head displacements, while being more sensitive to fine movements. Nevertheless, caution must be used when inferring the lumen of the ducts from that of the bony canals, as the amount of the SC cross-section occupied by perilymphatic space is variable depending on the species (Ramprashad et al., 1984; Spoor and Zonneveld, 1998).
Superimposed on the aforementioned hominoid and hominid synapomorphies (Table 8; Figure 12), there are also marked differences among hominid genera. Such differences mainly relate the relative size among the SCs (which varies particularly along bgPC2 and bgPC3), while hylobatid genera are less diverse in this regard. Chimpanzees and bonobos are characterized by equally developed SCs and a moderately short CC. Similarly, orangutans possess evenly proportioned SCs and can be distinguished from Pan by a shorter CC, more inflated canals, and a greater vertical compression of the anterior SC. Gorillas display the largest intrageneric variability in the studied sample, especially with regard to SCs slenderness, coupled with some other distinctive traits (obtuse angle of the CC apex, and longer lateral SC and CC).
Humans differ from apes in the enlarged vertical canals, a laterally protruding and inferiorly displaced posterior canal, and a moderately smaller lateral canal. Relatively enlarged vertical canals are also found in Theropithecus. The human morphology has been linked to bipedalism (Spoor et al., 1994; Spoor et al., 2003), as accelerations during bipedal walking mainly occur along the vertical axis, so that the broadly similar morphology of Theropithecus might be related to the bipedal shuffling characteristic of this genus during foraging, causing them spend an extremely large amount of time with an erect trunk posture (Wrangham, 1980).
Irrespective of the functional implications of the variation in vestibular morphology among anthropoids, our analyses show that this variation bears strong phylogenetic signal and, hence, has potential for reconstructing the evolutionary history of the group—particularly hominoids, which in spite of their extant decimated diversity are more variable in this regard (particularly when size differences among the SCs are considered) (Spoor and Zonneveld, 1998; Le Maître et al., 2017) than the taxonomically more diverse cercopithecoids. Although functional demands frequently lead to the independent evolution of similar morphologies (homoplasy), often function is not decoupled from—but superimposed on—phylogeny, with many clades being characterized by synapomorphic features linked to the adaptation for new functions. Therefore, to the extent that vestibular morphology appears to be linked to positional behavior (Spoor et al., 1994; Spoor et al., 2007; Walker et al., 2008; Silcox et al., 2009; Ryan et al., 2012; Le Maître et al., 2017), the higher variation of vestibular morphology displayed by hominoids compared to cercopithecoids agrees with the more diverse and varied locomotor repertoires of the former. This is because cercopithecoids as a whole are largely pronograde terrestrial quadrupeds that mostly differ in the degree of arboreality vs. terrestriality (Fleagle, 2013; Gosselin-Ildari, 2013), while crown hominoids are characterized by a derived and versatile orthograde body plan and associated adaptations throughout the body that enable very different and very specialized antipronograde behaviors—vertical climbing (all apes), ricochetal brachiation (hylobatids), arboreal quadrumanous suspension and clambering (orangutans), below-branch arm-swinging as well as semiterrestrial knuckle-walking (African apes), and terrestrial bipedalism (humans) (Hunt, 1991; Thorpe and Crompton, 2006).
The fact that, in terms of positional behavior, extant hominoid lineages have more significantly diverged in different directions from their last common ancestor with cercopithecoids explains why extant hominoids more strongly differ in vestibular features—even if the functional link of some vestibular features remains to be better determined. Determining the order in which these features evolved is therefore required to use them for inferring the phylogenetic placement of extinct hominoids. Besides proposing various potential synapomorphies for the hominoid and hominid clades, we further reconstruct the evolution of the vestibular apparatus in this group by estimating ancestral vestibular morphotypes by means of maximum likelihood and a molecular phylogeny. From the LCA of crown hominoids, great apes and humans appear derived in the opposite direction of hylobatids with regard to their volumetric proportions (stout vs. slender SCs, respectively). According to our ancestral state reconstruction for crown hominoids, in this regard hylobatids appear secondarily convergent with cercopithecoids. However, this hypothesis (and the alternate one, that hylobatids merely reflect more closely the primitive condition for crown catarrhines as a whole) should be tested by means of adding extinct stem cercopithecoids and fossil hominoids of less controversial affinities than Oreopithecus to the analysis. It is noteworthy that Hoolock, in agreement with the basal position of this genus among extant hylobatids, apparently retains a more primitive morphology than other hominoids for various features. While the characters related to an anteromedial displacement of the lateral SC appear synapomorphic for hominoids, the only incipient vertical compression of the anterior canal SC morphology of Hoolock might be plesiomorphic for hylobatids, in which case the marked vertical compression of this canal would be a synapomorphy of hominids only, with the remaining hylobatid genera having also evolved it in parallel.
While various great ape lineages and humans further diverged from one another from the more derived condition of the reconstructed crown hominoid LCA, as noted above gibbons and siamangs might have secondarily converged to some extent with cercopithecoids by evolving slenderer SCs, presumably as a result of similar evolutionary pressures posed by fast-moving types of locomotion. Based on our reconstructed ancestral morphotypes, bonobos and (to a lesser extent) chimpanzees would be closer to the LCAs of hominids, hominines, and the Pan-Homo clade, than either humans or the remaining extant great apes (gorillas and orangutans). The latter would have diverged from the hominid LCA in markedly different directions both in terms of SC configuration and stoutness. Our results therefore support the view that not only hominins, but also gorillas and, to a large extent, the orangutan lineage diverged from ancestors with a largely Pan-like vestibular morphology. It would be tempting to interpret this pattern in locomotor terms (e.g., by suggesting a semiterrestrial ancestry not only for hominines, but also for crown hominids as a whole). However, caution is required as other selection pressures and/or non-adaptive factors could have potentially played an equally, if not more significant, role in determining vestibular shape variation in this group.
With regard to the fossil hominoids analyzed here, Australopithecus not only displays the various hominoid synapomorphies mentioned above, but also hominid-like volumetric proportions of the SCs and, as expected in a bipedal hominin, human-like vestibular features such as large anterior and posterior canals. This is in agreement with the large amount of habitual bipedal behaviors inferred for the almost complete skeleton (StW 573) to which one of the analyzed specimens belongs (Heaton et al., 2019), as well as previous analyses of the inner ear as a whole (Beaudet et al., 2019b). Following our ancestral state reconstruction, Australopithecus appears derived in the same direction as humans, although it more closely resembles the morphology of the Australopithecus-Homo clade LCA, which is derived from the reconstructed LCA of the Pan-Homo clade in the opposite direction as chimpanzees and bonobos are. The fact that the two analyzed specimens of Australopithecus are classified as humans with a moderately high probability reflects the fact that their vestibular morphology already approximates the human condition, although maintaining plesiomorphic characters, particularly in the specimen that more closely resembles chimpanzees. On the other hand, their similar classification probability with other extant great ape genera is consistent with a more primitive vestibular morphology. It is noteworthy that, although the two analyzed individuals are very similar to one another (Beaudet et al., 2019b), they display noticeable differences in vestibular morphology, with StW 578 showing more human-like (even if stouter, approximating the orangutan morphology) canals and StW 573 retaining a more African great ape-like morphology. This could be related with diachronic changes within South African Australopithecus or to its previously noted heterogeneity (Clarke, 2013; Grine, 2013), and might help to discern, coupled with other features, the number of species represented among the current samples.
Unlike Australopithecus, Oreopithecus displays a mosaic vestibular morphology that defies a simple phylogenetic interpretation, as it does not fit well among the variation displayed by any extant hominoid genus. This is due to the combination of some hominoid and hominid synapomorphies with more plesiomorphic, cercopithecoid-like, or even platyrrhine-like, features. In particular, the hominid-like volumetric proportions of the Oreopithecus SCs would support the contention that this taxon is a great ape (Begun et al., 1997), or even a member of the European dryopithecine radiation, as previously argued by some authors (Köhler and Moyà-Solà, 1997; Harrison and Rook, 1997). However, other, apparently more plesiomorphic vestibular features, are at odds with such an interpretation, and even with the previous suggestion that Oreopithecus postcranium would be consistent with a giant hylobatid that emphasized cautious climbing (Sarmiento, 1987). In particular, the vestibular morphology of Oreopithecus does not overlap with that of extant hylobatids in any respect, particularly differing by the acute angle between the anterior and posterior canal planes and by the large posterior canal relative to the others. Oreopithecus also appears more primitive than crown hominoids in the shape of the lateral canal, which is flat (rather than displaying an upwards bent ampullar portion) and posteriorly displaced (especially in the junction between the ampulla and the vestibular recesses). In these regard, the Oreopithecus morphology more closely resembles that of cercopithecoids and platyrrhines, respectively, possibly reflecting a plesiomorphic condition that would be more consistent with a stem hominoid status, as recently supported by some other authors (Nengo et al., 2017).
If our interpretation above is correct, then the stout volumetric proportions of Oreopithecus would be homoplastic with those of great apes and humans, representing an independent acquisition that might be functionally related to the evolution of a slower mode of locomotion—in agreement with previous analysis of the inner ear of this taxon (Rook et al., 2004; Ryan et al., 2012) and the possession of an orthograde body plan with adaptations related to cautious vertical climbing and forelimb-dominated suspension (Sarmiento, 1987; Harrison and Rook, 1997). This is plausible given that hylobatids appear to some extent convergent in this regard to cercopithecoids, due to their agile locomotion. This suggests that volumetric proportions are quite labile in evolutionary terms, so that other features and functional considerations must also be considered when interpreting the vestibular morphology of extinct taxa in phylogenetic terms. Interestingly, Oreopithecus resembles australopiths and humans in the possession of a larger vertical and posterior canals relative to the lateral one, which apparently represents another homoplasy that would lend some support to the controversial claim that bipedalism featured prominently among the locomotor repertoire of this taxon (Köhler and Moyà-Solà, 1997; Rook et al., 1999). However, the SCs of Oreopithecus are clearly distinguishable from those present in Homo and Australopithecus regarding volumetric proportions, orientation, and shape. This rules out a hominin-like bipedalism for Oreopithecus—in further agreement with the lack of the lower torso features than in australopiths and humans are functionally linked to committed bipedalism (Hammond et al., 2020)—but would not be at odds with the possession of more varied orthograde positional behaviors combining climbing with a different type of bipedalism (more related to a stable bipedal stance and short distance shuffling instead of fast walking or running), as previously inferred based on the foot of this taxon (Köhler and Moyà-Solà, 1997).
In conclusion, our study provides new insight into the evolution of the vestibular apparatus in hominoids and confirms the potential of SC shape for investigating further the phylogenetic affinities of fossil apes, which are still controversial due to the inherent limitations of the fossil record. This is not to say that functional considerations must not be taken into account—rather the contrary, as several of the discussed vestibular features are arguably linked with the demands of particular positional behaviors, as noted by previous authors (Spoor et al., 2007; Walker et al., 2008; Silcox et al., 2009; Ryan et al., 2012; Perier et al., 2016; Le Maître et al., 2017). However, as exemplified by the analysis of the extinct hominin Australopithecus, the various characters identified as potentially synapomorphic for either crown hominoids or hominids offer the prospect of refining the phylogenetic placement of fossil apes for which their stem vs. crown hominoid status is controversial—as these features can be easily scored from CTs of the petrosal bone and incorporated into formal cladistic analyses including information from other anatomical areas. On the other hand, our ancestral state reconstructions rely mainly on living taxa, which is potentially problematic in the case of hominoids, which were much more diverse in the Miocene and appear quite prone to homoplasy, particularly with regard to the locomotor adaptations of the few surviving lineages. Even if the quantification of phylogenetic signal based on the phylogeny of extant taxa indicates that vestibular morphology overall is not significantly affected by homoplasy, the evolutionary history of vestibular morphology presented here on the basis of ancestral morphotypes should be treated with caution as a set of working hypotheses that require further testing based on the information provided by a larger fossil sample. In particular, given the relationship between vestibular morphology and positional behavior, and the fact that the locomotor apparatus of extinct hominoids frequently displays a mosaic of primitive and derived features unknown among the surviving lineages (e.g., Moyà-Solà et al., 2004; Alba, 2012; Alba et al., 2015; Böhme et al., 2019), it may be predicted that the vestibular morphology of extinct hominoids will similarly display unique combinations of features. This is illustrated here by the condition of Oreopithecus, which is nevertheless most consistent with that of a stem hominoid somewhat convergent with hominids in terms of locomotion. In any case, our conclusions should be subjected to further scrutiny in the future by means of the inclusion of additional fossil taxa, with emphasis on Miocene hominoids as well as stem cercopithecoids.
Materials and methods
Sample composition and acquisition
Request a detailed protocolThe analyzed sample includes microcomputed tomography (µCT) scans of 142 dried anthropoid crania belonging to 27 species and 25 genera, including all extant great ape genera and a selection of hylobatids, cercopithecoids, and platyrrhines (Table 9 and Supplementary file 1). A few specimens are juveniles instead of adults, but this should not affect their vestibular morphology as the bony labyrinth ossifies in early prenatal stages, bounding its shape and size (Jeffery and Spoor, 2004; Perier et al., 2016). The hominoid subsample consists of 48 individuals belonging to 8 species and seven genera (Supplementary file 1). For each specimen, the bony labyrinth was virtually extracted (from the left side when possible, or otherwise from the right side and mirrored) by segmenting the µCT image stacks (voxel size reported in Supplementary file 1). Virtual 3D models were generated using Avizo 9.0.1 software (FEI Visualization Sciences Group). The fossil sample consists of one left bony labyrinth belonging to the late Miocene stem-hominoid Oreopithecus bambolii (Rook et al., 2004) and of two right inner ear that have been virtually extracted from the Australopithecus specimens StW 573 and StW 578 from Sterkfontein (Beaudet et al., 2019b).The vestibular apparatus was separated from the cochlea by cutting right under the oval window and the saccule, and filling the resulting hole with a flat surface in Geomagic Studio 2014 software (3D Systems).
Shape analysis
Request a detailed protocolOur 3DGM approach is based on deformation methods (Durrleman et al., 2012b; Durrleman et al., 2012a), which do not rely on a priori defined landmarks but consider instead the geometrical correspondences between continuous surfaces, and are particularly convenient for comparing overall shape and complex 3D surface changes (Durrleman et al., 2012b; Dumoncel et al., 2014; Beaudet et al., 2016a). This method relies on the construction of a sample-average surface model (template) and its deformation to the investigated surfaces (Durrleman et al., 2012b; Durrleman et al., 2012a; Beaudet et al., 2016b). Unlike in classical landmark-based 3DGM analyses, the surfaces are represented by a set of oriented faces and the comparisons do not assume a point-to-point correspondence between samples (Durrleman et al., 2012b). Prior to the analysis, the unscaled vestibular surfaces were aligned and scaled using the ‘Align Surface’ module of Avizo 9.0. Then, deformations between surfaces were mathematically modeled as a diffeomorphism (i.e., a one-to-one deformation of the 3D space that is smooth, invertible, and with a smooth inverse), and a set of momenta (vectors representing the flow of deformations from the initial position of the control points on the template to the target shape) were estimated with Deformetrica 3 software. Due to its high-demanding computational power, analyses were run in the CALMIP supercomputing center (Toulouse, France).
We inspected interspecific major patterns of shape variation by means of between-group principal components analysis (bgPCA) of the deformation-based shape residuals, using major clades (i.e., platyrrhines, cercopithecoids, hylobatids and hominids) as grouping factor. The restricted platyrrhine sample included in this study is aimed to serve as an outgroup for catarrhines, since the description of the vestibular morphology variation among New World monkeys as a whole is beyond the scope of this paper. Each group has been designed to include a large number of individuals (>>10) in order to prevent spurious separation between the groups used in the analysis. We preferred bgPCA over linear discriminant analysis (LDA) because the latter produces overexaggerated separation among groups when the number of variables is close to the number of the analyzed individuals (Mitteroecker and Bookstein, 2011). Taking into account that some authors have recently recommended caution when interpreting bgPCA results, as they might present spurious grouping (Cardini et al., 2019), we compared our results with those of a landmark-based 3DGM analysis ran on the same sample and investigated the presence of a preexisting group structure (See ‘Exploration of a preexisting group structure in the tangent space of the vestibular shape’ section for further details).
Correlation between SC shape and size (allometry) was assessed using multivariate regression of the deformation fields against log-transformed cube root of the entire vestibular volume (ln Vol, in mm3), as well as bivariate regressions computed for each bgPC against ln Vol. We further inspected the correlation between log-transformed cube root of the SC volume (ln VolSC, in mm3; including the SCs and the CC cut at their connection with the vestibular recesses) and log-transformed length (ln L, in mm; measured along the streamline of the SCs and of the CC). All these regressions were performed for the whole anthropoid sample, as well as for hominids and non-hominids separately, as we detected several differences for the former. In the case of the ln VolSC vs. ln L regression, we also checked the homogeneity of slopes and intercepts between hominoids and the rest of the sample via analysis of covariance (ANCOVA). Regressions were performed by means of ordinary least-square linear regression (OLS) as well as by taking into account the phylogenetic non-independence between data, that is by fitting our linear models via phylogenetic generalized least squares (PGLS). The data used to compute the regressions are given in Figure 7—source data 1.
The hierarchical clustering analysis (HCA) of the deformation fields (used to identify preexisting group structure embedded in our shape data) has been computed caret v6.0–84 (Kuhn, 2008) and Factominer v1.34 (Le et al., 2008) R packages. The discrimination and the amount of overlap among the groups defined a priori for the bgPCA has been quantified by computing the number of correctly classified individuals with cross-validation using the Morpho v2.6 (Schlager, 2014) R package. Finally, we estimated the posterior probabilities of group membership for the fossil specimens based on the Mahalanobis distances between their bgPC scores and group centroids using the ‘typprobClass’ function of Morpho v2.6 (Schlager, 2014) R package. These probabilities, which were computed for both bgPCA (based on main anthropoid groups as well as on extant hominoid genera only) denote the likelihood of the specimens to belong to each group without assuming that it must belong to one of them (which is required when comparing the fossils with extant genera), so that the sum of probabilities for each specimen does not equal one. Probabilities < 0.05 indicate that the specimen falls outside the variability of that particular group. Further statistical analyses were carried out using different R packages in RStudio v.1.1.453 for R v.3.5.0: ape v5.1 (Paradis, 2012), phytools v0.6–60 (Revell, 2012), Morpho v2.6 (Schlager, 2014), caper v1.0.2 (Orme et al., 2013), and geomorph v3.1.1 (Adams et al., 2019).
Phylogenetic signal analysis
Request a detailed protocolTo assess phylogenetic signal, that is the degree to which related species resemble each other (Felsenstein, 1985; Harvey and Pagel, 1991), we used a phylogenetic tree (Figure 13) derived from a time-calibrated molecular phylogeny for the extant taxa (Springer et al., 2012). Estimated mean divergence dates for the various included extant clades are indicated in Figure 13 (see Table 1 and S1 in Springer et al. (2012) for 95% composite credibility intervals). The phylogenetic placement of Oreopithecus has been controversial; for our purposes here, we followed Nengo et al. (2017) in considering this taxon as a stem-hominoid, although other possibilities are discussed in the text. The node of the Oreopithecus-crown hominoid divergence has been placed 1 Ma older than the divergence of crown apes and humans and its tip corresponds to its last occurrence in the fossil record (7.0–6.5 Ma; Rook et al., 2000). For the South African Australopithecus sp., we used the published first appearance datum for Australopithecus africanus (4.02 Ma) that includes the Jacovec specimens into the taxon (Wood and Boyle, 2016).
We computed Pagel’s λ (Pagel, 1999), Blomberg’s K (Blomberg et al., 2003) and Kmult (Adams, 2014) using phytools v.0.6–44 (Revell, 2012) and geomorph v3.1.1 (Adams et al., 2019) R packages. These metrics compare the variance in the phylogenetic tree tips relative to that expected under a Brownian motion evolutionary model. Pagel’s λ is a scaling coefficient of the expected correlations between related species on the tree, and varies from 0 (no correlation due to absence of phylogenetic signal) to 1 (covariance proportional to phylogenetic distance, implying maximal phylogenetic signal). Blomberg’s K and its multivariate generalization Kmult inform on how precisely the variance-covariance patterns found in the data are matched by the phylogenetic tree and where variance accumulates: K ≈ one implies that the mode of evolution closely resembles that expected under Brownian motion; when K < 1, close relatives resemble each other less than expected (variance is accumulated within clades), implying an evolutionary pattern that departs from a purely stochastic model (which could be caused by adaptation uncorrelated with phylogeny, that is homoplasy); finally, K > 1 is found when close relatives are more similar than expected under Brownian motion (variance is among clades), which could indicate stabilizing selection.
Phylomorphospace and ancestral states analysis
Request a detailed protocolTo quantify major patterns of vestibular shape variation along the branches of the phylogeny we relied on a phylomorphospace approach (Sidlauskas, 2008), which allows us to intuitively visualize the extent and direction of the inferred shape change by means of branch length and orientation. This method projects the phylogenetic tree (Figure 13) onto the tangent space computed from the bgPCA and estimates the position in the morphospace of the internal nodes (ancestral morphologies) via a maximum likelihood (ML) method for continuous characters (Felsenstein, 1988; Schluter et al., 1997) using the ‘fastAnc’ function of phytools v0.6–60 R package (Revell, 2012). Subsequently, the bgPC scores of the ancestral states are rotated and translated from the shape data back into the configuration space for interpolation and 3D visualization using Deformetrica 3 software.
Data availability
All the results generated during this study are included as in the manuscript and as supporting files. Source data files are provided for Figures 2, 7, 8, and 9.
-
MorphoSourceA collection of non-human primate computed tomography scans housed in MorphoSource, a repository for 3D data.
References
-
Locomotor inferences in Pierolapithecus and Hispanopithecus: reply to Deane and begun (2008)Journal of Human Evolution 59:143–149.https://doi.org/10.1016/j.jhevol.2010.02.002
-
Fossil apes from the Vallès-Penedès basinEvolutionary Anthropology 21:254–269.https://doi.org/10.1002/evan.21312
-
Morphoarchitectural variation in South African fossil cercopithecoid endocastsJournal of Human Evolution 101:65–78.https://doi.org/10.1016/j.jhevol.2016.09.003
-
The inner ear of the Paranthropus specimen DNH 22 from Drimolen, South AfricaAmerican Journal of Physical Anthropology 170:439–446.https://doi.org/10.1002/ajpa.23901
-
Function, Phylogeny and Fossils: Miocene Hominoid Evolution and Adaptation389–415, Events in Hominoid Evolution, Function, Phylogeny and Fossils: Miocene Hominoid Evolution and Adaptation, New York, Plenum Press, 10.1007/978-1-4899-0075-3_18.
-
A Companion to Paleoanthropology398–416, The Miocene Hominoid Radiations, A Companion to Paleoanthropology, Oxford, Blackwell, 10.1002/9781118332344.ch21.
-
Handbook of Paleoanthropology1261–1332, Fossil Record of Miocene Hominoids, Handbook of Paleoanthropology, London, Springer, 10.1007/978-3-642-39979-4_32.
-
Comparative analysis of vestibular ecomorphology in birdsJournal of Anatomy 231:990–1018.https://doi.org/10.1111/joa.12726
-
Semicircular canal size and shape influence on disorientationAerospace Medicine and Human Performance 89:744–748.https://doi.org/10.3357/AMHP.5104.2018
-
The Paleobiology Of AustralopithecusAustralopithecus from Sterkfontein Caves, South Africa, The Paleobiology Of Australopithecus, Dordrecht, Springer, 10.1007/978-94-007-5919-0_7.
-
The bony labyrinth in the Aroeira 3 Middle Pleistocene craniumJournal of Human Evolution 124:105–116.https://doi.org/10.1016/j.jhevol.2018.08.003
-
The dentinoenamel junction in primatesInternational Journal of Primatology 8:99–114.https://doi.org/10.1007/BF02735159
-
Motion from the past. A new method to infer vestibular capacities of extinct speciesComptes Rendus Palevol 9:397–410.https://doi.org/10.1016/j.crpv.2010.07.012
-
Landmark-free 3D method for comparison of fossil hominins and hominids based on endocranium and EDJ shapesAmerican Journal of Physical Anthropology 153.
-
BookTopology Preserving Atlas Construction from Shape Data without Correspondence Using Sparse ParametersIn: Ayache N, Delingette H, Golland P, Mori K, editors. Proceedings of Medical Image Computing and Computer Aided Intervention. Berlin: Springer-Verlag. pp. 223–230.https://doi.org/10.1007/978-3-642-33454-2_28
-
Phylogenies and the comparative methodThe American Naturalist 125:1–15.https://doi.org/10.1086/284325
-
Phylogenies and quantitative charactersAnnual Review of Ecology and Systematics 19:445–471.https://doi.org/10.1146/annurev.es.19.110188.002305
-
The subarcuate Fossa and cerebellum of extant primates: comparative study of a skull-brain interfaceAmerican Journal of Physical Anthropology 77:143–164.https://doi.org/10.1002/ajpa.1330770202
-
BookTemplate Estimation Form Unlabeled Point Set Data and Surfaces for Computational AnatomyIn: Pennec X, editors. Proceedings of the International Workshop on the Mathematical Foundations of Computational Anatomy Sophia-Antipolis. INRIA. pp. 28–39.
-
BookThe Evolution of Cercopithecoid Locomotion: A Morphometric, Phylogenetic, and Character Mapping Approach. Ph.D DissertationStony Brook University.
-
BookThe Alpha Taxonomy of Australopithecus africanusIn: Delson E, Sargis EJ, editors. The Paleobiology of Australopithecus. Dordrecht: Springer. pp. 73–104.https://doi.org/10.1007/978-94-007-5919-0_6
-
Bony labyrinth shape variation in extant Carnivora: a case study of MusteloideaJournal of Anatomy 228:366–383.https://doi.org/10.1111/joa.12421
-
BookCatarrhine OriginsIn: Begun D. R, editors. A Companion to Paleoanthropology. Oxford: Blackwell Publishing. pp. 376–396.https://doi.org/10.1002/9781118332344
-
BookThe Fossil Record and Evolutionary History of Hylobatids. Evolution of Gibbons and SiamangIn: Reichard UH, Hirai H, Barelli C, editors. Phylogeny, Morphology, and Cognition. London: Springer. pp. 90–110.https://doi.org/10.1007/978-1-4939-5614-2_4
-
BookEnigmatic Anthropoid or Misunderstood Ape?In: Begun DR, Ward CV, editors. Function, Phylogeny, and Fossils. Boston: Springer. pp. 327–362.https://doi.org/10.1007/978-1-4899-0075-3_16
-
Positional behavior in the HominoideaInternational Journal of Primatology 12:95–118.https://doi.org/10.1007/BF02547576
-
BookCercopithecoideaIn: Werdelin L, Sanders WL, editors. Cenozoic Mammals of Africa. London: University of California Press. pp. 393–429.https://doi.org/10.1525/california/9780520257214.003.0023
-
Prenatal growth and development of the modern human labyrinthJournal of Anatomy 204:71–92.https://doi.org/10.1111/j.1469-7580.2004.00250.x
-
Analysis of vestibular labyrinthine geometry and variation in the human temporal boneFrontiers in Neuroscience 12:107.https://doi.org/10.3389/fnins.2018.00107
-
Building predictive models in R using the caret packageJournal of Statistical Software 28:1–26.https://doi.org/10.18637/jss.v028.i05
-
Parallel evolution in the hominoid trunk and forelimbEvolutionary Anthropology 6:87–99.https://doi.org/10.1002/(SICI)1520-6505(1998)6:3<87::AID-EVAN3>3.0.CO;2-T
-
FactoMineR: an R package for multivariate analysisJournal of Statistical Software 25:.https://doi.org/10.18637/jss.v025.i01
-
New data about semicircular canal morphology and locomotion in modern hominoidsJournal of Anatomy 231:95–109.https://doi.org/10.1111/joa.12619
-
Deep evolutionary roots of strepsirrhine primate labyrinthine morphologyJournal of Anatomy 216:368–380.https://doi.org/10.1111/j.1469-7580.2009.01177.x
-
The labyrinthine morphology of Pronycticebus gaudryi (Primates, Adapiformes)Palaeobiodiversity and Palaeoenvironments 92:527–537.https://doi.org/10.1007/s12549-012-0099-z
-
The phylogenetic relationships of Oreopithecus bambolii Gervais, 1872Comptes Rendus de l'Académie des Sciences Paris 324:141–148.
-
Semicircular duct dimensions and sensitivity of the vertebrate vestibular systemJournal of Theoretical Biology 167:239–256.https://doi.org/10.1006/jtbi.1994.1066
-
BookAnalysis of Phylogenetics and Evolution with RNew York: Springer.https://doi.org/10.1007/978-1-4614-1743-9
-
Different level of intraspecific variation of the bony labyrinth morphology in slow- versus fast-moving primatesJournal of Mammalian Evolution 23:353–368.https://doi.org/10.1007/s10914-016-9323-3
-
BookBiomechanics of the Semicircular Canals and Otolith OrgansIn: Highstein SM, Fay RR, Popper AN, editors. The Vestibular System. New York: Springer. pp. 153–201.https://doi.org/10.1007/0-387-21567-0_4
-
Semicircular canal size and locomotion in colobine monkeys: a cautionary taleFolia Primatologica 87:213–223.https://doi.org/10.1159/000449286
-
Dimensional analysis and dynamic response characterization of mammalian peripheral vestibular structuresAmerican Journal of Anatomy 169:295–313.https://doi.org/10.1002/aja.1001690306
-
Phytools: an R package for phylogenetic comparative biology (and other things)Methods in Ecology and Evolution 3:217–223.https://doi.org/10.1111/j.2041-210X.2011.00169.x
-
The bony labyrinth of Oreopithecus bamboliiJournal of Human Evolution 46:347–354.https://doi.org/10.1016/j.jhevol.2004.01.001
-
Evolution of locomotion in Anthropoidea: the semicircular canal evidenceProceedings of the Royal Society B: Biological Sciences 279:3467–3475.https://doi.org/10.1098/rspb.2012.0939
-
The phylogenetic position of Oreopithecus and its significance in the origin of the HominoideaAmerican Museum Novitates 2881:1–44.
-
Morpho: calculations and visualisations related to geometric morphometricsMorpho: calculations and visualisations related to geometric morphometrics, https://github.com/zarquon42b/Morpho.
-
Semicircular canal system in early primatesJournal of Human Evolution 56:315–327.https://doi.org/10.1016/j.jhevol.2008.10.007
-
Dental trait expression at the enamel-dentine junction of lower molars in extant and fossil hominoidsJournal of Human Evolution 54:173–186.https://doi.org/10.1016/j.jhevol.2007.09.012
-
The bony labyrinth of NeanderthalsJournal of Human Evolution 44:141–165.https://doi.org/10.1016/S0047-2484(02)00166-5
-
Absence of the subarcuate fossa in cercopithecidsJournal of Human Evolution 31:569–575.https://doi.org/10.1006/jhev.1996.0081
-
Comparative review of the human bony labyrinthAmerican Journal of Physical Anthropology 107:211–251.https://doi.org/10.1002/(SICI)1096-8644(1998)107:27+<211::AID-AJPA8>3.0.CO;2-V
-
Orangutan positional behavior and the nature of arboreal locomotion in JominoideaAmerican Journal of Physical Anthropology 131:384–401.https://doi.org/10.1002/ajpa.20422
-
The semicircular canal system and locomotion: the case of extinct lemuroids and lorisoidsEvolutionary Anthropology 17:135–145.https://doi.org/10.1002/evan.20165
-
Hominin taxic diversity: fact or fantasy?American Journal of Physical Anthropology 159:37–78.https://doi.org/10.1002/ajpa.22902
Article and author information
Author details
Funding
Agencia Estatal de Investigación (CGL2016-76431-P)
- David M Alba
Agencia Estatal de Investigación (CGL2017-82654-P)
- Alessandro Urciuoli
- Salvador Moyà-Solà
Generalitat de Catalunya (2017 SGR 86)
- Alessandro Urciuoli
- Salvador Moyà-Solà
Generalitat de Catalunya (2017 SGR 116)
- David M Alba
Agencia Estatal de Investigación (BES-2015–071318)
- Alessandro Urciuoli
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
Access to the HPC resources of CALMIP supercomputing center were granted under the allocation 2016-[P1440] attributed to the AMIS Laboratory, Toulouse (France). This work has been funded by the Agencia Estatal de Investigación and European Regional Development Fund of the European Union (projects CGL2016-76431-P and CGL2017-82654-P, AEI/FEDER, EU), by the Ministerio de Economía y Empresa (grant BES-2015–071318 to AU), the Generalitat de Catalunya (CERCA Programme, and consolidated research groups 2017 SGR 86 and 2017 SGR 116 GRC), and the French CNRS. We deeply thank Sergio Almécija for helpful advice and discussion during the elaboration of this paper and José Braga for valuable suggestions. We are grateful to Lorenzo Rook for giving access to the Oreopithecus bambolii petrosal bone originally published in Rook et al. (2004). Lynn Copes, Lynn Lucas, and the MCZ provided access to a part of the scans used in the study, originally appearing in Copes et al. (2016), funding for the collection of which was provided by NSF DDIG #0925793, and a Wenner-Gren Foundation Dissertation Grant #8102 (both to Lynn Copes). The Evolutionary Studies Institute provided access to the Australopithecus sp. data originally appearing in Beaudet (2019a), the collection of which was funded by PAST, DST-NRF, CoE-Pal, Claude Leon Foundation, IFAS, French embassy, Leakey Foundation, Wits University. The scans for fossil and extant species were downloaded from MorphoSource.org, a web-accessible archive for 3D digital data housed by Duke University.
Copyright
© 2020, Urciuoli et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 5,390
- views
-
- 549
- downloads
-
- 24
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Evolutionary Biology
Eyespot patterns have evolved in many prey species. These patterns were traditionally explained by the eye mimicry hypothesis, which proposes that eyespots resembling vertebrate eyes function as predator avoidance. However, it is possible that eyespots do not mimic eyes: according to the conspicuousness hypothesis, eyespots are just one form of vivid signals where only conspicuousness matters. They might work simply through neophobia or unfamiliarity, without necessarily implying aposematism or the unprofitability to potential predators. To test these hypotheses and explore factors influencing predators’ responses, we conducted a meta-analysis with 33 empirical papers that focused on bird responses to both real lepidopterans and artificial targets with conspicuous patterns (i.e. eyespots and non-eyespots). Supporting the latter hypothesis, the results showed no clear difference in predator avoidance efficacy between eyespots and non-eyespots. When comparing geometric pattern characteristics, bigger pattern sizes and smaller numbers of patterns were more effective in preventing avian predation. This finding indicates that single concentric patterns have stronger deterring effects than paired ones. Taken together, our study supports the conspicuousness hypothesis more than the eye mimicry hypothesis. Due to the number and species coverage of published studies so far, the generalisability of our conclusion may be limited. The findings highlight that pattern conspicuousness is key to eliciting avian avoidance responses, shedding a different light on this classic example of signal evolution.
-
- Evolutionary Biology
The rise of angiosperms to ecological dominance and the breakup of Gondwana during the Mesozoic marked major transitions in the evolutionary history of insect-plant interactions. To elucidate how contemporary trophic interactions were influenced by host plant shifts and palaeogeographical events, we integrated molecular data with information from the fossil record to construct a time tree for ancient phytophagous weevils of the beetle family Belidae. Our analyses indicate that crown-group Belidae originated approximately 138 Ma ago in Gondwana, associated with Pinopsida (conifer) host plants, with larvae likely developing in dead/decaying branches. Belids tracked their host plants as major plate movements occurred during Gondwana’s breakup, surviving on distant, disjunct landmasses. Some belids shifted to Angiospermae and Cycadopsida when and where conifers declined, evolving new trophic interactions, including brood-pollination mutualisms with cycads and associations with achlorophyllous parasitic angiosperms. Extant radiations of belids in the genera Rhinotia (Australian region) and Proterhinus (Hawaiian Islands) have relatively recent origins.