skip to main content


Title: Comparing and combining sliding semilandmarks and weighted spherical harmonics for shape analysis
Abstract

Quantifying morphological variation is critical for conducting anatomical research. Three‐dimensional geometric morphometric (3D GM) landmark analyses quantify shape using homologous Cartesian coordinates (landmarks). Setting up a high‐density landmark set and placing it on all specimens, however, can be a time‐consuming task. Weighted spherical harmonics (SPHARM) provides an alternative method for analyzing the shape of such objects. Here we compare sliding semilandmark and SPHARM analyses of the calcaneus ofGorilla gorilla gorilla(n = 20),Pan troglodytes troglodytes(n = 20), andHomo sapiens(n = 20) to determine whether the SPHARM and sliding semilandmark analyses capture comparable levels of shape variation. We also compare both the sliding semilandmark and SPHARM analyses to a novel combination of the two methods, here termed SPHARM–sliding. In SPHARM–sliding, the vertices of the surface models produced from the SPHARM analysis (that are the same in number and relative location) are used as the starting landmark positions for a sliding semilandmark analysis. Calcaneal shape variation quantified by all three analyses was summarized using separate principal components analyses. Results were compared using the root mean square (RMS) and maximum distance between surface models of species averages scaled (up) to centroid size created from each analysis. The average RMS was 0.23 mm between sliding semilandmark and SPHARM average surface models, 0.19 mm between SPHARM and SPHARM sliding average surface models, and 0.22 mm between sliding semilandmark and SPHARM sliding average surface models. Although results indicate that all three analyses are comparable methods for 3D shape analysis, there are advantages and disadvantages to each. While the SPHARM analysis is less time‐intensive, it is unable to capture the same level of detail around the sharp edges of articular facets on average surface models as the sliding semilandmark analysis. The SPHARM analysis also does not allow for individual articular facets to be analyzed in isolation. SPHARM–sliding, however, captures the same level of detail as the sliding semilandmark analysis, and (as in the sliding semilandmark analysis) allows for the evaluation of individual portions of bone. SPHARM is a comparable method to a 3D GM analysis for small, irregularly shaped bones, such as the calcaneus, and SPHARM–sliding allows for an expedited set up process for a sliding semilandmark analysis.

 
more » « less
NSF-PAR ID:
10364058
Author(s) / Creator(s):
 ;  ;  
Publisher / Repository:
Wiley-Blackwell
Date Published:
Journal Name:
Journal of Anatomy
Volume:
240
Issue:
4
ISSN:
0021-8782
Page Range / eLocation ID:
p. 678-687
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. The foot plays a prominent role in weight-bearing suggesting it may reflect locomotor variation. Despite the immense amount of foot research, the calcaneus has been relatively understudied. Here we analyzed the entire calcaneal shape of Gorilla gorilla gorilla (n=41), Gorilla beringei graueri (n=17) and Gorilla beringei beringei (n=8) to understand how morphology relates to locomotor behavior. Calcanei were surface scanned and external shape analyzed using a three-dimensional geometric morphometric sliding semilandmark analysis. Semilandmarks were slid to minimize the bending energy of the thin plate spline interpolation function relative to the updated Procrustes average. Generalized Procrustes Analysis was used to align landmark configurations and shape variation was summarized using a principal components analysis. Procrustes distances between species were calculated and resampling statistics were run to test for group differences. All subspecies demonstrate statistically different morphologies (p<0.005 for pairwise comparisons). G. b. graueri separates from other subspecies based on posterolateral morphology, with G. b. graueri demonstrating an elongated peroneal trochlea, and thus more bone superiorly than G. g. gorilla. Compared to G. b. beringei, G. b. graueri has less bone inferiorly near the tuberosity. Cuboid and posterior talar facet shapes correlate with arboreality. G. b. beringei (most terrestrial) has a flatter cuboid facet and a more transversely oriented/relatively smaller posterior talar facet than G. g. gorilla (most arboreal) and G. b. graueri represents an intermediate morphology. These differences demonstrate a relationship between calcaneal shape and locomotor behavior and suggest that G. b. graueri may load its foot differently from the other subspecies. This project was supported by NSF grant # BCS - 1824630. 
    more » « less
  2. Abstract Objectives

    Craniofacial fluctuating asymmetry (FA) refers to the random deviations from symmetry exhibited across the craniofacial complex and can be used as a measure of developmental instability for organisms with bilateral symmetry. This article addresses the lack of data on craniofacial FA in nonhuman primates by analyzing FA magnitude and variation in chimpanzees, gorillas, and macaques. We offer a preliminary investigation into how FA, as a proxy for developmental instability, varies within and among nonhuman primates.

    Materials and Methods

    We generated 3D surface models of 121 crania fromPan troglodytes troglodytes,Gorilla gorilla gorilla, andMacaca fascicularis fascicularis. Using geometric morphometric techniques, the magnitude of observed FA was calculated and compared for each individual, sex, and taxon, along with the variation of FA across cranial regions and for each bilateral landmark.

    Results

    Gorillas and macaques exhibited higher and more similar magnitudes of FA to each other than either taxon did to chimpanzees; variation in magnitude of FA followed this same trend. No significant differences were detected between sexes using pooled data across species, but sex did influence FA magnitude within taxa in gorillas. Further, variation in FA variance across cranial regions and by landmark was not distributed in any particular pattern.

    Conclusion

    Possible environmentally induced causes for these patterns of FA magnitude include differences in growth rate and physiological stress experienced during life. Developmental stability may be greatest in chimpanzees in this sample. Additionally, these results point to appropriate landmarks for future FA analyses and may help suggest more urgent candidate taxa for conservation efforts.

     
    more » « less
  3. Abstract The field of comparative morphology has entered a new phase with the rapid generation of high-resolution three-dimensional (3D) data. With freely available 3D data of thousands of species, methods for quantifying morphology that harness this rich phenotypic information are quickly emerging. Among these techniques, high-density geometric morphometric approaches provide a powerful and versatile framework to robustly characterize shape and phenotypic integration, the covariances among morphological traits. These methods are particularly useful for analyses of complex structures and across disparate taxa, which may share few landmarks of unambiguous homology. However, high-density geometric morphometrics also brings challenges, for example, with statistical, but not biological, covariances imposed by placement and sliding of semilandmarks and registration methods such as Procrustes superimposition. Here, we present simulations and case studies of high-density datasets for squamates, birds, and caecilians that exemplify the promise and challenges of high-dimensional analyses of phenotypic integration and modularity. We assess: (1) the relative merits of “big” high-density geometric morphometrics data over traditional shape data; (2) the impact of Procrustes superimposition on analyses of integration and modularity; and (3) differences in patterns of integration between analyses using high-density geometric morphometrics and those using discrete landmarks. We demonstrate that for many skull regions, 20–30 landmarks and/or semilandmarks are needed to accurately characterize their shape variation, and landmark-only analyses do a particularly poor job of capturing shape variation in vault and rostrum bones. Procrustes superimposition can mask modularity, especially when landmarks covary in parallel directions, but this effect decreases with more biologically complex covariance patterns. The directional effect of landmark variation on the position of the centroid affects recovery of covariance patterns more than landmark number does. Landmark-only and landmark-plus-sliding-semilandmark analyses of integration are generally congruent in overall pattern of integration, but landmark-only analyses tend to show higher integration between adjacent bones, especially when landmarks placed on the sutures between bones introduces a boundary bias. Allometry may be a stronger influence on patterns of integration in landmark-only analyses, which show stronger integration prior to removal of allometric effects compared to analyses including semilandmarks. High-density geometric morphometrics has its challenges and drawbacks, but our analyses of simulated and empirical datasets demonstrate that these potential issues are unlikely to obscure genuine biological signal. Rather, high-density geometric morphometric data exceed traditional landmark-based methods in characterization of morphology and allow more nuanced comparisons across disparate taxa. Combined with the rapid increases in 3D data availability, high-density morphometric approaches have immense potential to propel a new class of studies of comparative morphology and phenotypic integration. 
    more » « less
  4. ABSTRACT

    The third metacarpal has been a focus of study when examining questions surrounding early hominin locomotion since this bone is adapted to the diverse range of positional behaviors performed by extant hominoids. The shape of this bone is potentially under strong selective pressure related to the biomechanical demands of terrestrial knuckle‐walking, arboreal clambering, and brachiation performed by extant hominoids since the hand directly interacts with the substrate during the performance of these movements. The objective of the present study was to explore shape variation of the third metacarpal and examine how different parts of the bone discriminated between hominoid genera that perform these different locomotor behaviors. In addition to examining general interspecies variation, shape analysis was applied to testing the knuckle‐walking hypothesis for human evolution. Fourteen 3D landmark coordinates were collected on hominoid third metacarpals, and principal component analysis and Procrustes distances were used to examine metacarpal shape. Comparable measurements were collected on fossilized third metacarpals ofAustralopithecus afarensisas an early hominin test case for examining the knuckle‐walking hypothesis. Analyses that included landmarks collected on both ends of the bone distinguished humans from great apes and presented a strong functional signal related to suspensory locomotion among nonhuman hominoids, whereas the distal articular surface provided the strongest knuckle‐walking signal. The shapes ofAustralopithecus afarensismetacarpals examined in the current study did not provide evidence for a trajectory of shape change in early hominin evolution that started from a metacarpal adapted for terrestrial knuckle‐walking. Anat Rec, 302:983–998, 2019. © 2018 Wiley Periodicals, Inc.

     
    more » « less
  5. Abstract

    Current approaches to quantify phalangeal curvature assume that the long axis of the bone's diaphysis approximates the shape of a portion of a circle (included angle method) or a parabola (second‐degree polynomial method). Here we developed, tested, and employed an alternative geometric morphometrics‐based (GM) approach to quantify diaphysis shape of proximal phalanges in humans, apes and monkeys with diverse locomotor behaviors. One hundred landmarks of the central longitudinal axis were extracted from 3D surface models and analyzed using 2DGM methods, including generalized Procrustes analyses. Principal components analyses were performed and PC1 scores (>80% of variation) represented the dorsopalmar shape of the bone's central longitudinal axis and separated taxa consistently and in accord with known locomotor behavioral profiles. The most suspensory taxa, including orangutans, hylobatids and spider monkeys, had significantly lower PC1 scores reflecting the greatest amounts of phalangeal curvature. In contrast, bipedal humans and the quadrupedal cercopithecoid monkeys sampled (baboons, proboscis monkeys) exhibited significantly higher PC1 scores reflecting flatter phalanges. African ape (gorillas, chimpanzees and bonobos) phalanges fell between these two extremes and were not significantly different from each other. PC1 scores were significantly correlated with both included angle and theacoefficient of a second‐degree polynomial calculated from the same landmark dataset, but had a significantly higher correlation with included angles. Our alternative approach for quantifying diaphysis shape of proximal phalanges to investigate dorsopalmar curvature is replicable and does not assume a priori either a circle or parabola model of shape, making it an attractive alternative compared with existing methodologies.

     
    more » « less