A Study on the Asymmetry of the Human Left and Right Pubic Symphyseal Surfaces Using High-Definition Data Capture and Computational Shape Methods.

The pubic symphysis is among the most commonly used bilateral age indicators. Because of potential differences between right and left sides, it is necessary to investigate within-individual asymmetry, which can inflate age estimation error. This study uses 3D laser scans of paired pubic symphyses for 88 documented White males. Scan data are analyzed by numerical shape algorithms, proposed as an alternative to traditional visual assessment techniques. Results are used to quantify the within-individual asymmetry, evaluating if one side produces a better age-estimate. Relationships between the asymmetry and advanced age, weight, and stature are examined. This analysis indicates that the computational, shape-based techniques are robust to asymmetry (>80% of paired differences are within 10 years and >90% are within 15 years). For notably more asymmetric cases, differences in estimates are not associated with life history factors. Based on this study, either side can be used for age-at-death estimation by the computational methods.

Morphometric analysis of shape differences in Windover and Point Hope archaic human mandibles.

OBJECTIVES: The mandible can provide valuable information on both the life history and genetic makeup of Archaic human populations. The following analysis tests two hypotheses: (a) that there are significant differences in morphology in mandibular shape between the genders amongst Archaic North American Homo sapiens and (b) that there is a significant difference in variance in mandibular shape between Archaic Windover and Point Hope. MATERIALS AND METHODS: A sample made from mandible specimens taken from both populations is subjected to Principal Component Analyses (PCA). The component scores from the PCAs are subjected to both a Multivariate Analysis of Covariance (mancova) and a general Multivariate Analysis of Variance (manova) to determine whether significant differences in variance exist between the sexes and the populations. RESULTS: The mancova found that there are no significant interactions between the PC scores in population, sex, or size. Significant differences in variance were found between males and females and between the Windover and Point Hope populations. CONCLUSIONS: Differences in variance observed between the populations are suspected to be due to differences in subsistence strategies and possibly non-masticatory utilizations of teeth. Differences in variance between the genders are suspected to be genetic in origin.

Testing Reliability of the Computational Age-At-Death Estimation Methods between Five Observers Using Three-Dimensional Image Data of the Pubic Symphysis.

In an effort to standardize data collection and analysis in age estimation, a series of computational methods utilizing high-dimensional image data of the age indicator have recently been proposed as an alternative to subjective visual, trait-to-phase matching techniques. To systematically quantify the reproducibility of such methods, we investigate the intrascan variability and within- and between-observer reliability in initial scan data capturing and editing using 3D laser scans of the Suchey-Brooks pubic symphysis casts and five shape-based computational methods. Our results show that (i) five observers with various training background and experience levels edited the scans consistently for all three trials and the derived shape measures and age estimates were in excellent agreement among observers, and (ii) the computational methods are robust to a measured degree of scan trimming error. This study supports the application of computational methods to 3D laser scanned images for reliable age-at-death estimation, with reduced subjectivity.

Age-at-Death Estimation for Modern Populations in Mexico and Puerto Rico through the Use of 3D Laser Scans of the Pubic Symphysis.

Reliable age-at-death estimates from the adult skeleton are of fundamental importance in forensic anthropology, because it contributes to the identity parameters used in a medicolegal death investigation. However, reliable estimates are difficult because many traditional aging methods depend on a set of population-specific criteria derived from individuals of European and African descent. The absence of information on the potential differences in the aging patterns of underrepresented, especially Latinx, populations may hinder our efforts to produce useful age-at-death estimates. In response to these concerns, this study explores the utility of currently available aging techniques and whether population-specific aging methods among Latinx groups are needed. The authors obtained data from two skeletal collections representing modern individuals of Mexican and Puerto Rican origin. They examined five newly developed computational shape-based techniques using 3D laser scans of the pubic symphysis and one traditional bone-to-phase technique. A validation test of all computational and traditional methods was implemented, and new population-specific equations using the computational algorithms were generated and tested against a subsample. Results suggest that traditional and computational aging techniques applied to the pubic symphysis perform best with individuals within 35-45 years of age. Levels of bias and inaccuracy increase as chronological age increases, with overestimation of individuals younger than 35 years and underestimation of individuals older than 45 years. New regression models provided error rates comparable to, and in some occasions outperformed, the original computational models developed on white American males, but age estimates did not significantly improve. This study shows that population-specific models do not necessarily improve age estimates in Latinx samples. Results do suggest that computational methods can ultimately outperform the Suchey-Brooks method and provide improved objectivity when estimating age at death in Latinx samples.

The Development and Use of Computational Tools in Forensic Science.

Modern computational resources make available a rich tool kit of statistical methods that can be applied to forensic questions. This tool kit is built on the foundation of statistical developments dating back to the 19th century. To fully and effectively exploit these developments, both the makers and users of software must be keenly aware of the quality, that is, the accuracy and precision, of the data being modeled or analyzed, and end users must be sufficiently familiar with the underlying theory to understand the process and results of any analysis or software they use. This is especially important for medicolegal personnel who might be called upon to testify in a court of law and be subject to cross-examination. With respect to the development of computational tools, it is increasingly important that they be made available as open-source code to avoid the pitfalls of commercial software support and the potential dependence of end users on orphaned software.

A Computational Framework for Age-at-Death Estimation from the Skeleton: Surface and Outline Analysis of 3D Laser Scans of the Adult Pubic Symphysis.

In forensic anthropology, age-at-death estimation typically requires the macroscopic assessment of the skeletal indicator and its association with a phase or score. High subjectivity and error are the recognized disadvantages of this approach, creating a need for alternative tools that enable the objective and mathematically robust assessment of true chronological age. We describe, here, three fully computational, quantitative shape analysis methods and a combinatory approach that make use of three-dimensional laser scans of the pubic symphysis. We report a novel age-related shape measure, focusing on the changes observed in the ventral margin curvature, and refine two former methods, whose measures capture the flatness of the symphyseal surface. We show how we can decrease age-estimation error and improve prior results by combining these outline and surface measures in two multivariate regression models. The presented models produce objective age-estimates that are comparable to current practices with root-mean-square-errors between 13.7 and 16.5 years.

A Landmark-Free Method for Three-Dimensional Shape Analysis.

BACKGROUND: The tools and techniques used in morphometrics have always aimed to transform the physical shape of an object into a concise set of numerical data for mathematical analysis. The advent of landmark-based morphometrics opened new avenues of research, but these methods are not without drawbacks. The time investment required of trained individuals to accurately landmark a data set is significant, and the reliance on readily-identifiable physical features can hamper research efforts. This is especially true of those investigating smooth or featureless surfaces. METHODS: In this paper, we present a new method to perform this transformation for data obtained from high-resolution scanning technology. This method uses surface scans, instead of landmarks, to calculate a shape difference metric analogous to Procrustes distance and perform superimposition. This is accomplished by building upon and extending the Iterative Closest Point algorithm. We also explore some new ways this data can be used; for example, we can calculate an averaged surface directly and visualize point-wise shape information over this surface. Finally, we briefly demonstrate this method on a set of primate skulls and compare the results of the new methodology with traditional geometric morphometric analysis.

An enhanced computational method for age-at-death estimation based on the pubic symphysis using 3D laser scans and thin plate splines.

OBJECTIVES: The pubic symphysis is frequently used to estimate age-at-death from the adult skeleton. Assessment methods require the visual comparison of the bone morphology against age-informative characteristics that represent a series of phases. Age-at-death is then estimated from the age-range previously associated with the chosen phase. While easily executed, the "morphoscopic" process of feature-scoring and bone-to-phase-matching is known to be subjective. Studies of method and practitioner error demonstrate a need for alternative tools to quantify age-progressive change in the pubic symphysis. This article proposes a more objective, quantitative method that analyzes three-dimensional (3D) surface scans of the pubic symphysis using a thin plate spline algorithm (TPS). MATERIALS AND METHODS: This algorithm models the bending of a flat plane to approximately match the surface of the bone and minimizes the bending energy required for this transformation. Known age-at-death and bending energy were used to construct a linear model to predict age from observed bending energy. This approach is tested with scans from 44 documented white male skeletons and 12 casts. RESULTS: The results of the surface analysis show a significant association (regression p-value = 0.0002 and coefficient of determination = 0.2270) between the minimum bending energy and age-at-death, with a root mean square error of ≈19 years. DISCUSSION: This TPS method yields estimates comparable to established methods but offers a fully integrated, objective and quantitative framework of analysis and has potential for use in archaeological and forensic casework.

Modeling Bone Surface Morphology: A Fully Quantitative Method for Age-at-Death Estimation Using the Pubic Symphysis.

The pubic symphysis is widely used in age estimation for the adult skeleton. Standard practice requires the visual comparison of surface morphology against criteria representing predefined phases and the estimation of case-specific age from an age range associated with the chosen phase. Known problems of method and observer error necessitate alternative tools to quantify age-related change in pubic morphology. This paper presents an objective, fully quantitative method for estimating age-at-death from the skeleton, which exploits a variance-based score of surface complexity computed from vertices obtained from a scanner sampling the pubic symphysis. For laser scans from 41 modern American male skeletons, this method produces results that are significantly associated with known age-at-death (RMSE = 17.15 years). Chronological age is predicted, therefore, equally well, if not, better, with this robust, objective, and fully quantitative method than with prevailing phase-aging systems. This method contributes to forensic casework by responding to medico-legal expectations for evidence standards.

Biomechanical implications of intraspecific shape variation in chimpanzee crania: moving toward an integration of geometric morphometrics and finite element analysis.

In a broad range of evolutionary studies, an understanding of intraspecific variation is needed in order to contextualize and interpret the meaning of variation between species. However, mechanical analyses of primate crania using experimental or modeling methods typically encounter logistical constraints that force them to rely on data gathered from only one or a few individuals. This results in a lack of knowledge concerning the mechanical significance of intraspecific shape variation that limits our ability to infer the significance of interspecific differences. This study uses geometric morphometric methods (GM) and finite element analysis (FEA) to examine the biomechanical implications of shape variation in chimpanzee crania, thereby providing a comparative context in which to interpret shape-related mechanical variation between hominin species. Six finite element models (FEMs) of chimpanzee crania were constructed from CT scans following shape-space Principal Component Analysis (PCA) of a matrix of 709 Procrustes coordinates (digitized onto 21 specimens) to identify the individuals at the extremes of the first three principal components. The FEMs were assigned the material properties of bone and were loaded and constrained to simulate maximal bites on the P(3) and M(2) . Resulting strains indicate that intraspecific cranial variation in morphology is associated with quantitatively high levels of variation in strain magnitudes, but qualitatively little variation in the distribution of strain concentrations. Thus, interspecific comparisons should include considerations of the spatial patterning of strains rather than focus only on their magnitudes.

The feeding biomechanics and dietary ecology of Paranthropus boisei.

The African Plio-Pleistocene hominins known as australopiths evolved derived craniodental features frequently interpreted as adaptations for feeding on either hard, or compliant/tough foods. Among australopiths, Paranthropus boisei is the most robust form, exhibiting traits traditionally hypothesized to produce high bite forces efficiently and strengthen the face against feeding stresses. However, recent mechanical analyses imply that P. boisei may not have been an efficient producer of bite force and that robust morphology in primates is not necessarily strong. Here we use an engineering method, finite element analysis, to show that the facial skeleton of P. boisei is structurally strong, exhibits a strain pattern different from that in chimpanzees (Pan troglodytes) and Australopithecus africanus, and efficiently produces high bite force. It has been suggested that P. boisei consumed a diet of compliant/tough foods like grass blades and sedge pith. However, the blunt occlusal topography of this and other species suggests that australopiths are adapted to consume hard foods, perhaps including grass and sedge seeds. A consideration of evolutionary trends in morphology relating to feeding mechanics suggests that food processing behaviors in gracile australopiths evidently were disrupted by environmental change, perhaps contributing to the eventual evolution of Homo and Paranthropus.

Influence of edentulism on human orbit and zygomatic arch shape.

Edentulism, or tooth loss, seriously alters the appearance of the lower facial skeleton. The aim of this study was to determine if complete maxillary edentulism also impacts the curvature shape of the orbits and zygomatic arches in elderly adults. The study was conducted on 80 crania comprising two cross-sectional populations of elderly African- and European-Americans (60-80 years old). Forty of the crania possessed intact dentition; the remaining 40 exhibited complete edentulism with tooth socket resorption. Three-dimensional semilandmarks representing the curvature of the orbits and zygomatic arches were collected using a hand-held digitizer. Each craniofacial region's semilandmarks were aligned into a common coordinate system via generalized Procrustes superimposition. Regional variation in shape was explored via principal component analysis, multivariate analysis of variance, discriminant function analysis, cross-validation, and vector plots. Shape differences between the edentulous and dentate groups were detected in both the orbits (P = 0.0022) and zygomatic arches (P = 0.0026). Ancestry and sex differences were also identified in both regions. Orbit data correctly classified dentate crania 65% of the time and edentulous crania 72.5% of the time. Zygomatic arch data correctly classified 75% dentate and 60% of edentulous crania. The individual curves constituting each region also exhibited shape alteration with tooth loss, with the exception of the inferior zygomatic curve. Vector plots revealed patterns of superoinferior expansion, and medial and lateral recession depending on the region examined. These results suggest a relationship exists between maxillary edentulism and changes in the surrounding craniofacial structures.

A new method for the analysis of soft tissues with data acquired under field conditions.

Analyzing soft-tissue structures is particularly challenging due to the lack of homologous landmarks that can be reliably identified across time and specimens. This is particularly true when data are to be collected under field conditions. Here, we present a method that combines photogrammetric techniques and geometric morphometrics methods (GMM) to quantify soft tissues for their subsequent volumetric analysis. We combine previously developed methods for landmark data acquisition and processing with a custom program for volumetric computations. Photogrammetric methods are a particularly powerful tool for field studies as they allow for image acquisition with minimal equipment requirements and for the acquisition of the spatial coordinates of points (anatomical landmarks or others) from these images. For our method, a limited number of homologous landmarks, i.e., points that can be found on any specimen independent of space and time, and further distinctive points, which may vary over time, space and subject, are identified on two-dimensional photographs and their three-dimensional coordinates estimated using photogrammetric methods. The three-dimensional configurations are oriented by the spatial principal components (PCs) of the homologous points. Crucially, this last step orients the configuration such that x and y-information (PC1 and PC2 coordinates) constitute an anatomically-defined plane with the z-values (PC3 coordinate) in the direction of interest for volume computation. The z-coordinates are then used to estimate the volume of the tissue. We validate our method using a physical, geometric model of known dimensions and physical (wax) models designed to approximate perineal swellings in female macaques. To demonstrate the usefulness and potential of our method, we use it to estimate the volumes of Barbary macaque sexual swellings recorded in the field with video images. By analyzing both the artificial data and real monkey swellings, we validate our method's accuracy and illustrate its potential for application in important areas of biological research.

Viewpoints: diet and dietary adaptations in early hominins: the hard food perspective.

Recent biomechanical analyses examining the feeding adaptations of early hominins have yielded results consistent with the hypothesis that hard foods exerted a selection pressure that influenced the evolution of australopith morphology. However, this hypothesis appears inconsistent with recent reconstructions of early hominin diet based on dental microwear and stable isotopes. Thus, it is likely that either the diets of some australopiths included a high proportion of foods these taxa were poorly adapted to consume (i.e., foods that they would not have processed efficiently), or that aspects of what we thought we knew about the functional morphology of teeth must be wrong. Evaluation of these possibilities requires a recognition that analyses based on microwear, isotopes, finite element modeling, and enamel chips and cracks each test different types of hypotheses and allow different types of inferences. Microwear and isotopic analyses are best suited to reconstructing broad dietary patterns, but are limited in their ability to falsify specific hypotheses about morphological adaptation. Conversely, finite element analysis is a tool for evaluating the mechanical basis of form-function relationships, but says little about the frequency with which specific behaviors were performed or the particular types of food that were consumed. Enamel chip and crack analyses are means of both reconstructing diet and examining biomechanics. We suggest that current evidence is consistent with the hypothesis that certain derived australopith traits are adaptations for consuming hard foods, but that australopiths had generalized diets that could include high proportions of foods that were both compliant and tough.

Microwear, mechanics and the feeding adaptations of Australopithecus africanus.

Recent studies of dental microwear and craniofacial mechanics have yielded contradictory interpretations regarding the feeding ecology and adaptations of Australopithecus africanus. As part of this debate, the methods used in the mechanical studies have been criticized. In particular, it has been claimed that finite element analysis has been poorly applied to this research question. This paper responds to some of these mechanical criticisms, highlights limitations of dental microwear analysis, and identifies avenues of future research.

Regional shape change in adult facial bone curvature with age.

Life expectancies have increased dramatically over the last 100 years, affording greater opportunities to study the impact of age on adult craniofacial morphology. This article employs a novel application of established geometric morphometric methods to examine shape differences in adult regional facial bone curvature with age. Three-dimensional semilandmarks representing the curvature of the orbits, zygomatic arches, nasal aperture, and maxillary alveolar process were collected from a cross-sectional cranial sample of mixed sex and ancestry (male and female; African- and European-American), partitioned into three age groups (young adult = 18-39; middle-aged = 40-59 years; and elderly = 60+ years). Each facial region's semilandmarks were aligned into a common coordinate system via generalized Procrustes superimposition. Regional variation in shape was then explored via a battery of multivariate statistical techniques. Age-related shape differences were detected in the orbits, zygomatic arches, and maxillary alveolar process. Interactions between age, sex, and ancestry were also identified. Vector plots revealed patterns of superoinferior compression, lateral expansion, and posterior recession depending on the population/subpopulation, location, and age groups examined. These findings indicate that adult craniofacial curvature shape is not static throughout human life. Instead, age-related spatial modifications occur in various regions of the craniofacial skeleton. Moreover, these regional alterations vary not only through time, but across human populations and the sexes.

The structural rigidity of the cranium of Australopithecus africanus: implications for diet, dietary adaptations, and the allometry of feeding biomechanics.

Australopithecus africanus is an early hominin (i.e., human relative) believed to exhibit stress-reducing adaptations in its craniofacial skeleton that may be related to the consumption of resistant food items using its premolar teeth. Finite element analyses simulating molar and premolar biting were used to test the hypothesis that the cranium of A. africanus is structurally more rigid than that of Macaca fascicularis, an Old World monkey that lacks derived australopith facial features. Previously generated finite element models of crania of these species were subjected to isometrically scaled loads, permitting a direct comparison of strain magnitudes. Moreover, strain energy (SE) in the models was compared after results were scaled to account for differences in bone volume and muscle forces. Results indicate that strains in certain skeletal regions below the orbits are higher in M. fascicularis than in A. africanus. Moreover, although premolar bites produce von Mises strains in the rostrum that are elevated relative to those produced by molar biting in both species, rostral strains are much higher in the macaque than in the australopith. These data suggest that at least the midface of A. africanus is more rigid than that of M. fascicularis. Comparisons of SE reveal that the A. africanus cranium is, overall, more rigid than that of M. fascicularis during premolar biting. This is consistent with the hypothesis that this hominin may have periodically consumed large, hard food items. However, the SE data suggest that the A. africanus cranium is marginally less rigid than that of the macaque during molar biting. It is hypothesized that the SE results are being influenced by the allometric scaling of cranial cortical bone thickness.

The feeding biomechanics and dietary ecology of Australopithecus africanus.

The African Plio-Pleistocene hominins known as australopiths evolved a distinctive craniofacial morphology that traditionally has been viewed as a dietary adaptation for feeding on either small, hard objects or on large volumes of food. A historically influential interpretation of this morphology hypothesizes that loads applied to the premolars during feeding had a profound influence on the evolution of australopith craniofacial form. Here, we test this hypothesis using finite element analysis in conjunction with comparative, imaging, and experimental methods. We find that the facial skeleton of the Australopithecus type species, A. africanus, is well suited to withstand premolar loads. However, we suggest that the mastication of either small objects or large volumes of food is unlikely to fully explain the evolution of facial form in this species. Rather, key aspects of australopith craniofacial morphology are more likely to be related to the ingestion and initial preparation of large, mechanically protected food objects like large nuts and seeds. These foods may have broadened the diet of these hominins, possibly by being critical resources that australopiths relied on during periods when their preferred dietary items were in short supply. Our analysis reconciles apparent discrepancies between dietary reconstructions based on biomechanics, tooth morphology, and dental microwear.

Age and gender based biomechanical shape and size analysis of the pediatric brain.

Injuries caused by motor vehicle crashes (MVCs) are the leading cause of head injury and death for children in the United States. This study aims to describe the shape and size (morphologic) changes of the cerebrum, cerebellum, brainstem, and ventricles of the pediatric occupant to better predict injury and assess how these changes affect finite element model (FEM) response. To quantify morphologic differences in the brain, a Generalized Procrustes Analysis (GPA) with a sliding landmark method was conducted to isolate morphologic changes using magnetic resonance images of 63 normal subjects. This type of geometric morphometric analysis was selected for its ability to identify homologous landmarks on structures with few true landmarks and isolate the shape and size of the individuals studied. From the resulting landmark coordinates, the shape and size changes were regressed against age to develop a model describing morphologic changes in the pediatric brain as a function of age. The most statistically significant shape change was in the cerebrum with p-values of 0.00346 for males and 0.00829 for females. The age-based model explains over 80% of the variation in size in the cerebrum. Using size and shape models, affine transformations were applied to the SIMon FEM to determine differences in response given differences in size and size plus shape. The geometric centroid of the elements exceeding 15% strain was calculated and compared to the geometric centroid of the entire structure. Given the same Haversine pulse, the centroid location, a metric for the spatial distribution of the elements exceeding an injury threshold, varied based on which transformation was applied to the model. To assess the overall response of the model, three injury metrics were examined to determine the magnitude of the metrics each element sustained and the overall volume of elements that experienced that value. These results suggested that the overall response of the model was driven by the variation in size, with little variation due to changes in shape. This study demonstrates a new methodology to quantify the shape and size variation of the brain from infancy to adulthood. The use of the changes in shape and size when applied to a FEM suggests that there are differences in the spatial distribution of the elements that exceed a specific threshold based on shape but the overall volume of elements experiencing the specified magnitude was more dependent on the changes in the size of the model with little change due to shape.

Quantification of age-related shape change of the human rib cage through geometric morphometrics.

The aim of this study is to quantify patterns of age-related shape change in the human thorax using Procrustes superimposition. Landmarks (n=106) selected from anonymized computed tomography (CT) scans of 63 adult males free of skeletal pathology were used to describe the form of the rib cage. Multivariate linear regression was used to determine a relationship between landmark location and age. Linear and quadratic models were also investigated. A permutation test employing 1 x 10(5) random trials was used to assess the model significance for both model formulations. Linear relationships between the centroid size (CS) of a landmark set and the corresponding individual's height, weight, and BMI were conducted to enable scaling of the dimensionless results from the Procrustes analysis. A significance level of alpha=0.05 was used for all tests. The average age of the study subjects was 57.0+/-17.3 years. Complete landmark sets were obtained from most of the scans (44 of 63). The quadratic relationship between the age and landmark location was found to be significant (p=0.037), thereby establishing a relationship between the age and thoracic shape change. The linear relationship was mildly significant as well (p=0.073). Significant relationships between the centroid size of the dataset and subject weight, height and BMI were determined, with the best-correlated value being weight (p=0.002, R(2)=0.22). Landmark datasets calculated using the quadratic model exhibited shape change consistent with the clinical observations (increasing kyphosis and rounding of the thoracic cage). Procrustes superimposition represents a potential improvement in the approach used to generate computational models for injury biomechanics studies. The coefficients from the quadratic model are provided and can be used to generate the complete set of model landmark data points at a given age.