Reconstructing human population history from dental phenotypes

Peopling of the Americas: A new approach to assessing dental morphological variation in Asian and Native American populations

OBJECTIVES

Through biodistance analyses, anthropologists have used dental morphology to elucidate how people moved into and throughout the Americas. Here, we apply a method that focuses on individuals rather than sample frequencies through the application rASUDAS2, based on a naïve Bayes' algorithm.

MATERIALS AND METHODS

Using the database of C.G. Turner II, we calculated the probability that an individual could be assigned to one of seven biogeographic groups (American Arctic, North & South America, East Asia, Southeast Asia & Polynesia, Australo-Melanesia, Western Eurasia, & Sub-Saharan Africa) through rASUDAS2. The frequency of classifications for each biogeographic group was determined for 1418 individuals from six regions across Asia and the Americas.

RESULTS

Southeast Asians show mixed assignments but rarely to American Arctic or "American Indian." East Asians are assigned to East Asia half the time while 30% are assigned as Native American. People from the American Arctic and North & South America are assigned to Arctic America or non-Arctic America 75%-80% of the time, with 10%-15% classified as East Asian.

DISCUSSION

All Native American groups have a similar degree of morphological affinity to East Asia, as 10%-15% are classified as East Asian. East Asians are classified as Native American in 30% of cases. Individuals in the Western Hemisphere are decreasingly classified as Arctic the farther south they are located. Equivalent levels of classification as East Asian across all Native American groups suggests one divergence between East Asians and the population ancestral to all Native Americans. Non-arctic Native American groups are derived from the Arctic population, which represents the Native American founder group.

Human population dynamics in Upper Paleolithic Europe inferred from fossil dental phenotypes

Despite extensive archaeological research, our knowledge of the human population history of Upper Paleolithic Europe remains limited, primarily due to the scarce availability and poor molecular preservation of fossil remains. As teeth dominate the fossil record and preserve genetic signatures in their morphology, we compiled a large dataset of 450 dentitions dating between ~47 and 7 thousand years ago (ka), outnumbering existing skeletal and paleogenetic datasets. We tested a range of competing demographic scenarios using a coalescent-based machine learning Approximate Bayesian Computation (ABC) framework that we modified for use with phenotypic data. Mostly in agreement with but also challenging some of the hitherto available evidence, we identified a population turnover in western Europe at ~28 ka, isolates in western and eastern refugia between ~28 and 14.7 ka, and bottlenecks during the Last Glacial Maximum. Methodologically, this study marks the pioneering application of ABC to skeletal phenotypes, paving the way for exciting future research avenues.

Evaluating predictions of the patterning cascade model of crown morphogenesis in the human lower mixed and permanent dentition

OBJECTIVE

The patterning cascade model of crown morphogenesis has been studied extensively in a variety of organisms to elucidate the evolutionary history surrounding postcanine tooth form. The current research is the first to use a large modern human sample to examine whether the crown configuration of lower deciduous and permanent molars aligns with expectations derived from the model. This study has two main goals: 1) to determine if metameric and antimeric pairs significantly differ in size, accessory trait expression, and relative intercusp spacing, and 2) assess whether the relative distance among early-forming cusps accounts for observed variation in accessory cusp expression.

METHODS

Tooth size, intercusp distance, and morphological trait expression data were collected from 3D scans of mandibular dental casts representing participants of the Harvard Solomon Islands Project. Paired tests were utilized to compare tooth size, accessory trait expression, and relative intercusp distance between diphyodont metameres and permanent antimeres. Proportional odds logistic regression was implemented to investigate how the odds of greater accessory cusp expression vary as a function of the distance between early-developing cusps.

RESULTS/SIGNIFICANCE

Comparing paired molars, significant differences were identified for tooth size and cusp 5 expression. Several relative intercusp distances emerged as important predictors of cusp 6 expression, however, results for cusp 5 and cusp 7 did not match expected patterns. These findings support previous quantitative genetic results and suggest the development of neighboring crown structures represents a zero-sum partitioning of cellular territory and resources. As such, this study contributes to a better understanding of the foundations of deciduous and permanent molar crown variation in humans.

Inferring human neutral genetic variation from craniodental phenotypes

There is a growing consensus that global patterns of modern human cranial and dental variation are shaped largely by neutral evolutionary processes, suggesting that craniodental features can be used as reliable proxies for inferring population structure and history in bioarchaeological, forensic, and paleoanthropological contexts. However, there is disagreement on whether certain types of data preserve a neutral signature to a greater degree than others. Here, we address this unresolved question and systematically test the relative neutrality of four standard metric and nonmetric craniodental data types employing an extensive computational genotype-phenotype comparison across modern populations from around the world. Our computation draws on the largest existing data sets currently available, while accounting for geographically structured environmental variation, population sampling uncertainty, disparate numbers of phenotypic variables, and stochastic variation inherent to a neutral model of evolution. Our results reveal that the four data types differentially capture neutral genomic variation, with highest signals preserved in dental nonmetric and cranial metric data, followed by cranial nonmetric and dental metric data. Importantly, we demonstrate that combining all four data types together maximizes the neutral genetic signal compared with using them separately, even with a limited number of phenotypic variables. We hypothesize that this reflects a lower level of genetic integration through pleiotropy between, compared to within, the four data types, effectively forming four different modules associated with relatively independent sets of loci. Therefore, we recommend that future craniodental investigations adopt holistic combined data approaches, allowing for more robust inferences about underlying neutral genetic variation.

Dental size variation in admixed Latin Americans: Effects of age, sex and genomic ancestry

Dental size variation in modern humans has been assessed from regional to worldwide scales, especially under microevolutionary and forensic contexts. Despite this, populations of mixed continental ancestry such as contemporary Latin Americans remain unexplored. In the present study we investigated a large Latin American sample from Colombia (N = 804) and obtained buccolingual and mesiodistal diameters and three indices for maxillary and mandibular teeth (except third molars). We evaluated the correlation between 28 dental measurements (and three indices) with age, sex and genomic ancestry (estimated using genome-wide SNP data). In addition, we explored correlation patterns between dental measurements and the biological affinities, based on these measurements, between two Latin American samples (Colombians and Mexicans) and three putative parental populations: Central and South Native Americans, western Europeans and western Africans through PCA and DFA. Our results indicate that Latin Americans have high dental size diversity, overlapping the variation exhibited by the parental populations. Several dental dimensions and indices have significant correlations with sex and age. Western Europeans presented closer biological affinities with Colombians, and the European genomic ancestry exhibited the highest correlations with tooth size. Correlations between tooth measurements reveal distinct dental modules, as well as a higher integration of postcanine dentition. The effects on dental size of age, sex and genomic ancestry is of relevance for forensic, biohistorical and microevolutionary studies in Latin Americans.

An initial investigation of dental morphology variation among three southern Naga ethnic groups of Northeast India

OBJECTIVES

This study examines dental morphology trait prevalence among three southern Naga groups and compares them to 10 ethnic groups from other regions of South Asia to accomplish two objectives: assess the biological relationship of these Tibeto-Burman-speakers to speakers of non-Tibeto-Burman languages in other South Asian regions, and determine which traits distinguish northeast Indians from other South Asians.

METHODS

Dental morphology traits were scored with the Arizona State University Dental Anthropology System. Tooth-trait combinations were evaluated for significant inter-trait correlation and intra-trait correspondence within dental fields. Comparisons were based on simple trait prevalence and with Smith's MMD. Affinities based on the former were accomplished with correspondence analysis and principal components analysis. Affinities based on the latter were undertaken with neighbor-joining cluster analysis and multidimensional scaling.

RESULTS

After elimination due to inter-trait correlations and uniform prevalence, biodistances based on the remaining 17 tooth-trait combinations identify significant differences between northeast Indians and other South Asian ethnic groups due to high frequencies of shoveling on the maxillary incisors and Cusp 6 on the mandibular molars coupled with low frequencies of Carabelli's trait and Cusp 5 on UM1 and UM2, respectively.

CONCLUSIONS

Patterns of biodistances obtained from dental morphology are consilient with those obtained from DNA indicating statistically significant differences between northeast Indians from members of ethnic groups of other regions of South Asia. Researchers should explore the sex-specific patterns. Biodistances should not be limited to "key" teeth within dental fields, for in almost every case traits present on mesial and distal teeth yield non-redundant information.

Microevolutionary variation in molar morphology of Onychomys leucogaster decoupled from genetic structure

In neutral models of quantitative trait evolution, both genetic and phenotypic divergence scale as random walks, producing a correlation between the two measures. However, complexity in the genotype-phenotype map may alter the correlation between genotypic and phenotypic divergence, even when both are evolving neutrally or nearly so. Understanding this correlation between phenotypic and genetic variation is critical for accurately interpreting the fossil record. This study compares the geographic structure and scaling of morphological variation of the shape of the first lower molar of 77 individuals of the northern grasshopper mouse Onychomys leucogaster to genome-wide SNP variation in the same sample. We found strong genetic structure but weak or absent morphological structure indicating that the scaling of each type of variation is decoupled from one another. Low P_ST values relative to F_ST values are consistent with a lack of morphological divergence in contrast to genetic divergence between groups. This lack of phenotypic structure and the presence of notable within-sample phenotypic variance are consistent with uniform selection or constraints on molar shape across a wide geographic and environmental range. Over time, this kind of decoupling may result in patterns of phenotypic stasis masking underlying genetic patterns.

A novel system for classifying tooth root phenotypes

Human root and canal number and morphology are highly variable, and internal root canal form and count does not necessarily co-vary directly with external morphology. While several typologies and classifications have been developed to address individual components of teeth, there is a need for a comprehensive system, that captures internal and external root features across all teeth. Using CT scans, the external and internal root morphologies of a global sample of humans are analysed (n = 945). From this analysis a method of classification that captures external and internal root morphology in a way that is intuitive, reproducible, and defines the human phenotypic set is developed. Results provide a robust definition of modern human tooth root phenotypic diversity. The method is modular in nature, allowing for incorporation of past and future classification systems. Additionally, it provides a basis for analysing hominin root morphology in evolutionary, ecological, genetic, and developmental contexts.

Maxillary molar root and canal morphology of Neolithic and modern Chinese

OBJECTIVE

This study aimed to characterize Neolithic human maxillary molars from archeological remains at the Jiaojia site, Shandong, China, and compare their ultrastructural features with sex and age-matched modern locals.

DESIGN

Maxillary first (n = 86) and second (n = 80) molars in 5000-year-old individuals (n = 50) from the Jiaojia site were scanned by cone-beam computed tomography (CBCT). Sex and age-matched control groups were assigned from oral surgical patients at Shandong University. Images were analyzed for crown size, root length, root morphology, canal inter-orifice distances, mesiobuccal canal morphology, and second mesiobuccal (MB2) canal prevalence and location. Neolithic and modern values were compared statistically using Chi-squared and Mann-Whitney test at p < .05.

RESULTS

Crown and root size were smaller, and canal inter-orifice distances were shorter in Neolithic maxillary molars than their modern counterparts. For mesiobuccal roots, Weine's Type I single canals were the most prevalent in Neolithic and modern first and second molars. MB2 canal prevalence were not significantly different (p > .05) in Neolithic (53.3%) or modern (60.5%) first molars, and Neolithic (11.3%) or modern (21.3%) second molars. But, MB2 prevalence was significantly higher for modern than ancient male first (p = .032) and second (p = .005) molars. Additionally, MB2 were located more mesially and closer to MB1 in Neolithic than modern molars.

CONCLUSIONS

Maxillary molar root and canal morphology of ancient 5000-year-old remains at the Jiaojia site resemble that of local patients. A trend towards larger tooth size, and more dispersed MB2 canals over this short evolutionary period warrants additional investigation.

New Approach to Dental Morphometric Research Based on 3D Imaging Techniques

Recent progress in imaging and image processing techniques has provided for improvements in odontological research in a variety of aspects. Thus, the presented method has been developed precisely in order to assess metrically 3D reconstructions of teeth. Rapidly and accurately obtained data of a wide range and appropriate density are sufficient enough for morphometric studies rather than tooth size assessments which are inherent to conventional techniques. The main contributions providing for holistic and objective morphometric analysis of teeth are the following: (1) interpretation of basic dental morphological features; (2) automated of orientational coordinate system setup based on tooth surface analysis; (3) new tooth morphometric parameters which could not be obtained through application of conventional odontometric techniques; (4) methodological novelty for automated odontomorphometric analysis pipeline. Application of tomographic imaging, which has been used for obtaining 3D models, expands the proposed method potential further through providing detailed and comprehensive reconstructions of teeth. The current study was conducted on unique material from the archaeological site of Sunghir related to the Upper Palaeolithic period. Metric assessments of external and internal morphological layers of teeth were performed in common orientation and sectioning. The proposed technique allowed more profound analysis of Sunghirian teeth which date back to the times of modern human morphology formation.

Population history of Brazilian south and southeast shellmound builders inferred through dental morphology

OBJECTIVE

The Southeast and South Coast of Brazil was inhabited during most of the Holocene by shellmound builders. Although there are cultural differences in the archaeological record between regions, it is still debatable how these differences may relate to different population histories. Here, we contribute to this discussion by exploring dental morphological affinities between several regional series.

MATERIALS AND METHODS

Dental morphology of 385 individuals from 14 archaeological sites was analyzed using the Arizona State University Dental Anthropology System. Fifteen traits were used to explore morphological affinities among series through Euclidean distance, Mean Measure of Divergence, and Principal Component Analysis. Mantel matrix correlation and partial correlation tests were used to examine the association between biological, geographic, and chronological distances.

RESULTS

Morphological affinities show that ceramic and nonceramic South Coast groups cluster and differ from most Southeast series. In contrast, Southeast coastal and riverine groups display high morphological variance, showing less biological coherence among them. These biological distances between regions are partially explained by geography, but not by chronology.

CONCLUSIONS

The results support that these coastal populations were low-mobility groups. Although interactions between individuals of different regions likely existed, gene flow occurred mostly among individuals from local or adjacent areas. The introduction of ceramic in the South Coast is not associated with changes in dental morphology patterns, suggesting its adoption is not exclusively associated with the arrival of different biological groups. Southeast coastal and riverine groups show high phenotypic diversity, suggesting a different history of human occupation and cultural development than observed in the South Coast.

The genetic architecture of anterior tooth morphology in a longitudinal sample of Australian twins and families

OBJECTIVE

This study presents a quantitative genetic analysis of human anterior dental morphology in a longitudinal sample of known genealogy. The primary aim of this work is to generate a suite of genetic correlations within and between deciduous and permanent characters to access patterns of integration across the diphyodont dental complex.

DESIGN

Data were recorded from casted tooth crowns representing participants of a long-term Australian twin and family study (deciduous n = 290, permanent n = 339). Morphological trait expression was observed and scored following Arizona State University Dental Anthropology System standards. Bivariate genetic correlations were estimated using maximum likelihood variance decomposition models in SOLAR v.8.1.1.

RESULTS

Genetic correlation estimates indicate high levels of integration between antimeres but low to moderate levels among traits within a tooth row. Only 9% of deciduous model comparisons were significant, while pleiotropy was indicated for one third of permanent trait pairs. Canine characters stood out as strongly integrated, especially in the deciduous dentition. For homologous characters across dentitions (e.g., deciduous i¹ shoveling and permanent I¹ shoveling), ∼70% of model comparisons yielded significant genetic correlations.

CONCLUSIONS

Patterns of genetic correlation suggest a morphological canine module that spans the primary and secondary dentition. Results also point to the existence of a genetic mechanism conserving morphology across the diphyodont dental complex, such that paired deciduous and permanent traits are more strongly integrated than characters within individual tooth rows/teeth.

Uncloaking a Lost Cause: Decolonizing ancestry estimation in the United States

Since the professionalization of US‐based forensic anthropology in the 1970s, ancestry estimation has been included as a standard part of the biological profile, because practitioners have assumed it necessary to achieve identifications in medicolegal contexts. Simultaneously, forensic anthropologists have not fully considered the racist context of the criminal justice system in the United States related to the treatment of Black, Indigenous, and People of Color; nor have we considered that ancestry estimation might actually hinder identification efforts because of entrenched racial biases. Despite ongoing criticisms from mainstream biological anthropology that ancestry estimation perpetuates race science, forensic anthropologists have continued the practice. Recent years have seen the prolific development of retooled typological approaches with 21st century statistical prowess to include methods for estimating ancestry from cranial morphoscopic traits, despite no evidence that these traits reflect microevolutionary processes or are suitable genetic proxies for population structure; and such approaches have failed to critically evaluate the societal consequences for perpetuating the biological race concept. Around the country, these methods are enculturated in every aspect of the discipline ranging from university classrooms, to the board‐certification examination marking the culmination of training, to standard operating procedures adopted by forensic anthropology laboratories. Here, we use critical race theory to interrogate the approaches utilized to estimate ancestry to include a critique of the continued use of morphoscopic traits, and we assert that the practice of ancestry estimation contributes to white supremacy. Based on the lack of scientific support that these traits reflect evolutionary history, and the inability to disentangle skeletal‐based ancestry estimates from supporting the biological validity of race, we urge all forensic anthropologists to abolish the practice of ancestry estimation.

Neutral evolution of human enamel-dentine junction morphology

Teeth have been studied for decades and continue to reveal information relevant to human evolution. Studies have shown that many traits of the outer enamel surface evolve neutrally and can be used to infer human population structure. However, many of these traits are unavailable in archaeological and fossil individuals due to processes of wear and taphonomy. Enamel-dentine junction (EDJ) morphology, the shape of the junction between the enamel and the dentine within a tooth, captures important information about tooth development and vertebrate evolution and is informative because it is subject to less wear and thus preserves more anatomy in worn or damaged specimens, particularly in mammals with relatively thick enamel like hominids. This study looks at the molar EDJ across a large sample of human populations. We assessed EDJ morphological variation in a sample of late Holocene modern humans (n = 161) from archaeological populations using μ-CT biomedical imaging and geometric morphometric analyses. Global variation in human EDJ morphology was compared to the statistical expectations of neutral evolution and "Out of Africa" dispersal modeling of trait evolution. Significant correlations between phenetic variation and neutral genetic variation indicate that EDJ morphology has evolved neutrally in humans. While EDJ morphology reflects population history, its global distribution does not follow expectations of the Out of Africa dispersal model. This study increases our knowledge of human dental variation and contributes to our understanding of dental development more broadly, with important applications to the investigation of population history and human genetic structure.

Testing the utility of dental morphological trait combinations for inferring human neutral genetic variation

Researchers commonly rely on human dental morphological features in order to reconstruct genetic affinities among past individuals and populations, particularly since teeth are often the best preserved part of a human skeleton. Tooth form is considered to be highly heritable and selectively neutral and, therefore, to be an excellent proxy for DNA when none is available. However, until today, it remains poorly understood whether certain dental traits or trait combinations preserve neutral genomic signatures to a greater degree than others. Here, we address this long-standing research gap by systematically testing the utility of 27 common dental traits and >134 million possible trait combinations in reflecting neutral genomic variation in a worldwide sample of modern human populations. Our analyses reveal that not all traits are equally well-suited for reconstructing population affinities. Whereas some traits largely reflect neutral variation and therefore evolved primarily as a result of genetic drift, others can be linked to nonstochastic processes such as natural selection or hominin admixture. We also demonstrate that reconstructions of population affinity based on many traits are not necessarily more reliable than those based on only a few traits. Importantly, we find a set of highly diagnostic trait combinations that preserve neutral genetic signals best (up to [Formula: see text] _r = 0.580; 95% r range = 0.293 to 0.758; P = 0.001). We propose that these trait combinations should be prioritized in future research, as they allow for more accurate inferences about past human population dynamics when using dental morphology as a proxy for DNA.

Root and root canal diversity in human permanent maxillary first premolars and upper/lower first molars from a 14th-17th and 18th-19th century Radom population

OBJECTIVE

The aim of this study was to assess whether analyzed groups from two historical periods: Late Medieval (LMP), and Modern (MP) from Radom varied in the number of tooth roots and root canal system morphology.

METHODS

Root morphology of 229 permanent human teeth were analyzed using Cone Beam Computed Tomography. Additionally, the mitochondrial DNA (mtDNA) of 29 individuals from the LMP and 31 from the MP was analyzed.

RESULTS

In LMP, the maxillary first premolars were dominated by one root, while in MP second and third roots also appeared. Maxillary first molars in LMP presented three roots, while two-rooted forms occurred in MP. All mandibular first molars from the LMP and almost all (98%) from MP presented two roots. The greatest diversity in terms of root canal number occurred in one-rooted maxillary first premolars, the mesiobuccal root of the maxillary first molars, and the mesial and distal roots of the mandibular first molars in both groups. A few haplogroups from outside Europe (C, N, and R) were recorded in the MP Radom population. Moreover, this population had substantially higher haplogroup diversity compared with the LMP population.

CONCLUSION

Odontological research indicates an increase in the diversity in the number of roots and the shape of root canals in MP. This information corresponds to genetic research, which also indicates an increase in the diversity of haplogroups during the MP.

Variation of 3D outer and inner crown morphology in modern human mandibular premolars

Objectives

This study explores the outer and inner crown of lower third and fourth premolars (P₃, P₄) by analyzing the morphological variation among diverse modern human groups.

Materials and Methods

We studied three‐dimensional models of the outer enamel surface and the enamel–dentine junction (EDJ) from μCT datasets of 77 recent humans using both an assessment of seven nonmetric traits and a standard geometric morphometric (GM) analysis. For the latter, the dental crown was represented by four landmarks (dentine horns and fossae), 20 semilandmarks along the EDJ marginal ridge, and pseudolandmarks along the crown and cervical outlines.

Results

Certain discrete traits showed significantly different regional frequencies and sexual dimorphism. The GM analyses of both P₃s and P₄s showed extensive overlap in shape variation of the various populations (classification accuracy 15–69%). The first principal components explained about 40% of shape variance with a correlation between 0.59 and 0.87 of the features of P₃s and P₄s. Shape covariation between P₃s and P₄s expressed concordance of high and narrow or low and broad crowns.

Conclusions

Due to marked intragroup and intergroup variation in GM analyses of lower premolars, discrete traits such as the number of lingual cusps and mesiolingual groove expression provide better geographic separation of modern human populations. The greater variability of the lingual region suggests a dominance of functional constraints over geographic provenience or sex. Additional information about functionally relevant aspects of the crown surface and odontogenetic data are needed to unravel the factors underlying dental morphology in modern humans.

Dental evolutionary rates and its implications for the Neanderthal-modern human divergence

The origin of Neanderthal and modern human lineages is a matter of intense debate. DNA analyses have generally indicated that both lineages diverged during the middle period of the Middle Pleistocene, an inferred time that has strongly influenced interpretations of the hominin fossil record. This divergence time, however, is not compatible with the anatomical and genetic Neanderthal affinities observed in Middle Pleistocene hominins from Sima de los Huesos (Spain), which are dated to 430 thousand years (ka) ago. Drawing on quantitative analyses of dental evolutionary rates and Bayesian analyses of hominin phylogenetic relationships, I show that any divergence time between Neanderthals and modern humans younger than 800 ka ago would have entailed unexpectedly rapid dental evolution in early Neanderthals from Sima de los Huesos. These results support a pre-800 ka last common ancestor for Neanderthals and modern humans unless hitherto unexplained mechanisms sped up dental evolution in early Neanderthals.

Human bony labyrinth is an indicator of population history and dispersal from Africa

The cavity system of the inner ear—the so-called bony labyrinth—houses the senses of balance and hearing. This structure is embedded in dense petrous bone, fully formed by birth and generally well preserved in human skeletal remains, thus providing a rich source of morphological information about past populations. Here we show that labyrinthine morphology tracks genetic distances and geography in an isolation-by-distance model with dispersal from Africa. Because petrous bones have become prime targets of ancient DNA recovery, we propose that all destructive studies first acquire high-resolution 3D computed-tomography data prior to any invasive sampling. Such data will constitute an important archive of morphological variation in past and present populations, and will permit individual-based genotype–phenotype comparisons.

The dispersal of modern humans from Africa is now well documented with genetic data that track population history, as well as gene flow between populations. Phenetic skeletal data, such as cranial and pelvic morphologies, also exhibit a dispersal-from-Africa signal, which, however, tends to be blurred by the effects of local adaptation and in vivo phenotypic plasticity, and that is often deteriorated by postmortem damage to skeletal remains. These complexities raise the question of which skeletal structures most effectively track neutral population history. The cavity system of the inner ear (the so-called bony labyrinth) is a good candidate structure for such analyses. It is already fully formed by birth, which minimizes postnatal phenotypic plasticity, and it is generally well preserved in archaeological samples. Here we use morphometric data of the bony labyrinth to show that it is a surprisingly good marker of the global dispersal of modern humans from Africa. Labyrinthine morphology tracks genetic distances and geography in accordance with an isolation-by-distance model with dispersal from Africa. Our data further indicate that the neutral-like pattern of variation is compatible with stabilizing selection on labyrinth morphology. Given the increasingly important role of the petrous bone for ancient DNA recovery from archaeological specimens, we encourage researchers to acquire 3D morphological data of the inner ear structures before any invasive sampling. Such data will constitute an important archive of phenotypic variation in present and past populations, and will permit individual-based genotype–phenotype comparisons.