Patterns of cranial shape variation in the Papionini (Primates: Cercopithecinae)

Allometry and phylogenetic divergence: Correspondence or incongruence?

The potential connection between trends of within species variation, such as those of allometric change in morphology, and phylogenetic divergence has been a central topic in evolutionary biology for more than a century, including in the context of human evolution. In this study, I focus on size-related shape change in craniofacial proportions using a sample of more than 3200 adult Old World monkeys belonging to 78 species, of which 2942 specimens of 51 species are selected for the analysis. Using geometric morphometrics, I assess whether the divergence in the direction of static allometries increases in relation to phyletic differences. Because both small samples and taxonomic sampling may bias the results, I explore the sensitivity of the main analyses to the inclusion of more or less taxa depending on the choice of a threshold for the minimum sample size of a species. To better understand the impact of sampling error, I also use randomized subsampling experiments in the largest species samples. The study shows that static allometries vary broadly in directions without any evident phylogenetic signal. This variation is much larger than previously found in ontogenetic trajectories of Old World monkeys, but the conclusion of no congruence with phylogenetic divergence is the same. Yet, the effect of sampling error clearly contributes to inaccuracies and tends to magnify the differences in allometric change. Thus, morphometric research at the boundary between micro- and macro-evolution in primates, and more generally in mammals, critically needs very large and representative samples. Besides sampling error, I suggest other non-mutually exclusive explanations for the lack of correspondence between allometric and phylogenetic divergence in Old World monkeys, and also discuss why directions might be more variable in static compared to ontogenetic trajectories. Even if allometric variation may be a poor source of information in relation to phylogeny, the evolution of allometry is a fascinating subject and the study of size-related shape changes remains a fundamental piece of the puzzle to understand morphological variation within and between species in primates and other animals.

Description and taxonomic assessment of fossil Cercopithecidae from the Pliocene Galili Formation (Ethiopia)

The Mount Galili Formation in the Afar region, Ethiopia, samples a critical time in hominin evolution, 4.4 to 3.8 Ma, documenting the last appearance of Ardipithecus and the origin of Australopithecus. This period is also important in the evolution of cercopithecids, especially the origin of Theropithecus in general and Theropithecus oswaldi lineage in particular. Galili has provided a total of 655 cercopithecid specimens that include crania, mandibles, isolated teeth and postcrania. All the fossils were recovered from the Lasdanan (5.3-4.43 Ma), Dhidinley (4.43-3.9 Ma) and Shabeley Laag (∼3.92-3.8 Ma) Members. Here, we described and analyzed 362 fossils employing both qualitative and quantitative methods. Descriptions of the material were supplemented with dental metrics and cranial shape analysis using three-dimensional geometric morphometrics. Results indicate the presence of at least six cercopithecid taxa: Theropithecus oswaldi serengetensis (n = 28), Theropithecus sp. (n = 2), three non-Theropithecus papionin groups (n = 134) and one colobine-size group (n = 58). The T. o. serengetensis represents the earliest form of the lineage, documented from ∼3.9 Ma Galili sediments. The three Galili papionins include a smaller taxon, a medium-sized taxon comparable to Pliopapio alemui and a large papionin overlapping in size with Soromandrillus, Gorgopithecus and Dinopithecus. The majority of Galili colobines have closest affinities to Kuseracolobus aramisi and some overlap with other taxa. Papionins dominate the Galili cercopithecid collection, although colobines are still fairly common (approximately 25% of the sample). Thus, Galili sample is like Kanapoi (4.2-4.1 Ma) and Gona (5.2-3.9 Ma) localities but distinct from Aramis, suggesting paleoecological similarity to the former sites. On the other hand, Theropithecus is less abundant at Galili than geologically younger Hadar (3.4-3.2 Ma) and Woranso-Mille (3.8-3.6 Ma) sites. Whether this difference is due to sampling, time or landscape variation requires further investigation.

Pedomorphosis in the ancestry of marsupial mammals

Within mammals, different reproductive strategies (e.g., egg laying, live birth of extremely underdeveloped young, and live birth of well-developed young) have been linked to divergent evolutionary histories. How and when developmental variation across mammals arose is unclear. While egg laying is unquestionably considered the ancestral state for all mammals, many long-standing biases treat the extreme underdeveloped state of marsupial young as the ancestral state for therian mammals (clade including both marsupials and placentals), with the well-developed young of placentals often considered the derived mode of development. Here, we quantify mammalian cranial morphological development and estimate ancestral patterns of cranial shape development using geometric morphometric analysis of the largest comparative ontogenetic dataset of mammals to date (165 specimens, 22 species). We identify a conserved region of cranial morphospace for fetal specimens, after which cranial morphology diversified through ontogeny in a cone-shaped pattern. This cone-shaped pattern of development distinctively reflected the upper half of the developmental hourglass model. Moreover, cranial morphological variation was found to be significantly associated with the level of development (position on the altricial-precocial spectrum) exhibited at birth. Estimation of ancestral state allometry (size-related shape change) reconstructs marsupials as pedomorphic relative to the ancestral therian mammal. In contrast, the estimated allometries for the ancestral placental and ancestral therian were indistinguishable. Thus, from our results, we hypothesize that placental mammal cranial development most closely reflects that of the ancestral therian mammal, while marsupial cranial development represents a more derived mode of mammalian development, in stark contrast to many interpretations of mammalian evolution.

Morphology of the Bony Labyrinth Supports the Affinities of Paradolichopithecus with the Papionina

Discoveries in recent decades indicate that the large papionin monkeys Paradolipopithecus and Procynocephalus are key members of the Late Pliocene - Early Pleistocene mammalian faunas of Eurasia. However, their taxonomical status, phylogenetic relationships, and ecological profile remain unclear. Here we investigate the two latter aspects through the study of the inner ear anatomy, as revealed by applying micro-CT scan imaging techniques on the cranium LGPUT DFN3-150 of Paradolichopithecus from the lower Pleistocene (2.3 Ma) fossil site Dafnero-3 in Northwestern Greece. Using geometric morphometric methods, we quantified shape variation and the allometric and phylogenetic signals in extant cercopithecines (n = 80), and explored the morphological affinities of the fossil specimen with extant taxa. LGPUT DFN3-150 has a large centroid size similar to that of baboons and their relatives. It shares several shape features with Macacina and Cercopithecini, which we interpret as probable retention of a primitive morphology. Overall, its inner ear morphology is more consistent with a stem Papionini more closely related to Papionina than Macacina, or to a basal crown Papionina. Our results, along with morphometrical and ecological features from previous studies, call into question the traditional hypothesis of a Paradolichopithecus-Macacina clade, and provide alternative perspectives in the study of Eurasian primate evolution during the late Neogene-Quaternary.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10764-022-00329-4.

Divergent Genital Morphologies and Female-Male Covariation in Watersnakes

Genital evolution can be driven by diverse selective pressures. Across taxa we see evidence of covariation between males and females, as well as divergent genital morphologies between closely related species. Quantitative analyses of morphological changes in coevolving male and female genitalia have not yet been shown in vertebrates. This study uses 2D and 3D geometric morphometrics to quantitatively compare the complex shapes of vaginal pouches and hemipenes across three species of watersnakes (the sister taxa Nerodia fasciata, N. sipedon, and a close relative N. rhombifer) to address the relationship between genital morphology and divergence time in a system where sexual conflict may have driven sexually antagonistic coevolution of genital traits. Our pairwise comparisons of shape differences across species show that the sister species have male and female genitalia that are significantly different from each other, but more similar to each other than to N. rhombifer. We also determine that the main axes of shape variation are the same for males and females, with changes that relate to deeper bilobation of the vaginal pouch and hemipenes. In males, the protrusion of the region of spines at the base of the hemipene trades off with the degree of bilobation, suggesting amelioration of sexual conflict, perhaps driven by changes in the relative size of the entrance of the vaginal pouch that could have made spines less effective.

Methods for studying allometry in geometric morphometrics: a comparison of performance

Allometry has been the focus of growing interest in studies using geometric morphometric methods to address a wide range of research questions at the interface of ecology and evolution. This study uses computer simulations to compare four methods for estimating allometric vectors from landmark data: the multivariate regression of shape on a measure of size, the first principal component (PC1) of shape, the PC1 in conformation space, and a recently proposed method, the PC1 of Boas coordinates. Simulations with no residual variation around the allometric relationship showed that all four methods are logically consistent with one another, up to minor nonlinearities in the mapping between conformation space and shape tangent space. In simulations that included residual variation, either isotropic or with a pattern independent of allometry, regression of shape on size performed consistently better than the PC1 of shape. The PC1s of conformation and of Boas coordinates were very similar and very close to the simulated allometric vectors under all conditions. An extra series of simulations to elucidate the relation between conformation and Boas coordinates indicated that they are almost identical, with a marginal advantage for conformation. Empirical examples of ontogenetic allometry in rat skulls and rockfish body shape illustrate simple biological applications of the methods. The paper concludes with recommendations how these methods for estimating allometry can be used in studies of evolution and ecology.

Biomechanical correlates of zygomaxillary-surface shape in papionin primates and the effects of hard-object feeding on mangabey facial form

Extant African papioninans are distinguished from macaques by the presence of excavated facial fossae; however, facial excavation differs among taxa. Mangabeys (Cercocebus, Rungwecebus, and Lophocebus) exhibit fossae that invade the zygomatic forming pronounced suborbital fossae (SOFs). Larger-bodied Papio, Mandrillus, and Theropithecus have lateral rostral fossae with minimal/absent suborbital fossae. Because prior studies have shown that mangabeys exhibit adaptations to anterior dental loading (e.g., palatal retraction), it is plausible that mangabey SOFs represent structural accommodation to masticatory-system shape rather than facial allometry, as commonly hypothesized. We analyzed covariation between zygomaxillary-surface shape, masticatory-system shape, and facial size in 141 adult crania of Macaca fascicularis, Papio kindae, Cercocebus, and Lophocebus. These taxa represent the range of papionin SOF expression while minimizing size variation (narrow allometry). Masticatory-system landmarks (39) registered palate shape, bite points, masticatory muscle attachments, and the temporomandibular joint. Semilandmarks (450) captured zygomaxillary-surface shape. Following Procrustes superimposition with semilandmark sliding and principal components analyses, multivariate regression was used to explore allometry, and two-block partial least-squares analyses (within-configuration and separate-blocks) were used to examine covariation patterns. Scores on principal components 1-2 and the first partial least-square (PLS1) separate mangabeys from Macaca and Papio. Both zygomaxillary-surface shape and masticatory-system shape are correlated with size within taxa and facial morphotypes; however, regression distributions indicate morphotype shape differences are non-allometric. PLS1 accounts for ∼95% of shape covariance (p < 0.0001) and shows strong linear correlations (r-PLS = ∼0.95, p < 0.0001) between blocks. Negative PLS1 scores in mangabeys reflect deep excavation of the suborbital malar surface, palatal retraction, and anterior displacement of jaw adductor muscles and the temporomandibular joint. Neither PC1 nor PLS1 scores ordinate specimens by facial size. Taken together, these results fail to support the allometric hypothesis but suggest that mangabey zygomaxillary morphology is closely linked with adaptations to hard-object feeding.

Evolutionary drivers of morphological differentiation among three bottlenose dolphin lineages, Tursiops spp. (Delphinidae), in the northwest Indian Ocean utilizing linear and geometric morphometric techniques

Local adaptation and adaptive radiations are typically associated with phenotypic variation suited to alternative environments. In the marine environment, the nature of relevant ecological or environmental transitions is poorly understood, especially for highly mobile species. Here we compare three genetic lineages in the genus Tursiops (bottlenose dolphins), using linear measurements and geometric morphometric techniques, in the context of environmental variation in the northwest Indian Ocean. Cranial morphology was clearly differentiated comparing Tursiops truncatus and Tursiops aduncus, while a recently discovered genetic lineage, found in the Arabian Sea, was morphologically most similar to T. aduncus from the same region, but distinct for various measures, particularly metrics associated with the lateral dimension of the skull. The extent of divergence between T. truncatus and T. aduncus compared to differences between the T. aduncus lineages is consistent with the recent phylogeny for these species. Therefore, with the corroboration of genetic and morphological inference, we propose two conservation units of T. aduncus be recognized in the region at a sub-specific level so that their conservation can be managed effectively. We consider possible evolutionary mechanisms associated with regional habitat characteristics and the exploitation of distinct prey resources.

Accurate and semi-automated reassociation of intermixed human skeletal remains recovered from bioarchaeological and forensic contexts

Commingled remains describes the situation of intermixed skeletal elements, an extremely common occurrence in contemporary forensic cases, archaeological mass graves, as well as fossil hominin assemblages. Given that reliable identification is typically impossible for commingled contexts, a plethora of previous studies has focused on the development of refined methods for reassociating the bones of each individual skeleton. Here, a novel virtual approach for quantifying the degree of three-dimensional shape compatibility between two adjoining bone articular surfaces is put forth. Additionally, the integrability of this method with traditional osteometric techniques is evaluated. We focus on the paradigm of the hip joint, whose articulating bone elements (the femur and the innominate bone) are crucial for reconstructing the biological profile of unidentified human remains. The results demonstrate that this new semi-automated methodology is highly accurate both for large commingled assemblages (such as those resulting from mass disasters or burials) as well as smaller-scale contexts (such as those resulting from secondary burials).

Developmental instability in wild Nigerian olive baboons (Papio anubis)

Background

Developmental instability in archaeological samples can be detected through analysis of skeletal and dental remains. During life, disruptions to biological internal homeostasis that occur during growth and development redirect bodily resources to returning to homeostasis and away from normal processes such as symmetrical development. Because dental enamel does not remodel in life, any deviations from normal development are left behind. Even subtle disturbances to developmental trajectory may be detected in asymmetrical development of traits, specifically a random variation in sides termed fluctuating asymmetry. Human dental fluctuating asymmetry studies are common, but here we investigate the permanent dentition of a non-human primate Papio anubis, for potential fluctuating asymmetry relative to sex, weaning, and reproductive maturity. The sample stems from an outlier population that lives in the wettest and most humid habitat of any studied baboon group.

Methods

The skulls of adult baboons were collected after their natural death in Gashaka Gumti National Park, Nigeria. The permanent dentition of antimeric teeth (paired) were measured for maximum length and breadth using standard methods. The metrics were analyzed to assess the presence of fluctuating asymmetry in adult permanent mandibular and maxillary dentition. Measurement error and other forms of asymmetry (antisymmetry, directional asymmetry) were considered and dental measures expressing true fluctuating asymmetry were used to address three research questions.

Results

Males exhibit greater fluctuating asymmetry than females, suggesting that males experience greater overall instability during the developmental period. While weaning is not more stressful than other life history stages for males and females (using the first molar fluctuating asymmetry index as a proxy compared to other teeth), it is more stressful for females than males. The onset of reproduction is also not more stressful than other life history stages for males and females (using the third molar fluctuating asymmetry index as a proxy compared to other teeth), but it is more stressful for males than females. We explore possible explanations for these findings in the discussion.

Ontogenetic allometry and scaling in catarrhine crania

In studies of ontogenetic allometry, ontogenetic scaling has often been invoked to explain cranial morphological differences between smaller and larger forms of closely related taxa. These scaled variants in shape have been hypothesized to be the result of the extension or truncation of common growth allometries. In this scenario, change in size is the determining factor, perhaps under direct selection, and changes in cranial shapes are byproducts, not under direct selection themselves. However, many of these conclusions are based on studies that used bivariate generalizations of shape. Even among multivariate analyses of growth allometries, there are discrepancies as to the prevalence of ontogenetic scaling among primates, how shared the trajectories need to be, and which taxa evince properties of scaled variants. In this investigation, we use a large, comparative ontogenetic sample, geometric morphometric methods, and multivariate statistical tests to examine ontogenetic allometry and evaluate if differences in cranial shape among closely related catarrhines of varying sizes are primarily driven by size divergence, that is, ontogenetic scaling. We then evaluate the hypothesis of size as a line of least evolutionary resistance in catarrhine cranial evolution. We found that patterns of ontogenetic allometry vary among taxa, indicating that ontogenetic scaling sensu stricto does not often account for most morphological differences and that large and small taxa within clades are generally not scaled variants. The presence of a variety of ontogenetic pathways for the evolution of cranial shapes provides indirect evidence for selection acting directly on the cranial shape, rather than on size alone.

Three-dimensional surface scanning methods in osteology: A topographical and geometric morphometric comparison

Objectives

Three‐dimensional (3D) data collected by structured light scanners, photogrammetry, and computed tomography (CT) scans are increasingly combined in joint analyses, even though the scanning techniques and reconstruction software differ considerably. The aim of the present study was to compare the quality and accuracy of surface models and landmark data obtained from modern clinical CT scanning, 3D structured light scanner, photogrammetry, and MicroScribe digitizer.

Material and methods

We tested 13 different photogrammetric software tools and compared surface models obtained by different methods for four articulated human pelves in a topographical analysis. We also measured a set of 219 landmarks and semilandmarks twice on every surface as well as directly on the dry bones with a MicroScribe digitizer.

Results

Only one photogrammetric software package yielded surface models of the complete pelves that could be used for further analysis. Despite the complex pelvic anatomy, all three methods (CT scanning, 3D structured light scanning, photogrammetry) yielded similar surface representations with average deviations among the surface models between 100 and 200 μm. A geometric morphometric analysis of the measured landmarks showed that the different scanning methods yielded similar shape variables, but data acquisition via MicroScribe digitizer was most prone to error.

Discussion

We demonstrated that three‐dimensional models obtained by different methods can be combined in a single analysis. Photogrammetry proved to be a cheap, quick, and accurate method to generate 3D surface models at useful resolutions, but photogrammetry software packages differ enormously in quality.

Detecting Phylogenetic Signal and Adaptation in Papionin Cranial Shape by Decomposing Variation at Different Spatial Scales

Phylogenetic reconstruction based on morphometric data is hampered by homoplasies. For example, many similarities in cranial form between primate taxa more strongly reflect ecological similarities rather than phylogenetic relatedness. However, the way in which the different cranial bones constitute cranial form is, if at all, of less functional relevance and thus largely hidden from selection. We propose that these “constructional details” are better indicators of phylogenetic history than any large-scale shape feature or raw form variable. Within a geometric morphometric context, we show how to analyze the relative extent of bones independently of differences in overall shape. We also show how to decompose total shape variation into small-scale and large-scale shape variation. We apply both methods to the midsagittal cranial morphology of papionin monkeys, which are well known for the discrepancy between morphological similarities and phylogenetic relationships. We study phylogenetic signal and functional adaptation using a molecular phylogeny and contextual data on feeding ecology and locomotor behavior. As expected, total cranial shape, bone outline shape, and large-scale shape features were only weakly associated with phylogenetic distance. But the relative bone contributions and small-scale shape features were both highly correlated with phylogenetic distances. By contrast, the association with ecological and behavioral variables was strongest for the outline shape and large-scale shape features. Studies of morphological adaptation and phylogenetic history thus profit from a decomposition of shape variation into different spatial scales. [Adaptation; canalization; cranial shape; geometric morphometrics; papionini; partial warps; phylogeny.]

Australian Rodents Reveal Conserved Cranial Evolutionary Allometry across 10 Million Years of Murid Evolution

AbstractAmong vertebrates, placental mammals are particularly variable in the covariance between cranial shape and body size (allometry), with rodents being a major exception. Australian murid rodents allow an assessment of the cause of this anomaly because they radiated on an ecologically diverse continent notably lacking other terrestrial placentals. Here, we use 3D geometric morphometrics to quantify species-level and evolutionary allometries in 38 species (317 crania) from all Australian murid genera. We ask whether ecological opportunity resulted in greater allometric diversity compared with other rodents or whether conserved allometry suggests intrinsic constraints and/or stabilizing selection. We also assess whether cranial shape variation follows the proposed rule of craniofacial evolutionary allometry (CREA), whereby larger species have relatively longer snouts and smaller braincases. To ensure we could differentiate parallel versus nonparallel species-level allometric slopes, we compared the slopes of rarefied samples across all clades. We found exceedingly conserved allometry and CREA-like patterns across the 10-million-year split between Mus and Australian murids. This could support both intrinsic-constraint and stabilizing-selection hypotheses for conserved allometry. Large-bodied frugivores evolved faster than other species along the allometric trajectory, which could suggest stabilizing selection on the shape of the masticatory apparatus as body size changes.

Obtaining new resolutions in carnivore tooth pit morphological analyses: A methodological update for digital taphonomy

Modern day investigation in fields of archaeology and palaeontology can be greatly characterised by an exponential growth of integrated new technologies, nevertheless, while these advances are of great significance to multiple lines of research, their evaluation and update over time is equally as important. Here we present an application of inter and intra-observer analysis in taphonomy based geometric morphometrics, employing robust non-parametric statistical analyses for the study of experimental carnivore tooth pit morphologies. To fully understand the influence of measurement errors in the collection of this data, our statistical assessment was performed on fully superimposed, partially superimposed and raw landmark coordinates collected from 3D surface scanning. Experimental samples used to assess these errors includes wolf and dog tooth pits used in modern day ecological livestock predation analysis. Results obtained from this study highlight the importance of landmark type in the assessment of error, emphasising the value of semi-landmark models over the use of ambiguous Type III landmarks. In addition to this, data also reveals the importance of observer experience for the collection of data alongside an interesting increase in error when working with fully superimposed landmarks due to the "Pinocchio Effect". Through this study we are able to redefine the geometric morphometric models used for tooth pit morphological analyses. This final hybrid Type II fixed landmark and semi-landmark model presents a significant reduction in human induced error, generating a more metrically reliable and replicable method that can be used for data pooling in future inter-institutional research. These results can be considered a fundamental step forward for carnivore inspired studies, having an impact on archaeological, palaeontological, modern-day ecological research as well as applications in other forensic sciences.

The evolution of palate shape in the Lepilemur-Cheirogaleidae clade (Primates: Strepsirrhini)

OBJECTIVES

Phylogenies consistently group the folivorous Lepilemur species with the small-bodied insectivorous-frugivorous cheirogaleids. Juvenile lepilemurs and adult cheirogaleids share allometries in most aspects of skull morphology, except the palate. We investigated potential influences on palate shape in these taxa and several outgroups using geometric morphometrics.

MATERIALS AND METHODS

Our sample included representatives of four extant strepsirrhine families, Cheirogaleidae (including Lepilemurinae), Lemuridae, Indriidae, and Galagidae, and one subfossil Megaladapis. Our dataset comprised 32 landmarks collected from 397 specimens representing 15 genera and 28 species, and was analyzed using generalized procrustes analyses and between group principal component analysis. We explored the influence of size, phylogeny, diet, and the propagation of loud vocalizations on palate shape.

RESULTS

While congeneric species clustered within the morphospace, the phylomorphospace did not mirror molecular phylogenetic hypotheses of higher-order relationships. Four palate forms were distinguished within the Cheirogaleidae. Diet, strongly linked to body size, had the single greatest influence on palate shape. The production of long-distance advertisement calls was most often associated with positive scores on the PC1 axis.

DISCUSSION

Our results suggest that the extensive variation in palate shape among Cheirogaleidae is related to dietary shifts that accompanied changes in body size during the clade's radiation. Molecular phylogenies indicate that cheirogaleid diversification involved repeated dwarfing events, which in turn drove dietary shifts from ancestral folivory-frugivory to frugivory, gummivory, and faunivory in the descendant species. The elongated Lepilemur palate is probably related to accelerated eruption of the cheek teeth to render juveniles competent to shear leaves upon weaning.

Patterns and magnitudes of craniofacial covariation in extant cercopithecids

The cranium contains almost all of the vertebrate sensory organs and plays an essential role in vertebrate evolution. Research on the primate cranium has revealed that it is both highly integrated and modular, but studies have historically focused on covariance between the neurocranium and facial skeleton rather than on bones specific to special senses such as vision. The goal of this work is to investigate patterns and magnitudes of craniofacial covariation in extant cercopithecids with particular attention to the orbits. This study takes a quantitative approach using data collected from 38 homologous cranial landmarks across 11 genera of cercopithecid monkeys (Cercopithecidae, N = 291). These data demonstrate that both patterns and magnitudes of craniofacial covariation differ across Cercopithecidae at subfamily, tribe, and genus levels, with the strongest integration in the papionins (and specifically Papio) and significantly weaker covariation in the colobines, particularly Presbytis. Orbital height does not covary with other measurements of the cranium to the same degree as other cranial traits in Cercopithecidae and is highly constrained across the family. This study has important implications for our understanding of the evolution and development of morphological diversity in the cercopithecid cranium and evolution of the primate eye. This study also highlights the potential error of broad assumptions about generalizing patterns and magnitudes of modularity and integration across primates. Additionally, these findings reiterate the importance of trait selection for interpreting fossil taxonomy, as craniofacial covariation may impact phenotypes commonly used to differentiate fossil primate species.

Shape variation in the facial part of the cranium in macaques and African papionins using geometric morphometrics

Macaques are one of the most successful nonhuman primates, and morphological distinctions from their close relatives, African papionins, are easily detected by the naked eye. Nevertheless, evolutionary allometry often accounts for a large amount of the total variation and potentially hides and precludes the detection of morphological distinctions that exist between macaques and African papionins, thus distorting their phyletic comparison. Geometric morpgometric analyses were performed using landmark coordinates in cranial samples from macaques (N =  135) and African papionins (N =  152) to examine the variation in their facial shape. A common allometric trend was confirmed to represent a moderately long face in macaques as being small-to-moderate-bodied papionins. Macaques possessed many features that were distinct from those of African papionins, while they simultaneously showed a large intrageneric variation in every feature, which precluded the separation of some groups of macaques from African papionins. This study confirmed that a moderately smooth sagittal profile is present in non-Sulawesi macaques. It also confirmed that a well-developed anteorbital drop is distinct in Mandrillus and Theropithecus, but it showed that Papio resembles macaques regarding this feature. This finding showed that apparently equivalent features which can be detected by the naked eye were probably formed by different combinations of the principal patterns. It should be noted that the differences detected here between macaques and African papionins are revealed after appropriate adjustments are made to eliminate the allometric effects over the shape features. While landmark data sets still need to be customized for specific studies, the information provided by this article is expected to help such customization and to improve future phyletic evaluation of the fossil papionins.

Allometry and advancing age significantly structure craniofacial variation in adult female baboons

Primate craniofacial growth is traditionally assumed to cease upon maturation or at least be negligible, whereas bony remodeling is typically associated with advanced adult age and, in particular, tooth loss. Therefore, size and shape of the craniofacial skeleton of young and middle-aged adults should be stable. However, research on both modern and historic human samples suggests that portions of the CFS exhibit age-related changes in mature individuals, both related to and independent of tooth loss. These results demonstrate that the age-category 'adult' is heterogeneous, containing individuals demonstrating post-maturational age-related variation, but the topic remains understudied outside of humans and in the cranial vault and base. Our research quantifies variation in a sample of captive adult female baboons (n = 97) in an effort to understand how advancing age alters the mature CFS. Craniometric landmarks and sliding semilandmarks were collected from computed tomography (CT) scans of adult baboons aged 7-32 years old. To determine whether craniofacial morphology is sensitive to aging mechanisms and whether any such effects are differentially distributed throughout the cranium, geometric morphometric techniques were employed to compare the shapes of various cranial regions among individuals of increasing age. Unexpectedly, the biggest form differences were observed between young and middle-aged adults, rather than between adults with full dentitions and those with some degree of tooth loss. Shape variation was greatest in masticatory and nuchal musculature attachment areas. Our results indicate that the craniofacial skeleton changes form during adulthood in baboons, raising interesting questions about the molecular and biological mechanisms governing these changes.

A new three-dimensional geometric morphometrics analysis of the Ouranopithecus macedoniensis cranium (Late Miocene, Central Macedonia, Greece)

OBJECTIVES

This study aims to virtually reconstruct the deformed face (XIR-1) and maxilla (RPl-128) of the Late Miocene hominoid Ouranopithecus macedoniensis from Greece, through the application of mirror-imaging and segmentation. Additionally, analysis was conducted through 3D geometric morphometrics, utilizing a comparative sample of fossil hominoids, extant great apes (Gorilla, Pan, and Pongo) and humans, so as to explore shape variation and phenetic similarities between them.

MATERIALS AND METHODS

High-resolution computed tomography was used to create digital representations of the XIR-1 and RPl-128 specimens. The virtual reconstruction of the XIR-1 cranium was achieved by mirror-imaging, while the RPl-128 maxilla was virtually segmented and reattached in a correct anatomical position. Anatomical landmarks were registered in three dimensions on a comparative sample of adult crania of extant great apes, humans and fossil hominoids. The data were processed with Procrustes superimposition and analyzed using multivariate statistics methods.

RESULTS

Results show that Ouranopithecus macedoniensis falls within or close to the Gorilla convex hull in the principal component analyses, and it is closer to the mean Procrustes shape distance of primarily Gorilla. Both specimens, XIR-1 and RPl-128, are classified as Gorilla based on discriminant function analyses.

DISCUSSION

The results of our geometric morphometrics analyses indicate that Ouranopithecus macedoniensis is morphologically more similar to Gorilla than to Homo, Pan, or Pongo, results that can contribute to the evaluation of existing hypotheses about its phylogenetic position.

Parietal lobe variation in cercopithecid endocasts

In extant primates, the posterior parietal cortex is involved in visuospatial integration, attention, and eye-hand coordination, which are crucial functions for foraging and feeding behaviors. Paleoneurology studies brain evolution through the analysis of endocasts, that is molds of the inner surface of the braincase. These may preserve imprints of cortical structures, such as sulci, which might be of interest for locating the boundaries of major cortical regions. Old World monkeys (Cercopithecidae) represent an interesting zoological group for evolutionary studies, because of their diverse ecologies and locomotor behaviors. In this study, we quantify parietal lobe variation within the cercopithecid family, in a sample of 30 endocasts including 11 genera and 17 species, by combining landmark-based and landmark-free geometric morphometric analyses. More specifically, we quantitatively assess variation of the parietal proportions based on landmarks placed on reliable anatomical references and of parietal lobe surface morphology through deformation-based methods. The main feature associated with the cercopithecid endocranial variation regards the inverse proportions of parietal and occipital lobes, with colobines, Theropithecus, and Papio displaying relatively larger parietal lobes and smaller occipital lobes compared with cercopithecins. The parietal surface is anteroposteriorly longer and mediolaterally flatter in colobines, while longitudinally shorter but laterally bulging in baboons. Large parietal lobes in colobines and baboons are likely to be independent evolutionary traits, and not necessarily associated with analogous functions or morphogenetic mechanisms.

Measurement error in μCT-based three-dimensional geometric morphometrics introduced by surface generation and landmark data acquisition

Computed-tomography-derived (CT-derived) polymesh surfaces are widely used in geometric morphometric studies. This approach is inevitably associated with decisions on scanning parameters, resolution, and segmentation strategies. Although the underlying processing steps have been shown to potentially contribute artefactual variance to three-dimensional landmark coordinates, their effects on measurement error have rarely been assessed systematically in CT-based geometric morphometric studies. The present study systematically assessed artefactual variance in landmark data introduced by the use of different voxel sizes, segmentation strategies, surface simplification degrees, and by inter- and intra-observer differences, and compared their magnitude to true biological variation. Multiple CT-derived surface variants of the anuran (Amphibia: Anura) pectoral girdle were generated by systematic changes in the factors that potentially influence the surface geometries. Twenty-four landmarks were repeatedly acquired by different observers. The contribution of all factors to the total variance in the landmark data was assessed using random-factor nested permanovas. Selected sets of Euclidean distances between landmark sets served further to compare the variance among factor levels. Landmark precision was assessed by landmark standard deviation and compared among observers and days. Results showed that all factors, except for voxel size, significantly contributed to measurement error in at least some of the analyses performed. In total, 6.75% of the variance in landmark data that mimicked a realistic biological study was caused by measurement error. In this landmark dataset, intra-observer error was the major source of artefactual variance followed by inter-observer error; the factor segmentation contributed < 1% and slight surface simplification had no significant effect. Inter-observer error clearly exceeded intra-observer error in a different landmark dataset acquired by six partly inexperienced observers. The results suggest that intra-observer error can potentially be reduced by including a training period prior to the actual landmark acquisition task and by acquiring landmarks in as few sessions as possible. Additionally, the application of moderate and careful surface simplification and, potentially, also the use of case-specific optimal combinations of automatic local thresholding algorithms and parameters for segmentation can help reduce intra-observer error. If landmark data are to be acquired by several observers, it is important to ensure that all observers are consistent in landmark identification. Despite the significant amount of artefactual variance, we have shown that landmark data acquired from microCT-derived surfaces are precise enough to study the shape of anuran pectoral girdles. Yet, a systematic assessment of measurement error is advisable for all geometric morphometric studies.

A missing piece of the Papio puzzle: Gorongosa baboon phenostructure and intrageneric relationships

Most authors recognize six baboon species: hamadryas (Papio hamadryas), Guinea (Papio papio), olive (Papio anubis), yellow (Papio cynocephalus), chacma (Papio ursinus), and Kinda (Papio kindae). However, there is still debate regarding the taxonomic status, phylogenetic relationships, and the amount of gene flow occurring between species. Here, we present ongoing research on baboon morphological diversity in Gorongosa National Park (GNP), located in central Mozambique, south of the Zambezi River, at the southern end of the East African Rift System. The park exhibits outstanding ecological diversity and hosts more than 200 baboon troops. Gorongosa National Park baboons have previously been classified as chacma baboons (P. ursinus). In accordance with this, two mtDNA samples from the park have been placed in the same mtDNA clade as the northern chacma baboons. However, GNP baboons exhibit morphological features common in yellow baboons (e.g., yellow fur color), suggesting that parapatric gene flow between chacma and yellow baboons might have occurred in the past or could be ongoing. We investigated the phenostructure of the Gorongosa baboons using two approaches: 1) description of external phenotypic features, such as coloration and body size, and 2) 3D geometric morphometric analysis of 43 craniofacial landmarks on 11 specimens from Gorongosa compared to a pan-African sample of 352 baboons. The results show that Gorongosa baboons exhibit a mosaic of features shared with southern P. cynocephalus and P. ursinus griseipes. The GNP baboon phenotype fits within a geographic clinal pattern of replacing allotaxa. We put forward the hypothesis of either past and/or ongoing hybridization between the gray-footed chacma and southern yellow baboons in Gorongosa or an isolation-by-distance scenario in which the GNP baboons are geographically and morphologically intermediate. These two scenarios are not mutually exclusive. We highlight the potential of baboons as a useful model to understand speciation and hybridization in early human evolution.

Muscle-induced loading as an important source of variation in craniofacial skeletal shape

The shape of the craniofacial skeleton is constantly changing through ontogeny and reflects a balance between developmental patterning and mechanical-load-induced remodeling. Muscles are a major contributor to producing the mechanical environment that is crucial for "normal" skull development. Here, we use an F₅ hybrid population of Lake Malawi cichlids to characterize the strength and types of associations between craniofacial bones and muscles. We focus on four bones/bone complexes, with different developmental origins, alongside four muscles with distinct functions. We used micro-computed tomography to extract 3D information on bones and muscles. 3D geometric morphometrics and volumetric measurements were used to characterize bone and muscle shape, respectively. Linear regressions were performed to test for associations between bone shape and muscle volume. We identified three types of associations between muscles and bones: weak, strong direct (i.e., muscles insert directly onto bone), and strong indirect (i.e., bone is influenced by muscles without a direct connection). In addition, we show that although the shape of some bones is relatively robust to muscle-induced mechanical stimulus, others appear to be highly sensitive to muscular input. Our results imply that the roles for muscular input on skeletal shape extend beyond specific points of origin or insertion and hold significant potential to influence broader patterns of craniofacial geometry. Thus, changes in the loading environment, either as a normal course of ontogeny or if an organism is exposed to a novel environment, may have pronounced effects on skeletal shape via near and far-ranging effects of muscular loading.

Differential evolution of cerebral and cerebellar fossae in recent Homo: A new methodological approach

The endocranium shows the influence of the shape and development of brain tissues and overall brain modifications. During the late Upper Pleistocene and Holocene smaller brains appeared and the higher position of endinion relative to inion might indicate changes in cerebellar and occipital lobes. In previous studies, the depths of the cerebral and cerebellar fossae were not specifically considered; new tools for quantitatively measuring these irregular, problematic curved areas need to be developed. This paper's main objective is to investigate to what degree changes in the fossae's depths of extant humans have occurred with respect to fossil anatomically modern humans (AMH) and older Homo species. The proportions of the occipital and nuchal planes are compared measuring the inner and outer surfaces of the bone. Additionally, this paper proposes a quantitative geometric methodology based on endocranial landmarks that create a plane with which to measure the position of the deepest part of the fossa: it represents a curvature maxima - concavity - associated with local structures. The four points thus obtained could be framed in Bookstein's Type II landmarks but without biomechanical implication. Through univariate, bivariate and multivariate analyses (principal components analysis) of raw and size-corrected data we study the differential evolution in recent Homo species, which presents a more vertical occipital area than ancient fossils. Our results corroborate this derived trait; additionally, we have observed a tendency towards a relative decrease in the depth of the cerebral fossae and maintenance of the cerebellar fossae.

Creating diversity in mammalian facial morphology: a review of potential developmental mechanisms

Mammals (class Mammalia) have evolved diverse craniofacial morphology to adapt to a wide range of ecological niches. However, the genetic and developmental mechanisms underlying the diversification of mammalian craniofacial morphology remain largely unknown. In this paper, we focus on the facial length and orofacial clefts of mammals and deduce potential mechanisms that produced diversity in mammalian facial morphology. Small-scale changes in facial morphology from the common ancestor, such as slight changes in facial length and the evolution of the midline cleft in some lineages of bats, could be attributed to heterochrony in facial bone ossification. In contrast, large-scale changes of facial morphology from the common ancestor, such as a truncated, widened face as well as the evolution of the bilateral cleft possessed by some bat species, could be brought about by changes in growth and patterning of the facial primordium (the facial processes) at the early stages of embryogenesis.

Phylogenetic relationship of a fossil macaque (Macaca cf. robusta) from the Korean Peninsula to extant species of macaques based on zygomaxillary morphology

Little is known about the biogeographical and evolutionary histories of macaques (Macaca spp.) in East Asia because the phylogenetic positions of fossil species remain unclear. Here we examined the zygomaxillary remains of a fossil macaque (M. cf. robusta) from the Durubong Cave Complex, South Korea, that dates back to the late Middle to Late Pleistocene, to infer its phylogenetic relationship to extant species. We took 195 fixed- and semi-landmarks from the zygomaxillary regions of the fossil specimen and from 147 specimens belonging to 14 extant species. We then conducted a generalized Procrustes analysis followed by a multivariate statistical analysis to evaluate the phenetic affinities of the fossil to the extant species and reconstructed the most parsimonious phylogenetic tree using a phylogenetic morphometric approach. We found that the fossil was most similar to Macaca fuscata (Japanese macaque) in the zygomaxillary morphospace although it was at the limit of the range of variation for this species. The second closest in the morphospace was the continental Macaca mulatta (rhesus macaque). Parsimonious reconstruction confirmed that the fossil was most closely related to M. fuscata, even after controlling for the effects of allometry. These findings suggest that in the late Middle to Late Pleistocene, close relatives of M. fuscata that looked like the extant species were distributed on the Korean Peninsula, where no species of macaques are found today. Thus, some morphological characteristics of M. fuscata may have developed before its ancestor dispersed into the Japanese archipelago.

The evolution of genital shape variation in female cetaceans

Male genital diversification is likely the result of sexual selection. Female genital diversification may also result from sexual selection, although it is less well studied and understood. Female genitalia are complex among whales, dolphins, and porpoises, especially compared to other vertebrates. The evolutionary factors affecting the diversity of vaginal complexity could include ontogeny, allometry, phylogeny, sexual selection, and natural selection. We quantified shape variation in female genitalia using 2D geometric morphometric analysis, and validated the application of this method to study soft tissues. We explored patterns of variation in the shape of the cervix and vagina of 24 cetacean species (n = 61 specimens), and found that genital shape varies primarily in the relative vaginal length and overall aspect ratio of the reproductive tract. Extensive genital shape variation was partly explained by ontogenetic changes and evolutionary allometry among sexually mature cetaceans, whereas phylogenetic signal, relative testis size, and neonate size were not significantly associated with genital shape. Female genital shape is diverse and evolves rapidly even among closely related species, consistent with predictions of sexual selection models and with findings in invertebrate and vertebrate taxa. Future research exploring genital shape variation in 3D will offer new insights into evolutionary mechanisms because internal vaginal structures are variable and can form complex spirals.

Additive genetic variation in the craniofacial skeleton of baboons (genus Papio) and its relationship to body and cranial size

OBJECTIVES

Determining the genetic architecture of quantitative traits and genetic correlations among them is important for understanding morphological evolution patterns. We address two questions regarding papionin evolution: (1) what effect do body and cranial size, age, and sex have on phenotypic (V_P ) and additive genetic (V_A ) variation in baboon crania, and (2) how might additive genetic correlations between craniofacial traits and body mass affect morphological evolution?

MATERIALS AND METHODS

We use a large captive pedigreed baboon sample to estimate quantitative genetic parameters for craniofacial dimensions (EIDs). Our models include nested combinations of the covariates listed above. We also simulate the correlated response of a given EID due to selection on body mass alone.

RESULTS

Covariates account for 1.2-91% of craniofacial V_P . EID V_A decreases across models as more covariates are included. The median genetic correlation estimate between each EID and body mass is 0.33. Analysis of the multivariate response to selection reveals that observed patterns of craniofacial variation in extant baboons cannot be attributed solely to correlated response to selection on body mass, particularly in males.

DISCUSSION

Because a relatively large proportion of EID V_A is shared with body mass variation, different methods of correcting for allometry by statistically controlling for size can alter residual V_P patterns. This may conflate direct selection effects on craniofacial variation with those resulting from a correlated response to body mass selection. This shared genetic variation may partially explain how selection for increased body mass in two different papionin lineages produced remarkably similar craniofacial phenotypes.

Cranial sexual dimorphism in the Kinda baboon (Papio hamadryas kindae)

OBJECTIVES

The smallest extant member of genus Papio, the Kinda baboon exhibits low sexual dimorphism and a distinctive cranial shape. Ontogenetic scaling accounts for most cranial-shape differences within Papio, but studies have shown that the Kinda follows a separate ontogenetic trajectory. If so, its cranial-dimorphism pattern should differ from other subspecies. To evaluate this hypothesis, morphometric analysis was used to investigate cranial dimorphism in Papio.

MATERIALS AND METHODS

Three-dimensional landmarks were digitized on 434 adult crania representing six Papio subspecies. Size- and shape-dimorphism magnitudes were quantified using centroid size and Procrustes distances. Patterns of sex- and size-related variation were explored using MAN(C)OVA, multivariate regression, and form-space PCA. Canine dimorphism was investigated using dental metrics.

RESULTS

Kinda size and shape dimorphism are significantly lower than in other Papio subspecies. The relative magnitude of Kinda shape dimorphism is similar to other southern baboons; Kinda canine dimorphism is unremarkable. MAN(C)OVA results support subspecies differences in cranial dimorphism and scaling. Allometric and dimorphism vectors differ significantly in some subspecies, and their vector-angle matrices are strongly correlated. The Kinda's allometric vector angles are divergent. Form-space PC3, summarizing size-independent dimorphism, separates the Kinda from other subspecies.

DISCUSSION

The Kinda baboon exhibits significantly lower size and shape dimorphism than other baboons, but its relative dimorphism levels are unexceptional. The Kinda differs from other subspecies in patterns of allometry, size-related shape dimorphism, and residual shape dimorphism. Kinda facial shape is "masculinized" relative to size, especially in females, suggesting female sexual selection contributed to the evolution of Kinda dimorphism.

New cranium of the large cercopithecid primate Theropithecus oswaldi leakeyi (Hopwood, 1934) from the paleoanthropological site of Makuyuni, Tanzania

The Pleistocene hominin site of Makuyuni, near Lake Manyara, Tanzania, is known for fossils attributable to Homo and Acheulean artifacts (Ring et al., 2005; Kaiser et al., 2010; Frost et al., 2012). Here we describe the fossil primate material from the Manyara Beds, which includes the first nearly complete female cranium of Theropithecus oswaldi leakeyi and a proximal tibia from the same taxon. The cranium is dated to between 633 and 780 Ka and the tibia to the Pleistocene. The T. oswaldi lineage is one of the most important among Neogene mammals of Africa: it is both widespread and abundant. The size of the dentition, cranium, and tibia all confirm the previously recognized trend of increasing body size in this lineage and make their taxonomic assignments secure. The morphology of this specimen provides new insights into the evolution of this lineage through time, as well as its geographic variation and sexual dimorphism. The cranium also shows damage consistent with a mammalian carnivore, most likely a felid. The identification of this material as representing T. o. leakeyi agrees with the Middle Pleistocene age estimates for the MK4 locality in particular and the Manyara Beds in general.

Resolving relationships between several Neolithic and Mesolithic populations in Northern Eurasia using geometric morphometrics

OBJECTIVES

Remains from several Eastern European and Siberian Mesolithic and Neolithic sites are analysed to clarify their biological relationships. We assume that groups' geographical distances correlate with genetic and, therefore, morphological distances between them.

MATERIALS AND METHODS

Material includes complete male crania from several Mesolithic and Neolithic burial sites across Northern Eurasia and from several modern populations. Geometric morphometrics and multivariate statistical techniques are applied to explore morphological trends, group distances, and correlations with their geographical position, climate, and the time of origin.

RESULTS

Despite an overlap in the morphology among the modern and archeological groups, some of them show significant morphological distances. Geographical parameters account for only a small proportion of cranial variation in the sample, with larger variance explained by geography and age together. Expectations of isolation by distance are met in some but not in all cases. Climate accounts for a large proportion of autocorrelation with geography. Nearest-neighbor joining trees demonstrate group relationships predicted by the regression on geography and on climate.

DISCUSSION

The obtained results are discussed in application to relationships between particular groups. Unlike the Ukrainian Mesolithic, the Yuzhny Oleni Ostrov Mesolithic displays a high morphological affinity with several groups from Northern Eurasia of both European and Asian origin. A possibility of a common substrate for the Yuzhny Oleni Ostrov Mesolithic and Siberian Neolithic groups is reviewed. The Siberian Neolithic is shown to have morphological connection with both modern Siberian groups and the Native North Americans.

Socioeconomic Status Is Not Related with Facial Fluctuating Asymmetry: Evidence from Latin-American Populations

The expression of facial asymmetries has been recurrently related with poverty and/or disadvantaged socioeconomic status. Departing from the developmental instability theory, previous approaches attempted to test the statistical relationship between the stress experienced by individuals grown in poor conditions and an increase in facial and corporal asymmetry. Here we aim to further evaluate such hypothesis on a large sample of admixed Latin Americans individuals by exploring if low socioeconomic status individuals tend to exhibit greater facial fluctuating asymmetry values. To do so, we implement Procrustes analysis of variance and Hierarchical Linear Modelling (HLM) to estimate potential associations between facial fluctuating asymmetry values and socioeconomic status. We report significant relationships between facial fluctuating asymmetry values and age, sex, and genetic ancestry, while socioeconomic status failed to exhibit any strong statistical relationship with facial asymmetry. These results are persistent after the effect of heterozygosity (a proxy for genetic ancestry) is controlled in the model. Our results indicate that, at least on the studied sample, there is no relationship between socioeconomic stress (as intended as low socioeconomic status) and facial asymmetries.