Main

Humans exhibit several aspects of life history different from great apes, such as helpless newborns, early weaning and allomaternal infant care, a prolonged postnatal brain growth period, delayed somatic maturation, late onset of reproduction, transgenerational co-operation and a long post-reproductive life1. In particular, humans have an extended period of immaturity called childhood—the time after weaning and before the attainment of nutritional independence—that is characterized by (allo)parental care and ample opportunities for social learning8. These traits do not fossilize, however, so inferences on the evolution of human life history must be drawn from developmental proxies preserved in hominin fossils4,5,9,10. Teeth are of special relevance because their crowns and roots grow by the daily incremental deposition of enamel and dentin layers, respectively, thus conserving detailed information about the maturation schedule of an individual’s masticatory system5,11. Humans, compared to great apes, show a marked slow-down in dental maturation, especially of the permanent molars, which correlates well with the pace of cerebral development and somatic maturation across primates3,5. Dental incremental growth data from fossil hominins have yielded reliable estimates of the age at death of immature specimens12, permitting researchers to use dental growth chronologies to calibrate cerebral, cranial and somatic development. For example, microstructural dental evidence from the immature Australopithecus afarensis specimen from Dikika, Ethiopia13, obtained using propagation phase-contrast synchrotron microtomography (PPC-SRµCT), suggests that this individual essentially followed an ape-like dental ontogenetic schedule. However, its endocranial volume (ECV) at the time of death (2.4 years) had reached a smaller proportion of adult ECV than in great apes, which has been interpreted as the first evidence of prolonged postnatal brain growth. There is still relatively little information available on the exact timing of dental ontogeny in early Homo. The Nariokotome individual (KNM-WT 15 000), dated to approximately 1.6–1.5 million years ago (Ma), is thought to have followed a great-ape-like rather than human-like dental ontogenetic schedule14,15,16,17, but these inferences are based on estimates rather than direct measurements of its dental microstructural growth. When exactly the characteristic life-history slow-down occurred in the evolution of our genus Homo is therefore still an open question.

Here we present a comprehensive and non-destructive microstructural analysis of dental development in an early Homo individual from Dmanisi and draw inferences on the life-history evolution of early Homo. The Dmanisi site, dated to 1.85–1.77 Ma (ref. 18), has yielded five well-preserved crania, associated mandibles and various postcranial elements19,20,21,22,23. Together with a rich record of faunal remains and lithic artefacts, the Dmanisi ensemble has provided unique insights into the palaeobiology of an early Homo population outside Africa, its life history, social organization and cognitive capabilities24,25. Here we focus on a subadult individual that is represented by cranium D2700, mandible D2735 and various postcranial elements20,22.

The individual’s permanent dentition is well preserved, and each tooth type is represented by at least one antimere, except for the lower first incisors (Fig. 1). The first molar teeth (M1s) exhibit extensive occlusal wear, resulting in cuspal dentin exposure. The second molars (M2s) only exhibit moderate wear. The lower right third molar (LRM3) was found in isolation, but fits well in its alveolar cavity. The corresponding alveolus for the lower left M3 (LLM3) is not present, indicating LLM3 agenesis23,26. The upper right M3 (URM3) exhibits a supernumerary cusp (Extended Data Fig. 1) and was excluded from further microstructural analyses to ensure data consistency, its antimere having a normal morphology. All M3s exhibit well-developed, but not fully formed, roots, while already exhibiting some wear facets. The individual therefore died before reaching dental maturity, which allows an accurate determination of the individual’s age at death through quantification of the incremental growth structures preserved in the tooth enamel and dentin.

Fig. 1: Dentition of Dmanisi skull D2700/D2735.
figure 1

a–c, Computed tomography (CT)-based visualization of the skull and associated teeth showing right lateral (a), frontal (b) and left lateral (c) radiographic projections. d,f, Occlusal views of reconstructed mandibular (d) and maxillary (f) dentition (green, high-resolution PPC-SRµCT of isolated teeth used for microstructural analyses; orange, medium-resolution PPC-SRµCT of teeth in place, not included in microstructural analyses; red, low-resolution scans (medical CT, not included in analyses); light blue, missing teeth mirrored from preserved antimeres; dark blue, missing mandibular first incisors, reconstructed from second incisors). e, Frontal view of the skull. Scale bar, 5 cm.

Full size image

We used PPC-SRµCT virtual histology6,27 to track the development of the individual’s best-preserved maxillary and mandibular teeth (Fig. 1, Methods, Supplementary Data 1 and Extended Data Fig. 2). To obtain a complete record of dental development, the worn crown cusps were reconstructed using comparative data from better-preserved teeth (Extended Data Fig. 3; Methods). The cusps of the M1 typically comprise the neonatal line (a perinatal stress marker) in their apical-most region. Because this information was lost through dental wear, all following analyses are based on the assumption that the timing of M1 crown initiation was similar to that in living humans and great apes, that is, shortly before, to around, birth (Methods). Following this preparatory work, we reconstructed the developmental profiles of all teeth, as follows (Methods). (1) We identified major stress markers (indicating discrete systemic stress events during development), which are visible as hypoplastic discontinuities in the dental growth line sequences (Fig. 2a and Extended Data Figs. 4–6). Stress lines were matched across all teeth, resulting in a relative chronology of dental development (Fig. 2b). (2) High-resolution (0.72 μm voxel size) virtual cross-sections through the best-preserved enamel and dentin regions revealed a long-period (Retzius/Andresen) growth line timing of 6 days (which is at the fast end of variation in fossil hominins, living humans and great apes28) (Supplementary Data 2), and an average dentin secretion rate of 4.28 μm per day (Supplementary Data 3). (3) We counted the well-preserved long-period growth lines of the dentin (Andresen lines), and combined this information with the long-period line spacing and the dentin secretion rate data to establish the absolute timing of the major stress events, and of root formation from initiation to completion. (4) These data were combined with stress marker information from the enamel and calibrated with the Andresen line counts in the dentin to reconstruct crown formation from initiation to completion (Extended Data Figs. 4–6). (5) Finally, the crown and root formation schedules of all the teeth were aggregated into a composite record of dental ontogeny from birth to death (Fig. 2b).

Fig. 2: Dental development of the D2700/D2735 individual.
figure 2

a, PPC-SRµCT cross-sectional images through the developmental planes of all analysed teeth (for abbreviations, see text). Coloured lines indicate the six major stress markers (hypoplasias) used for cross-matching the teeth. b, Composite chart showing dental development. Tooth abbreviations as in a; individual cusps: mb, mesiobuccal; ml, mesiolingual; db, distobuccal; dl, distolingual. Grey bars, crown formation periods from initiation to termination; orange bars, root formation periods (light bands mark maximum root extension rate as a proxy for dental eruption). The vertical coloured lines correspond to the stress events shown in a. Scale bar, 1 cm.

Full size image

Our data indicate that individual D2700/D2735 died at an age of 11.36 ± 0.65 years (the error margin results from interdental variation in the distances between the stress markers, variation in the reconstruction of worn cusps, variation in the assumed M1 crown initiation and variation in the daily dentin secretion rate) (Supplementary Data 1–3). The time to root apex closure of the third molars was estimated to be 1.3–1.5 years (Methods). Thus, D2700/D2735 died shortly before reaching full dental maturity (12.0–13.5 years).

The almost complete incremental growth record of this individual’s permanent teeth allowed us to reconstruct the longitudinal ontogenetic trajectory of its dentition from birth to dental maturity (Figs. 3 and 4 and Supplementary Videos 1 and 2). We recorded the ontogeny of D2700/D2735’s dentition at 6-month intervals from birth to death and compared this with corresponding data from living humans, great apes and fossil hominins (Methods and Supplementary Data 4). To characterize taxon-specific modes of dental ontogeny, we distinguished between the pattern and rate of permanent dentition formation. The pattern describes how the eight tooth types (I1, I2, C, P3, P4, M1, M2, M3) matured relative to each other (Fig. 3), while the rate describes how the dentition, as a whole, matures with increasing age (Fig. 4). In quantitative terms, the pattern is given by the dental maturation scores (DMSs) of all tooth types of an individual’s dentition (with the tooth types indicated by the subscript labels DMSI1, DMSI2, DMSC, DMSP3, DMSP4, DMSM1, DMSM2 and DMSM3), while the rate is given by the DMS of the entire dentition as a function of individual age (Methods).

Fig. 3: Comparative analysis of patterns of dental ontogeny.
figure 3

a, Dental ontogenetic schedule by tooth type (Dmanisi, red; humans, green; chimpanzees (captive individuals), blue; dark/light hues indicate crown/root formation periods; small squares indicate estimated eruption time). Error bars (±1 s.d.) drawn at crown initiation and completion, tooth emergence and root completion (data represent n = 25 independent human population-level studies and n = 10 independent studies on chimpanzees; see source data for details). The symbols below the abscissa indicate estimated mean ages at M1 emergence in other great apes (captive gorillas, orange triangle; captive orangutans, orange square) and in fossil hominins (A. afarensis, circle; A. africanus, + sign; P. robustus, Y sign; early Homo, × sign). b, PCA of dental maturation patterns (DMS-PCA) (Dmanisi, red trajectory; humans, green dots and trajectory; captive known-age chimpanzees, blue dots and trajectory; numbers on trajectories indicate age in years; wild-living chimpanzees, blue circles; wild-living bonobos, blue diamonds; wild/captive gorillas, open/filled orange triangles; wild/captive orangutans, open/filled orange squares; A. afarensis, black circles (Dikika, AL333-105); A. africanus, + signs (Sts24, Taung); P. robustus, Y signs (DNH107, SK63); early Homo, × signs (KNM-ER1590, ZKD B-I4, KNM-WT 15000)). Each data point summarizes the maturation states of teeth I1–M3 of a given dentition. Insets show maturation patterns at selected locations in PC space. At comparable ages (6 years), the posterior dentition in humans is delayed compared to chimpanzees. At similar stages of advancement along the ontogenetic trajectory (chimpanzees: 6 years; humans: 9 years), the anterior dentition of chimpanzees appears relatively delayed compared to that of humans48. For data sources, see Supplementary Data 4.

Source Data

Full size image
Fig. 4: Comparative analysis of rates of dental ontogeny.
figure 4

a, Total DMS versus chronological age. Dmanisi data (red) per year. Chimpanzee data (blue): mean ± 1 s.d. of 1-year bins (n ≥ 10 individuals per bin; Supplementary Data 4). Human data (green): mean ± 1 s.d. of population-level 1-year bins (n = 13 independent population samples; Supplementary Data 4). Orange symbols indicate known-age individuals of Gorilla gorilla (triangles) and Pongo spp. (squares). Black/grey symbols indicate fossil hominin specimens of known/estimated individual age. Circles, A. afarensis (Dikika 2.4 years and AL333-105 approximately 2.4 years); + signs, A. africanus (Taung approximately 3.5 years and Sts24 4.35 years); Y signs, P. robustus (SK63 approximately 3 years and DNH107 4.8 years); × signs, early Homo (KNM-WT 15000, approximately 8 years). b, Inferred instantaneous ontogenetic rate of the dentition (change in total DMS per year) as a function of chronological age. Note early, middle and late timing of DGS (arrows) in chimpanzees (early infancy), Dmanisi (approximately 5 years) and humans (7 years; interpopulation range: 6.0–8.5 years).

Source Data

Full size image

Figure 3a summarizes the timing of the key dental ontogeny events in humans, chimpanzees and Dmanisi. The intraspecific variation in crown initiation and completion, tooth emergence and root completion is substantial29,30,31,32,33,34. However, a principal component analysis (PCA) of the DMSI1…M3 data (DMS-PCA (ref. 35)) (Methods) showed that humans have a distinctly different pattern of dental ontogeny compared to all great apes (chimpanzees, bonobos, gorillas and orangutans) (Fig. 3b). Notably, the Dmanisi pattern is within the human range of ontogenetic variation and outside the great ape range of variation (Fig. 3b).

Figure 4a shows taxon-specific maturation trajectories of the entire dentition. Chimpanzee dentition matures rapidly during the first few years of postnatal life, reaching full maturity at around 11–13 years of age36. Data from gorilla and orangutan specimens of known age suggest trajectories in the upper and lower ranges of the chimpanzee variation, respectively37,38. By contrast, human dentition is characterized by a slow early postnatal maturation, reaching full maturity around 18–22 years of age. The dentition of the Dmanisi individual initially followed a slow human-like rate trajectory until approximately 4 years of age, then catching up to reach a chimpanzee-like trajectory around 8 years, and finally reaching dental maturity at 12.0–13.5 years.

Figure 4b shows the inferred instantaneous maturation rate (total DMS per year) along the trajectories. In the chimpanzee dentition, the rate is highest in the first 2 years of postnatal life and then declines towards adulthood. By contrast, in living humans, the rate is comparatively low in the early postnatal period, but then increases to a peak at around 7.0 years of age, after which it declines steadily towards maturity. We refer to the peak rate as the dentition growth spurt (DGS). Notably, the Dmanisi individual also exhibits a DGS, but this occurred at approximately 5.3 years of age, which is significantly later than in chimpanzees (shortly after birth) and close to the lower range of the modern human DGS (in the range 6.0–8.5 years) (Fig. 4b and Extended Data Table 1).

Discussion

Dmanisi’s dental ontogeny appears as a combination of derived human-like and primitive great-ape-like features, on both relative and absolute time scales. Its pattern of ontogeny is within the range of human variation (Fig. 3b), with a delayed initiation of the posterior teeth (P3, P4, M2, M3, except M1) compared to the anterior teeth (I1, I2, C), and a wide spacing between the M1, M2 and M3 initiation times compared to great apes (Fig. 3a). As an effect, Dmanisi’s dentition, as a whole, matured slowly during the first approximately 5 years of life and exhibited a late DGS (Fig. 4). However, Dmanisi’s tooth-specific maturation rates differed from those of humans in several ways. Once the crowns of the premolars, M2 and M3, were initiated, their formation rates were exceptionally high, resulting in Dmanisi catching up to a chimpanzee-like maturation trajectory (Fig. 4a), and the estimated eruption times of the molars (M1: around 3.8 years; M2: around 7.6 years; M3: around 10.6 years) being more similar to those observed in wild chimpanzees31,39 than in humans40. Also, the predicted age of dental maturity (around 12.0–13.5 years) is more similar to chimpanzees and other great apes (11–13 years)36 than to modern humans (18–22 years).

How did Dmanisi’s mode of dental ontogeny evolve? Our data suggest a diversity of dental ontogenetic patterns and rates among both fossil hominin and extant great ape taxa. Great apes show a common pattern of ontogeny (Fig. 3b), but taxon-specific differences in rate38,41 (Fig. 4a). The data for A. afarensis (specimens AL 333-105, Dikika; geological age around 3.3 Ma) and Australopithecus africanus (Taung, Sts 24; around 2.8 Ma) suggest that the dental ontogenies of these species were probably within the range of ontogenetic pattern and rate variation seen in extant great apes, whereas Paranthropus (DNH 107, SK 63; around 2.0–1.8 Ma) demonstrates a mode of dental ontogeny not found in extant taxa41,42 (Figs. 3b and 4a). Fossils representing African early Homo (KNM-WT 15 000 Nariokotome, KNM-ER 1590; 1.8–1.6 Ma) conform with the Dmanisi/human pattern of ontogeny (Figs. 3b and 4a), supporting the notion that this pattern is an evolutionary hallmark of the genus Homo.

What are the implications for the evolution of hominin life history? Because teeth are central to masticatory function, evolutionary changes in their maturation schedules tend to reflect changes in masticatory biomechanics43 and in the timing of dietary transitions during an individual’s life, from breastfeeding to solid foods and from nutritional dependence to independence. These changes, in turn, are thought to reflect evolutionary changes in the allocation of metabolic resources to cerebral, somatic and reproductive maturation, and thus changes in life history.

In Dmanisi, compared to great apes, the characteristically delayed maturation of the posterior permanent dentition (note the relatively late eruption of P3 and P4 in Fig. 3a) and the late DGS imply a prolonged use of the deciduous molars. This probably indicates an extended period of nutritional dependence.

The long childhood period in humans8 is characterized by delayed dental and somatic maturation relative to brain growth and maturation3,5,10,44. This pattern is typically thought to have evolved because the large human brain requires considerably more metabolic energy to grow to adult size than in great apes2,8,44. The resulting delay in somatic maturation, in turn, would have provided new opportunities for immature individuals to acquire advanced cognitive skills in a complex social environment while the brain was still growing and maturing8.

However, this evolutionary scenario is not applicable to early Homo, as the advanced dental and skeletal maturity of Dmanisi and Nariokotome at their respective ages of death are more similar to those of great apes than living humans22,45,46,47 (Extended Data Table 2). Furthermore, the brain size of early Homo at Dmanisi (ECV range: 546–730 cm3) was only moderately larger than that of great apes, and therefore probably had a moderate impact on the metabolic trade-off between brain and body growth.

We therefore propose an alternative scenario for early Homo that assumes a diversity of life-history strategies among Plio-Pleistocene hominins, as reflected by the diversity of dental ontogenies42. Australopithecus probably exhibited a primitive ape-like mode of dental ontogeny48 (Figs. 3 and 4). However, immature Australopithecus individuals show evidence of a prolonged phase of postnatal brain growth and maturation compared to great apes13,49, suggesting that this pattern emerged before a marked evolutionary expansion of the brain.

In early Homo, according to this scenario, the extended growth phase did not evolve as a consequence of brain expansion and the associated trade-offs in resource allocation8,10. Rather, the selective pressures underlying the observed changes must be sought elsewhere. The palaeoanthropological and archaeological context of Dmanisi provides valuable information in this regard. The presence of Mode I stone tools and cut marks on animal bones indicates tool-assisted meat consumption at Dmanisi18,50. Therefore, it remains to be investigated whether the evolution of the early Homo pattern of dental ontogeny reflects a shift towards tool-assisted meat consumption, implying a reduced overall masticatory load51. Furthermore, the sociocultural context and its potential effects on life history and reproduction must be considered8,52. The four mature individuals found at Dmanisi19,21,22,23 span an age range from early to late adulthood, including one senile individual who had survived several years without teeth21,24. The possible presence of post-reproductive individuals in the Dmanisi palaeopopulation24 can be viewed in the context of biocultural reproduction7, which implies that intergenerational co-operation facilitated an extended period of growth sustained by alloparental provisioning52, and ultimately contributed to increased reproductive success.

This raises the question of whether the inferred life-history characteristics of early Homo reflect an evolutionary change in reproduction strategies7,8—before substantial brain expansion and reorganization25, and before a general slow-down in life history. To test these hypotheses, further dental microstructural analyses of key early Homo fossils, such as the Nariokotome individual, and of small-brained fossil hominin taxa, such as Australopithecus sediba53 and Homo naledi54, are required. Overall, our analyses add to the growing body of evidence that the evolution of hominin life histories did not follow a single linear trajectory, but was a branched and modular process, involving diverse selective pressures and adaptive responses.

Methods

Propagation phase-contrast synchrotron microtomography

Cranium D2700 and associated mandible D2735 were scanned at the ID19 beamline of the European Synchrotron Radiation Facility, Grenoble, France. Volume scans were performed using PPC-SRµCT. For general dental morphological assessment, as well as for the counting of Andresen lines, two configurations were used in two successive series of experiments. The first used voxel sizes of 4.89 µm with a monochromatic beam set at 51 keV to take advantage of the gadolinium K-edge of the 5-μm-thick gadox scintillator. The second series took place a few years later, after substantial technical evolution of the ID19 beamline, and was performed using a 5.06-µm voxel size using a polychromatic beam whose spectrum was shaped with a tungsten filter to isolate a peak between the K-edge of the lutecium from the lutecium–aluminium–garnet–(cerium) scintillator (63.3 keV) and the K-edge of the tungsten filter (69.5 keV). This combination enabled much faster scans than with the monochromatic setup due to the higher flux and higher energy, and with perfect beam stability and improved coherence, resulting in much better phase-contrast imaging. In both cases, the propagation distance was 5,000 mm and the scans were performed in half-acquisition to almost double the field of view of the detector. In addition, scans for daily periodicity assessment were performed at 0.75 μm at 65 KeV, also in half-acquisition. All scanning parameters are listed in Supplementary Data 1, and a comparison of the monochromatic and polychromatic results is presented in Extended Data Fig. 7.

Acquisition of dental developmental data

Determination of developmental planes

The dental development was tracked based on the PPC-SRµCT data of at least one tooth of each type (Supplementary Data 1). Using the software VGStudio Max v.3.4 and later (Volume Graphics GmbH), virtual cross-sections through the teeth were generated to optimally visualize the microstructural growth lines in the enamel (Retzius lines) and dentin (Andresen lines). The optimum cross-sectional plane, named the developmental plane55, was defined as the plane passing through a given dentin horn and the top of the underlying pulp cavity. The plane’s x axis was oriented in the bucco-lingual direction and its y axis along the tooth’s longitudinal axis. The plane’s thickness (z-axis depth) was set to values between 100 and 150 µm. The use of these ‘thick sections’ enhanced the visibility of the developmental microstructures in the tooth6.

Reconstruction of worn cusp tips

Several teeth of the Dmanisi individual had undergone substantial wear since eruption, as evidenced by missing enamel cusp tips and dentin horn tips (Extended Data Fig. 3). Because the most relevant microstructural developmental information in the teeth of this specimen came from the dentin (see ‘Daily secretion rate’), we focused on reconstructing the relatively small missing apical volumes of the dentin horns (Extended Data Fig. 3), using better-preserved teeth as guides. This process was repeated in several independent sessions. The reconstruction uncertainty inherent in this inference procedure is reflected in the calculations of the dental developmental timing and the age at death (Supplementary Data 1).

Matching of stress markers across the entire dentition

The general sequence of tooth development is defined by the order in which teeth initiate, emerge and complete their growth. Episodes of major stress during an individual’s development tend to affect the dentition as a whole, resulting in growth discontinuities, which become manifest as stress lines on the outer crown surface and in the microscopic record of enamel and dentin growth lines. A given stress event leaves its mark on all teeth that are still in the enamel and/or dentin deposition phase at this moment. This property is an important prerequisite for placing the development of the individual teeth in a common time frame (Extended Data Fig. 2). In the Dmanisi dentition, six major stress events could be identified (numbered 1 to 6 in Supplementary Data 1 and assigned different colours in Fig. 2). The matching of these stress lines across all teeth was ensured by cross-correlation between the Retzius and Andresen lines, and was further verified by comparison with the three-dimensional topology of the enamelo-dentin junction as well as the hypoplasias of the root surface (Extended Data Fig. 4) and lateral enamel surface (Extended Data Figs. 5 and 6), following the methodology proposed in ref. 56.

Long-period duration

The deposition of dental enamel and dentin exhibit both short-term (daily) and long-term (several days) periodicities. In virtual PPC-SRµCT cross-sections, the daily incremental growth structures are visible as cross-striations (or laminations) in the enamel (with a size of about 1–3 μm)6. Typically, the cross-striations are most visible in the relatively flat areas of the mid-lateral cuspal enamel. In enamel, the long-period structures are visible as Retzius lines (size of about 30–40 μm, measured along the enamel prisms). The number of daily increments between consecutive Retzius lines gives the value of the long-period line periodicity. This is a property that varies between individuals of a given species (for example, in humans, between 5 and 12 days28,57), but can also sometimes vary within a given individual, despite general consensus to the contrary. For example, in a modern human sample of n = 40, 18 individuals (45%) showed variation in long periodicity, with longer duration in the anterior than the posterior dentition57.

To determine the actual long-periodicity of the D2700/D2735 dentition, scans of the enamel at a voxel resolution of 0.72 µm were performed on the URP4, URC, ULM3 and LRM3, and the cross-striations were directly traced and counted between the long-period lines (Supplementary Data 1 and 2 and Extended Data Fig. 2). Each tooth was evaluated by two independent observers, and 178 independent measurements were taken all over the teeth. Owing to strong artefacts in the data and the non-perfect preservation of the specimens, periodicities of between 3 and 9 days were observed. The 6-day periodicity represented 41% of the measurements, the 5-day periodicity 31% and the other periodicities 28%. Given the non-normal distribution of these values (Supplementary Data 2) (Kolmogorov–Smirnov test: F = 0.2465, P < 9.1693 × 10−10), a Wilcoxon test was performed for comparison with expected periodicities of 5, 6 and 7 days, respectively. The only match between the expected and actual data was observed for the 6-day periodicity (P < 0.23), while the results for 5 (8.51 × 10−9) and 7 days (7.31 ×  10−10) suggest significant differences between the expected and observed distributions. Our data show no evidence of an antero-posterior gradient in the long-period duration57. We thus used the 6-day periodicity to reconstruct the developmental profiles of all teeth by counting the Retzius and/or Andresen lines, and multiplying the total number of lines by the long-period duration.

Daily secretion rate

In the D2700/D2735 specimen, the long-period lines in the cuspal enamel (Retzius lines) were barely visible at the available voxel resolution (4.89 and 5.06 µm), which prevented direct measurements of the daily enamel secretion rate. However, the D2700/D2735 teeth exhibit exceptionally well-preserved dentin microstructures, which permitted an inference of the daily secretion rates by counting the Andresen lines in those areas where they were clearly visible. Specifically, the daily secretion rate of dentin was evaluated as follows. The length (L, in microns) of a given well-preserved primary dentin tubule was measured (secondary dentin tends to show compression of the Andresen line patterns), and the number, A, of the Andresen line intervals along the tubule was counted. The daily secretion rate (S, in microns per day) was then evaluated as S = L/(A × P), where P is the periodicity in days (Extended Data Fig. 2). The resulting value of 4.28 μm was high compared to humans and chimpanzees30,58, but was cross-validated by two observers, and consistent results were obtained when the different sources of information were combined. Finally, this information was combined with the Andresen line counts and the cross-matched stress markers—visible in both the enamel and the dentin—to reconstruct the absolute developmental chronology of the whole dentition (Supplementary Data 1).

Dental eruption time

The time of dental eruption was estimated from the maximum extension rate of the root, termed the root growth spurt32,58 (Extended Data Fig. 2). The method used to determine the root extension rates was adapted from previous studies58,59 and proceeded as follows. Along the cemento-dentin junction, a first 200-μm-diameter circle was placed at the cervix and a second was aligned at the junction between the cemento-dentin junction and the Andresen line, passing tangentially on the side of the first circle (Extended Data Fig. 2). This procedure was repeated all along the length of the root. The greatest distance along the external surface of the root between the two circles marked the maximum extension rate, that is, the root growth spurt.

In great apes, the root growth spurt coincides with the period of dental eruption, whereas in modern humans, the growth spurt precedes eruption32,58 (Supplementary Data 4). In the Dmanisi individual, therefore, the maximum rates of root growth should be regarded as minimum estimates of tooth eruption times. Furthermore, tooth eruption is a process that takes place over several months. We therefore represent the estimated eruption period as the time point of the maximum root extension rate ± 3 months (Fig. 2 and Extended Data Fig. 2). It should be noted that such a root growth spurt is visible in the third molars on the last third of the developing roots, while these teeth erupted quite some time before the death of the individual, as indicated by the light wear facets. This supports our basic assumption that the root growth spurt in this individual tends to coincide with the dental eruption.

Age at death and dental ontogenetic trajectory

Synopsis of dental development

The time-calibrated microstructural growth data of the tooth crowns and roots permitted reconstruction of the complete sequence of dental development throughout the life of D2700/D2735, from the initiation to completion of each tooth. The resulting diagram is visualized in Fig. 2 and the underlying data are presented in Supplementary Data 2. These data were used to determine the age at death of the individual, as well as its longitudinal dental ontogenetic trajectory.

Age at death

The age at death corresponds to the total time of dental development, from the crown initiation of the first molars (estimated from reconstructions of their worn cusp tips) to the last deposited dentin surface in the roots of the third molars. Because the M1 neonatal line was lost through dental wear, Dmanisi’s M1 crown initiation age had to be estimated from comparative data. In humans, the M1 initiation range is from approximately 0.15 years before to approximately 0.03 years after birth, that is, −0.06 ± 0.09 years (refs. 60,61). With the M1 initiation set to birth (0 years), two independent analyses, performed by P.T. and V.B., yielded an age at death of 11.42 ± 0.56 years (Supplementary Data 1). Adding to this number the M1 initiation range yielded the following upper and lower bounds for the age at death of Dmanisi: −0.15 + 11.42 − 0.56 = 10.71 and 0.03 + 11.42 + 0.56 = 12.01 years, which together yield 11.36 ± 0.65 years.

Estimating the age at dental maturity

Dental maturity is reached with the root closure of the third molars. The time to root closure in the D2700/D2735 M3s was estimated by comparison with the homologous complete M2 roots. The M3 roots were scaled down to match the size of the M2 roots. Following the M2 dentin tubules, the distance from the open end of the M3 root to the closed M2 root tip was measured. Scaling this distance up to the original M3 size, and assuming an average dentin secretion rate of 4.28 μm, gives 469–561 days (about 1.3–1.5 years) to M3 root completion. The lower and upper bounds for the age at dental maturity are 10.71 + 1.3 = 12.0 years and 12.01 + 1.5 = 13.5 years.

Visualization of the ontogenetic trajectory

The ontogenetic trajectory was visualized using a video capturing the reconstructed developmental stages of each tooth at 6-month intervals, from birth to death, in the form of developmental plane sections (Supplementary Video 1). A second video (Supplementary Video 2) shows the same data in the form of virtual pseudo-orthopantomograms, which can be compared, to some extent, with clinical radiographic orthopantomograms.

The videos were constructed from consecutive virtual snapshots taken at specific time points along the ontogenetic trajectory as follows. Each tooth’s developmental plane (as defined above) was used as a base image. On each image, we first annotated the stress lines (Fig. 2 and Extended Data Figs. 4–6), which served as absolute time marks. Using the dentin and enamel extension rates as relative measurements of time, we then interpolated the locations of the Andresen and Retzius lines every 6 months, from birth to death. Then, working backwards from the complete developmental plane, dental material was removed virtually, and the resulting images were saved as snapshots for each 6-month step. The videos were recomposed from the single-tooth image sequences as follows. For each series (upper and lower dentition), two reference lines were defined to determine the position of the teeth before and after emergence, respectively. Pre-emergence teeth were positioned with their distal-most part at the pre-emergence line, while post-emergence teeth were positioned with their proximal-most part (crown cusps) at the post-emergence line. Finally, the images of all the teeth for each time step were merged together, and the timed snapshots were combined to form a video. The effects of tooth wear were not simulated in these videos because it would have been too difficult to make accurate interpolations of this process, so the crowns are shown with reconstructed cusp tips (recognizable by them having higher density than the preserved enamel), and the actual tooth wear is only shown at the time of death.

Ontogenetic trajectory of the dentition

The long dental development record of D2700/D2735 up to the subadult stage (well-developed, but open, M3 roots) allowed an almost complete dental ontogenetic trajectory to be reconstructed. The recovery of an individual longitudinal trajectory of 11.4 years is unique for fossil hominins, and even for individuals of living species, where trajectories are typically obtained by combining data from ontogenetic series rather than longitudinal ontogenetic samples. A scoring system developed previously13 that makes it possible to compare virtual sections with radiographic surveys was used to determine the overall maturation state of the entire dentition at 6-month intervals, from birth to death. Specifically, the maturation state of each tooth (I1, …, M3) was scored at each time step, following the procedures described in ref. 13. The resulting DMSs (DMSI1, …, DMSM3) obtained from the virtual histology sections were then converted into equivalent clinical radiographic scores13 to permit direct comparisons with the radiographic data from the literature. Finally, at each time step, the DMSs of all the teeth were summed to obtain a total DMS (standardized to 100%). The resulting longitudinal ontogenetic trajectories of D2700/D2735 are visualized in Fig. 3 (DMSI1, …, DMSM3) and Fig. 4 (total DMS).

Comparative data

Comparative data on the developing dentition of known-age chimpanzees and humans were collected from the literature (Supplementary Data 4). The data for chimpanzees were from refs. 35,36,62 and comprise all eight permanent tooth types (I1–M3). The known-age human dataset was combined from studies examining all eight teeth (I1–M3) (refs. 40,63,64,65) or seven teeth (I1–M2) (refs. 63,65,66,67,68,69,70,71,72,73). Comparative data for Pan paniscus came from ref. 74. The data for Gorilla and Pongo were acquired from medical CT scans of collection specimens (Supplementary Data 4 and ref. 38). Data on the dental ontogeny of fossil hominin specimens representing Australopithecus, Paranthropus and early Homo were collected from the literature9,12,13,48,75,76,77. Comparative data on M1 eruption were from refs. 32,34,38,78.

The primary data used in our comparative analyses were the DMSs of each tooth type (DMSI1, …, DMSM3), and of the dentition as a whole (total DMS = sum of DMSI1, …, DMSM3). The included studies were based on different dental-imaging methods (dental radiography, CT, synchrotron-based tomography) and used different DMS systems, which required prior calibration and standardization13,79. As a common reference, we used the Demirjian DMS system66, which subdivides dental maturation into eight stages (four stages each for crown and root formation). The calibration schemes used to convert the 10-, 12- and 14-stage DMS systems to the 8-stage Demirjian system are listed in Supplementary Data 4. All calibrated data were standardized to total DMS = 100% for a fully formed dentition. Consistency checks to ensure the comparability of the data are described in detail further below. These checks show that, after calibration, the potential residual bias due to the use of different imaging methods and different scoring systems did not exceed the natural variation found in 1-year age bins of chimpanzees and/or humans.

Data analysis

We analysed the ontogeny of the dentition as a whole (hereafter referred to as the dental ontogeny) and focused on two aspects, the ontogenetic pattern and rate.

The ontogenetic pattern of an individual’s dentition was given by the set of DMSs of the eight tooth types (DMSI1, …, DMSM3). For a sample of n specimens, which resulted in an n × 8 DMS matrix. To analyse these multivariate data, we used the approach proposed in ref. 35. Using the software JMP v.15.2 (SAS Institute Inc.), a PCA was performed on the n × 8 DMS matrix, resulting in eight PCs that accounted for the largest to smallest proportions of the total variance in the sample. Typically, the first few components of a PCA comprise a large proportion of the total variance in the sample, which permits visualization of the relevant patterns of variation in the sample in a low-dimensional subspace.

The DMS-PCA method, as we call it here, made it possible to document variation along and across ontogenetic trajectories and to identify taxon-specific differences between the maturation states of the different tooth types relative to each other (Fig. 3b). Ontogenetic trajectories start with no teeth (0,0,0,0,0,0,0,0) and end with all eight fully formed teeth (1,1,1,1,1,1,1,1), so that, in multivariate space, all trajectories have the same starting point, then diverge, finally converging again to a single end point. Therefore, taxon discrimination is best in the middle of the trajectory.

In our analysis of Fig. 3b, PC1 accounted for variation along the ontogenetic trajectories and was essentially equivalent to the total DMS (PC1 = −16.105 + 0.284 × DMS; R2 = 0.9999; P < 0.0001), as already noted in ref. 35, with PC2 accounting for taxon-specific differences between trajectories. It is important to note that individual age was not part of the DMS-PCA, but served as an external variable to time-stamp taxon-specific trajectories through multivariate space. Therefore, specimens of unknown individual age, such as wild-lived great apes and many hominin fossils, can also be included in the analysis.

While the DMS-PCA method provided a detailed and comprehensive picture of dental ontogenetic patterns, we noted that alternative statistical methods have been used to reach similar conclusions regarding the differences between great ape, human and fossil hominin patterns of dental ontogeny48,80.

Ontogenetic rate

To quantify the rate at which the dentition of a given individual or taxon developed, the total DMS (sum of DMSI1, …, DMSM3) was plotted against age. The resulting total DMS trajectories are visualized in Fig. 4a. Furthermore, to summarize the differences between taxon-specific trajectories, we evaluated instantaneous ontogenetic rates = (DMS(ti) − DMS(ti−Δt))/Δt, with Δt = 0.5 years (Fig. 4b). Using MS Excel v.16.16.27 software (Microsoft Office), a moving-average function, comprising time steps ti−Δt, ti and ti+Δt, with Δt = 0.5 years, was then applied to level out, to some extent, the discontinuities resulting from the discrete nature of the scoring procedure (Supplementary Data 4). The instantaneous ontogenetic rates trajectory of chimpanzees was evaluated from the data in refs. 35,36,62. The human instantaneous ontogenetic rate trajectories were evaluated per population from the maturation score data for eight teeth40,63,64,65,81 or seven teeth63,65,67,68,69,70,72,73. The human DMS trajectories show a DGS, the timing of which varies between populations (6.3–8.5 years) (see main text and Extended Data Table 1). The upper range (8.5 years) largely coincides with M3 crown initiation, such that seven-teeth scoring schemes are unlikely to bias DGS estimates compared to eight-teeth scoring schemes.

Consistency checks

Influence of imaging methods and scoring schemes on dental maturation trajectories

To assess how different imaging methods and different scoring schemes affected the evaluation of the dental ontogenetic rate and pattern, we used data from two independent studies of dental ontogeny in Finnish children, here named Haavikko70 (ref. 40) and Nyström07 (ref. 82). The Haavikko70 data were based on a ten-stage DMS scheme, whereas the Nyström07 data used the eight-stage Demirjian scheme66. To assess the effects of imaging methodology (traditional dental radiography versus tomography-based virtual histology), we recoded the Haavikko70 data from the original radiographic to histological (tomographic) scores, using the conversion scheme presented in ref. 13. To assess the effects of data standardization, we recoded the Haavikko70 data with the ten-to-eight-stage recoding scheme presented in Supplementary Data 4. The analyses presented in Extended Data Fig. 8a,b thus compare four datasets: Nyström07 eight-stage, Haavikko70 ten-stage, Haavikko70 eight-stage radiography and Haavikko70 eight-stage virtual histology. We evaluated the DMS-versus-age profiles (Extended Data Fig. 8a), and used a DMS-PCA to visualize the ontogenetic trajectories through multivariate space (Extended Data Fig. 8b). Then results show that the different scoring schemes yielded different DMS-versus-age profiles, but that standardization led to good correspondence between the two datasets. As noted earlier12,13,83, the uncalibrated data from the virtual histology yielded total DMS-versus-age profiles that were consistently more advanced than the radiography-based profiles. Interestingly, however, in the DMS-PCA space, all four datasets (calibrated and uncalibrated) showed dental maturation trajectories that largely coincided, indicating that the DMS-PCA was robust to differences in imaging methods and scoring systems.

Dental maturation profiles and evaluation of the dentition growth spurt

To assess the influence of different DMS schemes on the evaluation of the DGS, we compared dental maturation data from three independent studies on Finnish children40,67,82. These studies used different scoring schemes (ten-stage versus eight-stage), and different tooth arrays (including/excluding M3) (Extended Data Fig. 8c and Supplementary Data 4). Nevertheless, the DMS rate peaks were at similar locations along the age axis (Extended Data Fig. 8d), indicating that the evaluation of DGS was robust to differences in DMS scoring schemes.

Assessment of age-related variation in dental maturation patterns along taxon-specific ontogenetic trajectories

DMS-PCA was performed on known-age humans and chimpanzees. As shown in Extended Data Fig. 8e, there is intra-taxon overlap between the 1 s.d.-density ellipses around consecutive 1-year age groups. The natural variation within 1-year age groups along the trajectories tended to be greater than the potential bias remaining after calibration for different data acquisition and scoring methods (Extended Data Fig. 8b).

Reporting summary

Further information on research design is available in the Nature Portfolio Reporting Summary linked to this article.