Less tautology, more biology? A comment on “high-density” morphometrics

A computerized facial approximation method for Homo sapiens based on facial soft tissue thickness depths and geometric morphometrics

Facial approximation (FA) provides a promising means of generating the possible facial appearance of a deceased person. It facilitates exploration of the evolutionary forces driving anatomical changes in ancestral humans and can capture public attention. Despite the recent progress made toward improving the performance of FA methods, a limited understanding of detailed quantitative craniofacial relationships between facial bone and soft tissue morphology may hinder their accuracy, and hence subjective experience and artistic interpretation are required. In this study, we explored craniofacial relationships among human populations based upon average facial soft tissue thickness depths (FSTDs) and covariations between hard and soft tissues of the nose and mouth using geometric morphometrics. Furthermore, we proposed a computerized method to assign the learned craniofacial relationships to generate a probable facial appearance of Homo sapiens, reducing human intervention. A smaller resemblance comparison (an average Procrustes distance was 0.0258 and an average Euclidean distance was 1.79 mm) between approximated and actual faces and a greater recognition rate (91.67%) tested by a face pool indicated that average dense FSTDs contributed to raising the accuracy of approximated faces. Results of partial least squares (PLS) analysis showed that nasal and oral hard tissues have an effect on their soft tissues separately. However, relatively weaker RV correlations (<0.4) and greater approximation errors suggested that we need to be cautious about the accuracy of the approximated nose and mouth soft tissue shapes from bony structures. Overall, the proposed method can facilitate investigations of craniofacial relationships and potentially improve the reliability of the approximated faces for use in numerous applications in forensic science, archaeology, and anthropology.

A landmarking protocol for geometric morphometric analysis of squamate endocasts

Landmark-based geometric morphometrics is widely used to study the morphology of the endocast, or internal mold of the braincase, and the diversity associated with this structure across vertebrates. Landmarks, as the basic unit of such methods, are intended to be points of correspondence, selected depending on the question at hand, whose proper definition is essential to guarantee robustness and reproducibility of results. In this study, 20 landmarks are defined to provide a framework to analyze the morphological variability in squamate endocasts. Ten species representing a cross-section of the diversity of Squamata from both phylogenetic and ecological (i.e., habitat) perspectives were considered, to select landmarks replicable throughout the entire clade, regardless of the degree of neuroanatomical resolution of the endocast. To assess the precision, accuracy, and repeatability of these newly defined landmarks, both intraobserver and interobserver error were investigated. Estimates of measurement error show that most of the landmarks established here are highly replicable, and preliminary results suggest that they capture aspects of endocast shape related to both phylogenetic and ecologic signals. This study provides a basis for further examinations of squamate endocast disparity using landmark-based geometric morphometrics.

Normal human craniofacial growth and development from 0 to 4 years

Knowledge of human craniofacial growth (increase in size) and development (change in shape) is important in the clinical treatment of a range of conditions that affects it. This study uses an extensive collection of clinical CT scans to investigate craniofacial growth and development over the first 48 months of life, detail how the cranium changes in form (size and shape) in each sex and how these changes are associated with the growth and development of various soft tissues such as the brain, eyes and tongue and the expansion of the nasal cavity. This is achieved through multivariate analyses of cranial form based on 3D landmarks and semi-landmarks and by analyses of linear dimensions, and cranial volumes. The results highlight accelerations and decelerations in cranial form changes throughout early childhood. They show that from 0 to 12 months, the cranium undergoes greater changes in form than from 12 to 48 months. However, in terms of the development of overall cranial shape, there is no significant sexual dimorphism in the age range considered in this study. In consequence a single model of human craniofacial growth and development is presented for future studies to examine the physio-mechanical interactions of the craniofacial growth.

A Comparison of Semilandmarking Approaches in the Analysis of Size and Shape

Often, few landmarks can be reliably identified in analyses of form variation and covariation. Thus, ‘semilandmarking’ algorithms have increasingly been applied to surfaces and curves. However, the locations of semilandmarks depend on the investigator’s choice of algorithm and their density. In consequence, to the extent that different semilandmarking approaches and densities result in different locations of semilandmarks, they can be expected to yield different results concerning patterns of variation and co-variation. The extent of such differences due to methodology is, as yet, unclear and often ignored. In this study, the performance of three landmark-driven semilandmarking approaches is assessed, using two different surface mesh datasets (ape crania and human heads) with different degrees of variation and complexity, by comparing the results of morphometric analyses. These approaches produce different semilandmark locations, which, in turn, lead to differences in statistical results, although the non-rigid semilandmarking approaches are consistent. Morphometric analyses using semilandmarks must be interpreted with due caution, recognising that error is inevitable and that results are approximations. Further work is needed to investigate the effects of using different landmark and semilandmark templates and to understand the limitations and advantages of different semilandmarking approaches.

A Comparison of Semilandmarking Approaches in the Visualisation of Shape Differences

In landmark-based analyses of size and shape variation and covariation among biological structures, regions lacking clearly identifiable homologous landmarks are commonly described by semilandmarks. Different algorithms may be used to apply semilandmarks, but little is known about the consequences of analytical results. Here, we assess how different approaches and semilandmarking densities affect the estimates and visualisations of mean and allometrically scaled surfaces. The performance of three landmark-driven semilandmarking approaches is assessed using two different surface mesh datasets with different degrees of variation and complexity: adult human head and ape cranial surfaces. Surfaces fitted to estimates of the mean and allometrically scaled landmark and semilandmark configurations arising from geometric morphometric analyses of these datasets are compared between semilandmarking approaches and different densities, as well as with those from warping to landmarks alone. We find that estimates of surface mesh shape (i.e., after re-semilandmarking and then re-warping) made with varying numbers of semilandmarks are generally consistent, while the warping of surfaces using landmarks alone yields surfaces that can be quite different to those based on semilandmarks, depending on landmark coverage and choice of template surface for warping. The extent to which these differences are important depends on the particular study context and aims.

Ecological signal in the size and shape of marine amniote teeth

Amniotes have been a major component of marine trophic chains from the beginning of the Triassic to present day, with hundreds of species. However, inferences of their (palaeo)ecology have mostly been qualitative, making it difficult to track how dietary niches have changed through time and across clades. Here, we tackle this issue by applying a novel geometric morphometric protocol to three-dimensional models of tooth crowns across a wide range of raptorial marine amniotes. Our results highlight the phenomenon of dental simplification and widespread convergence in marine amniotes, limiting the range of tooth crown morphologies. Importantly, we quantitatively demonstrate that tooth crown shape and size are strongly associated with diet, whereas crown surface complexity is not. The maximal range of tooth shapes in both mammals and reptiles is seen in medium-sized taxa; large crowns are simple and restricted to a fraction of the morphospace. We recognize four principal raptorial guilds within toothed marine amniotes (durophages, generalists, flesh cutters and flesh piercers). Moreover, even though all these feeding guilds have been convergently colonized over the last 200 Myr, a series of dental morphologies are unique to the Mesozoic period, probably reflecting a distinct ecosystem structure.

Determination of Morphogeometric Patterns in Individuals with Total Mandibular Edentulism in the Interforaminal Region from Cone Beam Computed Tomography (CBCT) Scans: A Pilot Study

The aim of this study was to determine the morphogeometric patterns of the interforaminal region from cone beam computed tomography (CBCT) scans of individuals with total mandibular edentulism. CBCT images were obtained from 40 patients with total edentulism who are older (12 men and 28 women; average age of 69.5 ± 9.4 years) and who wore a non-implant-supported, lower, removable, total prosthesis. We conducted a two-dimensional (2D) morphogeometric analysis of the Digital Imaging and Communication in Medicine (DICOM) files from the CBCT scans, and five equidistant cross sections were planned. For the three-dimensional (3D) morphogeometric analysis, standard triangular language (STL) files were obtained after segmentation of the interforaminal mandibular region, and four anatomical landmarks and their respective curves were digitized. The patterns among the shapes were determined using principal component analysis (PCA) on MorphoJ software (version 1.07a). The results of the 2D morphogeometric analyses for PCA of the interforaminal mandibular paramedian region were PC1 or elongated drop shape, 54.78%; PC2 or wineskin shape, 17.65%; PC3 or pear shape, 11.77%; and PC4 or eggplant shape, 5.71%, and those for PCA of the symphyseal region were PC1 or elongated drop shape, 62.13%; PC2 or ovoid shape, 11.64%; PC3 or triangular shape, 9.71%; and PC4 or tuber shape, 4.96%. The results of the 3D morphogeometric analyses for the interforaminal hemimandibular region were PC1, 59.83%; PC2, 10.39%; PC3, 7.67%; and PC4, 5.09%. This study provides relevant information for future clinical guidelines on prosthetics and implants, in addition to proposing the use of new technologies that support diagnosis and treatment in patients with edentulism.

Evolution of the bovid cranium: morphological diversification under allometric constraint

The role of environmental selection in generating novel morphology is often taken for granted, and morphology is generally assumed to be adaptive. Bovids (antelopes and relatives) are widely differentiated in their dietary and climatic preferences, and presumably their cranial morphologies are the result of adaptation to different environmental pressures. In order to test these ideas, we performed 3D geometric morphometric analyses on 141 crania representing 96 bovid species in order to assess the influence of both extrinsic (e.g. diet, habitat) and intrinsic (size, modularity) factors on cranial shape. Surprisingly, we find that bovid crania are highly clumped in morphospace, with a large number of ecologically disparate species occupying a very similar range of morphology clustered around the mean shape. Differences in shape among dietary, habitat, and net primary productivity categories are largely non-significant, but we found a strong interaction between size and diet in explaining shape. We furthermore found no evidence for modularity having played a role in the generation of cranial differences across the bovid tree. Rather, the distribution of bovid cranial morphospace appears to be mainly the result of constraints imposed by a deeply conserved size-shape allometry, and dietary diversification the result of adaptation of existing allometric pathways.

Shape-changing chains for morphometric analysis of 2D and 3D, open or closed outlines

Morphometrics is a multivariate technique for shape analysis widely employed in biological, medical, and paleoanthropological applications. Commonly used morphometric methods require analyzing a huge amount of variables for problems involving a large number of specimens or complex shapes. Moreover, the analysis results are sometimes difficult to interpret and assess. This paper presents a methodology to synthesize a shape-changing chain for 2D or 3D curve fitting and to employ the chain parameters in stepwise discriminant analysis (DA). The shape-changing chain is comprised of three types of segments, including rigid segments that have fixed length and shape, scalable segments with a fixed shape, and extendible segments with constant curvature and torsion. Three examples are presented, including 2D mandible profiles of fossil hominin, 2D leaf outlines, and 3D suture curves on infant skulls. The results demonstrate that the shape-changing chain has several advantages over common morphometric methods. Specifically, it can be applied to a wide range of 2D or 3D profiles, including open or closed curves, and smooth or serrated curves. Additionally, the segmentation of profiles is a flexible and automatic protocol that can consider both biological and geometric features, the number of variables obtained from the fitting results for statistical analysis is modest, and the chain parameters that characterize the profiles can have physical meaning.

Can morphotaxa be assessed with photographs? Estimating the accuracy of two-dimensional cranial geometric morphometrics for the study of threatened populations of African monkeys

The classification of most mammalian orders and families is under debate and the number of species is likely greater than currently recognized. Improving taxonomic knowledge is crucial, as biodiversity is in rapid decline. Morphology is a source of taxonomic knowledge, and geometric morphometrics applied to two dimensional (2D) photographs of anatomical structures is commonly employed for quantifying differences within and among lineages. Photographs are informative, easy to obtain, and low cost. 2D analyses, however, introduce a large source of measurement error when applied to crania and other highly three dimensional (3D) structures. To explore the potential of 2D analyses for assessing taxonomic diversity, we use patas monkeys (Erythrocebus), a genus of large, semi-terrestrial, African guenons, as a case study. By applying a range of tests to compare ventral views of adult crania measured both in 2D and 3D, we show that, despite inaccuracies accounting for up to one-fourth of individual shape differences, results in 2D almost perfectly mirror those in 3D. This apparent paradox might be explained by the small strength of covariation in the component of shape variance related to measurement error. A rigorous standardization of photographic settings and the choice of almost coplanar landmarks are likely to further improve the correspondence of 2D to 3D shapes. 2D geometric morphometrics is, thus, appropriate for taxonomic comparisons of patas ventral crania. Although it is too early to generalize, our results corroborate similar findings from previous research in mammals, and suggest that 2D shape analyses are an effective heuristic tool for morphological investigation of small differences.

Cranial shape diversification in horses: variation and covariation patterns under the impact of artificial selection

The potential of artificial selection to dramatically impact phenotypic diversity is well known. Large-scale morphological changes in domestic species, emerging over short timescales, offer an accelerated perspective on evolutionary processes. The domestic horse (Equus caballus) provides a striking example of rapid evolution, with major changes in morphology and size likely stemming from artificial selection. However, the microevolutionary mechanisms allowing to generate this variation in a short time interval remain little known. Here, we use 3D geometric morphometrics to quantify skull morphological diversity in the horse, and investigate modularity and integration patterns to understand how morphological associations contribute to cranial evolvability in this taxon. We find that changes in the magnitude of cranial integration contribute to the diversification of the skull morphology in horse breeds. Our results demonstrate that a conserved pattern of modularity does not constrain large-scale morphological variations in horses and that artificial selection has impacted mechanisms underlying phenotypic diversity to facilitate rapid shape changes. More broadly, this study demonstrates that studying microevolutionary processes in domestic species produces important insights into extant phenotypic diversity.

Quantifying bite force in coexisting tayassuids and feral suids: a comparison between morphometric functional proxies and in vivo measurements

Background

Measuring mammals’ bite force in laboratory conditions is not a simple task, let alone on wild medium-sized mammals in the field. Thus, morphometric-proxies are usually used to infer morphofunctional properties of musculoskeletal features. For instance, the study of bite force-indexes suggests that different capacities to crack food items reduce the competition between coexistent collared and white-lipped peccaries (Pecari tajacu and Tayassu pecari). The presence of exotic feral hogs (Sus scrofa) in peccaries’ endemic areas gives rise to new ecological interactions between them. An example is the Brazilian Pantanal wetland, where ecomorphological mechanisms may play a role in their ecological relations. Taking this scenario as a case of study, we aimed to verify if the morphometric-proxies are de facto reliable tools, by comparing bite forces-indexes with the in vivo bite forces of these species.

Methods

We captured 21 collared and white-lipped peccaries and feral hogs in the Brazilian Pantanal to assess their bite force at first molar. The Bite Force Measuring Tube (BiTu) is a robust and simple mechanical device designed to be used in field conditions. Only 11 individuals successfully bit the BiTu before being released. Their body measurements were compared and correlated with their bite force. The in vivo bite forces were compared with bite force-indexes of two papers based on independent morphometric methods and datasets: Sicuro & Oliveira (2002) used classic morphometrics to infer the bite forces of these three species in the Brazilian Pantanal, and Hendges et al. (2019) used geometric morphometrics to compare bite forces-indexes and feeding habits of the extant peccary species. The results of all species were standardized (Z-curves) according to each method. Doing so, we obtained comparable dimensionless comparable values but maintaining the differences between them.

Results

The morphometric-proxies-based studies presented similar results: collared peccaries present weaker bites than white-lipped peccaries and feral hogs, while these two species presented no significant differences in their bite force-indexes. The in vivo bite force results suggest the same relations predicted by the morphometric models, including the high variation among the feral hogs. We found a significant correlation between the individuals’ weight (kg) and their actual bite force (N) but no significant correlations with the head length.

Conclusions

The BiTu proved to be a functional and low-cost tool to measure bite force in field conditions. The in vivo results presented a good correspondence with the predictions based on morphometric-proxies by Sicuro & Oliveira (2002) and Hendges et al. (2019). The results denote that these studies succeed in capturing the biomechanical signal of the three species’ skull-jaw systems. This empirical validation confirms that these morphometric-proxies analyses are reliable methods to ecomorphological and evolutionary inferences.

Automated description of the mandible shape by deep learning

PURPOSE

The shape of the mandible has been analyzed in a variety of fields, whether to diagnose conditions like osteoporosis or osteomyelitis, in forensics, to estimate biological information such as age, gender, and race or in orthognathic surgery. Although the methods employed produce encouraging results, most rely on the dry bone analyses or complex imaging techniques that, ultimately, hamper sample collection and, as a consequence, the development of large-scale studies. Thus, we proposed an objective, repeatable, and fully automatic approach to provide a quantitative description of the mandible in orthopantomographies (OPGs).

METHODS

We proposed the use of a deep convolutional neural network (CNN) to localize a set of landmarks of the mandible contour automatically from OPGs. Furthermore, we detailed four different descriptors for the mandible shape to be used for a variety of purposes. This includes a set of linear distances and angles calculated from eight anatomical landmarks of the mandible, the centroid size, the shape variations from the mean shape, and a group of shape parameters extracted with a point distribution model.

RESULTS

The fully automatic digitization of the mandible contour was very accurate, with a mean point to the curve error of 0.21 mm and a standard deviation comparable to that of a trained expert. The combination of the CNN and the four shape descriptors was validated in the well-known problems of forensic sex and age estimation, obtaining 87.8% of accuracy and a mean absolute error of 1.57 years, respectively.

CONCLUSION

The methodology proposed, including the shape model, can be valuable in any field that requires a quantitative description of the mandible shape and a visual representation of its changes such as clinical practice, surgery management, dental research, or legal medicine.

Sexual differences in human cranial morphology: Is one sex more variable or one region more dimorphic?

The quantification of cranial sexual dimorphism (CSD) among modern humans is relevant in evolutionary studies of morphological variation and in a forensic context. Despite the abundance of quantitative studies of CSD, few have specifically examined intra-sex variability. Here we quantify CSD in a geographically homogeneous sample of adult crania, which includes Italian individuals from the 19th and 20th centuries. Cranial morphology is described with 92 3D landmarks analyzed using Procrustean geometric morphometrics (PGMM). Size and shape variables are used to compare morphological variance between sexes in the whole cranium and four individual regions. The same variables, plus Procrustes form, are used to quantify average sex differences and explore classification accuracy. Our results indicate that: (a) as predicted by Wainer's rule, males present overall more variance in size and shape, albeit this is statistically significant only for total cranial size; (b) differences between sexes are dominated by size and to a lesser extent by Procrustes form; (c) shape only accounts for a minor proportion of variance; (d) the cranial base shows almost no dimorphism for shape; and (e) facial Procrustes form is the most accurate predictor of skeletal sex. Overall, this study suggests developmental factors underlying differences in CSD among cranial regions; stresses the need for population-specific models that describe craniofacial variation as the basis for models that facilitate the estimation of sex in unidentified skeletal remains; and provides one of the first confirmations of "Wainer's rule" in relation to sexual dimorphism in mammals specific to the human cranium.

Optimizing digitalization effort in morphometrics

Quantifying phenotypes is a common practice for addressing questions regarding morphological variation. The time dedicated to data acquisition can vary greatly depending on methods and on the required quantity of information. Optimizing digitization effort can be done either by pooling datasets among users, by automatizing data collection, or by reducing the number of measurements. Pooling datasets among users is not without risk since potential errors arising from multiple operators in data acquisition prevent combining morphometric datasets. We present an analytical workflow to estimate within and among operator biases and to assess whether morphometric datasets can be pooled. We show that pooling and sharing data requires careful examination of the errors occurring during data acquisition, that the choice of morphometric approach influences amount of error, and that in some cases pooling data should be avoided. The demonstration is based on a worked example (Sus scrofa teeth) using a combinations of 18 morphometric approaches and datasets for which we identified and quantified several potential sources of errors in the workflow. We show that it is possible to estimate the analytical power of a study using a small subset of data to select the best morphometric protocol and to optimize the number of variables necessary for analysis. In particular, we focus on semi-landmarks, which often produce an inflation of variables in contrast to the number of available observations use in statistical testing. We show how the workflow can be used for optimizing digitization efforts and provide recommendations for best practices in error management.