Skip to main content
Advertisement
Browse Subject Areas
?

Click through the PLOS taxonomy to find articles in your field.

For more information about PLOS Subject Areas, click here.

  • Loading metrics

Shape variation in modern human upper premolars

  • Petra G. Šimková ,

    Roles Formal analysis, Investigation, Project administration, Visualization, Writing – original draft, Writing – review & editing

    petra.simkova@univie.ac.at

    Affiliations Department of Evolutionary Anthropology, University of Vienna, Vienna, Austria, Human Evolution and Archaeological Sciences HEAS, University of Vienna, Vienna, Austria

  • Lisa Wurm,

    Roles Investigation, Writing – review & editing

    Affiliation Medical Technology Cluster, Business Upper Austria–OÖ Wirtschaftsagentur GmbH, Linz, Austria

  • Cinzia Fornai,

    Roles Methodology, Writing – review & editing

    Affiliations Department of Evolutionary Anthropology, University of Vienna, Vienna, Austria, Human Evolution and Archaeological Sciences HEAS, University of Vienna, Vienna, Austria, Department of Research in Occlusion Medicine, Vienna School of Interdisciplinary Dentistry–VieSID, Klosterneuburg, Austria, Center for Clinical Research, University Clinic of Dentistry Vienna, Medical University of Vienna, Vienna, Austria, Institute of Evolutionary Medicine, University of Zurich, Zürich, Switzerland

  • Viktoria A. Krenn,

    Roles Methodology, Writing – review & editing

    Affiliations Department of Evolutionary Anthropology, University of Vienna, Vienna, Austria, Institute of Evolutionary Medicine, University of Zurich, Zürich, Switzerland, Fraunhofer Austria Research GmbH, Graz, Austria

  • Gerhard W. Weber

    Roles Conceptualization, Methodology, Project administration, Supervision, Writing – review & editing

    Affiliations Department of Evolutionary Anthropology, University of Vienna, Vienna, Austria, Human Evolution and Archaeological Sciences HEAS, University of Vienna, Vienna, Austria, Core Facility for Micro-Computed Tomography, University of Vienna, Vienna, Austria

Abstract

Morphological variation in modern human dentition is still an open field of study. The understanding of dental shape and metrics is relevant for the advancement of human biology and evolution and is thus of interest in the fields of dental anthropology, as well as human anatomy and medicine. Of concern is also the variation of the inner aspects of the crown which can be investigated using the tools and methods of virtual anthropology. In this study, we explored inter- and intra-population morphometric variation of modern humans’ upper third and fourth premolars (P3s and P4s, respectively) considering both the inner and outer aspects of the crown, and discrete traits. We worked by means of geometric morphometrics on 3D image data from a geographically balanced sample of human populations from five continents, to analyse the shape of the dentinal crown, and the crown outline in 78 P3s and 76 P4s from 85 individuals. For the study of dental traits, we referred to the Arizona State University Dental Anthropology System integrated with more recent classification systems. The 3D shape variation of upper premolar crowns varied between short and mesio-distally broad, and tall and mesio-distally narrow. The observed shape variation was independent from the geographical origin of the populations, and resulted in extensive overlap. We noted a high pairwise correlation (r1 = 0.83) between upper P3s and P4s. We did not find any significant geographic differences in the analysed non-metric traits. Our outcomes thus suggest that geographical provenance does not play a determinant role in the shaping of the dental crown, whose genesis is under strict genetic control.

Introduction

Dental morphology has been investigated widely in order to study the variation within and between species based on general shape, metric and non-metric traits [114]. The crown of the premolars is usually composed of two main cusps (buccal and lingual) each continuing into a root cone [15]. These cones can form two separate roots or be fused into one root. Rarely, a three-rooted premolar with one lingual and two buccal roots can be observed [11,15]. By position within the dental row and function, maxillary premolars are intermediate between front and molar teeth. In this respect, third premolars possess some functional characteristics of a canine, while fourth premolars are more engaged in grinding functions comparable to molars [6].

3D morphological variation of upper premolar crowns among modern human populations is poorly investigated since most of the studies focused on the roots and root canals [16,17] also in clinical settings [1820]. Previous investigations of the human maxillary premolar variation using traditional metric approaches [21] and landmark-based multivariate techniques [22], but on the outer enamel surface (OES), revealed morphological differences that help to distinguish between modern human populations. On the other hand, upper premolar morphology delivered varying results for the purposes of taxonomic classification in paleoanthropological studies using non-metric trait analyses [2,23,24], and geometric morphometric approach delivered significant results only for the P4s [23]. Metric measurements such as bucco-lingual and mesio-distal crown dimensions or crown height have also been used to study premolar size variation between and within populations [3,25,26]. Other studies focused on the occurrence and expressions of non-metric traits such as the maxillary premolar accessory ridges [11,27,28], the premolar mesial and distal accessory cusp [11,29,30], or the essential crests [30] among various modern human populations.

Compared to the OES, the examination of the enamel-dentine junction (EDJ) presents advantages because it represents the primary developmental structure of the crown and is consequently the precursor of the enamel cap’s morphology [3140], naturally resulting in a strong correlation between these two surfaces [31,33,41,42]. Furthermore, the overlying enamel protects the EDJ from early deterioration processes such as abrasion, erosion or chipping. The EDJ thus represents an appropriate alternative to the OES since it also bears relevant taxonomic information [31,37,39,43,44]. Morphometric studies of the EDJ of both permanent and deciduous molars are well represented in the literature [31,33,4446], while for maxillary premolars, extensive comparative studies are very rare [47], or only focused on the descriptive assessment of the EDJ and OES [4851]. Currently, 3D geometric morphometric (GM) analyses of the inner tooth crown variation in geographically diverse modern human samples have been published for permanent and deciduous molars [23,33,44,46,47,52], as well as lower premolars [53,54], and central incisors [55]. These studies did not detect geographically-dependent dental shape variation [33,54], thus, in the current study, we applied quantitative approaches to explore inter- and intra-population shape variation and covariation of modern human P3s and P4s, expecting to find similar outcomes. We want to test the hypothesis that the ranges of shape variation of upper premolars in different world populations overlap, meaning that recent human populations cannot be distinguished based on the shape of their upper premolar’s dentinal crown. We used high-resolution 3D image data and GMM to investigate a geographically heterogeneous sample consisting of modern human populations from Africa, South America, Europe, the Near East, Southeast Asia and Oceania. Size variation was investigated based on the natural logarithm of Centroid Size. Moreover, we analysed the expression of the non-metric trait on both the EDJ and the OES. Comprehensive knowledge of these factors is essential in biology for taxonomic and evolutionary research, including human genotype-phenotype associations, and is also relevant in the various fields of medicine including forensics and dentistry.

Materials

We included 78 maxillary P3s and 76 maxillary P4s from 85 recent modern humans (Table 1). In particular, our sample comprised individuals from Oceania (n = 15), including the indigenous people of Papua New Guinea (n = 6), the former diverging from Eurasian populations around 62 to 75 thousand years ago [5658], South America (n = 19), including the native population of Selk’nam from Tierra del Fuego (n = 7), Europe (n = 17), including Avars from the 8th century (n = 7), sub-Saharan Africa (n = 15), Southeast Asia (n = 6) and the Near East (n = 14) including Bedouins (n = 4) and Natufians (n = 10). All teeth used for the EDJ analyses were free of decay or fillings and did not exceed Molnar’s [59] wear stage 3 (slight to moderate wear). Accordingly, populations from various geographical regions, climatic zones, and environments, featuring different life styles, were represented in our sample.

thumbnail
Table 1. List of third and fourth premolars (P3 and P4, respectively) used in this study.

(Wear recorded according to Molnar, 1977; f = female; m = male; + = present;— = absent; N.A. = not available; EDJ = enamel-dentine junction).

https://doi.org/10.1371/journal.pone.0301482.t001

Methods

Scanning & segmentation

All specimens were scanned at the Vienna μCT Lab, Austria, with an industrial Viscom X8060 NDT scanner (scanning parameters: 110–140 kV, 280–410 mA, 1400–2000 ms, 0.75 mm copper filter with a voxel size ranging between 20 and 50 μm). Virtual segmentation of the μCT data was performed in Amira software (www.thermofisher.com) to obtain triangulated surface models of the enamel and dentine. Models were extracted from the left premolars but, if absent or unusable, the right-side premolars were considered (n = 49) since directional asymmetry is not expected in teeth [60]. If this was the case, the right-side datasets were flipped before surface extraction. Slightly worn dentinal horn tips were virtually reconstructed in Amira software by extending the contours of the existing part of the dentinal horn tips (n = 38).

Reorientation of the models and data collection on the cervical and crown outlines

The segmented dental surface models were reoriented in Geomagic Design X 64 (www.3dsystems.com) according to an established protocol [61]. Thus, the best-fit plane of the cervical margin, or cervical plane, was computed and was used to reorient the dataset to set the cervical plane parallel to the x-y plane of the virtual working space. Afterwards, the crown was rotated to align the mesial ridge to the y-axis. Cervical and crown outlines of the reoriented surfaces were collected and projected onto the cervical plane. The outlines were imported into Rhinoceros 6 (www.rhino3d.com) and split into 24 curve segments by equiangular radial vectors that originated from the centroid of the outline. The intersection points between the radii and the outline marked 24 pseudo-landmarks.

Landmark collection on the EDJ

The reoriented surfaces of the dentinal crowns were imported into the free-to-download EVAN-Toolbox software (www.evan-society.org) for further landmark collection. A total of 24 landmarks were placed on the EDJ of each of the premolars. On all premolars, two fixed landmarks were placed on the paracone and protocone horn tips. On the P3s, the remaining two fixed landmarks were placed at the deepest mesial and distal points of the central groove. Since the central groove is not as well defined in P4s, the landmarks were instead placed on the deepest points of the distal and mesial ridges, which are always well visible. Additionally, the marginal ridge of both tooth types was represented by 20 curve semilandmarks. The sliding of the semilandmarks was performed using the bending energy technique [6264]. Our approach matched that performed on lower premolars in Krenn et al. [54], who validated this technique by assessing the observer error. We do not expect the validity of these protocols to change in the upper premolars.

Geometric morphometric analyses

We performed GM analyses using the EVAN-Toolbox for each set of landmarks (namely, cervical outline, crown outline, EDJ, and dentinal crown) for P3s and P4s, separately. The dentinal crown dataset was created by merging the cervical outline and the EDJ landmark configurations. A General Procrustes Analysis (GPA [65]) was performed for all data sets to normalize the landmark configurations. The resulting Procrustes shape coordinates were analyzed via principal component analysis (PCA). The Thin-Plate Spline technique was used to visualize the shape changes along the principal components (PCs) by means of warping [62,66]. To explore the shape covariation between P3s and P4s, as well as between different aspects within the same dental type (i.e. between the cervical outline and the EDJ), we performed a 2-block partial least squares (2B-PLS) analysis. To analyze the effect of allometry, multivariate regression was performed. The size was investigated by using the natural logarithm of Centroid Size (lnCS [67]). By using the lnCS derived from the dentinal crown dataset (combining the EDJ with cervical outline) we obtained a measure of size representing crown height as well. For testing tooth size differences between populations, we performed Kruskal-Wallis and Mann-Whitney-U tests (in PAST 4.03; www.softpedia.com), since the obtained data was not normally distributed. Permutational Multivariate Analysis of Variance (PERMANOVA) was performed using PAST 4.03 to test for shape differences between the groups (using all PC-scores), allowing for a distribution-free setting using permutational algorithms.

Non-metric traits

The following non-metric traits were considered:

  1. Premolar mesial and distal accessory cusps
    Accessory mesial or distal cusps originate from a bifurcation on the margin of the sagittal sulcus that is separating the buccal and lingual cusps, forming a free-standing cusp on the OES [11,13,68]. Sakai et al. [69] describe accessory cusps as indistinct on the EDJ, with a protrusion or swelling of the marginal ridge where the cusp could be. We used the following classifications to evaluate the presence of the accessory cusps: absent (0), mesial cusp (1), distal cusp (2), and mesial and distal cusp (3).
  2. Maxillary premolar accessory ridges (MxPAR)
    Accessory ridges may or may not be present on the lingual lobe segment of the buccal cusp, mesially and distally from the essential crest of both P3s and P4s [11], and can be extended into the central groove. These features are prevalent in both the OES and the EDJ, and their distribution was found to be heterogeneous among populations [28,70]. To represent the expression of these features, we used the following classification: absent (0), mesial accessory ridge (1), distal accessory ridge (2), and mesial and distal accessory ridges (3).
  3. Essential crest
    The essential crests of the central lobe of the paracone and protocone are usually separated by the sagittal sulcus and central groove. However, if this is not the case, i.e., both essential crests are connected, they create a transverse crest. Additionally, the separated essential crest of the buccal cusp can also be bifurcated. Several studies were able to detect these features on both EDJ and OES [4851]. In this study, we classified every kind of essential crest expression separately, namely the separated buccal and lingual crests (0), transverse crest (1), and the bifurcated buccal essential crest (2).

Traits 1 and 2 were derived from the Arizona State University Dental Anthropological System [11,12]. For trait 3 (essential crest) we referred to Bailey’s [2] buccal essential crest and traverse crest grading system. All traits were scored both on the OES and EDJ of each specimen and evaluated according to their presence and expression (Fig 1), performing a Chí2 test in PAST 4.03.

thumbnail
Fig 1. Illustration of scored non-metric traits and their expression categories.

https://doi.org/10.1371/journal.pone.0301482.g001

Ethics statement

In this study, we used non-invasive μCT (micro computed tomography) data (a three-dimensional X-ray method) from human remains included in the collections of the Department of Evolutionary Anthropology, University of Vienna and the Natural History Museum Vienna, Austria. Since only osteological material was used, no ethical approval or guidance was required as the research did not involve present-day human samples. The specimens were exported to Europe between the 19th century and the beginning of the 20th century, being curated in Austrian institutions for over a century. They were regarded and treated as precious parts of natural history collections and used for research since then. Today, their acquisition raises concerns according to current ethical principles. During the last 15 years, skulls were μCT-scanned at the Vienna μCT Lab of the University of Vienna for the purpose of preservation and fundamental research. Specimens were handled only by highly trained personnel in the most respectful way. No damages or alterations of any kind were caused to the remains.

Results

Shape variation in P3s

Results obtained from the PCA analysis of the P3 dentinal crown (Fig 2A) revealed an extensive overlap of the populations. Variation along PC1 (27% of the total variance explained) was mainly driven by the expansion of the distal fossa relative to the mesial fossa and the relative mesio-distal broadening of the cervical outline. The shape of the P3 dentinal crown varied from relatively short with a broad base, and rounded lingual occlusal aspect, to tall and mesiodistally narrow crowns. The horn tips shifted towards each other or apart, relative to the expansion and reduction of the cervical aspect of the crown, with broad crowns expressing horn tips closer to each other and narrow crowns the other way around. Shape variation along PC2 (14%) reflected the relative expansion of the mesio-buccal and mesio-lingual areas of the dentinal crowns with respect to each other. With the expansion of the mesio-lingual side, the relative position of the lingual cusp shifted mesially. Along PC3 (10%), the relative shape changes reflected the expansion of the distal aspect of the tooth as well as the relative position of the lingual cusp. The lingual cusp was more mesially positioned in disto-buccally expanded dentinal crowns, and disto-lingually positioned in disto-lingually expanded crowns.

thumbnail
Fig 2.

PC1—PC2 plot for third premolars (a) and fourth premolars (b) in shape space for the dentinal crown. The warpings show the real variation in occlusal and distal view at values PC1±0.15/ PC2 ±0.10; c) Landmark configuration capturing the dentinal crown of the P3s; and d) the P4s (orange = landmarks, blue = semi- and pseudolandmarks).

https://doi.org/10.1371/journal.pone.0301482.g002

Shape variation in P4s

The shape of the P4 dentinal crowns (Fig 2B) along PC1 (29%) varied between short with mesio-distally expanded lingual cusp and bucco-lingually elongated crown base to tall with a mesio-distally narrow lingual cusp and bucco-lingually shorter crown base. A narrower occlusal aspect corresponds with a taller buccal cusp, positioned further from the lingual cusp relative to the crown base, while in broader occlusal aspects the cusps are more even and shifted towards each other. We observed a separation of Southeast Asian individuals from South Americans and Oceanians. However, this was probably caused by the small size of the Southeast Asian sample. On the other hand, Near Easterners and South Americans, although partially overlapping, tended to show the opposite shape variation along PC1. Shape variation along PC2 (18%) was mainly driven by the relative bucco-lingual position of the mesial fossa. A taller buccal cusp and relatively buccolingually elongated cervical outline were associated with a lingually shifted mesial fossa and vice versa. The variation along PC3 (14%) reflected the relative buccolingual position of the distal fossa as well as the height of the crown and relative expansion of the cervical outline.

The PCA analyses of the cervical outline, crown outline and EDJ dataset in both P3s and P4s (Supplementary Fig A, B, C, D in S1 File) showed that the populations in our sample overlapped widely and could not be differentiated based on the morphology of the upper premolars crown. Supplementary Table A in S1 File contains percentages of variance explained by the first five PCs in every dataset. Furthermore, our observations were confirmed by insignificant results obtained from the PERMANOVA analyses for P3s’ (p = 0.572), as well as P4s’ (p = 0.613) dentinal crowns, EDJ (P3, p = 0.53; p4, p = 0.87), cervical outlines (P3, p = 0.45; P4 = 0.06) and the P3s’ crown outlines (p = 0.17) (Supplementary Tables B and C in S1 File). We have found significant results only in the crown outline analysis for P4s (p = 0.03).

Covariation

The 2B-PLS analysis demonstrated a strong pairwise correlation (r1 = 0.83) between P3 and P4 dentinal crowns (Fig 3), which showed concurrent trends of dentinal crown shape variation. For the Singular Warp scores 1 of the dentinal crown analysis (explaining 81% of total covariance), both tooth types covaried between short-crowned with mesio-distally broad base or tall-crowned with mesio-distally narrow base. The occlusal aspect associated with the short and broad crown type was mesio-distally expanded and squared with a broad central groove. Tall- and slender-crowned teeth, on the other hand, showed mesio-distally reduced occlusal aspects and were lingually narrow. Furthermore, higher dentinal horn tips were observed in short-crowned P4s with respect to tall-crowned premolars.

thumbnail
Fig 3. 2B-PLS plot capturing the covariation of upper third and fourth premolar dentinal crowns (dataset combining enamel-dentine junction and cervical outline; the warping shows the real shape variation in occlusal and distal view at the extremities of the range of distribution).

https://doi.org/10.1371/journal.pone.0301482.g003

The covariation between the EDJ occlusal aspect and the crown base (represented by the cervical outline) within tooth types was not as high as the covariation between the dentinal crowns of the two different tooth types. Nevertheless, pairwise correlation within the P4s parameters was higher (r1 = 0.56) than within P3s (r1 = 0.46) (Table 2).

thumbnail
Table 2. Results of the 2B-PLS (single warp score 1) for the upper third and fourth premolars (P3 and P4, respectively).

(EDJ = enamel-dentine junction).

https://doi.org/10.1371/journal.pone.0301482.t002

Size

The lnCS was compared across the sample separately for the dentinal crown, cervical and crown outlines, and EDJ (Fig 4, Supplementary Fig E in S1 File). Kruskal-Wallis test showed significant differences for lnCS between all populations (Table 3). Comparing the groups against each other, the Mann-Whitney-U test showed that the Oceanian teeth were significantly larger than any other group for all features except the EDJ analysis. Similarly, Europeans and Africans were consistently smaller. Supplementary Table D in S1 File shows size differences between sexes in upper premolars, assessed via Mann-Whitney-U test.

thumbnail
Fig 4. Boxplots of the natural logarithm of centroid sizes from the dentinal crown.

https://doi.org/10.1371/journal.pone.0301482.g004

thumbnail
Table 3. Results of the Kruskal–Wallis test.

Differences in dental sizes (expressed by the logarithm of Centroid Sizes) for the upper third and fourth premolars (P3 and P4, respectively) among all populations across the sample. (EDJ = enamel-dentine junction).

https://doi.org/10.1371/journal.pone.0301482.t003

Multivariate regression of size on shape for dentinal crowns was performed to analyse allometry and showed that only 1.9% and 2.1% of the total P3, respectively P4 shape variance could be explained by size.

Non-metric traits

Table 4 reports the prevalence of the premolar mesial and distal accessory cusps, accessory ridges and essential crest expressions on both the OES and the EDJ. Mesial and distal accessory cusps were present in only 10% of our sample. None of the Southeast Asian individuals exhibited this feature, and within the European, Near Eastern, South American and African samples, it could only be detected in a few individuals either on the OES or the EDJ of the P3 or P4. Various expressions of the accessory cusps on the OES were three times more frequent in P4s (12%) than in P3s (4%). Only five of these specimens expressed also an elevation of the marginal ridge on the EDJ. Interestingly, three accessory cusps were detected on the EDJ of the P3s and another three in the P4s, without any sign of expression on the OES.

thumbnail
Table 4. List of non-metric traits scored for the third and fourth premolars (P3 and P4, respectively).

(Numbers 0,1,2 and 3 refer to the manifestation (see Methods); n = Number of individuals).

https://doi.org/10.1371/journal.pone.0301482.t004

MxPARs were present in 55% of P3s and 65% of P4s. In P3s, a distally located accessory ridge occurred commonly on both OES and EDJ. Conversely, most of the P4 specimens expressed accessory ridges in both mesial and distal positions. The highest frequency of expression of accessory ridges in P3s was found in South Americans (50% on the OES, 75% on the EDJ), while in P4s, Europeans possessed most MxPARs (76% on the OES, 70% on the EDJ).

Essential crest expression was homogenously distributed on both surfaces in both tooth types. Approximately 90% of the individuals possessed a simple essential crest, separated by the sagittal sulcus. The prevalence of the transverse crests and bifurcated crest expressions were comparable among the samples in both tooth types on both OES and EDJ.The results of the Chí2 test performed for every non-metric trait in both tooth types for both surfaces showed no significant difference in trait occurrence.

Discussion

Dental development is generally under strong genetic influence resulting in constrained morphological variation of the dentition [7173]. The applied 3D methodology delivered fine-scale results about the morphometric variation of premolars’ crowns. Maxillary premolars varied mainly in the length-to-breadth ratio of the occlusal aspect as well as the height of the dentinal crown and horn tips. A broader tooth shape is associated with short crowns but higher horn tips positioned closer to each other. On the contrary, narrow teeth displayed a tall crown but lower horn tips more apart from each other.

Based on our analyses of the dentinal crown, it was evident that the populations in our sample could not be characterized or distinguished by their upper premolars crown morphology, despite their different population history. A subtle divergence could be detected visually in the PCA plots for the P4 dentinal crown, separating Southeast Asians from Oceanians and South Americans, but a statistical shape difference between these groups could not be confirmed by the PERMANOVA. Nonetheless, it is interesting that similar observations were made for lower premolars [54]. However, this outcome should be confirmed using a larger sample since our results are based on 3 specimens only. Additionally, the overlap between South American and Near Eastern populations was minor. Prevailing short and broad posterior dental crown types including its extreme expressions, like in the Near Easterners, were already described for African populations [46,54]. As mentioned above, small sample sizes of various populations have to be acknowledged, as well as the different number of sub-samples within geographical populations. Our goal for this study was to create a sample containing as much diversity as possible, thus including more sub-samples from geographical regions. In this sense, our sample composition works against its harmonization, nevertheless we find that these world populations have very similar upper premolar shape variation. Since our research requires dental specimens with minimal signs of wear and abrasion, no dental treatments, and high-resolution scanning, living or recently deceased individuals are not usable for our study. While these are restrictions not in our favour, our approach allowed for an accurate 3D shape analyses of dentinal crowns, which are not accessible with traditional methods.

Compared to the dentinal crown and EDJ analyses, PERMANOVA analysis of the 2D outline datasets, revealed a couple of significant differences in group centroids in five out of sixty pairs of populations (Supplementary Table B in S1 File). Since the cervical and crown outline analyses do not contain 3D information, and thus are more similar to traditional linear measurements that have been proven nonreliable in population variation analyses [14], we suggest using them only in the absence of combined datasets rather than alone.

We interpret the strong covariation between P3 and P4 shapes (r1 = 0.83), both either broad and short, or narrow and tall, as a sign of a strong genetic component during tooth development in an environment with shared genetic influence [72,73]. This high correlation might be surprising, considering that P3s show a rather prominent buccal cusp and share the tearing function with the canines, while the P4s’ more even cusps engage in grinding activities, like in molars [8]. These are the first 3D outcomes on upper premolars and they are in agreement with those obtained for the lower premolars using a similar protocol [54], showing the same general pattern of variation. Thus, combining the evidence from our study and that of Krenn et al. [54], the dentinal crowns of both upper and lower P3s and P4s showed an extensive overlap between diverse populations that vary from short and broad to tall and narrow. European and Southeast Asian upper and lower premolars expressed extreme variations of the tall crown variant in both analyses. However, we found the relative cusp position of both cusps in upper P4s rather stable, with changes observable mainly in horn height, while in lower P4s, significant variation of the mesio-distal position of the lingual cusp was observed. The covariation expressed by upper and lower premolars can be explained by the need for morphological compatibility to guarantee stable occlusion and efficient function [6]. Since only the buccal cusp of the lower premolars is engaged at maximum intercuspation of the dentitions, the lingual cusp is less morphologically constrained [6,9,74]. Differently, in upper premolars both cusps are directly involved in occlusion [6], thus a lower degree of morphological variation might be expected. Consequently, the lingual aspect of the buccal cusp (paracone) and both the lingual and buccal aspects of the lingual cusp (protocone) should be more stable than other parts of the crown [3,6] because they are directly involved in occlusion.

Such constraints in shape variation could also be reflected in more local features. In our non-metric trait analysis, we observed a higher rate of trait expression on the paracone than on the protocone, in the central groove and on the marginal ridges, which would be in line with the predictions by Kraus et al. [6]. However, we also found considerable variation in the lingual facet of the paracone. Therefore, not all of our results support the assumption of higher morphological conservation due to functional constraints. On the other hand, non-metric traits are understood to be under strong genetic control and vary neutrally in comparison to metric traits [14,75], therefore we should not expect them to be strictly related to function only. Our findings were mostly consistent with those of the previous studies on non-metric traits focusing on the OES [2,3,13,2123,27,28,68,7678] as well as studies based on OES and EDJ [4851,69]. However, in contrast to our results, Sakai et al. [69] reported that mesial and distal accessory cusps occurred more often in P3s than in P4s, working, however, only on Japanese dentitions. Overall, the rate of expression of accessory cusps in different human populations is heterogeneous, most of them observed in Oceanian individuals, while Southeast Asians and Europeans showed nearly none. Our findings are in agreement with Xing et al. [49], who found that accessory cusps on the OES are not necessarily associated with elevations on the EDJ. Turner et al. [68] also stated that there is no dentine involvement necessary in the formation of mesial and distal accessory cusps. Uneven enamel disposition can also correspond to accessory cusps, as well as other features displayed on the OES [38,79,80]. Similar patterns can be found in the derived forms of the essential crest that are more common in the African sample. We detected all different manifestations of the essential crest on the EDJ and OES of both tooth types; however, more often in P3s. Accessory ridges are also equally present on the EDJ and OES, where more than half of the P3s in our sample and nearly two-thirds of the P4s possessed an accessory ridge. This is consistent with the results of Burnett et al. [70] and Mihailidis et al. [28], which were, however, solely based on the OES. In agreement with these studies, accessory ridges in our sample were primarily present on the distal aspect of the tooth. In P3s, we observed the highest occurrence of accessory ridges in Europeans, followed by Africans and Near Easterners. In contrast to our work, Burnett et al. [70] found the highest frequency of accessory ridges in North-East Asian or Asian-derived populations and lower frequencies in Africans and Europeans. For the P4s, again, the distribution of accessory ridges in our sample was unequal, while Burnett et al. [70] described a relatively even distribution across the whole sample.

Regarding premolar crown size, Townsend’s [81] findings showed a significant sexual dimorphism for both maxillary premolars in the mesiodistal and buccolingual crown dimensions, with females possessing more broad/circular tooth crowns and males more narrow/ellipsoid crowns. We could not verify these results in our study due to the lack of reliable information about sex for most of our individuals. However, with regard to dental size, as represented by the lnCS, we did not find significant differences between males and females for the upper P4s, which parallels the results in Krenn et al. [54] for the lower premolars from a partially overlapping human sample. Conversely, we found significant differences for the upper P3s, but since our sample included more males (n = 26) than females (n = 9) our results should be interpreted with caution because of the uneven composition of the sample. According to Brace et al. [82] and Brace and Mahler [83], there is a gradient in increasing tooth size from north to south of the world. We indeed observed the smallest dental size in Europeans, while the largest teeth of our sample belonged to the Oceanian population. However, we also found African teeth of small size, which would not support Brace’s theory, but is in agreement with Hanihara and Ishida [24], who found sub-Saharan African dental sizes close to the average. Regardless of the unverifiable impact of sex on our results, we can state that static allometry in maxillary premolars is generally quite small in our sample. The effect of size was also small, explaining about 2% of the shape variation within our sample.

Our study generally demonstrated a generally large morphological variation of modern human upper premolars, yet a high correlation between first and second upper premolar. This variation does not seem to be the result of geographical separation in the course of modern human evolution, since most studied populations from different continents, climatic zones and environments did not differ from each other strong enough to be discriminated based on their upper premolar morphology. Our results instead show that the upper premolar morphological variation of recent modern humans is very similar among world populations, and thus, based on our data, we accept the tested hypothesis that the ranges of upper premolar shape variation of various human populations overlap and they cannot be distinguished based on their upper premolar’s dentinal crown morphology only. For evolutionary studies, comparing modern humans with archaic hominins or non-human primates, this means that the geographic composition of the comparative modern human sample is of lesser importance. More work is needed to unravel the mechanisms behind this concordance of premolar shape independent of geographical origin.

Supporting information

S1 File. Šimková et al.

Supporting information.

https://doi.org/10.1371/journal.pone.0301482.s001

(DOCX)

Acknowledgments

We are grateful to Sabine Eggers, Karin Wiltschke-Schrotta, Eduard Winter and Maria Teschler Nicola (Natural History Museum, Vienna, Austria), and Israel Hershkovitz (Department of Anatomy and Anthropology, Sackler Faculty of Medicine, Tel Aviv University, Israel) who granted access to the materials. Katarina Matiasek provided background information on collections and Dominika Teplanova assisted with scanning. Special gratitude goes to Martin Dockner (Vienna μCT Lab) for his technical support.

References

  1. 1. Alt KW, Rösing FW, Teschler-Nicola M. Dental anthropology: fundamentals, limits and prospects. Springer Science & Business Media; 2012.
  2. 2. Bailey SE. Neandertal dental morphology: implications for modern human origins. Arizona State University; 2002.
  3. 3. Hillson S. Dental anthropology. Cambridge University Press; 1996.
  4. 4. Irish JD, Nelson GC. Technique and application in dental anthropology. Cambridge University Press; 2008.
  5. 5. Kelley MA, Larsen CS. Advances in dental anthropology. Wiley; 1991.
  6. 6. Kraus BS, Jordan RE, Abrams L. Dental anatomy and occlusion: A study of the masticatory system. Baltimore: Baltimore, MD: Williams and Wilkins; 1969.
  7. 7. Lucas PW. Dental functional morphology: how teeth work. Cambridge University Press; 2004.
  8. 8. Nelson SJ, Ash M. Wheeler’s dental anatomy, physiology and occlusion. St. Louis: Elsevier; 2010.
  9. 9. Scott GR, Turner CG, Townsend GC, Martinón-Torres M. The anthropology of modern human teeth: Dental morphology and its variation in recent and fossil Homo sapiens. Cambridge: Cambridge University Press; 2018.
  10. 10. Scott GR, Irish JD. A companion to dental anthropology. John Wiley & Sons; 2015.
  11. 11. Scott GR, Irish JD. Human tooth crown and root morphology. Cambridge University Press; 2017.
  12. 12. Scott GR, Turner CG. Dental Anthropology. Annu Rev Anthropol. 1988;17: 99–126.
  13. 13. Scott GR, Turner CG. Anthropology of modern human teeth. Cambridge University Press Cambridge; 1997.
  14. 14. Rathmann H, Perretti S, Porcu V, Hanihara T, Scott GR, Irish JD, et al. Inferring human neutral genetic variation from craniodental phenotypes. PNAS Nexus. 2023;2: pgad217. pmid:37457893
  15. 15. Scott GR, Maier C, Heim K. Identifying and Recording Key Morphological (Nonmetric) Crown and Root Traits. A Companion to Dental Anthropology. 2015. pp. 245–264.
  16. 16. Ahmad IA, Alenezi MA. Root and Root Canal Morphology of Maxillary First Premolars: A Literature Review and Clinical Considerations. J Endod. 2016;42: 861–872. pmid:27106718
  17. 17. de Lima CO, de Souza LC, Devito KL, do Prado M, Campos CN. Evaluation of root canal morphology of maxillary premolars: a cone-beam computed tomography study. Aust Endod J. 2019;45: 196–201. pmid:30230115
  18. 18. Neelakantan P, Subbarao C, Ahuja R, Subbarao CV. Root and canal morphology of Indian maxillary premolars by a modified root canal staining technique. Odontology. 2011;99: 18–21. pmid:21271321
  19. 19. Özcan E, Çolak H, Hamidi MM. Root and canal morphology of maxillary first premolars in a Turkish population. J Dent Sci. 2012;7: 390–394.
  20. 20. Saber SEDM Ahmed MHM, Obeid M Ahmed HMA. Root and canal morphology of maxillary premolar teeth in an Egyptian subpopulation using two classification systems: a cone beam computed tomography study. Int Endod J. 2019;52: 267–278. pmid:30225932
  21. 21. Morris DH. Maxillary first premolar angular differences between North American Indians and non-North American Indians. Am J Phys Anthropol. 1981;54: 431–433. pmid:7234988
  22. 22. Lavelle CLB. A metrical comparison of maxillary first premolar form. Am J Phys Anthropol. 1984;63: 397–403. pmid:6731607
  23. 23. Gómez-Robles A, Martinón-Torres M, Bermúdez de Castro JM, Prado-Simón L, Arsuaga JL. A geometric morphometric analysis of hominin upper premolars. Shape variation and morphological integration. J Hum Evol. 2011;61: 688–702. pmid:22047673
  24. 24. Wood BA, Engleman CA. Analysis of the dental morphology of Plio-Pleistocene hominids. V. Maxillary postcanine tooth morphology. J Anat. 1988;161: 1–35. pmid:3254883
  25. 25. Hanihara T, Ishida H. Metric dental variation of major human populations. Am J Phys Anthropol Off Publ Am Assoc Phys Anthropol. 2005;128: 287–298. pmid:15838862
  26. 26. Pilloud MA, Hefner JT, Hanihara T, Hayashi A. The Use of Tooth Crown Measurements in the Assessment of Ancestry. J Forensic Sci. 2014;59: 1493–1501. pmid:25060236
  27. 27. Burnett SE. A New Look at Premolar Trait Variation: Maxillary Premolar Accessory Ridges (MxPAR). Dent Anthopology Publ Dent Anthropol Assoc. 1996.
  28. 28. Mihailidis S, Scriven G, Khamis M, Townsend G. Prevalence and patterning of maxillary premolar accessory ridges (MxPARs) in several human populations. Am J Phys Anthropol. 2013;152: 19–30. pmid:23913636
  29. 29. Adams D, Swenson VM, Scott GR. Global Distribution of Marginal Accessory Tubercles of the Maxillary Premolars. Dent Anthropol J. 2019;32: 8–15.
  30. 30. Liao W, Xing S, Li D, Martinon-Torres M, Wu X, Soligo C, et al. Mosaic dental morphology in a terminal Pleistocene hominin from Dushan Cave in southern China. Sci Rep. 2019;9: 1–14.
  31. 31. Fornai C, Bookstein FL, Weber GW. Variability of Australopithecus second maxillary molars from Sterkfontein Member 4. J Hum Evol. 2015;85: 181–192. pmid:26163295
  32. 32. Korenhof C. Evolutionary trends of the inner enamel anatomy of deciduous molars from Sangiran (Java, Indonesia). Teeth Funct Evol. 1982; 350–365.
  33. 33. Monson TA, Fecker D, Scherrer M. Neutral evolution of human enamel–dentine junction morphology. Proc Natl Acad Sci. 2020;117: 26183–26189. pmid:33020281
  34. 34. Nager G. Der Vergleich zwischen dem räumlichen Verhalten des Dentin-Kronenreliefs und dem Schmelzrelief der Zahnkrone. Cells Tissues Organs. 1960;42: 226–250.
  35. 35. Olejniczak A J, Martin L B, Ulhaas L. Quantification of dentine shape in anthropoid primates. Ann Anat—Anat Anz. 2004;186: 479–485. pmid:15646281
  36. 36. Sasaki K, Makiguchi H, Yoshino K, Kanazawa E. Morphological Study of the Dentino-enamel Junction of the Upper First Deciduous Molar. Anthropol Sci. 2002;110: 349–363.
  37. 37. Skinner MM, Wood BA, Hublin J-J. Protostylid expression at the enamel-dentine junction and enamel surface of mandibular molars of Paranthropus robustus and Australopithecus africanus. J Hum Evol. 2009;56: 76–85. pmid:18986683
  38. 38. Skinner MM, Evans A, Smith T, Jernvall J, Tafforeau P, Kupczik K, et al. Brief communication: Contributions of enamel-dentine junction shape and enamel deposition to primate molar crown complexity. Am J Phys Anthropol. 2010;142: 157–163. pmid:20091839
  39. 39. Skinner MM, Gunz P, Wood BA, Hublin J-J. Enamel-dentine junction (EDJ) morphology distinguishes the lower molars of Australopithecus africanus and Paranthropus robustus. J Hum Evol. 2008;55: 979–988. pmid:18824253
  40. 40. Skinner MM, Wood BA, Boesch C, Olejniczak AJ, Rosas A, Smith TM, et al. Dental trait expression at the enamel-dentine junction of lower molars in extant and fossil hominoids. Dent Tissue Stud 2D 3D Insights Hum Evol. 2008;54: 173–186. pmid:18048081
  41. 41. Bailey SE, Skinner MM, Hublin J-J. What lies beneath? An evaluation of lower molar trigonid crest patterns based on both dentine and enamel expression. Am J Phys Anthropol. 2011;145: 505–518. pmid:21312178
  42. 42. Guy F, Lazzari V, Gilissen E, Thiery G. To What Extent is Primate Second Molar Enamel Occlusal Morphology Shaped by the Enamel-Dentine Junction? PLOS ONE. 2015;10: e0138802. pmid:26406597
  43. 43. Gómez-Robles A, Bermúdez de Castro JM, Arsuaga J-L, Carbonell E, Polly PD. No known hominin species matches the expected dental morphology of the last common ancestor of Neanderthals and modern humans. Proc Natl Acad Sci. 2013;110: 18196. pmid:24145426
  44. 44. Weber GW, Fornai C, Gopher A, Barkai R, Sarig R, Hershkovitz I. The Qesem Cave hominin material (part 1): A morphometric analysis of the mandibular premolars and molar. Quat Int. 2016;398: 159–174.
  45. 45. Benazzi S, Coquerelle M, Fiorenza L, Bookstein F, Katina S, Kullmer O. Comparison of dental measurement systems for taxonomic assignment of first molars. Am J Phys Anthropol. 2011;144: 342–354. pmid:21302262
  46. 46. Šimková PG, Weber GW, Ramirez Rozzi FV, Slimani L, Sadoine J, Fornai C. Morphological variation of the deciduous second molars in the Baka Pygmies. Sci Rep. 2021;11: 16480. pmid:34389746
  47. 47. Hershkovitz I, Weber GW, Quam R, Duval M, Grün R, Kinsley L, et al. The earliest modern humans outside Africa. Science. 2018;359: 456–459. pmid:29371468
  48. 48. Liu W, Schepartz LA, Xing S, Miller-Antonio S, Wu X, Trinkaus E, et al. Late Middle Pleistocene hominin teeth from Panxian Dadong, South China. J Hum Evol. 2013;64: 337–355. pmid:23465337
  49. 49. Xing S, Sun C, Martinón-Torres M, Bermúdez de Castro JM, Han F, Zhang Y, et al. Hominin teeth from the Middle Pleistocene site of Yiyuan, Eastern China. J Hum Evol. 2016;95: 33–54. pmid:27260173
  50. 50. Xing S, Martinón-Torres M, Bermúdez de Castro JM, Wu X, Liu W. Hominin teeth from the early Late Pleistocene site of Xujiayao, Northern China. Am J Phys Anthropol. 2015;156: 224–240. pmid:25329008
  51. 51. Xing S, Martinón-Torres M, Bermúdez de Castro JM, Zhang Y, Fan X, Zheng L, et al. Middle Pleistocene Hominin Teeth from Longtan Cave, Hexian, China. PLOS ONE. 2015;9: e114265. pmid:25551383
  52. 52. Fornai C, Benazzi S, Gopher A, Barkai R, Sarig R, Bookstein FL, et al. The Qesem Cave hominin material (part 2): A morphometric analysis of dm2-QC2 deciduous lower second molar. State Art Multidiscip Res Middle Pleistocene Qesem Cave Isr 2015. 2016;398: 175–189.
  53. 53. Davies TW, Delezene LK, Gunz P, Hublin J-J, Skinner MM. Endostructural morphology in hominoid mandibular third premolars: Geometric morphometric analysis of dentine crown shape. J Hum Evol. 2019;133: 198–213. pmid:31358180
  54. 54. Krenn VA, Fornai C, Wurm L, Bookstein FL, Haeusler M, Weber GW. Variation of 3D outer and inner crown morphology in modern human mandibular premolars. Am J Phys Anthropol. 2019;169: 646–663. pmid:31099892
  55. 55. Kato A, Kouchi M, Mochimaru M, Isomura A, Ohno N. A geometric morphometric analysis of the crown form of the maxillary central incisor in humans. Dent Anthropol J. 2011;24: 1–10.
  56. 56. Malaspinas A-S, Westaway MC, Muller C, Sousa VC, Lao O, Alves I, et al. A genomic history of Aboriginal Australia. Nature. 2016;538: 207–214. pmid:27654914
  57. 57. Mallick S, Li H, Lipson M, Mathieson I, Gymrek M, Racimo F, et al. The Simons Genome Diversity Project: 300 genomes from 142 diverse populations. Nature. 2016;538: 201–206. pmid:27654912
  58. 58. Rasmussen M, Guo Xiaosen, Wang Yong, Lohmueller Kirk E., Rasmussen Simon, Albrechtsen Anders, et al. An Aboriginal Australian Genome Reveals Separate Human Dispersals into Asia. Science. 2011;334: 94–98. pmid:21940856
  59. 59. Molnar S. Human tooth wear, tooth function and cultural variability. Am J Phys Anthropol. 1971;34: 175–189. pmid:5572602
  60. 60. Frederick M, Gallup G. Fluctuating dental asymmetry in great apes, fossil hominins, and modern humans: Implications for changing stressors during human evolution. 2007.
  61. 61. Benazzi S, Panetta D, Fornai C, Toussaint M, Gruppioni G, Hublin J-J. Technical Note: Guidelines for the digital computation of 2D and 3D enamel thickness in hominoid teeth. Am J Phys Anthropol. 2014;153: 305–313. pmid:24242830
  62. 62. Bookstein FL. Principal warps: Thin-plate splines and the decomposition of deformations. IEEE Trans Pattern Anal Mach Intell. 1989;11: 567–585.
  63. 63. Gunz P, Mitteroecker P, Bookstein FL. Semilandmarks in Three Dimensions. In: Slice DE, editor. Modern Morphometrics in Physical Anthropology. Boston, MA: Springer US; 2005. pp. 73–98. https://doi.org/10.1007/0-387-27614-9_3
  64. 64. Gunz P, Mitteroecker P. Semilandmarks: a method for quantifying curves and surfaces. Hystrix Ital J Mammal. 2013;24: 103–109.
  65. 65. Gower JC. Generalized procrustes analysis. Psychometrika. 1975;40: 33–51.
  66. 66. Bookstein FL. Thin-plate splines and the atlas problem for biomedical images. Springer; 1991. pp. 326–342.
  67. 67. Mitteroecker P, Gunz P. Advances in geometric morphometrics. Evol Biol. 2009;36: 235–247.
  68. 68. Turner CGI, Nichol CR, Scott GR. Scoring produces for key morphological traits of the permanent dentition: The Arizona State University dental anthropology system. Adv Dent Anthropol. 1991; 13–31.
  69. 69. Sakai T, Sasaki I, Hanamura H. A Morphological Study of Enamel-Dentin Border on the Japanese Dentition. J Anthropol Soc Nippon. 1967;75: 207–223.
  70. 70. Burnett SE, Hawkey DE, Turner CG II. Brief communication: Population variation in human maxillary premolar accessory ridges (MxPAR). Am J Phys Anthropol. 2010;141: 319–324. pmid:19953530
  71. 71. Brook AH. Multilevel complex interactions between genetic, epigenetic and environmental factors in the aetiology of anomalies of dental development. Int Workshop Oral Growth Dev. 2009;54: S3–S17. pmid:19913215
  72. 72. Jernvall J, Thesleff I. Reiterative signaling and patterning during mammalian tooth morphogenesis. Mech Dev. 2000;92: 19–29. pmid:10704885
  73. 73. Tucker A, Sharpe P. The cutting-edge of mammalian development; how the embryo makes teeth. Nat Rev Genet. 2004;5: 499–508. pmid:15211352
  74. 74. Alt KW, Rösing FW, Teschler-Nicola M. Dental Anthropology—An Introduction. In: Alt KW, Rösing FW, Teschler-Nicola M, editors. Dental Anthropology: Fundamentals, Limits and Prospects. Vienna: Springer Vienna; 1998. pp. 1–3. https://doi.org/10.1007/978-3-7091-7496-8_1
  75. 75. Hanihara T. Morphological variation of major human populations based on nonmetric dental traits. Am J Phys Anthropol Off Publ Am Assoc Phys Anthropol. 2008;136: 169–182. pmid:18257017
  76. 76. Delgado-Burbano ME, Scott GR, Turner CG II. The Uto-Aztecan premolar among North and South Amerindians: Geographic variation and genetics. Am J Phys Anthropol. 2010;143: 570–578. pmid:20623681
  77. 77. Johnson KM, Stojanowski CM, Miyar KO, Doran GH, Ricklis RA. New evidence on the spatiotemporal distribution and evolution of the Uto‐Aztecan premolar. Am J Phys Anthropol. 2011;146: 474–480. pmid:21994018
  78. 78. Martinón-Torres M, Bastir M, Bermúdez de Castro JM, Gómez A, Sarmiento S, Muela A, et al. Hominin lower second premolar morphology: evolutionary inferences through geometric morphometric analysis. J Hum Evol. 2006;50: 523–533. pmid:16472839
  79. 79. Corruccini RS. The dentinoenamel junction in primates. Int J Primatol. 1987;8: 99.
  80. 80. Kono RT, Suwa G, Tanijiri T. A three-dimensional analysis of enamel distribution patterns in human permanent first molars. Arch Oral Biol. 2002;47: 867–875. pmid:12450518
  81. 81. Townsend GC. Intercuspal Distances of Maxillary Pre-molar Teeth in Australian Aboriginals. J Dent Res. 1985;64: 443–446. pmid:3855896
  82. 82. Brace CL, Rosenberg KR, Hunt KD. Gradual change in human tooth size in the late Pleistocene and post-Pleistocene. Evolution. 1987;41: 705–720. pmid:28564360
  83. 83. Brace CL, Mahler PE. Post-Pleistocene changes in the human dentition. Am J Phys Anthropol. 1971;34: 191–203. pmid:5572603