Secular change of sexually dimorphic cranial variables in Euro-Americans and Germans

The influence of biological relatedness on sexual dimorphism and sex classification based on external morphology of the frontal bone

The most significant sexual differences in the human skull are located in the upper third of the face (the frontal bone), which is a useful research object, mainly in combination with virtual anthropology methods. However, the influence of biological relatedness on sexual dimorphism and frontal bone variability remains unknown. This study was directed at sexual difference description and sex classification using the form and shape of the external surface of the frontal bones from a genealogically documented Central European osteological sample (nineteenth to twentieth centuries). The study sample consisted of 47 cranial CT images of the adult members of several branches of one family group over 4 generations. Three-dimensional virtual models of the frontal bones were analyzed using geometric morphometrics and multidimensional statistics. Almost the entire external frontal surface was significantly different between males and females, especially in form. Significant differences were also found between this related sample and an unrelated one. Sex estimation of the biologically related individuals was performed using the classification models developed on a sample of unrelated individuals from the recent Czech population (Čechová et al. in Int J Legal Med 133: 1285 1294, 2019), with a result of 74.46% and 63.83% in form and shape, respectively. Failure of this classifier was caused by the existence of typical traits found in the biologically related sample different from the usual manifestation of sexual dimorphism. This can be explained as due to the increased degree of similarity and the reduction of variability in biologically related individuals. The results show the importance of testing previously published methods on genealogical data.

TestDimorph: An R package for analysis of interpopulation sexual dimorphism differences using summary statistics

OBJECTIVE

The degree of sexual dimorphism in certain traits between males and females differ from one sample to another. Although trait differences by sex are often reported in bioanthropological research, few studies test for statistical significance or make raw data available. TestDimorph is the first R package dedicated to testing and comparing the degree of sexual dimorphism among different samples by leveraging summary statistics.

MATERIALS AND METHODS

We provide two approaches of analysis of inter-sample differences in degree of sexual dimorphism: univariate and multivariate for two or more samples. The methods follow upon publications primarily from the AJBA. Within-sex size variability between samples is compared using one-way ANOVA followed by control for multiple pairwise comparisons. In addition, we compute the overlapping area between the density functions of two normal distributions from the mixture intersection index or the non-overlapping area using the dissimilarity index as well as Hedges' g with inferential support using the 95% confidence interval. Finally, we use a multivariate analysis of differences in patterning of sexual dimorphism between samples.

RESULTS

We demonstrate various results from applying TestDimorph functions to data supplied with the package.

DISCUSSION

The package has many features including functionality for working with summary statistics, simulating data from summary statistics, and the extraction of summary statistics from raw data, so that the entire analysis can be performed through the package.

A Geometric Morphometric Study on Sexual Dimorphism in Viscerocranium

The level of sexual dimorphism manifested by human bones is an important factor for development of effective sex estimation methods. The aim of the study was to investigate the sexual dimorphism in the size and shape of the viscerocranium using geometric morphometric techniques. It also aimed to explore the sex differences in distinct viscerocranial regions and to establish the most dimorphic region with regard to size and shape. Computed tomography images of 156 males and 184 females were used in the study. Three-dimensional coordinates of 31 landmarks were acquired. Five landmark configurations were constructed from the viscerocranium and its orbital, nasal, maxillary, and zygomatic region. Generalized Procrustes superimposition, principal component analysis, and discriminant analysis were applied to each configuration. The significance of the sex differences in size and shape was assessed and significant differences were found in all configurations. The highest accuracy was obtained from both shape and size of the whole viscerocranium. Based on size only, the highest accuracy was achieved by the nasal region. The accuracy based on shape was generally low for all configurations, but the highest result was attained by the orbital region. Hence, size is a better sex discriminator than shape.

Age-related differences in cranial sexual dimorphism in contemporary Europe

Biomechanical load and hormonal levels tended to change just like the soft and skeletal tissue of the elderly with age. Although aging in both sexes shared common traits, it was assumed that there would be a reduction of sexual dimorphism in aged individuals. The main goals of this study were (1) to evaluate age-related differences in cranial sexual dimorphism during senescence, (2) to determine age-related differences in female and male skulls separately, and (3) to compare skull senescence in Czech and French adult samples as discussed by Musilová et al. (Forensic Sci Int 269:70-77, 2016). The cranial surface was analyzed using coherent point drift-dense correspondence analysis. The study sample consisted of 245 CT scans of heads from recent Czech (83 males and 59 females) and French (52 males and 51 females) individuals. Virtual scans in the age range from 18 to 92 years were analyzed using geometric morphometrics. The cranial form was significantly greater in males in all age categories. After size normalization, sexual dimorphism of the frontal, occipital, and zygomatic regions tended to diminish in the elderly. Its development during aging was caused by morphological changes in both female and male skulls but secular changes must also be taken into account. The most notable aging changes were the widening of the neurocranium and the retrusion of the face, including the forehead, especially after the age of 60 in both sexes. Sexual dimorphism was similar between the Czech and French samples but its age-related differences were partially different because of the population specificity. Cranial senescence was found to degrade the accuracy of sex classification (92-94%) in the range of 2-3%.

Sex estimation of the humerus: A geometric morphometric analysis in an adult sample

Sex estimation is the keystone for positive identification when an unidentified human body is recovered in forensic contexts. However, in complex death scenes such as mass disasters, the remains are often fleshed, mutilated, burned, and/or commingled. In situations such as these where it is not possible to analyze pelvis and/or cranium data, traditional metric and qualitative morphological methods on postcranial bones can yield unsatisfactory results. In such cases, geometric morphometric techniques offer an alternative to the analysis of both shape and size components of morphological variation that can be of great utility for sex estimation in forensic investigations. The study population consisted of 72 well-preserved adult humeri (40 males and 32 females; mean age of 62 years) that were photographed in standardized positions with landmarks located in four two-dimensional views of the humerus (anterior surface of the proximal epiphysis, and anterior, posterior and inferior surface of distal epiphysis). Principal components analysis, canonical variates analysis and discriminant analysis were applied. The data indicated that males and females were classified with low levels of accuracy (54.95-77.92% for males; 56.87-71.78% for females) based on shape variables. However, when the shape variable was combined with the centroid size, the levels of accuracy increased (81.86-94.92% for males; 84.08-94.88% for females). To obtain larger differences between males and females, it is necessary the combination of centroid size with shape variables; the shape of the humerus is insufficient to discriminate sex with accuracy.

Data mining for sex estimation based on cranial measurements

The aim of the present study is to develop effective and understandable classification models for sex estimation and to identify the most dimorphic linear measurements in adult crania by means of data mining techniques. Furthermore, machine learning models and models developed through logistic regression analysis are compared in terms of performance. Computed tomography scans of 393 adult individuals were used in the study. A landmark-based approach was applied to collect the metric data. The three-dimensional coordinates of 47 landmarks were acquired and used for calculation of linear measurements. Two datasets of cranial measurements were assembled, including 37standard measurements and 1081 interlandmark distances, respectively. Three data mining algorithms were applied: the rule induction algorithms JRIP and Ridor, and the decision tree algorithm J48. Two advanced attribute selection methods (Weka BestFirst and Weka GeneticSearch) were also used. The best accuracy result (91.9 %) was achieved by a set of rules learnt by the JRIP algorithm from the dataset constructed by application of the GeneticSearch selection algorithm to the dataset of standard cranial measurements. The set consisted of five rules including seven cranial measurements. Its accuracy was even better than the classification rates achieved by the logistic regression models. Concerning the second dataset of nonstandard measurements, the best accuracy (88.3 %) was obtained by using classification models learnt by two algorithms - JRIP with a dataset preprocessed by the BestFirst selection algorithm and Ridor with preprocessing by the GeneticSearch selection algorithm. Our experiments show that for the two datasets mentioned above the rule-based models contain smaller sets of rules with shorter lists of measurements and achieve better classification accuracy results in comparison with decision tree-based models.

Effects of Multi-Generational Soft Diet Consumption on Mouse Craniofacial Morphology

Variations in craniofacial morphology may arise as a result of adaptation to different environmental factors such as soft diet (SD), which lessens functional masticatory load. Prior studies have shown that changes in the masticatory muscle function associated with a switch to short-term SD led to changes in craniofacial morphology and alveolar bone architecture. However, the long-term effects of SD and the associated adaptive changes in craniofacial shape are unclear. Our novel study set out to profile prospective skull changes in mice fed with SDs over multiple generations using three-dimensional (3D) geometric morphometric analysis (GMA). Our results revealed that short-term SD consumption led to a significant decrease in craniofacial size, along with numerous shape changes. Long-term SD consumption over 15 continuous generations was not associated with changes in craniofacial size; however, shape analysis revealed mice with shortened crania and mandibles in the anteroposterior dimension, as well as relative widening in the transverse dimension compared to the average shape of all mice analyzed in our study. Moreover, changes in shape and size associated with different functional loads appeared to be independent - shape changes persisted after diets were switched for one generation, whereas size decreased after one generation and then returned to baseline size. Our study is the first to study the role of prolonged, multi-generational SD consumption in the determination of craniofacial size and shape.

Population specific data improves Fordisc®'s performance in Italians

Using discriminant functions based on morphometric data is one of the most approved methods for sex and ancestry estimation on skeletons. Nevertheless, population data from various parts of the world is needed in order to reliably classify an individual into a group. Due to population variation even sex estimation is biased when there is a lack of adequate data. Software that computes discriminant functions based on morphometric data is Fordisc^®. Unfortunately, the above mentioned effects reduce its applicability in countries other than the US. For improvement of this situation data collection is currently performed extensively. The present paper shows a comparison of an Italian sample from the identified modern skeletal collection of CAL (Collezione Antropologica Labanof) [1] (Cattaneo, 2018) at the Institute of Legal Medicine Milan, Italy with a Euro-American sample from the Forensic Data Bank at the University of Tennessee Knoxville, USA. Fordisc^® 3.1 was used to study population differences and sexual dimorphism. The analyzes were performed on a selection of 19 highly influential measurements that are present in most individuals of both groups respectively. Italian crania show in relation to Euro-Americans wider and lower vaults with shorter cranial bases and wider faces. The degree of sexual dimorphism is similar in both groups. Yet there is a shift in the absolute value range for males and females that biases sex estimation by almost 25% when an individual is classified on the respective other discriminant function. Our results provide explanations for Fordisc^®'s unsatisfying performance on non-US individuals. At the same time they show that significant improvement is achieved by adding more population samples to its dataset.

Cranial secular change from the nineteenth to the twentieth century in modern German individuals compared to modern Euro-American individuals

Studying secular changes on human skulls is a central issue in anthropological research, which is however insufficiently investigated for modern German populations. With our study, we focus on morphological cranial variations within Germans during the nineteenth and twentieth centuries. To study this, we recorded different facial landmarks from a cohort study of about 540 German individuals of different age and sex by calculating their cranial size, shape dimensions, and cranial module and cranial capacity to get information about variations occurring during the decades. According to this, measured variables for Germans and Americans, to which we compared our results, were maximum cranial length (glabello-occipital length), basion-bregma height (BBH), basion-nasion length (BNL), maximum cranial breadth (XCB), and cranial base breadth (AUB). Cranial size was calculated as the geometric mean of GOL, BBH, and XCB. Samples were organized into quarter century birth cohorts, with birth years ranging from 1800 to 1950. One-way ANOVA was used to test for variation among cohorts. Over the past 150 years, Americans and Germans showed significant parallel changes, but the American cranium remained relatively higher, with a longer cranial base, as well as narrower than the German cranium. Our results should also lead to the extension of the range of populations listed and investigated for Fordisc®, a forensic software to identify unknown individuals as from their skeletal remains or just parts of them. Fordisc cannot provide a satisfying identification of European individuals yet because the database is missing enough European reference samples.

Investigating the sex-related geometric variation of the human cranium

Accurate sexing methods are of great importance in forensic anthropology since sex assessment is among the principal tasks when examining human skeletal remains. The present study explores a novel approach in assessing the most accurate metric traits of the human cranium for sex estimation based on 80 ectocranial landmarks from 176 modern individuals of known age and sex from the Athens Collection. The purpose of the study is to identify those distance and angle measurements that can be most effectively used in sex assessment. Three-dimensional landmark coordinates were digitized with a Microscribe 3DX and analyzed in GNU Octave. An iterative linear discriminant analysis of all possible combinations of landmarks was performed for each unique set of the 3160 distances and 246,480 angles. Cross-validated correct classification as well as multivariate DFA on top performing variables reported 13 craniometric distances with over 85% classification accuracy, 7 angles over 78%, as well as certain multivariate combinations yielding over 95%. Linear regression of these variables with the centroid size was used to assess their relation to the size of the cranium. In contrast to the use of generalized procrustes analysis (GPA) and principal component analysis (PCA), which constitute the common analytical work flow for such data, our method, although computational intensive, produced easily applicable discriminant functions of high accuracy, while at the same time explored the maximum of cranial variability.