Untangling the relationship between developmental and evolutionary integration.

Patterns of integration and modularity among organismal traits are prevalent across the tree of life, and at multiple scales of biological organization. Over the past several decades, researchers have studied these patterns at the developmental, and evolutionary levels. While their work has identified the potential drivers of these patterns at different scales, there appears to be a lack of consensus on the relationship between developmental and evolutionary integration. Here, we review and summarize key studies and build a framework to describe the conceptual relationship between these patterns across organismal scales and illustrate how, and why some of these studies may have yielded seemingly conflicting outcomes. We find that among studies that analyze patterns of integration and modularity using morphological data, the lack of consensus may stem in part from the difficulty of fully disentangling the developmental and functional causes of integration. Nonetheless, in some empirical systems, patterns of evolutionary modularity have been found to coincide with expectations based on developmental processes, suggesting that in some circumstances, developmental modularity may translate to evolutionary modularity. We also advance an extension to Hallgrímsson et al.'s palimpsest model to describe how patterns of trait modularity may shift across different evolutionary scales. Finally, we also propose some directions for future research which will hopefully be useful for investigators interested in these issues.

Deep learning to capture leaf shape in plant images: Validation by geometric morphometrics.

Plant leaves play a pivotal role in automated species identification using deep learning (DL). However, achieving reproducible capture of leaf variation remains challenging due to the inherent "black box" problem of DL models. To evaluate the effectiveness of DL in capturing leaf shape, we used geometric morphometrics (GM), an emerging component of eXplainable Artificial Intelligence (XAI) toolkits. We photographed Ranunculus auricomus leaves directly in situ and after herbarization. From these corresponding leaf images, we automatically extracted DL features using a neural network and digitized leaf shapes using GM. The association between the extracted DL features and GM shapes was then evaluated using dimension reduction and covariation models. DL features facilitated the clustering of leaf images by source populations in both in situ and herbarized leaf image datasets, and certain DL features were significantly associated with biological leaf shape variation as inferred by GM. DL features also enabled leaf classification into morpho-phylogenomic groups within the intricate R. auricomus species complex. We demonstrated that simple in situ leaf imaging and DL reproducibly captured leaf shape variation at the population level, while combining this approach with GM provided key insights into the shape information extracted from images by computer vision, a necessary prerequisite for reliable automated plant phenotyping.

Dental data challenge the ubiquitous presence of Homo in the Cradle of Humankind.

The origins of Homo, as well as the diversity and biogeographic distribution of early Homo species, remain critical outstanding issues in paleoanthropology. Debates about the recognition of early Homo, first appearance dates, and taxonomic diversity within Homo are particularly important for determining the role that southern African taxa may have played in the origins of the genus. The correct identification of Homo remains also has implications for reconstructing phylogenetic relationships between species of Australopithecus and Paranthropus, and the links between early Homo species and Homo erectus. We use microcomputed tomography and landmark-free deformation-based three-dimensional geometric morphometrics to extract taxonomically informative data from the internal structure of postcanine teeth attributed to Early Pleistocene Homo in the southern African hominin-bearing sites of Sterkfontein, Swartkrans, Drimolen, and Kromdraai B. Our results indicate that, from our sample of 23 specimens, only 4 are unambiguously attributed to Homo, 3 of them coming from Swartkrans member 1 (SK 27, SK 847, and SKX 21204) and 1 from Sterkfontein (Sts 9). Three other specimens from Sterkfontein (StW 80 and 81, SE 1508, and StW 669) approximate the Homo condition in terms of overall enamel-dentine junction shape, but retain Australopithecus-like dental traits, and their generic status remains unclear. The other specimens, including SK 15, present a dominant australopith dental signature. In light of these results, previous dietary and ecological interpretations can be reevaluated, showing that the geochemical signal of one tooth from Kromdraai (KB 5223) and two from Swartkrans (SK 96 and SKX 268) is consistent with that of australopiths.

Got milkweed? Genetic assimilation as potential source for the evolution of nonmigratory monarch butterfly wing shape.

Monarch butterflies (Danaus plexippus) are well studied for their annual long-distance migration from as far north as Canada to their overwintering grounds in Central Mexico. At the end of the cold season, monarchs start to repopulate North America through short-distance migration over the course of multiple generations. Interestingly, some populations in various tropical and subtropical islands do not migrate and exhibit heritable differences in wing shape and size, most likely an adaptation to island life. Less is known about forewing differences between long- and short-distance migrants in relation to island populations. Given their different migratory behaviors, we hypothesized that these differences would be reflected in wing morphology. To test this, we analyzed forewing shape and size of three different groups: nonmigratory, lesser migratory (migrate short-distances), and migratory (migrate long-distances) individuals. Significant differences in shape appear in all groups using geometric morphometrics. As variation found between migratory and lesser migrants has been shown to be caused by phenotypic plasticity, and lesser migrants develop intermediate forewing shapes between migratory and nonmigratory individuals, we suggest that genetic assimilation might be an important mechanism to explain the heritable variation found between migratory and nonmigratory populations. Additionally, our research confirms previous studies which show that forewing size is significantly smaller in nonmigratory populations when compared to both migratory phenotypes. Finally, we found sexual dimorphism in forewing shape in all three groups, but for size in nonmigratory populations only. This might have been caused by reduced constraints on forewing size in nonmigratory populations.

Utilizing geometric morphometrics to investigate gene function during organ growth: Insights through the study of beetle horn shape allometry.

Static allometry is a major component of morphological variation. Much of the literature on the development of allometry investigates how functional perturbations of diverse pathways affect the relationship between trait size and body size. Often, this is done with the explicit objective to identify developmental mechanisms that enable the sensing of organ size and the regulation of relative growth. However, changes in relative trait size can also be brought about by a range of other distinctly different developmental processes, such as changes in patterning or tissue folding, yet standard univariate biometric approaches are usually unable to distinguish among alternative explanations. Here, we utilize geometric morphometrics to investigate the degree to which functional genetic manipulations known to affect the size of dung beetle horns also recapitulate the effect of horn shape allometry. We reasoned that the knockdown phenotypes of pathways governing relative growth should closely resemble shape variation induced by natural allometric variation. In contrast, we predicted that if genes primarily affect alternative developmental processes, knockdown effects should align poorly with shape allometry. We find that the knockdown effects of several genes (e.g., doublesex, Foxo) indeed closely aligned with shape allometry, indicating that their corresponding pathways may indeed function primarily in the regulation of relative trait growth. In contrast, other knockdown effects (e.g., Distal-less, dachs) failed to align with allometry, implicating these pathways in potentially scaling-independent processes. Our findings moderate the interpretation of studies focusing on trait length and highlight the usefulness of multivariate approaches to study allometry and phenotypic plasticity.

Replicated Functional Evolution in Cichlid Adaptive Radiations.

AbstractAdaptive radiations highlight the mechanisms by which species and traits diversify and the extent to which these patterns are predictable. We used 1,110 high-speed videos of suction feeding to study functional and morphological diversification in 300 cichlid species from three African Great Lake radiations of varying ages (Victoria, Malawi, and Tanganyika) and an older, spatially dispersed continental radiation in the Neotropics. Among African radiations, standing diversity was reflective of time. Morphological and functional variance in Lake Victoria, the youngest radiation, was a subset of that within Lake Malawi, which itself was nested within the older Tanganyikan radiation. However, functional diversity in Neotropical cichlids was often lower than that in Lake Tanganyika, despite being much older. These two radiations broadly overlapped, but each diversified into novel trait spaces not found in the youngest lake radiations. Evolutionary rates across radiations were inversely related to age, suggesting extremely rapid trait evolution at early stages, particularly in lake radiations. Despite this support for early bursts, other patterns of trait diversity were inconsistent with expectations of adaptive radiations. This work suggests that cichlid functional evolution has played out in strikingly similar fashion in different radiations, with contingencies eventually resulting in lineage-specific novelties.

Freshwater Habitats Promote Rapid Rates of Phenotypic Evolution in Sculpin Fishes (Perciformes: Cottoidea).

AbstractInvasions of freshwater habitats by marine fishes provide exceptional cases of habitat-driven biological diversification. Freshwater habitats make up less than 1% of aquatic habitats but contain ∼50% of fish species. However, while the dominant group of freshwater fishes (Otophysi) is older than that of most marine fishes (Percomorphaceae), it is less morphologically diverse. Classically, scientists have invoked differences in the tempo and/or mode of evolution to explain such cases of unequal morphological diversification. We tested for evidence of these phenomena in the superfamily Cottoidea (sculpins), which contains substantial radiations of marine and freshwater fishes. We find that the morphology of freshwater sculpins evolves faster but under higher constraint than that of marine sculpins, causing widespread convergence in freshwater sculpins and more morphological disparity in marine sculpins. The endemic freshwater sculpins of Lake Baikal, Siberia, are exceptions that demonstrate elevated novelty akin to that of marine sculpins. Several tantalizing factors may explain these findings, such as differences in habitat stability and/or habitat connectivity between marine and freshwater systems.

A shared pattern of midfacial bone modelling in hominids suggests deep evolutionary roots for human facial morphogenesis.

Midfacial morphology varies between hominoids, in particular between great apes and humans for which the face is small and retracted. The underlying developmental processes for these morphological differences are still largely unknown. Here, we investigate the cellular mechanism of maxillary development (bone modelling, BM), and how potential changes in this process may have shaped facial evolution. We analysed cross-sectional developmental series of gibbons, orangutans, gorillas, chimpanzees and present-day humans (n = 183). Individuals were organized into five age groups according to their dental development. To visualize each species's BM pattern and corresponding morphology during ontogeny, maps based on microscopic data were mapped onto species-specific age group average shapes obtained using geometric morphometrics. The amount of bone resorption was quantified and compared between species. Great apes share a highly similar BM pattern, whereas gibbons have a distinctive resorption pattern. This suggests a change in cellular activity on the hominid branch. Humans possess most of the great ape pattern, but bone resorption is high in the canine area from birth on, suggesting a key role of canine reduction in facial evolution. We also observed that humans have high levels of bone resorption during childhood, a feature not shared with other apes.

Intra-leaf modeling of Cannabis leaflet shape produces leaf models that predict genetic and developmental identities.

The iconic, palmately compound leaves of Cannabis have attracted significant attention in the past. However, investigations into the genetic basis of leaf shape or its connections to phytochemical composition have yielded inconclusive results. This is partly due to prominent changes in leaflet number within a single plant during development, which has so far prevented the proper use of common morphometric techniques. Here, we present a new method that overcomes the challenge of nonhomologous landmarks in palmate, pinnate, and lobed leaves, using Cannabis as an example. We model corresponding pseudo-landmarks for each leaflet as angle-radius coordinates and model them as a function of leaflet to create continuous polynomial models, bypassing the problems associated with variable number of leaflets between leaves. We analyze 341 leaves from 24 individuals from nine Cannabis accessions. Using 3591 pseudo-landmarks in modeled leaves, we accurately predict accession identity, leaflet number, and relative node number. Intra-leaf modeling offers a rapid, cost-effective means of identifying Cannabis accessions, making it a valuable tool for future taxonomic studies, cultivar recognition, and possibly chemical content analysis and sex identification, in addition to permitting the morphometric analysis of leaves in any species with variable numbers of leaflets or lobes.

Phenotypic selection patterns in a hybrid zone between two Calceolaria species with contrasting pollinators: insights from field surveys and fitness assessments.

Hybrid zones provide natural experimental settings to test hypotheses about species divergence. We concentrated on a hybrid swarm in which oil-collecting bees and flower-pecking birds act as pollinators of two Calceolaria species. We asked whether both pollinators contributed to flower divergence by differentially promoting prezygotic fitness at the phenotypic extremes that represent parentals. We studied pollinator-mediated selection on phenotypic traits critical in plant-pollinator mechanical interaction, namely plant height, reward-to-stigma distance, and flower shape. We utilised the quantity and quality of pollen deposited as fitness measures and distinguished between the contribution of the two pollinator types. Results showed uni- and bivariate disruptive selection for most traits through pollen grains deposited by both pollinators. Bird-mediated fitness favoured low plants with a long reward-to-stigma distance and a straight corolla, while bee-mediated fitness favoured tall plants with a short reward-to-stigma distance and curved corolla. In addition, stabilising selection at one end of the phenotypic range showed a bird-mediated reproductive asymmetry within the swarm. The disruptive pattern was countered, albeit weakly, by hybrids receiving higher-quality pollen on the stigmas. Results suggest that pollinator-mediated selection promotes divergence of integrated flower phenotypes mechanically adjusted either to bees or birds underscoring the importance of pollinator specialisation in diversification.

Signatures of convergence in Neotropical cichlid fish.

Convergent evolution of similar phenotypes suggests some predictability in the evolutionary trajectories of organisms, due to strong and repeated selective pressures, and/or developmental constraints. In adaptive radiations, particularly in cichlid fish radiations, convergent phenotypes are commonly found within and across geographical settings. Cichlids show major repeated axes of morphological diversification. Recurrent changes in body patterns reveal adaption to alternative habitats, and modifications of the trophic apparatus respond to the exploitation of different food resources. Here we compare morphologically and genetically two Neotropical cichlid assemblages, the Mexican desert cichlid and the Nicaraguan Midas cichlid, with similar polymorphic body and trophic adaptations despite their independent evolution. We found a common morphological axis of differentiation in trophic structures in both cichlid radiations, but two different axes of differentiation in body shape, defining two alternative limnetic body patterns. Adaptation to limnetic habitats implied regulation of immune functions in the Midas cichlid, while morphogenesis and metabolic functions in the desert cichlid. Convergent phenotypic adaptions could be associated to divergent gene regulation.

Charting facial growth and development for Bantu Africans: Central tendencies, variational properties and sexual dimorphisms.

Current studies on facial growth and development have been largely based on European populations. Less studied are African populations, who because of their distinct genetic makeup and environmental conditions, provide deeper insights into patterns of facial development. Patterns of facial shape development in African populations remain largely uncharacterised. Our study aimed to establish facial growth and development trajectories based on a cohort of 2874 Bantu Africans from Tanzania aged 6-18 years, with particular focus on identifying morphogenetic processes that lead to observed developmental shape changes. Procrustes ANCOVA suggested sexually dimorphic patterns of facial shape development (p = 0.0036). The forehead was relatively contracted during development in both sexes. The glabella region was more anteriorly displaced in females due to expansion in the region laterosuperior to the eyes. Nasal protrusion increased with development, which was found to arise from local expansion in the nasal alae and columella. Local expansion in the upper and lower labial regions resulted in forward displaced lips in both sexes, with the effect more pronounced in males. The mentum was displaced more anteriorly in females due to comparatively more expanded mental regions with development. The lateral facial region corresponding to the underlying body of the mandible were developmentally expanded but were posteriorly positioned due to protrusive growth of surrounding structures. Generalised additive modelling of Procrustes variance suggested that facial variation decreased non-linearly with age (p < 0.05). Relative principal component analysis suggested that variations in facial outline shape were developmentally constrained, whereas nasolabial and mental regions, where developmental changes were significant, became morphologically diversified with development. In contrast to simple descriptive illustration of facial shape development, we gained transformative insights into patterns of facial shape development by analysing morphogenetic processes and variational properties. Our analytical framework is broadly applicable to morphometric studies on ontogenetic shape changes.

A 3D geometric morphometric analysis of the bovid distal humerus, with special reference to Rusingoryx atopocranion (Pleistocene, Eastern Africa).

The family Bovidae [Mammalia: Artiodactyla] is speciose and has extant representatives on every continent, forming key components of mammal communities. For these reasons, bovids are ideal candidates for studies of ecomorphology. In particular, the morphology of the bovid humerus has been identified as highly related to functional variables such as body mass and habitat. This study investigates the functional morphology of the bovid distal humerus in isolation due to its increased likelihood of preservation in the fossil record, and the resulting opportunity for a better understanding of the ecomorphology of extinct bovids. A landmark scheme of 30 landmarks was used to capture the 3D distal humerus morphology in 111 extant bovid specimens. We find that the distal humerus has identifiable morphologies associated with body mass, habitat preference and tribe affiliation and that some characteristics are shared between high body mass bovids and those living on hard, flat terrain which is likely due to the high stress on the bone in both cases. We directly apply our findings regarding extant bovids to the extinct alcelaphine bovid, Rusingoryx atopocranion from the mid to late Pleistocene (>33-45 ka) Lake Victoria region of Kenya. This species is known for some peculiar morphologies including a domed cranium with hollow nasal crests, and having small hooves for a bovid of its size. Another interesting aspect of Rusingoryx's skeletal morphology which has not been addressed is an unusual protrusion on the lateral epicondyle of the distal humerus. Despite considerable individual variation in the Rusingoryx specimens, we find evidence to support its historical assignment to the tribe Alcelaphini, and that it likely preferred open grassland habitats, which is consistent with independent reconstructions of the palaeoenvironment. We also provide the most accurate body mass estimate for Rusingoryx to date, based on distal humerus centroid size. Overall, we are able to conclude that the distal humerus in extant bovids is highly informative regarding body mass, habitat preference and tribe, and that this can be applied directly to a fossil taxon with promising results.

Geometric growth of the normal human craniocervical junction from 0 to 18 years old.

The craniocervical junction (CCJ) forms the bridge between the skull and the spine, a highly mobile group of joints that allows the mobility of the head in every direction. The CCJ plays a major role in protecting the inferior brainstem (bulb) and spinal cord, therefore also requiring some stability. Children are subjected to multiple constitutive or acquired diseases involving the CCJ: primary bone diseases such as in FGFR-related craniosynostoses or acquired conditions such as congenital torticollis, cervical spine luxation, and neurological disorders. To design efficient treatment plans, it is crucial to understand the relationship between abnormalities of the craniofacial region and abnormalities of the CCJ. This can be approached by the study of control and abnormal growth patterns. Here we report a model of normal skull base growth by compiling a collection of geometric models in control children. Focused analyses highlighted specific developmental patterns for each CCJ bone, emphasizing rapid growth during infancy, followed by varying rates of growth and maturation during childhood and adolescence until reaching stability by 18 years of age. The focus was on the closure patterns of synchondroses and sutures in the occipital bone, revealing distinct closure trajectories for the anterior intra-occipital synchondroses and the occipitomastoid suture. The findings, although based on a limited dataset, showcased specific age-related changes in width and closure percentages, providing valuable insights into growth dynamics within the first 2 years of life. Integration analyses revealed intricate relationships between skull and neck structures, emphasizing coordinated growth at different stages. Specific bone covariation patterns, as found between the first and second cervical vertebrae (C1 and C2), indicated synchronized morphological changes. Our results provide initial data for designing inclusive CCJ geometric models to predict normal and abnormal growth dynamics.

Maxillary morphology of chimpanzees: Captive versus wild environments.

Morphological studies typically avoid using osteological samples that derive from captive animals because it is assumed that their morphology is not representative of wild populations. Rearing environments indeed differ between wild and captive individuals. For example, mechanical properties of the diets provided to captive animals can be drastically different from the food present in their natural habitats, which could impact cranial morphology and dental health. Here, we examine morphological differences in the maxillae of wild versus captive chimpanzees (Pan troglodytes) given the prominence of this species in comparative samples used in human evolution research and the key role of the maxilla in such studies. Size and shape were analysed using three-dimensional geometric morphometric methods based on computed tomography scans of 94 wild and 30 captive specimens. Captive individuals have on average larger and more asymmetrical maxillae than wild chimpanzees, and significant differences are present in their maxillary shapes. A large proportion of these shape differences are attributable to static allometry, but wild and captive specimens still differ significantly from each other after allometric size adjustment of the shape data. Levels of shape variation are higher in the captive group, while the degree of size variation is likely similar in our two samples. Results are discussed in the context of ontogenetic growth trajectories, changes in dietary texture, an altered social environment, and generational differences. Additionally, sample simulations show that size and shape differences between chimpanzees and bonobos (Pan paniscus) are exaggerated when part of the wild sample is replaced with captive chimpanzees. Overall, this study confirms that maxillae of captive chimpanzees should not be included in morphological or taxonomic analyses when the objective is to characterise the species.

How to make a digital reconstruction of the human ribcage.

Up to now, there have been no publication standardizing the digital reconstruction of the modern human ribcage from commingled costo-vertebral material. Consequently, we designed a validated protocol based on anatomical features observed in the literature and the CT scanned ribcages of 10 adult European individuals. After quantifying the shape of these ribcages using 3D geometric morphometrics, we split each vertebra and rib within their corresponding (semi)landmarks. Subsequently, individual bones + (semi)landmarks were imported to LhpFusionBox, commingled and 3D reconstructed. To validate the accuracy of the protocol, we first reconstructed a randomly chosen ribcage three times and then compared these reconstructions to the rest of the sample. Since these reconstructions were closer to their original counterpart than to the others, the remaining sample was reconstructed once. Next, we tested the intra-observer error during reconstructing using the Procrustes distances among the original ribcages and the reconstructions. We observed that first each ribcage reconstruction was clustered to its original counterpart and second there was a learning curve showing an improvement in the reconstruction process over time. Subsequently, we explored general size and shape differences among the original and reconstructed ribcages through a study of centroid size and a permutation test on the Procrustes distances (10,000 permutations), respectively. Specific shape differences between both groups were further examined through a principal component analysis in shape space. None of these analyses found statistical differences between the original and reconstructed ribcages (p > 0.05). Eventually, we extracted the mean shapes of the original ribcages and the reconstructions in order to visualize potential deviations caused by the anatomical considerations of the researcher. These results demonstrate that the protocol is accurate enough to be used when reconstructing a disarticulated human ribcage.

Sexual dimorphism and allometry in human scapula shape.

Scapula shape is highly variable across humans and appears to be sexually dimorphic-differing significantly between biological males and females. However, previous investigations of sexual dimorphism in scapula shape have not considered the effects of allometry (the relationship between size and shape). Disentangling allometry from sexual dimorphism is necessary because apparent sex-based differences in shape could be due to inherent differences in body size. This study aimed to investigate sexual dimorphism in scapula shape and examine the role of allometry in sex-based variation. We used three-dimensional geometric morphometrics with Procrustes ANOVA to quantify scapula shape variation associated with sex and size in 125 scapulae. Scapula shape significantly differed between males and females, and males tended to have larger scapulae than females for the same body height. We found that males and females exhibited distinct allometric relationships, and sexually dimorphic shape changes did not align with male- or female-specific allometry. A secondary test revealed that sexual dimorphism in scapula shape persisted between males and females of similar body heights. Overall, our findings indicate that there are sex-based differences in scapula shape that cannot be attributed to size-shape relationships. Our results shed light on the potential role of sexual selection in human shoulder evolution, present new hypotheses for biomechanical differences in shoulder function between sexes, and identify relevant traits for improving sex classification accuracy in forensic analyses.

Inner ear morphology in wild versus laboratory house mice.

The semicircular canals of the inner ear are involved in balance and velocity control. Being crucial to ensure efficient mobility, their morphology exhibits an evolutionary conservatism attributed to stabilizing selection. Release of selection in slow-moving animals has been argued to lead to morphological divergence and increased inter-individual variation. In its natural habitat, the house mouse Mus musculus moves in a tridimensional space where efficient balance is required. In contrast, laboratory mice in standard cages are severely restricted in their ability to move, which possibly reduces selection on the inner ear morphology. This effect was tested by comparing four groups of mice: several populations of wild mice trapped in commensal habitats in France; their second-generation laboratory offspring, to assess plastic effects related to breeding conditions; a standard laboratory strain (Swiss) that evolved for many generations in a regime of mobility reduction; and hybrids between wild offspring and Swiss mice. The morphology of the semicircular canals was quantified using a set of 3D landmarks and semi-landmarks analyzed using geometric morphometric protocols. Levels of inter-population, inter-individual (disparity) and intra-individual (asymmetry) variation were compared. All wild mice shared a similar inner ear morphology, in contrast to the important divergence of the Swiss strain. The release of selection in the laboratory strain obviously allowed for an important and rapid drift in the otherwise conserved structure. Shared traits between the inner ear of the lab strain and domestic pigs suggested a common response to mobility reduction in captivity. The lab-bred offspring of wild mice also differed from their wild relatives, suggesting plastic response related to maternal locomotory behavior, since inner ear morphology matures before birth in mammals. The signature observed in lab-bred wild mice and the lab strain was however not congruent, suggesting that plasticity did not participate to the divergence of the laboratory strain. However, contrary to the expectation, wild mice displayed slightly higher levels of inter-individual variation than laboratory mice, possibly due to the higher levels of genetic variance within and among wild populations compared to the lab strain. Differences in fluctuating asymmetry levels were detected, with the laboratory strain occasionally displaying higher asymmetry scores than its wild relatives. This suggests that there may indeed be a release of selection and/or a decrease in developmental stability in the laboratory strain.

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.

Molar enamel-dentine junction shape of Pliobates cataloniae and other Iberian pliopithecoids.

The phylogenetic relationships of the small-bodied catarrhine Pliobates cataloniae (∼11.6 Ma, NE Iberian Peninsula) have been controversial since its original description. However, the recent report of additional dentognathic remains has supported its crouzeliid pliopithecoid status. Based on the available hypodigm, the molar enamel-dentine junction (EDJ) shape of P. cataloniae is compared with that of other pliopithecoids from the same basin as well as both extinct and extant hominoids to further evaluate its pliopithecoid affinities. We also quantitatively compare the EDJ shape among these taxa by means of landmark-based three-dimensional geometric morphometrics using principal component analysis (PCA), canonical variate analysis, and between-group PCA. Permutation tests are performed to test whether Pliobates variation exceeds that of extant hominoid genera. Results indicate that Pliobates is similar in molar EDJ shape to other pliopithecoids, particularly crouzeliids. The variation displayed by Pliobates upper molars is less marked at the EDJ level than at the outer enamel surface, probably owing to differential enamel wear and intraspecific differences in enamel thickness. Multivariate analyses of EDJ shape show that all pliopithecoids (including Pliobates) cluster together in the PCAs, canonical variate analyses, and between-group PCAs and occupy a different portion of the morphospaces from extinct and extant hominoids. Posterior and typicality probabilities strongly support the classification of Pliobates as a pliopithecoid, wheras permutation tests fail to reject the single-genus hypothesis for the P. cataloniae hypodigm. We conclude that P. cataloniae is a crouzeliid pliopithecoid, as recently supported by cladistic analyses of craniodental characters, and that previous cladistic results that supported a stem hominoid status are attributable to postcranial convergences with crown hominoids. Our results further highlight the potential of three-dimensional geometric morphometrics analyses of the EDJ shape for better informing fossil primate alpha-taxonomy by means of quantitatively testing hypotheses about tooth shape variation.