Disentangling isolated dental remains of Asian Pleistocene hominins and pongines

Scholars have debated the taxonomic identity of isolated primate teeth from the Asian Pleistocene for over a century, which is complicated by morphological and metric convergence between orangutan (Pongo) and hominin (Homo) molariform teeth. Like Homo erectus, Pongo once showed considerable dental variation and a wide distribution throughout mainland and insular Asia. In order to clarify the utility of isolated dental remains to document the presence of hominins during Asian prehistory, we examined enamel thickness, enamel-dentine junction shape, and crown development in 33 molars from G. H. R. von Koenigswald's Chinese Apothecary collection (11 Sinanthropus officinalis [= Homo erectus], 21 “Hemanthropus peii,” and 1 “Hemanthropus peii” or Pongo) and 7 molars from Sangiran dome (either Homo erectus or Pongo). All fossil teeth were imaged with non-destructive conventional and/or synchrotron micro-computed tomography. These were compared to H. erectus teeth from Zhoukoudian, Sangiran and Trinil, and a large comparative sample of fossil Pongo, recent Pongo, and recent human teeth. We find that Homo and Pongo molars overlap substantially in relative enamel thickness; molar enamel-dentine junction shape is more distinctive, with Pongo showing relatively shorter dentine horns and wider crowns than Homo. Long-period line periodicity values are significantly greater in Pongo than in H. erectus, leading to longer crown formation times in the former. Most of the sample originally assigned to S. officinalis and H. erectus shows greater affinity to Pongo than to the hominin comparative sample. Moreover, enamel thickness, enamel-dentine junction shape, and a long-period line periodicity value in the “Hemanthropus peii” sample are indistinguishable from fossil Pongo. These results underscore the need for additional recovery and study of associated dentitions prior to erecting new taxa from isolated teeth.


Introduction
Fossil remains attributed to Homo erectus have been recovered from numerous localities in mainland Asia and island Southeast Asia driven, in part, by intense scholarly interest during the last century in a potential Asian origin for the genus Homo. These recoveries derive from more than a dozen sites with well-preserved hominin cranial material (reviewed in [1][2][3]). More problematic are localities where only isolated teeth or fragmentary jaws have been discovered (e.g., Longgupo, China: [4]; most Vietnamese localities) ( Table A in S1 File), as H. erectus likely coexisted with large-bodied apes [5,6]. Morphological and metrical convergence between Asian hominin and pongine post-canine dental material has confounded paleoanthropologists for decades [5,[7][8][9][10][11][12]. This is exemplified by the phenomenon of variably wrinkled (crenulated) molar occlusal surfaces with low cuspal relief.
Taxonomic assessments of mixed primate faunal collections, such as those from Chinese Apothecaries [13,14], the Sangiran dome [5,15], and Mohui Cave [12,[16][17][18] have yielded varied opinions about the affinities of these isolated dental remains. An example is von Koenigswald's identification of a new hominin "Hemanthropus peii" from a subset of his Chinese Apothecary collection [14,19], which has since been considered to be an australopithecine ( [20,21], but see [22,23]), a member of the genus Homo [24], an assortment of fossil orangutans [9,12,25], and/or a new type of hominoid [9,12,18]. Von Koenigswald [14] originally reported that this taxon was "difficult to define" and distinguished several isolated molar teeth from orangutans solely on the lack of "fine wrinkles" and resulting "better defined cusps," promising a more detailed description that was never published. Yet both orangutans and humans show marked intraspecific variation in the expression of crenulations and the underlying cuspal morphology [7,26]. Moreover crenulations are compromised by wear, making this a tenuous morphological feature on which to erect a new genus and species.
This study quantifies aspects of internal tooth structure and development using nondestructive imaging to examine the taxonomic affiliation of isolated teeth from the Asian Pleistocene. Recent studies have explored aspects of tooth structure and development with propagation phase contrast synchrotron micro-computed tomography (SR microCT) and conventional microCT for better resolution of Pleistocene Asian primates than is possible from traditional assessments of tooth size and shape [30,33,36,[41][42][43]. We have expanded this approach in the current study by considering a suite of dental characteristics in two historic collections of teeth, which are compared to a large sample of recent and fossil orangutans and hominins.

Imaging
We initially employed a Skyscan 1172 microCT housed in the Department of Human Evolution at the Max Planck Institute for Evolutionary Anthropology to scan 33 teeth from the Chinese Apothecary collections and 57 teeth from the Sangiran dome at 100 kV and 100 mA, with a metallic filter (0.04 mm of copper plus 0.5 mm of aluminum) and an isometric voxel size of 14-27 μm. Unfortunately, due to diagenetic modification [58,59], it was not possible to virtually distinguish the interface between enamel and dentine in the great majority of the Sangiran dome dental sample. Seven teeth from the Sangiran dome collection, along with 11 Sinanthropus officinalis (= Homo erectus) teeth and the "Hemanthropus peii" type (CA 673) from the Chinese Apothecary collection, were subsequently scanned on beamline ID 19 at the European Synchrotron Radiation Facility (ESRF, Grenoble, France) following the protocol detailed below. Four H. erectus teeth from Zhoukoudian and one Chinese fossil orangutan tooth (CA 260) were also similarly imaged as part of the comparative sample.
We used two optical configurations to reveal overall tooth crown structure and aspects of the microstructure. Scans of entire teeth were performed with a FReLoN 2K CCD camera coupled to a LuAG:Ce scintillator of 200 μm thick with an optical system yielding a pixel size of 4.96 μm. The beam was set at an average energy of 60.8 keV by filtering the white beam of the W150 wiggler set at a gap of 65 mm with 2 mm of aluminum and 0.25 mm of tungsten. We used a propagation distance of 4 m in order to reveal the general growth pattern and the enamel-dentine junction. The scans were performed in half-acquisition mode with 5000 projections of 0.3 seconds each over 360 degrees. Most of the high-resolution propagation phase contrast scans used to investigate the longperiod line periodicity were performed using a FReLoN E2V CCD camera coupled to a 10 μm thick GGG:Eu scintillator through a microscope optic yielding a pixel size of 0.6 μm. The beam was set at an average energy of 59.7 keV by filtering the white beam of the two U32 undulators closed respectively at 11.88 and 11.72 mm gaps with 2 mm of aluminum and 0.25 mm of tungsten. We used a propagation distance of 150 mm to reveal incremental lines with scans in half-acquisition mode using 5000 projections of 0.9 seconds each over 360 degrees. The H. erectus comparative sample was scanned using a different microscope optic adapted for white beam, producing a pixel size of 0.73 μm. In this case we used a FReLoN 2K14 CCD camera in frame transfer mode, a 24 μm thick GGG:eu scintillator, and 5000 projections of 0.3 seconds each in half-acquisition mode. The average energy of the beam (68.3 keV) was set by filtering the white beam of the ID19 W150 wiggler set at a gap of 35 mm by 3 mm of aluminum, 0.25 mm of copper and 0.06 mm of tungsten. A suitable phase contrast effect was obtained with a propagation distance of 150 mm.
For some teeth, the high resolution scans led to significant local darkening of the enamel. The samples were restored to their original color by illumination with low energy UV light (black light neon tubes) for a few hours [60]. Recent refinements subsequent to the acquisition of these data have reduced the radiation dose to avoid this darkening effect [61]. Volumes were reconstructed using filtered-backprojection algorithm (PyHST software, ESRF) in edge detection mode. Residual ring artefacts were corrected on reconstructed slices with a custom Matlab code [62]. The 3D renderings, virtual slices, and raw data employed in this study are available in the open access database for paleontology hosted by the ESRF (http://paleo.esrf.fr/ picture.php?/3229/category/2223).

Enamel structure analysis
Due to diagenesis and wear of the molars it was not possible to conduct three-dimensional (3D) measurements of enamel thickness [63][64][65], nor 3D enamel-dentine junction (EDJ) morphology [33,37,64,66], so a cross-sectional two-dimensional (2D) approach was taken. Virtual 2D section planes of the mesial cusps were generated from 3D models with VG Studio MAX 2.0 software (Volume Graphics, Inc.) according to published protocols [36,57,65,67]. Several variables were quantified on 2D section planes with Sigma Scan Pro software interfaced to a Wacom digitizing tablet following Martin [28,68]: enamel cap area (c), EDJ length (e), and coronal dentine area enclosed by the enamel cap (b) (Fig 1). Average enamel thickness (AET) was calculated as [c/e], yielding the average straight-line distance (in mm), or thickness, from the enamel-dentine junction to the outer enamel surface. AET was scaled for comparisons between taxa of different size by calculation of relative enamel thickness (RET; unitless index): Comparative molar enamel thickness data for recent Pongo (n = 135), fossil Pongo (n = 139), H. erectus (n = 4), and recent human molars (n = 271) was taken from [53,55,57,67].
Molar EDJ morphology was quantified by collecting nine landmarks and semi-landmarks from each intact and lightly worn cross-section [34,35] (Fig 1), and calculating a series of nine relative distances from these landmarks. These distances were combined with a database of homologous measurements of five H. erectus molars, 90 recent orangutans, 141 fossil orangutans, and 258 recent humans [30,35,53,55,57]. The relative distances were subjected to discriminant function analysis (DFA) using SPSS software (v. 21.0, IBM, Inc.). Variation in mesial cross-section EDJ shape from first to third molars is minimal and does not overwhelm the ability of this technique to distinguish taxa [35]. Moreover, as a number of these molars are of uncertain position within the dental arcade, molars from all three maxillary or mandibular positions were combined in the analyses. Sixteen putative hominin or pongine molars from the Chinese Apothecary and Sangiran dome collections were treated as ungrouped, or nonclassified, cases in the DFA, and the sample size of each taxon was not used to adjust the probability of molars belonging to any group. A separate DFA of molar EDJ shape in 19 "Hemanthropus peii" molars was run, which also included specimens CA 770 and 771 in the H. erectus sample since their maxillary molar EDJ shapes and long-period line periodicity values strongly support their hominin affinities [30].

Dental development analysis
Casts and photographs of the fossil sample revealed that very few of the tooth surfaces were sufficiently preserved for developmental assessment due to extreme postmortem pitting, etching and/or extensive tooth wear prior to death. Phase contrast SR microCT was employed to assess long-period line periodicity and crown formation times from virtual histological slices and 3D models of the original material (illustrated in [69]: Fig 1; [70]: Fig 8). Virtual histological slices (typically 30-100 μm thick) were produced following 3D optimized orientation to capture the full growth axis (dentine and pulpal horn tips). The periodicity, or number of daily lines between successive long-period incremental features (Fig 2), was examined in 23 fossil teeth. The Mann-Whitney U test was performed for comparisons of periodicity values between H. erectus and fossil Pongo teeth in the comparative sample.
Crown formation time was calculated for the few teeth that were sufficiently preserved as the sum of the cuspal and lateral enamel formation times, save for one instance when a The area of the enamel cap is represented as c, the area of the dentine under the enamel cap is represented as b, and the length of the enameldentine junction is represented as e. Landmarks are defined as follows: 1) tip of the lingual enamel cervix; 2) lingual intersection of the EDJ and a line parallel to the cervical diameter and bisecting the length between the cervical diameter and landmark 5; 3) lingual intersection of the EDJ and a line parallel to the cervical diameter and running through landmark 5; 4) mesiolingual dentine horn tip; 5) lowest point of the EDJ between the mesiolingual and mesiobuccal cusp tips; 6) mesiobuccal dentine horn tip; 7) buccal intersection of the EDJ and a line parallel to the cervical diameter and running through landmark 5; 8) buccal intersection of the EDJ and a line parallel to the cervical diameter and bisecting the length between the cervical diameter and landmark 5; and, 9) tip of the buccal enamel cervix. Landmark 1 was made to lie at x, y coordinate (0, 0) and landmark 9 at (0, 100) in every specimen examined in order to account for differences in tooth size. Reproduced from [30]. https://doi.org/10.1371/journal.pone.0204737.g001 Disentangling isolated dental remains of Asian Pleistocene hominins and pongines complete series of long-period lines was identified and counted from the coronal dentine and then multiplied by the long-period line periodicity (CA 796). Cuspal formation time was estimated by dividing the linear enamel thickness of each tooth cusp by 4.2 μm/day, which is the average cuspal daily secretion rate for fossil Pongo (4.19 μm/day, n = 9 teeth: [54]) and H. erectus (4.23 μm/day, n = 1-2 teeth: [40]). Direct measurements of cuspal daily secretion rates were not made due to the time-and data-intensive nature of quantifying cross-striation spacing virtually in the full thickness of cuspal enamel.
Lateral enamel formation time was calculated by multiplying the number of enamel longperiod lines (internal Retzius lines or external perikymata) by the long-period line periodicity (or periodicity range when a single integer was not discernible). Minor estimates of longperiod line numbers were made for worn or chipped enamel; specimens were excluded when these could not be counted for at least 90% of the total crown height. Comparative data on long-period line periodicity and crown formation times were chosen to minimize potential bias from methodological differences or interobserver error; data sources are provided in corresponding tables and figures. Because crown formation times vary among molar cusps and between upper and lower molars [56], comparisons among samples and taxa were limited to specific cusp types.  Disentangling isolated dental remains of Asian Pleistocene hominins and pongines fossil Pongo and Homo. The "Hemanthropus peii" sample (n = 19) shows the greatest similarity to fossil Pongo values, particularly for maxillary molars, with only a narrow overlap with Homo. The type specimen of "Hemanthropus peii" (CA 673) has a maxillary molar RET value of 13.1, which is similar to the fossil Pongo mean (13.8). Individual average and relative enamel thickness values are given in Table 2, along with range data for H. erectus and fossil Pongo.

Enamel-dentine junction morphology
Discriminant function analysis of molar EDJ landmark relative distances demonstrates that the Homo and Pongo comparative sample are grouped reliably at the generic level based on the nine distance ratios (91% of 247 maxillary teeth correctly classified; 96% of 247 mandibular teeth correctly classified). The main maxillary molar analysis had three significant functions with a combined Χ 2 (24) = 387.4, Wilk's λ = 0.199 (p < 0.001). After removal of the first function, there was still a strong association between groups and predictors: Χ 2 (14) = 96.8, Wilk's λ = 0.668 (p < 0.001). After removal of the second function, the significant relationship between groups and predictors persisted: Χ 2 (6) = 19.7, Wilk's λ = 0.921 (p = 0.003). A plot of the first two discriminant functions is shown in Fig 4. The first discriminant function accounts for 83.5% of the variance and has positive relationships with the buccal horn height relative to the center of the occlusal basin (landmark 5), the lingual horn height relative to the center of the occlusal basin, and the buccal dentine horn height. The second function accounts for 13.4% of the variance and has positive relationships with the width of the dentine midway between the crown base and the occlusal basin, the width of the dentine below the occlusal basin, the width of the dentine horns, and the lingual dentine horn height.
The discriminant function analysis of mandibular molars revealed two significant functions with a combined Χ 2 (16) = 449.9, Wilk's λ = 0.154 (p < 0.001). After removal of the first function, there was still a strong association between groups and predictors: Χ 2 (7) = 65.2, Wilk's λ = 0.763 (p < 0.001). The first discriminant function accounts for 92.7% of the variance and has positive relationships with the lingual dentine horn height, the buccal dentine horn height, the buccal horn height relative to the center of the occlusal basin, and the lingual horn height relative to the center of the occlusal basin. The second function accounts for 7.3% of the variance and has a positive relationship with the width of the dentine midway between the crown base and the occlusal basin, the width of the dentine below the occlusal basin, and the width of the dentine horns.
Only three of 16 putative hominin or pongine molars from the Chinese Apothecary and Sangiran dome collections were classified as Homo when left as ungrouped cases (CA 770, CA 771, S7-9), falling within the range of Homo sapiens molars and outside the Pongo range in two of three instances (Fig 4). Two mandibular molars fell within the range of H. sapiens but were classified as Pongo as their most likely group (CA 804, CA 816). A similar DFA analysis of the "Hemanthropus peii" sample resulted in the classification of 14 molars as fossil Pongo, four molars as extant Pongo (including the type specimen, CA 673), and one maxillary molar (CA 705) as H. sapiens (Fig 5).

Dental development
Long-period line periodicities in the Chinese Apothecary and Sangiran dome teeth range from 7 to at least 12 days (possibly 13 days) ( Table 3). The periodicity of the H. erectus comparative sample ranges from 6-8 days (n = 7), which is significantly lower than the fossil Pongo sample (range 9-12 days, n = 19) (Z = -3.934, p < 0.001; calculated using the highest possible values for those H. erectus teeth that are uncertain). The Sinanthropus officinalis type (CA 770) has a periodicity of 8 days, which is within the range of the H. erectus sample, and below the fossil  Tooth types as in Table 1. Values for CA 770 and 771 have changed slightly since their original publication [30] due to an update in sectioning protocols [63].
https://doi.org/10.1371/journal.pone.0204737.t002 Disentangling isolated dental remains of Asian Pleistocene hominins and pongines Pongo range. The value of the "Hemanthropus peii" type (CA 673) was somewhat indistinct, falling between 9-11 days, which exceeds the known range of H. erectus, but is within the fossil Pongo range. It was not possible to estimate the long-period line periodicity in two teeth from the study sample due to poor incremental feature visibility (S7-53, S7-76). Crown formation times were estimated for nine molar cusps of seven teeth in the Chinese Apothecary and Sangiran dome sample (Table 4), although three of these are minimum values (indicated by ">") due to natural tooth wear. Four crown formation times are given as multiple values due to uncertainty in the long-period line periodicity (CA 796 mesiolingual cusp, CA 804 mesiobuccal and mesiolingual cusps, S7-65 mesiolingual cusp). Comparative samples are limited to one H. erectus upper first mesiobuccal cusp, which formed in approximately 2.4 years [71], and ten fossil Pongo molar cusps, which range from 3.2-5.4 years ( [54] and this study]. Crown formation times for CA 770, CA 796, and S7-65 are similar to values for Homo (and extant Pongo), while CA 804 and CA 816 are similar to fossil Pongo (Fig 6).

Enamel structure
Enamel thickness appears to be of limited value for sorting hominin from pongine molars in mixed Asian faunas, as Pongo and Homo show overlapping ranges (see expanded extant data in [55,57,67]). This has been recognized since Miller's [27] radiographic study of the "Piltdown Man" jaw, which was ultimately determined to be a filed-down orangutan mandible [73]. Further complicating the issue is the presence of the Asian Pleistocene ape Gigantopithecus blacki; its relative enamel thickness values overlap thick-enameled fossil Pongo molars and most Homo Disentangling isolated dental remains of Asian Pleistocene hominins and pongines  [71]; H. sapiens: [72]; extant and fossil Pongo: [54] and this study.
https://doi.org/10.1371/journal.pone.0204737.g006 specimens [36,74] although its molar size and morphology are distinctive [14,44]. Studies of the three-dimensional distribution of enamel throughout the crown may provide better taxonomic resolution, as Pongo molar enamel appears to be more uniformly distributed than enamel in Homo molars [36,63,64]. This approach will be less useful for worn teeth or those that show poor tissue contrast with microCT scanning, as was the case for much of the fossil sample studied here.
While the presence of moderately thick enamel (sensu Martin [28]) does not distinguish Pongo and Homo molars, enamel thickness at the low end of the Pongo range-such as majority of the "Hemanthropus peii" sample-is more consistent with non-hominin hominoid primates. This enigmatic fossil group appears to have relatively thinner enamel than Neanderthals, the thinnest-enameled members of the genus Homo [65,67]. Evidence from molar enamel-dentine junction (EDJ) shape further supports their pongine status, as discriminant function analyses aligned 18 of 19 "Hemanthropus peii" molars with Pongo rather than Homo.
Our EDJ shape analysis confirms Martin's [28] observation that Pongo molars have shorter dentine horns than Homo, creating a flatter enamel-dentine junction and occlusal surface with less relief in the former taxon. Further study may clarify if there are other useful taxon-specific aspects of molar EDJ shape for sorting teeth in this region, as has been suggested by 3D analysis of maxillary molars of extant and fossil Pongo [75]. The small H. erectus maxillary molar comparative sample largely overlaps with both Pongo and recent H. sapiens, similar to the results of conventional metric and morphological analyses (reviewed in [30]). One tooth from our study sample showed greater affinities to H. erectus than to other groups (S7-9), and two were most like H. sapiens (CA 770, CA 771). These findings are consistent with information from hominin molar enamel thickness values, long-period line periodicities, and crown formation time (in the case of CA 770) ( Table 5), which we discuss further below.

Dental development
Developmental characters have been employed in taxonomic studies of hominins and hominoids for decades, but these studies have been complicated by methodological differences and a limited understanding of biological variation within and among serial tooth positions [39,56]. Daily secretion rates, which vary among certain hominins [40], are similar in fossil Pongo, H. erectus, and recent H. sapiens (cuspal average = 4.11 μm/day: [56]). Long-period line periodicity may be more distinct, as fossil Pongo shows the highest average periodicity value of any primate, although a single G. blacki specimen has a reported periodicity of 11 days [74]. The average periodicity in our H. erectus comparative sample is similar to early Homo specimens from South Africa and also Neanderthals, which fall between 7-8 days [69,72]. While expanded samples may lead to overlapping ranges for these fossil groups, periodicity values appear to be useful for distinguishing certain isolated Pongo and Homo teeth.
Crown formation times also vary between hominins and other hominoids (reviewed in [54,69]). Fossil orangutans show the longest molar crown formation times, along with G. blacki, due in part to their elevated long-period line periodicities and large crown size [54,74]. Limited data are available on early members of the genus Homo; post-canine formation times typically fall in the lower end of modern human ranges ( Table 6). The sole exception to this is CA 770, which has a mesiobuccal cusp formation time that exceeds the known range of upper first molar values in modern humans.

Presence and distribution of Asian hominins
Considerable efforts over the past 125 years have been dedicated to identifying hominins in Asian fossil assemblages, and this continues to the present day [16,18,30,33,41,[76][77][78][79][80][81][82]. It is now well established that H. erectus reached both mainland Asia and island Southeast Asia during the Early Pleistocene [83,84]. The earliest recovered teeth of H. erectus (originally termed Pithecanthropus erectus and Sinanthropus pekinensis) from Trinil, Java and Zhoukoudian, China have been thoroughly discussed and debated (e.g., [85][86][87]). However, we concur with Ciochon [18] that the frequency of occurrence of H. erectus in Asia may have been overestimated due to misidentification of isolated pongine dental remains.
We have shown here and in [30] that consideration of internal structural and developmental characters reveal subtle differences between H. erectus and fossil Pongo, providing additional evidence that may be used in concert with examinations of molar crown size and Tooth: types as in Table 1; Cusp: b-buccal; mb-mesiobuccal, ml-mesiolingual; Early Homo data from South Africa (SA): [69]; H. erectus: [71] and this study (M3549 LP3); recent H. sapiens: [69,72]. https://doi.org/10.1371/journal.pone.0204737.t006 Disentangling isolated dental remains of Asian Pleistocene hominins and pongines external shape (also see [33,41]). This approach could aid in the reassessment of dental remains from Tham Khuyen Cave, Vietnam [18,76,88] and Longgupo, China [4,12,89,90], which were originally thought to belong to early Homo but have been recently reclassified as pongines by Ciochon [18]. Similarly, we suggest that von Koenigswald's Chinese Apothecary collection includes more pongine remains than has been previously appreciated, as information on long-period line periodicities and crown formation times is broadly consistent with results from EDJ shape and enamel thickness analyses. This is particularly apparent when the "Hemanthropus peii" sample is compared with fossil Pongo. Wolpoff ([9]: p. 507) described this sample as "worn versions of the thick-enameled dental remains that have been attributed to "giant orangs."" Early Pleistocene fossil orangutans from China are especially large in comparison to modern-day orangutans, while Late Pleistocene forms are smaller [25]. Enamel thickness values and EDJ shapes that are nearly identical to Pongo call into question the proposition that the small "Hemanthropus peii" molar teeth belong to new Pleistocene hominoid (contra Ciochon [12,18]). Moreover, the "Hemanthropus peii" type specimen shows a long-period periodicity value that is consistent with our fossil Pongo sample. While it is possible that additional pongine genera may be included in the Chinese Apothecary collection, or in other sites from mainland Asia [91], we concur with Harrison et al. [25] that the "Hemanthropus peii" teeth are most likely fossil orangutans.
The identification of pongine material from Sangiran is somewhat less secure. Grine and Franzen [5] were unable to agree on the presence of fossil orangutans in the surface finds made by G.H.R. von Koenigswald. Our analysis suggests that two of the teeth originally attributed to Pongo by Grine are more likely to be hominins (S7-20; S7-65), while three of the teeth originally attributed to H. erectus (S7-53; S7-64; S7-76) are more similar to orangutans, particularly in terms of their EDJ shape. This conclusion could be strengthened with an expansion of our H. erectus molar comparative sample, as the current one is limited to maxillary molars. Mandibular molar EDJ shapes of extant humans differ from those of extant and fossil orangutans, but it is possible that H. erectus mandibular molar EDJ shapes were less distinctive than those of extant humans. Faunal assemblages from this region have been suggested to indicate an open woodland environment that predates the arrival of orangutans on Java [92,93]. Yet Drawhorn [94] notes that teeth in the Sangiran surface collection are unlikely to derive from a single biological population, reaching the same conclusion as Grine that some of the teeth derive from fossil orangutans. A recent analysis of enamel thickness and EDJ shape in a deciduous molar from Sangiran concluded that it belongs to Pongo [41], lending further support to the idea that von Koenigswald's Javanese collection contains both hominin and pongine material.
While structural and developmental integration may be useful for generic distinction of primate teeth, it remains to be seen whether these characters will aide in the distinction of Asian hominins, as is the case for comparisons of modern humans and Neanderthals (e.g., [65,67,72]). Smith et al. [67] documented considerable variation in enamel thickness within members of the genus Homo, with Neanderthals being an especially distinct group. Similarly, Neanderthals show lower long-period line periodicities and faster extension rates than recent humans, leading to shorter crown formation times [72]. Developmental analyses of Asian hominin teeth are limited to studies of the H. erectus Sangiran maxillary fragment S7-37 [71] and the Trinil molars [30], as well as an enigmatic child's maxilla from Xujiayao, China [66,79,95]. Fossils from this latter individual are described as possessing a combination of structural traits that are more similar to early Pleistocene Homo and Neanderthal fossils than to modern humans, yet the distribution and number of long-period growth lines are most similar to recent humans. It is possible that at least two and as many as four Homo species overlapped in Southeast Asia during the Pleistocene. Multiple lines of evidence place H. sapiens in China before 70,000 years ago [96,97], which may be concurrent with the enigmatic hominin material from Xujiayao [66]. Homo sapiens is now known to have crossed island Southeast Asia 63-73,000 years ago, reaching Australia by at least 65,000 years ago [81,98]. Moreover, a hominin toe bone from the Philippines at 67,000 years ago has been provisionally attributed to a small-bodied H. sapiens individual [99]. These groups may have encountered other hominins such as Homo floresiensis, which survived on Flores until 50,000-60,000 years ago [100]. Even if H. erectus, H. floresiensis, H. sapiens or other hominins did not occupy the same areas at the same time, uncertainties in dating, geological and taphonomic processes, and the provenience of mixed material such as von Koenigswald's Chinese Apothecary collection underscore the need to reassess the identification of other isolated, large-bodied, primate dental remains from the Pleistocene of Asia.
Supporting information S1 File. Contains supporting tables, figures, and references. Table A. Asian Pleistocene sites with isolated dental remains reported to belong to Hominins. Fig A. Chinese Apothecary dental material attributed to "Hemanthropus peii" by von Koenigswald. Fig B. Chinese Apothecary dental material attributed to Sinanthropus officinalis (= Homo erectus) by von Koenigswald.