Anne-Marie Bacon
email: anne-marie.bacon@parisdescartes.fr

Pierre-Olivier Antoine
email: pierre-olivier.antoine@umontpellier.fr

Thi Mai

Huong Nguyen

Kira Westaway
email: kira.westaway@mq.edu.au

Jian-Xin Zhao
email: j.zhao@uq.edu.au

Anh Tuan Nguyen

Philippe Duringer
email: duringer@unistra.fr

Jean-Luc Ponche
email: ponche@unistra.fr

Dung Canh
email: samcanhdung@gmail.com

Huu Nghia Sam

Truong

Antoine Pierre-Olivier

Thi Mai Nguyen

Huong

Westaway Kira

Tuan

Duringer Philippe

Ponche Jean-Luc

Sam Canh

Huu Truong

Nghia

Tran Thi Minh
email: tranthiminh86@gmail.com

Thanh Pham

Son

Nguyen Thi

Kim Thuy

Demeter Fabrice
email: demeter@mnhn.fr

Linear enamel hypoplasia in large-bodied mammals of Pleistocene northern Vietnam, with a special focus on Pongo

Keywords: Orangutans, large-bodied Ungulates, Life-history, Weaning, Southeast Asia. 1. Introduction

ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

make them vulnerable throughout their period of development (Galdikas and Wood, 1990;[START_REF] Delgado | The Behavioral Ecology and Conservation of the Orangutan (Pongo pygmaeus): A Tale of Two Islands[END_REF][START_REF] Jablonski | The influence of life history and diet on the distribution of catarrhine primates during the Pleistocene in eastern Asia[END_REF][START_REF] Knott | Female reproductive ecology of the apes: implications for human evolution[END_REF][START_REF] Wich | Determinants of orangutan density in the dryland forests of the Leuser Ecosystem[END_REF][START_REF] Wich | Orangutan life history variation[END_REF][START_REF] Van Noordwijk | Development of ecological competence in Sumatran orangutans[END_REF].

Studies show that, among wild populations of great apes, there is a great interspecific variability, but it should be noted that the analyses are based on different sample sizes, and different methods of recording LEH-affected individuals [START_REF] Guatelli-Steinberg | Linear enamel hypoplasia as an indicator of physiological stress in great apes: Reviewing the evidence in light of enamel growth variation[END_REF]; Table 1). Orangutans (Pongo abelii and P. pygmaeus) have the highest prevalence of hypoplasia, with ~79 to 100 % of individuals [START_REF] Guatelli-Steinberg | Linear enamel hypoplasia in gibbons (Hylobates lar carpenteri)[END_REF][START_REF] Hannibal | Linear Enamel hypoplasia in the great apes: analysis by genus and locality[END_REF][START_REF] Guatelli-Steinberg | Linear enamel hypoplasia as an indicator of physiological stress in great apes: Reviewing the evidence in light of enamel growth variation[END_REF]. Proportions vary significantly between Gorilla species, with ~5.3 to 11 % (G. beringei), ~39 to 95 % (G. gorilla), ~75 % (G. grauei), and between Pan species but to a lesser extent, with ~69.7 -95.7 % (P. troglodytes) and ~95.7 -98 % (P. paniscus) (Skinner, 1986a;[START_REF] Eckhardt | Tooth crown development: non-human primate perspectives on the interpretation of linear enamel hypoplasia frequencies in present and past hominid populations[END_REF][START_REF] Stottlemire | Presence of enamel hypoplasia in wild-shot chimpanzees (genus Pan) and gorillas (genus Gorilla)[END_REF][START_REF] Guatelli-Steinberg | Linear enamel hypoplasia in gibbons (Hylobates lar carpenteri)[END_REF][START_REF] Hürner | Indicateurs de stress squelettiques et dentaires : qu'en est-il chez les grands singes africains ?[END_REF][START_REF] Hannibal | Linear Enamel hypoplasia in the great apes: analysis by genus and locality[END_REF][START_REF] Guatelli-Steinberg | Linear enamel hypoplasia as an indicator of physiological stress in great apes: Reviewing the evidence in light of enamel growth variation[END_REF]. All three taxa show the repetitive occurrence of LEH marks on incisors and canines, suggesting repeated stress events, but they differ in the mean number of recorded episodes, with ranges from 4.5 -7.4 in Pongo, 3.1 -4.7 in Pan, and 0.1 -5.4 in Gorilla [START_REF] Guatelli-Steinberg | Linear enamel hypoplasia as an indicator of physiological stress in great apes: Reviewing the evidence in light of enamel growth variation[END_REF]). An array of factors have been explored to attempt to explain the significance of these differences between great apes: biological factors such as enamel formation time, the level of abrasion of the surface of dental crowns, the angles between striae of Retzius and the enamel surface [START_REF] Hannibal | Linear Enamel hypoplasia in the great apes: analysis by genus and locality[END_REF][START_REF] Guatelli-Steinberg | Linear enamel hypoplasia as an indicator of physiological stress in great apes: Reviewing the evidence in light of enamel growth variation[END_REF], or ecological factors related to the environmental context such as habitat, climate and seasonality, or again related to food availability, particularly the abundance of fruit (versus leaves) in ape diets. Overall, studies indicate that linear enamel hypoplasia is most likely influenced by a complex combination of factors [START_REF] Dean | Retrieving chronological age from dental remains of early fossil hominins to reconstruct human growth in the past[END_REF], and that there is evidently a great variability in the way individuals feel stress.

In relation to wild populations of orangutans, it has been proposed that marked repetitive sequences of linear hypoplasia might be explained by a significant stressor, with an average periodicity of six months (or a multiple of six months), lasting 6 -8 weeks. The nature of this stressor, however, remains unknown, although it has been suggested that seasonal changes are deeply involved [START_REF] Skinner | Periodicity of repetitive linear enamel hypoplasia in Asian and African apes[END_REF][START_REF] Skinner | Hypothesis for the causes and periodicity of repetitive linear enamel hypoplasia in large, wild African (Pan troglodytes and Gorilla gorilla) and Asian (Pongo pygmaeus) Apes[END_REF][START_REF] Skinner | Orangutans, enamel defects, and developmental health: A comparison of Borneo and Sumatra[END_REF]. Temperature and humidity greatly influence the degree of virulence of parasites responsible for diseases and malnutrition in immature individuals, but some other factors may be important: the density of animals, which increases infection rates [START_REF] Foitová | Parasites and their impacts on orangutan health[END_REF], genetics [START_REF] Goodman | Assessment of systemic physiological perturbations from dental enamel hypoplasias and associated histological structures[END_REF][START_REF] Miles | Colyer's variations and diseases of the teeth of animals[END_REF], and annual or supra-annual variation in fruit cycling [START_REF] Knott | Changes in orangutan caloric intake, energy balance, and ketones in response to fluctuating fruit availability[END_REF].

In this paper, we present the results of a macroscopic analysis of enamel hypoplastic defects in three fossil collections of Pongo teeth: the late Middle Pleistocene Coc Muoi Pongo cf. devosi (148 -117 ka; Bacon et al., 2018a); the Late Pleistocene Duoi U'Oi P. pygmaeus (70 -60 ka; [START_REF] Bacon | The Late Pleistocene Duoi U'Oi cave in northern Vietnam: palaeontology, sedimentology, taphonomy, palaeoenvironments[END_REF][START_REF] Bacon | Late Pleistocene mammalian assemblages of Southeast Asia: new dating, mortality profiles and evolution of the predator-prey relationships in an environmental context[END_REF]Bacon et al., , 2018b)). The Middle Pleistocene Tham Khuyen series of isolated teeth have been allocated to various species, P. pygmaeus [START_REF] Cuong | Fossile Menschenfunde aus Nordvietnam[END_REF], P. pygmaeus and P. hooijeri [START_REF] Schwartz | A review of the Pleistocene hominoid fauna of the Socialist Republic of Vietnam[END_REF], or P. weidenreichi [START_REF] Harrison | Fossil Pongo from the Early Pleistocene Gigantopithecus fauna of Chongzuo, Guangxi, southern China[END_REF], the age of which is however not precisely known (>475 ka; [START_REF] Ciochon | Dated co-occurrence of Homo erectus and Gigantopithecus from Tham Khuyen Cave, Vietnam[END_REF]. Despite multiple biases due to the nature of the fossil collections, our analysis seeks to test whether the Pleistocene Pongo experienced a level of developmental stress comparable to that of infants and juveniles (weaned immature individuals) from extant orangutan populations. Finally, the analysis of hypoplasia in the Coc Muoi large-bodied mammals, allows us to highlight the health status and level of stress of a late Middle Pleistocene fauna in a tropical ecosystem.

Previous results on Coc Muoi and Duoi U'Oi large-bodied ungulates

2.1.

Dating and faunal context Both faunal assemblages, almost essentially composed of isolated teeth of medium-to large-sized mammals, were recovered from karstic breccias in cave sites in northern Vietnam (Figure 1SI). Their analysis revealed that they share comparable taphonomic pathways, due to the action of biotic (bone accumulator, either rodents or carnivores) and abiotic (water flows) agents, through processes of deposition [START_REF] Bacon | Late Pleistocene mammalian assemblages of Southeast Asia: new dating, mortality profiles and evolution of the predator-prey relationships in an environmental context[END_REF]2018a).

Duoi U'Oi yielded 871 teeth (NISP), representing a fully modern fauna [START_REF] Bacon | The Late Pleistocene Duoi U'Oi cave in northern Vietnam: palaeontology, sedimentology, taphonomy, palaeoenvironments[END_REF]. The estimated age of the fauna (70 -60 ka) is based on U-series dating of the calcitic floor (66 ± 3 ka) that crosses the fossiliferous breccia, and luminescence dating of the underlying sandy and silty clay unit [START_REF] Bacon | Late Pleistocene mammalian assemblages of Southeast Asia: new dating, mortality profiles and evolution of the predator-prey relationships in an environmental context[END_REF]. Rhinocerotid teeth have been attributed to three species, Rhinoceros sondaicus, R. unicornis and Dicerorhinus sumatrensis.

Coc Muoi yielded 1323 isolated teeth (NISP) (Bacon et al., 2018a). The fauna is a mix of modern and archaic taxa, characteristic of the mammalian communities of the Middle Pleistocene to the first half of the Late Pleistocene in the region. The dating of the cave sediments (OSL and pIR-IRSL) indicate a potential age range of 148 -117 ka. Rhinocerotid teeth have been attributed to the same three Asian species, whereas the large bovid teeth have been tentatively assigned to Bos cf. sauveli, the kouprey.

LEH data in rhinocerotids and bovids

In this first macroscopic analysis of ungulates, we described the frequency and location of linear marks of hypoplasia on deciduous and permanent teeth. In Coc Muoi rhinocerotids, 73.1 % of teeth showed no signs of LEH, whereas the percentage of Duoi U'Oi teeth was much higher (90.3 %), with a notable difference in the numbers of teeth between samples (264 versus 52, respectively) (Table 1). The analyses had highlighted the small percentage of milk molars (4/144, 2.7 %) affected by linear hypoplasia in the Coc Muoi sample with respect to permanent teeth, and only 2 out of 5 milk molars in the Duoi U'Oi assemblage.

Enamel hypoplasia in the Coc Muoi rhinocerotids is most expressed on the last premolars (P4 and p4; 11/30, 36.6 %) and molars (M3 and m3; 27/64, 42.1 %) (Bacon et al., 2018a), i.e. the last erupting teeth. These teeth recorded between one to four LEH episodes (vs one to two LEH episodes for Duoi U'Oi but with a smaller dataset) (Table 1).

In Coc Muoi bovids, 84 % of the teeth showed no signs of LEH (Table 1). No hypoplasia has been observed on the four milk teeth (4 out of 69) (Bacon et al., 2018a). Among the 65 permanent teeth, eleven molars recorded one (10/11, 14.4 %,) to two LEH (1/11, 1.4 %) episodes (Table 1), with proportions of 6.6 % (1/15) of m1, 31.2 % (5/16) of m2, and 5.5 % (5/9) of m3 (Bacon et al., 2018a).

Considering the chronology of tooth development and the location of the LEH defect on the surface of the crowns, it was possible to deduce within what age range the stress event occurred for each individual (Bacon et al., 2018a). In the studied ungulates, crown formation is a process starting in utero (deciduous dentition and the first permanent teeth begin to form) and continuing around and after birth (permanent dentition). In rhinoceroses, P4/p4 and M3/m3 are the last teeth to form, erupting at around ~7 -9 years old in Ceratotherium simum and Diceros bicornis [START_REF] Goddard | Age criteria and vital statistics of a black rhinoceros population[END_REF][START_REF] Hillman-Smith | Age estimation of the white rhinoceros (Ceratotherium simum)[END_REF], whereas in large bovids, all three premolars and M3/m3 are the last teeth to form, erupting after circa 2 years of age (Bison bison; [START_REF] Niven | An inter-site comparison of enamel hypoplasia in bison: implications for paleoecology and modeling Late Plains Archaic subsistence[END_REF]) (Table 2). Results highlighted perinatal events and the interval coinciding with abandonment of calf and sexual maturity in rhinocerotids, particularly at Coc Muoi (~4 -6 years; 10 out of 25 individuals), as physiologically stressful periods. Bovids were exposed to a relatively high degree of undernutrition and disease during weaning (~9 months), and when they reached sexual maturity (~2 years).

LEH data in recent and fossil orangutans

3.1.

Material We used a sample of seventeen adult individuals (twelve males, four females, and one of unknown sex), to assess the prevalence of hypoplasia and its range of expression on the permanent teeth of recent orangutans. Ten additional immature individuals (four males, one female, and five of unknown sex) with full deciduous dentition and erupted M1/m1 were also selected to evaluate the prevalence of defects on deciduous teeth and the first permanent molars. Within this overall sample, sixteen adult and eight immature individuals had been caught in the wild and spent the rest of their lives in the zoo, while three were born at the zoo (Ménagerie du Jardin des Plantes, Paris). Only three specimens come from the island of Borneo, and can therefore be attributed to Pongo pygmaeus [START_REF] Brandon-Jones | Asian primate classification[END_REF], but their populations of origin are unknown. The island origin of the other 21 specimens is unknown. Therefore, the recent Pongo collection used here does not reflect a natural population. The collection is housed in the Laboratoire Mammifères et Oiseaux, Anatomie comparée of the Muséum national d'Histoire naturelle, in Paris (Table 1SI).

Due to the poor preservation of some fossil teeth, damaged or heavily worn, we selected 20 out of 25 permanent teeth in the Coc Muoi assemblage, and 60 out of 66 permanent teeth in the Tham Khuyen assemblage. The Duoi U'Oi sample consists of 30 permanent teeth and 2 milk teeth (out of 34 teeth) (Table 3). The orangutans from Coc Muoi have been assigned to Pongo cf. devosi (Bacon et al., 2018a), those from Duoi U'Oi to P. pygmaeus [START_REF] Bacon | The Late Pleistocene Duoi U'Oi cave in northern Vietnam: palaeontology, sedimentology, taphonomy, palaeoenvironments[END_REF], whereas those from Tham Khuyen are now considered to belong to P. weidenreichi [START_REF] Harrison | Fossil Pongo from the Early Pleistocene Gigantopithecus fauna of Chongzuo, Guangxi, southern China[END_REF]. However, tooth sizes between the three fossil samples show significant overlap (Table 2SI). Moreover, as mentioned by [START_REF] Harrison | Fossil Pongo from the Early Pleistocene Gigantopithecus fauna of Chongzuo, Guangxi, southern China[END_REF], taxonomic assignment of species based on isolated teeth must be cautiously considered. The data were collected in Vietnam on the original specimens at the Lang Son museum, Lang Son province (Coc Muoi), and at the Institute of Archaeology in Hanoi (Tham Khuyen and Duoi U'Oi).

3.2.

Method Given the precious nature of fossils, histological analyses [START_REF] O'hara | Applying standard perikymata profiles to Pongo pygmaeus canines to estimate perikymata counts between linear enamel hypoplasias[END_REF][START_REF] Mcgrath | Quantifying linear eamel hypoplasia in Viruga Mountain gorillas and other great apes[END_REF] were excluded. Moreover, some teeth presented worn crown surfaces with smoothened incremental marks, preventing the counting of perikymata. Thus, to address this issue, we conducted a macroscopic analysis of enamel hypoplasia. Furthermore, we deliberately limited our data collection to linear enamel hypoplasia [LEH], and ignored all other defects randomly distributed on crown surfaces. Some teeth were quite difficult to analyze due to the cratering of the enamel, especially those from Tham Khuyen and Coc Muoi. The causes of these marks (accidents, feeding behavior, damages due to transport, deposition and burial in the karsts, and/or recovery) could not however be established.

Such study of isolated teeth inevitably introduces important limitations and constraints: the different compositions of the three assemblages regarding tooth types, the fact that some molars cannot be assigned to a specific position on dental rows, and the fact that a set of isolated teeth can comprise several teeth belonging to a single individual. Furthermore, the calculation of the minimum number of individuals [MNI], based on the most frequent type of tooth, reveals the small datasets: MNI = 2 at Coc Muoi (4 upper right M1/M2), MNI = 3 at Duoi U'Oi (3 lower right p3), and MNI = 6 at Tham Khuyen (6 upper right P3).

LEH expression has been recorded from the original specimens, either modern or fossil, in the museums, with the naked eye. The scoring process has then been repeated in the laboratory using the images of each specimen, under low magnification (x10) until replicability was attained. Data were collected on the labial side of teeth of modern and fossil specimens by a single observer (AMB) to minimize inter-observer variability. As Figure 1 shows, LEH has multiple expressions. However, we did not observe in both fossil and recent samples, pit-type defects scattered over the surface of tooth crowns, as described by Witzel et al., (2006;Fig.1C). The two criteria chosen to describe the morphology are those commonly used for mammals [START_REF] Hillson | The relationship of enamel hypoplasia to the pattern of tooth crown growth: a discussion[END_REF][START_REF] Hillson | Teeth[END_REF]: (1) a line of pits, thin or wide, visible or partly visible on the contour of the crown; (2) a furrow, thin or wide, visible or partly visible on the contour of the crown. The pit-form defects appear when clusters of ameloblasts stop forming enamel matrix, whereas ameloblasts in between the disrupted clusters continue to produce normal enamel matrix [START_REF] Hillson | The relationship of enamel hypoplasia to the pattern of tooth crown growth: a discussion[END_REF]. The line displays either small and shallow pits, or deep pits located in a marked groove (Figure 1A-B). In the furrow-form defects, ameloblasts ceased simultaneously matrix production, giving them a more regular aspect (Figure 1C-D). We considered only the lines and furrows that can be matched on teeth forming at the same time, and which have similar systemic causes [START_REF] Hillson | The relationship of enamel hypoplasia to the pattern of tooth crown growth: a discussion[END_REF] (Figure 2). The criteria used for the location of the defects with regards to crown height are those used in previous analyses (Bacon et al., 2018a): (1) basal third; (2) middle third, and (3) occlusal third, estimated between the cervix and the tip of cusps.

Table 3 shows a clear discrepancy in the number of canines between samples. This is why we have divided the data analysis into two parts: (1) the first analysis of hypoplasia on the crown surfaces of the incisors, premolars, and molars from the three sites (Coc Muoi, Duoi U'Oi, and Tham Khuyen), and (2) the second analysis of hypoplasia on the set of canines from Tham Khyuen.

In order to estimate the biological age range of fossil individuals on isolated teeth at the time of physiological defects, we used a reference Table based on our observations of LEH on the seventeen skulls of recent adult Pongo individuals (Table 4). Incisors, canines and premolars develop overall concurrently, with however a great variability in the sequence and the period of development, initiation of crown mineralization, crown completion time, root formation, and eruption time [START_REF] Swindler | Primate dentition: an introduction to the teeth of non-human primates[END_REF][START_REF] Dean | Developing pongid dentition and its use for ageing individual crania in comparative cross-sectional growth studies[END_REF][START_REF] Dean | Tooth crown heights, tooth wear, sexual dimorphism and jaw growth in hominoids[END_REF][START_REF] Beynon | Histological study on the chronology of the developing dentition in gorilla and orangutan[END_REF][START_REF] Winkler | Aspects of dental development in the orangutan prior to the eruption of the permanent dentition[END_REF][START_REF] Winkler | Development of the orangutan permanent dentition: assessing patterns and variation in tooth development[END_REF][START_REF] Schwartz | Ontogeny of canine dimorphism in extant hominoids[END_REF][START_REF] Kelley | Dental development and life history in living African and Asian apes[END_REF][START_REF] Smith | Dental development in living and fossil orangutans[END_REF]. In contrast, the three molars develop separately with little temporal overlap in crown formation time (the first molars M1/m1 are the first teeth to form around birth; the initiation of calcification of M2/m2 occurs before -or at least simultaneous -the end of the crown completion of M1/m1; likewise, the initiation of calcification of M3/m3 occurs before -or at least simultaneous -the end of the crown completion of M2/m2) (Figure 2SI). Therefore, we used crown formation time span of the permanent molars and defined three age groups: ~birth -3 years (M1/m1), ~2 -5 years (M2/m2), and ~4 -8 years (M3/m3) (Figure 2SI). Given the large variability observed in extant individuals, the intervals used here are large enough to encompass mean estimates of molar crown formation times in Asian apes [START_REF] Dean | Developing pongid dentition and its use for ageing individual crania in comparative cross-sectional growth studies[END_REF][START_REF] Anemone | Dental development of known-age chimpanzees, Pan troglodytes (Primates, Pongidae)[END_REF][START_REF] Kuydendall | Dental development in chimpanzees (Pan troglodytes): The timing of tooth calcification stages[END_REF][START_REF] Winkler | Development of the orangutan permanent dentition: assessing patterns and variation in tooth development[END_REF][START_REF] Reid | A histological reconstruction of dental development in the common chimpanzee, Pan troglodytes[END_REF][START_REF] Smith | Dental development in living and fossil orangutans[END_REF][START_REF] Dean | Retrieving chronological age from dental remains of early fossil hominins to reconstruct human growth in the past[END_REF].

In the second analysis, considering the differential representation of tooth types in sites (Table 3), only seven lower canines of Tham Khuyen could be addressed in a comparative study with those of recent Pongo. Canines have the longest crown formation time of all tooth classes, with crown mineralization beginning shortly after birth until around ten years of age [START_REF] Beynon | Histological study on the chronology of the developing dentition in gorilla and orangutan[END_REF][START_REF] Dean | A longitudinal study of tooth growth in a single individual based on long-and short-period incremental markings in dentine and enamel[END_REF][START_REF] Winkler | Development of the orangutan permanent dentition: assessing patterns and variation in tooth development[END_REF][START_REF] Schwartz | Ontogeny of canine dimorphism in extant hominoids[END_REF]. In Pongo, like in other highly dimorphic large-bodied apes (Gorilla), the difference in canine crown height between males and females [START_REF] Kelley | Sexual dimorphism in canine shape among great apes[END_REF][START_REF] Waddle | Sexing isolated hominoid canines with special reference to the Middle Miocene specimens from Paṣalar, Turkey[END_REF], is due to distinct growth durations. With a range of 7.2 -9.6 years (8.83 years on average), males complete enamel formation later, after a longer time than females (4.5 -6.3 years; 5.47 years on average; [START_REF] Schwartz | Ontogeny of canine dimorphism in extant hominoids[END_REF]. Therefore, this long duration of crown formation combined with a high susceptibility to record physiological stresses [START_REF] Goodman | Factors affecting the distribution of enamel hypoplasias within the human permanent dentition[END_REF]Skinner, 1986a), makes canines display the largest number of LEH episodes of all types of teeth. Moreover, in males, canines record hypoplastic episodes when enamel formation of the M3/m3 has ceased and can then reveal additional information on the stress events that followed.

Each Tham Khuyen canine has been sexed based on measurements of crown height taken in specimens of recent Pongo pygmaeus (P. p. pygmaeus and P. p. abelii): males (range, 17.2 -28.2 mm, n = 20) and females (range, 14.1 -23.1 mm, n = 20) [START_REF] Waddle | Sexing isolated hominoid canines with special reference to the Middle Miocene specimens from Paṣalar, Turkey[END_REF]. The Tham Khuyen collection is composed of six canines from males, and one canine from a female. Taking into account that two canines (TK138 and TK157) could belong to a single individual, we obtained a MNI of 4 (three males and one female).

In addition to the LEH frequency and expression, we used some parameters defined by [START_REF] Skinner | Hypothesis for the causes and periodicity of repetitive linear enamel hypoplasia in large, wild African (Pan troglodytes and Gorilla gorilla) and Asian (Pongo pygmaeus) Apes[END_REF], such as the relative spacing between two adjacent episodes of LEH (the interval measured between the end of one LEH episode and the beginning of the next), the timing of the apparent onset of stress (the distance from the location of the first LEH to the cervix), and the duration of stress inferred from the span between the first and the last episode of LEH), all standardized against crown height. All the parameters used here have been defined using spatial measurements, and not the count of perikymata, for a comparable precision, and thus an accurate comparison. A Mitutoyo digital caliper was used to measure distances on canine crowns.

In the present intrageneric analysis which includes past and present Pongo species, we made the assumption that the enamel of all specimens shares common intrinsic characteristics. These include histological properties, the susceptibility of teeth to hypoplasia, the spacing of perikymata over the height of the crown, striae angles, the crown formation times, and the enamel thickness (see [START_REF] Hannibal | Linear Enamel hypoplasia in the great apes: analysis by genus and locality[END_REF], for further details on the intergeneric variability of these factors in great apes).

Results

LEH prevalence on Pongo milk teeth

The observation on recent immature individuals (N = 10) having complete deciduous dentition and permanent M1/m1, shows the low prevalence of hypoplasia on milk teeth. Only one (out of 10 individuals) has affected milk dentition. In this individual, one LEH defect is located on the basal third of M1/m1, and another one is visible as a thin depression of enamel on the mid-portion of deciduous canines (in a second individual, one defect is also located on the basal third of M1/m1, and there is no observable LEH on milk teeth). We did not observe localized hypoplasia in this Pongo sample, as described by several authors (Skinner, 1986b;[START_REF] Lukacs | Enamel hypoplasia in deciduous teeth of great apes: do differences in defect prevalence imply differential levels of stress?[END_REF][START_REF] Skinner | Localized hypoplasia of the primary canine in bonobos, orangutans and gibbons[END_REF].

Only the Duoi U'Oi collection contains two milk teeth, one upper central incisor and one lower canine (MNI = 1), both without LEH or localized hypoplasia.

LEH prevalence on Pongo permanent teeth

In our recent Pongo sample (N = 17), 94.1 % of adult individuals (16/17) were affected by hypoplasia (individuals showing at least one LEH episode). Interestingly, the only specimen which shows no LEH marks is the individual born in the zoo.

In the fossil samples, Coc Muoi has the highest percentage of teeth with LEH marks (70 %), Duoi U'Oi the lowest (43.3 %), and Tham Khuyen is in between (46.6 %) (Table 5).

At Coc Muoi (MNI = 2), both individuals were affected by hypoplasia, 2 out of 3 individuals at Duoi U'Oi, and 4 out of 6 individuals at Tham Khuyen.

Number of LEH episodes (incisors, premolars, and molars)

In recent Pongo, results indicate that 29.4 % of individuals (5/17) experienced one LEH episode, 23.5 % (4/17) two LEH episodes, and 29.4 % three LEH episodes. Cases with four or five LEH episodes were few (2/17) (Table 6).

Using the three tooth types (incisors, premolars, and molars), results show that the Coc Muoi teeth recorded one (47.3 %) to two (26.3 %) LEH episodes, one (9.4 %) up to three (9.4 %) at Tham Khuyen, and one (27.5 %) up to four (3.4 %) at Duoi U'Oi (Table 6). If we consider overall fossil data, teeth with more than three LEH episodes are few, following the pattern observed in the recent Pongo sample.

4.4.

Number of LEH episodes per period Table 7 illustrates the number of LEH episodes observed in recent Pongo (N = 17), using the three periods defined on crown formation time of the permanent molars. Most individuals (16/17) show no LEH episode between birth and ~3 years of age (i.e. strong enough to have been recorded during crown formation of the first molars, and visible under low magnification). Only one specimen displays LEH on M1/m1 (a faint furrow only visible on the lingual side of crown in one female individual). The observation of the additional immature individuals (N = 10) confirmed the low prevalence of hypoplasia in the first three years of life with only 2 out of 10 individuals showing a line of hypoplasia on the permanent first molars. All other LEH episodes occurred in the subsequent periods, i.e. during crown formation times of the second and third molars, approximately between ~2 and 8 years of age, with great variability between individuals.

Table 8 illustrates the number of LEH episodes (range and mean) per period in recent Pongo (only individuals showing at least one episode have been selected here; N = 16). Results indicate that the period ~2 -5 years was critical for immature individuals. Males had two LEH episodes on average (8 out of 11 individuals), and females 1.5 episodes (all 4 individuals). The frequency of episodes decreased in the ~4 -8 years period. Males had 1.5 LEH episodes on average (6 out of 11 individuals), and females only one episode (2 out of 4 individuals). The amount of data is however insufficient for a statistical analysis of the prevalence in males versus females. When considering the overall sample, the ~2 -5 years period evidently appears to be the period recording the greatest number of episodes with 1.92 on average (13/16 individuals), compared to 1.37 for the following period ~4 -8 years (8/16 individuals).

Table 9 indicates the number and range of LEH episodes observed on fossil teeth in each period (see Table 4SI for details). Using teeth that can be placed easily on dental rows (incisors, premolars, and third molars), the results show a greater record of LEH episodes during the periods ~2 -5 years (up to 3) and ~4 -8 years (up to 2), than during the first three years of life, like in their modern counterparts.

Although based on samples of small size, data also suggest that the ~2 -5 years interval was a critical period for immatures, with one individual at Coc Muoi, and two individuals at Duoi U'Oi and Tham Khuyen. Only Tham Khuyen shows however more affected individuals (all four individuals) in the ~4 -8 years period, than is in the previous period (~2 -5 years) (two individuals), in contrast to what is observed in recent Pongo. But the results are most likely biased due to the greater number of identified M3/m3 in this collection (Table 4SI).

In Table 10, we considered two distinct intervals, ~0 -5 years of age (in which we gathered incisors, premolars and all undistinguished first and second molars), and ~4 -8 years (third molars) (Table 4SI). The results support a greater record of LEH episodes before ~5 years (up to 4 at Duoi U'Oi), than between ~4 -8 years (up to 2 at Coc Muoi and Tham Khuyen). Results also indicate, as previously noted, a relatively greater number of affected individuals in the last period (~4 -8 years) at Tham Khuyen. In this site, three individuals (out of the 4 impacted) suffered stresses over the two periods.

4.5.

Type of defects Assessment of the type of the defects (furrow vs line of pits) in our recent Pongo sample reveals that furrows are the most frequent expression of hypoplasia, in all teeth analyzed (incisors, premolars, and molars: 62.5 to 100 % of individuals), compared to lines of pits with 11.1 to 25 % of individuals (Table 11). Contrary to incisors and premolars, molars rarely exhibit both types of hypoplasia (8.3 %).

Four individuals exhibit a line of small and shallow pits visible on one or few tooth types, while the same stress event generated a furrow on adjacent teeth, most likely due to different susceptibility to developmental hypoplasia between tooth types (Table 5SI). In other three individuals, the line which has the appearance of a groove with deep pits can be matched on several or all teeth, that suggests a greater severity of the stressor (Figure 1B).

The results also highlight that among the various tooth types, the P3/p3 are the most frequently impacted (15 out of 16 individuals, 93.7 %; Table 5SI). We observed on jaws of recent Pongo, that the lower p3 is more likely to display LEH with sometimes deep horizontal troughs that follow lines of hypoplasia, visible on the labial side of the crowns (Figure 1A). These deep troughs seem to be due to a greater fragility of enamel in areas already affected by hypoplasia, a phenomenon most likely accentuated by the position of the tooth on the arcade, its particular labial profile and elongated shape.

Fossil data indicate the predominance of furrows as the commonest type of hypoplasia: all teeth at Duoi U'Oi; 13 out of 14 teeth (92.8 %) at Coc Muoi; 20 out of 21 teeth (95.2 %) at Tham Khuyen) (Table 6SI). One out of the 2 individuals at Coc Muoi displays one line with small and shallow pits on a M1/M2, and 1 out of the 4 individuals at Tham Khuyen shows two lines of small pits on a lower i1. We did not observe any fossil tooth with both lines of pits and furrows.

Distribution of LEH episodes over canine crowns

In the comparison between the Tham Khuyen lower canines and those of recent Pongo (Figure 3), using three parameters (relative spacing between two adjacent episodes, timing of apparent onset of stress, and duration of stress), the Tham Khuyen data fall within ranges of variability observed in two recent Pongo samples (Table 7SI).

It is also the case for the number of LEH episodes visible throughout the crown height. Results in Table 12 show that this number in the Tham Khuyen sample is close to that observed in recent Pongo, with a range of 1 -10, and a mean of 6.4.

According to the median (Table 12), the value obtained for Tham Khuyen is high (x = 8), and close to that of recent P. abelii (n = 7) (data of [START_REF] Guatelli-Steinberg | Linear enamel hypoplasia as an indicator of physiological stress in great apes: Reviewing the evidence in light of enamel growth variation[END_REF]. Results are most likely biased due to a greater number of male individuals (vs females) in both samples, 3 males and one female at Tham Khuyen, and 7 males and 2 females in the P. abelii population.

Location, type and age range estimate of the most accentuated LEH episode on canines

The assessment of the pattern of LEH expression on canines shows that furrows are more common than line of pits: only 2 out of 9 affected individuals exhibit one line of pits (vs 3 to 8 furrows) in our recent Pongo sample, and 2 out of 4 individuals exhibit two lines of pits in Tham Khuyen (vs 1 to 10 furrows).

In order to identify a specific stress event during the development of immature individuals in recent and fossil samples, we assessed the location of the most accentuated linear marks (either furrow or line of pits) visible on the labial surface of the canines, and estimated the age range of individuals based on crown deciles from Skinner and Skinner (2017) (Table 13).

In recent Pongo, 3 out of 9 individuals have similar furrows over the crown height. In the other six individuals, results show that the expression of the most accentuated LEH episode can vary (one wide line of pits, one wide furrow, two close furrows). Among these six individuals, five present a defect which occurred when around two-thirds of the crowns were formed (0.45 -0.23). The corresponding age range can be estimated for the males between 4.84 and 7.25 years, and for the female between 3.45 and 4.77 years [START_REF] Skinner | Orangutans, enamel defects, and developmental health: A comparison of Borneo and Sumatra[END_REF] (Table 13). The last individual was the only one showing an accentuated wide line of pits (0.66 -0.60) that developed early during crown formation (Figure 1A), broadly before 4 years of age.

The observation of two canines from Tham Khuyen belonging to two male individuals, shows that the most accentuated defects were recorded when more than two-thirds of crowns were formed (0.30 -0.17), corresponding to a biological age range of 6.46 -8.09 years (Table 13). Furthermore, two canines (TK138 and TK157) which could belong to one male individual, display a large zone of disturbed enamel, deep enough to have possibly involved a defect of the underlying dentine (which cannot be ascertained here without microCT-scan) (Figure 3C). A comparable depression-type hypoplasia has been also observed in pigs [START_REF] Dobney | A protocol for recording linear enamel hypoplasia on archaeological pig teeth[END_REF]Witzel et al., 2006). It has also the characteristics of "coronal waisting" described in a chimpanzee population [START_REF] Skinner | Developmental defects of the dental crown in chimpanzees from the Taï National Park, Côte d'Ivoire: Coronal waisting[END_REF]: macroscopic constriction located on the mid-canine crown, well visible on the labial side. On one canine (TK138), the width of the depression could be estimated in the range ~0.20 -0.48, that corresponds broadly to a developmental stress that lasted between 4.84 and 7.25 years (Table 13).

Discussion

Hypoplasia and life history

As shown in several publications [START_REF] Guatelli-Steinberg | Linear enamel hypoplasia in gibbons (Hylobates lar carpenteri)[END_REF][START_REF] Skinner | Hypothesis for the causes and periodicity of repetitive linear enamel hypoplasia in large, wild African (Pan troglodytes and Gorilla gorilla) and Asian (Pongo pygmaeus) Apes[END_REF][START_REF] Hannibal | Linear Enamel hypoplasia in the great apes: analysis by genus and locality[END_REF][START_REF] Guatelli-Steinberg | Linear enamel hypoplasia as an indicator of physiological stress in great apes: Reviewing the evidence in light of enamel growth variation[END_REF], LEH is a common phenomenon in extant Pongo species (94.1 % of individuals in our reference sample, 100 % when excluding the zoo-born individual). Although based on small datasets, results on incisors, premolars, and molars show that LEH frequencies have different distributions in each period of development: rare episodes from birth to ~3 years of age, frequent episodes from ~ 2 to 5 years, and comparatively fewer episodes from ~4 to 8 years. Our observations on the seven tooth types match those obtained on lower incisors and canines, all of which show that a specific period during the development of young individuals contains the greatest number of LEH episodes. The period corresponds to the formation of lateral enamel in the mid section of canine crowns [START_REF] Skinner | Hypothesis for the causes and periodicity of repetitive linear enamel hypoplasia in large, wild African (Pan troglodytes and Gorilla gorilla) and Asian (Pongo pygmaeus) Apes[END_REF][START_REF] Guatelli-Steinberg | Linear enamel hypoplasia as an indicator of physiological stress in great apes: Reviewing the evidence in light of enamel growth variation[END_REF]. Regarding fossil Pongo, despite the small datasets, the present analysis suggests that hypoplasia might also have been frequent in past populations, with the two Coc Muoi individuals affected, 2 out of 3 individuals at Duoi U'Oi, and 4 out of 6 individuals at Tham Khuyen.

With regard to the morphology of LEH defects, it seems that both recent and fossil samples have a comparable pattern, with furrows more frequent than line of pits. Furrows and linear pits have similar systemic causes, and are associated with stresses due to undernutrition and disease [START_REF] Goodman | Assessment of systemic physiological perturbations from dental enamel hypoplasias and associated histological structures[END_REF][START_REF] Hillson | The relationship of enamel hypoplasia to the pattern of tooth crown growth: a discussion[END_REF] (differently, scattered of pits over the surface of tooth crowns are associated with postnatal deficiencies of vitamins A and D, and of calcium; [START_REF] Lovell | Patterns of Dental Enamel Defects at Ancient Mendes, Egypt[END_REF]Witzel et al., 2006).

The distribution of LEH episodes in fossil Pongo also follows a comparable pattern to that of the modern species. A combination of factors might explain the low frequencies of hypoplasia in the period "birth -3 years": first, the small number of involved teeth and the fact that teeth are less likely to display LEH. With regard to incisors, only the appositional zone of crowns is concerned (Table 3SI). In this zone, enamel layers do not appear on the surface and, therefore, don't show evidence of hypoplasia [START_REF] Hillson | The relationship of enamel hypoplasia to the pattern of tooth crown growth: a discussion[END_REF]. The first molars, which begin to form around birth, is also few affected by hypoplasia, most likely due to the low susceptibility to record LEH events (3.4 -4.8 %, Hannibal and Guatelli-Steinberg, 2005;6.2 % in our reference Pongo sample, Table 5SI). Immature individuals also have a particular immunity against stressors, obtained from the mother by the nursing, as previously noted by [START_REF] Skinner | Hypothesis for the causes and periodicity of repetitive linear enamel hypoplasia in large, wild African (Pan troglodytes and Gorilla gorilla) and Asian (Pongo pygmaeus) Apes[END_REF] who estimated the timing of the onset of stress in apes at around 2.6 years of age. The first three years of life constitute the period of most frequent body contacts between mother and offspring. During the first three months, infants are always in contact with the mother (100% of the time). This percentage decreases continuously to reach 50 % ca. 2 years of age, and 16 % ca. 2.5 years of age [START_REF] Van Noordwijk | Development of ecological competence in Sumatran orangutans[END_REF][START_REF] Van Noordwijk | Development of independence[END_REF][START_REF] Van Adrichem | The development of wild immature Sumatran orangutans (Pongo abelii) at Ketambe[END_REF]. During the first stage of their lives, young orangutans keep this proximity with their mothers during most activities (feeding, transportation, resting, nesting, etc.), and are less confronted to disturbances due to fluctuations in fruit availability since they feed from their mother's milk, or are less confronted to infections due to parasites since they still have few contacts with the vegetation.

The following period, from around ~2 to ~5 years of age, is obviously a period of vulnerability for immature individuals, as results obtained on our recent Pongo data show, with 13 out of 17 adult individuals who suffered up to 5 stress events during that time. In fossil Pongo, results show that 2 out of 3 individuals at Duoi U'Oi, and 2 out of 6 individuals at Tham Khuyen suffered up to 3 stresses. During this period, a great number of tooth types form concomitantly, and defects are found either in the whole imbricational zone of crowns (I1/i1, P4/p4, M2/m2), or in a portion of the imbricational zone (I2/i2, P3), that could explain in part high hypoplastic rates. Some teeth also have the higher hypoplastic susceptibility: I1/i1 (44.1 -50 %), I2/i2 (37.5 -63 %), P3/p3 (47.2 -62.7 %), P4/p4 (34.4 -35.7 %), M2/m2 (9.6 -19.3 %) [START_REF] Hannibal | Linear Enamel hypoplasia in the great apes: analysis by genus and locality[END_REF], and between 56.2 to 93.7 % in our reference Pongo sample (Table 5SI). Yet, important changes in the behavior of infants also characterize this period and could provide further explanation for why they react more intensively to stressors. Although still largely dependent on their mothers (for nutrition, protection, warmth, etc.), they begin to explore their immediate environment (less than 10 meters away), acquire locomotor skills, become able to process food items, and learn to build nests [START_REF] Van Noordwijk | Development of independence[END_REF]. This acquisition of foraging independence is also highly correlated with the onset of the weaning process [START_REF] Delgado | The Behavioral Ecology and Conservation of the Orangutan (Pongo pygmaeus): A Tale of Two Islands[END_REF], with all subsequent psychological and physiological consequences on the health status of individuals [START_REF] Skinner | Developmental defects of the dental crown in chimpanzees from the Taï National Park, Côte d'Ivoire: Coronal waisting[END_REF]. The etiology of LEH defects for this period is however impossible to reveal since it can be linked to multiple causes: diseases due to the virulence of parasites and other pathogens while infants are not yet resistant, malnutrition due to difficulties in processing food or fluctuations in fruit availability, psychologically stressful situations linked to first rejections by the mother, interaction with the other members of the group, or possible first encounters with arboreal predators, such as clouded leopards or pythons [START_REF] Rijksen | A field study on Sumatran orang utans (Pongo pygmaeus abelii Lesson 1827): Ecology, behavior, and conservation[END_REF].

The ~4 -8 years period shows lower LEH frequencies in the recent Pongo sample with 8 out of 17 adult individuals who suffered up to 3 stress events. At Duoi U'Oi, one individual suffered one stress event and, at Tham Khuyen, the 4 individuals suffered up to 2 stress events. The period coincides with the end of the crown formation of some teeth, I2/i2, p3 (only distal thirds of crowns), and with the crown formation of M3/m3, also less susceptible to record LEH episodes (15.3 -36.1 %;[START_REF] Hannibal | Linear Enamel hypoplasia in the great apes: analysis by genus and locality[END_REF]56.2 % in our reference Pongo sample, Table 5SI). That transition phase corresponds to the progressive learning of independence for young individuals who then explore their environment far from their mothers (>50 m) [START_REF] Van Noordwijk | Development of independence[END_REF], and progression towards complete weaning.

In recent orangutans, weaning is a long and gradual process, the aspects of which, social, nutritional and psychological, are difficult to describe by observers in the wild. In a Sumatran population, [START_REF] Van Noordwijk | Development of ecological competence in Sumatran orangutans[END_REF] observed that mothers began to reject suckling visits of the youngest infant around 5 years old, while the last observed suckling visit was that of a 7-year-old individual. The complete weaning age was estimated by these authors at around 8 years. However, it seems to occur earlier, around 5 -6 years old, in some individuals of Bornean populations [START_REF] Van Noordwijk | Development of independence[END_REF]. The interbirth interval also provides data on the age of juveniles at the time of weaning, a new pregnancy occurring when they are completely weaned (but see [START_REF] Smith | Cyclical nursing patterns in wild orangutans[END_REF]. Based on data compiled by Wish et al. (2009;Table 5.1) and by [START_REF] Knott | The ecology of female reproduction in wild orangutans[END_REF]Table 11.1), mean intervals evaluated in two females of a Sumatran population (8.2 -9.3 years) also appear significantly longer than those evaluated in three females of a Bornean population (6.1 -7.7 years).

In relation to wild Sumatran populations, this ecological background of a late weaning has been be related to the high-quality of the food available, with more fruit and insects (vs more bark, unripe fruit and seeds in the diet of Bornean populations), and to fewer fluctuations in fruit availability over time (Galdikas and Wood, 1990;[START_REF] Delgado | The Behavioral Ecology and Conservation of the Orangutan (Pongo pygmaeus): A Tale of Two Islands[END_REF][START_REF] Jablonski | The influence of life history and diet on the distribution of catarrhine primates during the Pleistocene in eastern Asia[END_REF][START_REF] Wich | Determinants of orangutan density in the dryland forests of the Leuser Ecosystem[END_REF][START_REF] Wich | Dietary and energetic responses of Pongo abelii to fruit availability fluctuations[END_REF][START_REF] Wich | Orangutan life history variation[END_REF][START_REF] Van Noordwijk | Development of ecological competence in Sumatran orangutans[END_REF][START_REF] Knott | The ecology of female reproduction in wild orangutans[END_REF].

The location of the most pronounced LEH defects on canine crowns in our recent Pongo sample corresponds to stress events which occurred in this ~4.8 -7.2 years interval (four males), or earlier ~3.4 -4.7 years (one female). Such morphology of enamel hypoplasia (in contrast with the morphology of adjacent lines within the same tooth), suggests a greater intensity of the stressor, or the interaction of two or more factors [START_REF] Goodman | Assessment of systemic physiological perturbations from dental enamel hypoplasias and associated histological structures[END_REF]. That might be consistent with the long dietary and behavioral transition of the weaning, known to be one of the most stressful periods in various large-bodied mammals [START_REF] Mead | Enamel hypoplasia in Miocene rhinoceroses (Teleoceras) from Nebraska: Evidence of severe physiological stress[END_REF][START_REF] Dobney | Interpreting developmental stress in archaeological pigs: the chronology of linear enamel hypoplasia[END_REF][START_REF] Franz-Odendaal | Insights from stable light isotopes on enamel defects and weaning in Pliocene herbivores[END_REF][START_REF] Franz-Odendaal | Enamel hypoplasia provides insights into the early systemic stress in wild and captive giraffes (Giraffa cameloparalis)[END_REF][START_REF] Niven | An inter-site comparison of enamel hypoplasia in bison: implications for paleoecology and modeling Late Plains Archaic subsistence[END_REF]Witzel et al., 2006).

The depression-type defect observed in one male individual (TK138; Table 13) from Tham Khuyen might correspond to the duration of weaning, estimated to have lasted around 2.5 years (~4.8 -7.5 years old). For Witzel et al. (2006), this type of defect reflects a prolonged period of undernutrition. The most accentuated episodes recorded in two other male individuals (TK46, TK118; Table 13) between ~6.5 and 8 years, also suggest stressful events during the weaning period.

5.2.

Paleoecology of large-bodied mammals The Coc Muoi data indicated that rhinoceroses and wild cattle were exposed to multiple physiological and psychological stress events specific to age, at around birth (rhinoceros), weaning (wild cattle, orangutans), abandonment of calf (rhinoceros), or when juveniles reached sexual maturity (both groups).

The Coc Muoi fauna gives a good representation of a tropical megaherbivore-bearing ecosystem that occurred in this region of Southeast Asia, in the late Middle Pleistocene, during the MIS6 -5, considering the time frame for breccia deposition between 148 and 117 ka (Bacon et al., 2018a). During the temperate MIS 6 phase, the climate remained particularly warm and wet at these latitudes. The age interval also includes the Eemian episode (5e: 124 -119 ka), which was marked by the strengthening of the East Asian summer monsoon [START_REF] Li | The stable isotope record in cervid tooth enamel from Tantang cave, Guangxi: Implications for the Quaternary East Asian monsoon[END_REF][START_REF] Wang | Millenial-and orbital-scale changes in the East Asian monsoon over the past 224,000 years[END_REF]. Under its influence, the climate was particularly warm and humid with most likely a great variation in seasonal rainfalls (Yuan et al., 2004;[START_REF] Ao | Different orbital rythms in the Asian summer monsoon records from the North and South China during the Pleistocene[END_REF]. The climatic factors of temperature and humidity, which generate damp and flooded soils, are considered to be the most important environmental factors influencing the proliferation and virulence of parasites, which are sources of diseases for immature largebodied mammals [START_REF] Skinner | Hypothesis for the causes and periodicity of repetitive linear enamel hypoplasia in large, wild African (Pan troglodytes and Gorilla gorilla) and Asian (Pongo pygmaeus) Apes[END_REF][START_REF] Foitová | Parasites and their impacts on orangutan health[END_REF].

In terms of paleoecology, Pleistocene orangutans from the Asian mainland could survive in different environmental conditions than those they occupy today, in nonanalog mammalian communities [START_REF] Stewart | The evolutionary consequence of the individualistic response to climate change[END_REF]. They occupied dense forests [START_REF] Bocherens | Flexibility of diet and habitat in Pleistocene South Asian mammals: Implications for the fate of the giant fossil ape Gigantopithecus[END_REF]Bacon et al., 2018b). Considering the location of fossiliferous deposits and the geological context of the area in northern Vietnam [START_REF] Duringer | Karst development, breccias history, and mammalian assemblages in Southeast Asia: A brief review[END_REF], it is highly plausible that they lived in the limestone forests on hills and tower karsts. Today, limestone karst forests are inhabited by wild populations of low density, due to relatively low tree species diversity and limited productivity (vs forested habitats such as peat-swamp forests, lowland Dipterocarp forests or dry forests) [START_REF] Marshall | Use of limestone karst forests by Bornean orangutans (Pongo pygmaeus morio) in the Sangkulirang Peninsula, East Kalimantan, Indonesia[END_REF][START_REF] Husson | Orangutan distribution, density, abundance and impacts of disturbance[END_REF].

Body weight estimates made on the only complete skeleton of a female Pongo to be found in the mainland (with a range of 81.6 -89.3 kg based on the crown areas of M1/m1, and 58.8 -83.1 kg based on long bones; [START_REF] Bacon | L'orang-outan de Hoà Binh (Viêt-Nam) : poids corporel et hypothèse locomotrice[END_REF], support a greater body size in the fossil species compared to the body size of females in extant populations (Eckhart, 1975;[START_REF] Markham | Brief communication: Weights of wild orangutans[END_REF][START_REF] Smith | Body mass in comparative primatology[END_REF]. Considering the influence of body mass on life-histories in large mammals (Owen [START_REF] Owen-Smith | Megaherbivores. The influence of very large body size on ecology[END_REF]Guadelli-Steinberg et al., 2000;[START_REF] Jablonski | The influence of life history and diet on the distribution of catarrhine primates during the Pleistocene in eastern Asia[END_REF], and the different ecological constraints due to past climatic conditions, the issue will be whether Pleistocene Pongo exhibited a different pattern of growth and development in comparison with extant orangutan populations.

Conclusion

In this study, we explored the prevalence and expression (furrow vs line of pits) of linear enamel hypoplasia (LEH) in Pongo from the Middle Pleistocene (Tham Khuyen, Coc Muoi) and the Late Pleistocene (Duoi U'Oi) karstic sites in northern Vietnam.

The results show that the LEH analysis based on isolated teeth is constrained by numerous biases including: the number of teeth; the differential representation of tooth types; the difficulty in distinguishing first from second molars; the small number of individuals [MNI]; the differential representation of males versus females within each Pongo sample.

The results indicate, nevertheless, that Pleistocene Pongo individuals experienced multiple LEH episodes over their long period of growth and development, with the two Coc Muoi individuals affected (MNI = 2), 2 out of 3 individuals at Duoi U'Oi (MNI = 3), and 4 out of 6 individuals at Tham Khuyen (MNI = 6). They had a comparable pattern of defects to those of their modern counterparts, most likely linked to the same causes (undernutrition, disease). Further analyses on larger samples from Southeast Asian sites will allow for a better understanding of the evolution of ecological variables in the Pongo lineage. Witzel, C., Kierdorf, U., Dobney, K., Ervynck, A., Vanpoucke, S., Kierdorf, H., 2006. Reconstructing impairment of secretory ameloblast function in porcine teeth by analysis of morphological alterations in dental enamel. Journal of Anatomy, 209, 93-110. Yuan, D., Cheng, H., Edwards, R.L., Dykoski, C.A., Kelly, M.J., Zhang, M., Qing, J., Lin, Y., Wang, Y., Wu, J., Dorale, J.A., An, Z., Cai, Y., 2004. Timing, duration, and transitions of the last interglacial Asian monsoon. Science 304, 575-578. 4SI). The distinction between first and second molars (M1/m1 and M2/m2) using only external morphology is often impossible (close morphology, presence of two facets of contact), whereas the distinction of M3/m3 from the two first molars is possible (morphology of the hypoconulid on m3 and absence of a facet of contact on the distal side of the crown on M3/m3). Figure in brackets refers to the number of teeth. 2012) are based on the scoring of accentuated lines or grooves under 20x magnification. m, male; f, female.

Table 13: Number of LEH episodes over the crown height (N), type (furrow vs line of pits) and location (distance from the cervix, standardized against crown height) of the most accentuated LEH episode visible on lower and (*) upper canines of nine individuals of our recent Pongo sample, and seven canines from Tham Khuyen (MNI = 4; three males and one female; TK138 and TK157 could belong to a single individual). Age ranges are based on age estimates of canine crown deciles from Skinner and Skinner (2017; Figure 5).

Figure caption

  (adult individuals). Location of LEH episodes on the crown surface of molars, incisors and premolars (occlusal third, middle third, and basal third). The three periods of deficiency are delimited using crown formation times of M1/m1 (~birth -3 years), M2/m2 (~2 -5 years), and M3/m3 (~4 -8 years). * not observed in the sample. Table 5: Number and percentage of permanent teeth showing LEH episodes in the three fossil collections. (1) Two M1; (2) Three m1 and four m2.Table 6: Number and percentage of individuals in the recent Pongo sample (N = 17), and number and percentage of teeth in the fossil collections, showing 0 to 5 LEH episodes on incisors, premolars, and molars. Table 7: Number of LEH episodes in recent Pongo individuals (N = 17) using incisors, premolars and molars, by period of development based on crown formation times of M1/m1 (~birth -3 years), M2/m2 (~2 -5 years), and M3/m3 (~4 -8 years). See Tables 3SI for the location of defects on teeth. Data indicate a great variability (1 out of 17 individuals shows no LEH episode; 8 out of 17 show LEH episodes only in the ~2 -5 years interval; 3 out of 17 only in the ~4 -8 years interval; 4 out of 17 individuals in the ~2 -8 years interval; and 1 out of 17 from birth to 8 years of age). The catalog number of skulls is noted in the second row of the Table.

Figure 1 :

 1 Figure captionFigure 1: Photos of permanent teeth of recent Pongo skulls showing the two morphologies (arrows) used in the paper for scoring LEH: A pit-form defect displays either small and shallow pits (A), or deep pits located in a marked groove (B); A furrow-form defect displays either a faint groove (C), or a deep groove (D).

Figure 2 :

 2 Figure 2: Photos of skulls of recent Pongo showing different types of LEH. We considered only the lines and furrows that can be matched on teeth forming at the same time. A. Young male adult mandible (1872-547) with one LEH episode; B. Male adult mandible (1927-98) with multiple LEH episodes on incisors, canines, and premolars (deep furrows and troughs are visible on lower p3 and p4 crowns); C. Young adult of unknown sex (AC1990) showing both lines of pits and furrows on upper incisors; D. Young female adult maxilla (A10683) showing faint furrows on crowns of upper premolars and molars (lingual view). Collections

Table caption Table 1 :

 caption1 Number and percentage of the teeth of large-sized bovids and rhinocerotids in the Coc Muoi and Duoi U'Oi collections showing no hypoplasia and one to four LEH episodes. In contrast to linear hypoplasia, permanent hypoplasia stands for long-lasting defects, observable on the entire height of the crown.

Table 2 :

 2 LEH in large ungulates from Coc Muoi and Duoi U'Oi. NIhyp, number of individuals with enamel hypoplasia; MNI, minimum number of individuals; Tooth type refers to the tooth used to calculate the MNI (the most frequent type of tooth, either right or left); Time range of crown formation (bovids,[START_REF] Niven | An inter-site comparison of enamel hypoplasia in bison: implications for paleoecology and modeling Late Plains Archaic subsistence[END_REF] rhinocerotids, Hillman-Smith et al., 1986;[START_REF] Tafforeau | Nature of laminations and mineralization in rhinoceros enamel using histology and X-ray synchrotron microtomography: potential implications for palaeoenvironmental isotopic studies[END_REF].

(1)

All m3s assigned to R. sondaicus;

(2)

One right m3 and one left m2 assigned to R. sondaicus and D. sumatrensis, respectively. No deciduous bovid teeth from Coc Muoi display LEH defects. Data on mean body weight of animals are from

[START_REF] Francis | A field guide to the mammals of South-East Asia[END_REF]

: Bos sauveli (~800 kg), Rhinoceros sondaicus (~1750 kg), Rhinoceros unicornis (~2250 kg), Dicerorhinus sumatrensis (~950 kg).

Table 3 :

 3 Composition of the three fossil Pongo collections, by tooth type, and position on jaws (upper vs lower).

Table 4 :

 4 Reference Table based on the observations of the 17 skulls of recent Pongo

Table 8 :

 8 Number (range and mean) of LEH episodes by period in males, females, and all recent Pongo individuals. Only individuals showing at least one episode have been considered (N = 16). n, number of affected individuals; N, number of individuals per group.

Table 9 :

 9 Number of LEH episodes (range) on incisors (I1/1, I2/i2), premolars (P3/p3, P4/p4), and molars (M3/m3) in the fossil Pongo samples. n, number of individuals; N, number of individuals per group (see Table

Table 10 :

 10 Number of LEH episodes (range) on all fossil Pongo teeth. The period ~0 -5 years allows for gathering all incisors, premolars and M1-2/m1-2 in a single subsample. n, number of individuals; N, number of individuals per group. Figure in brackets refers to the number of teeth.

Table 11 :

 11 Number and percentage of individuals showing different types of defects (either furrow or line of pits, or both) by tooth type, in the recent Pongo sample. n, number of individuals; N, number of individuals per group. We consider the individuals showing at least one LEH episode (N = 16). See Table 5SI for details.

Table 12 :

 12 MNI, sex of individuals, and number of LEH episodes (range, mean, and median) on lower canines of Tham Khuyen, and canines of two recent Pongo samples (our reference sample with 3 upper and 4 lower canines; and another sample of lower canines from Guatelli-Steinberg et al. (2012). Data from Guatelli-Steinberg et al. (

Acknowledgments

We thank all members of the Institute of Archaeology in Hanoi, and particularly Ass. Prof. Nguyen Giang Hai, Director of the Institute, for the agreement to study the Duoi and Tham Khuyen collections. We thank also Mr. Nong Duc Kien, present Director of the Lang Son Museum, and Mr. Nguyen Gia Quyen and the other members of the Lang Son museum, regarding the study of the Coc Muoi collection.

Our thanks also go to Christine Lefèvre, Joséphine Lesur and Aurélie Verguin of the Laboratoire Mammifères et Oiseaux, Anatomie comparée, Muséum national d'Histoire naturelle, in Paris, for agreement and access to the Pongo collection ; Clément Zanolli, laboratoire PACEA, UMR 5199, Université de Bordeaux, for his helpful comments, Amandine Blin, Muséum national d'Histoire naturelle, Centre national de la Recherche scientifique, UMS 2700 2AD, in Paris, for her advice on possible statistical tests; Vincent Zvénigorosky for editing the text in English, and Professor Nguyen Lan Cuong from the Insitute of Archaeology in Hanoi, for his kind help to take photos of teeth of the Tham Khuyen collection. Finally, the authors thank the anonymous reviewers for their helpful critiques and suggestions.