# Dental macrowear reveals ecological diversity of Gorilla spp.

Teagan Harty, Michael Berthaume, Eugenio Bortolini, Alistair Evans, Jordi Galbany, Franck Guy, Ottmar Kullmer, Vincent Lazzari, Alejandro Romero, Luca Fiorenza

► To cite this version:

Teagan Harty, Michael Berthaume, Eugenio Bortolini, Alistair Evans, Jordi Galbany, et al., Dental macrowear reveals ecological diversity of Gorilla spp.. Scientific Reports, 2022, 12 (1), pp.9203. 10.1038/s41598-022-12488-8. http://hal-03874666

HAL Id: hal-03874666

https://hal.science/hal-03874666

Submitted on 14 Nov 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

OPEN

# Dental macrowear reveals ecological diversity of *Gorilla* spp.

Teagan Harty<fnref n="1" />, Michael A. Berthaume², Eugenio Bortolini³,⁴, Alistair R. Evans⁵,⁶, Jordi Galbany⁷,⁸, Franck Guy⁹, Ottmar Kullmer¹⁰,¹¹, Vincent Lazzari⁹, Alejandro Romero¹²,¹³ & Luca Fiorenza<fnref n="1" />

Size and shape variation of molar crowns in primates plays an important role in understanding how species adapted to their environment. Gorillas are commonly considered to be folivorous primates because they possess sharp cusped molars which are adapted to process fibrous leafy foods. However, the proportion of fruit in their diet can vary significantly depending on their habitats. While tooth morphology can tell us what a tooth is capable of processing, tooth wear can help us to understand how teeth have been used during mastication. The objective of this study is to explore if differences in diet at the subspecies level can be detected by the analysis of molar macrowear. We analysed a large sample of second lower molars of Grauer's, mountain and western lowland gorilla by combining the Occlusal Fingerprint Analysis method with other dental measurements. We found that Grauer's and western lowland gorillas are characterised by a macrowear pattern indicating a larger intake of fruit in their diet, while mountain gorilla's macrowear is associated with the consumption of more folivorous foods. We also found that the consumption of herbaceous foods is generally associated with an increase in dentine and enamel wear, confirming the results of previous studies.

Gorillas have been traditionally described as specialised folivorous primates based on the early behavioural studies from populations of the Virunga Mountains<fnref n="1" />,². However, additional ecological studies have shown that the most common gorillas, who inhabit western lowland rainforests, eat a considerable variety of foods including large amounts of fruits³,⁴ (Fig. 1).

Gorillas are found in Central and Western Equatorial Africa, and they are separated by the Congo River into western and eastern gorillas⁶. The western gorillas are divided into two subspecies, the western lowland gorillas (*Gorilla gorilla gorilla*; Savage, 1847) and the Cross River gorillas (*Gorilla gorilla diehli*; Matschie, 1904)⁷. The western lowland gorillas are characterised by a coarse and short brown hair that covers most of their body, often red on crown, by flared nostrils, and by a prominent lip above the nasal septum⁶,⁷. They are geographically and genetically separated from the eastern gorillas, and they are distributed in southern Cameroon, Equatorial Guinea, Gabon, Republic of Congo, northern Angola, and western Democratic Republic of Congo⁶,⁸. The Cross River gorillas are smaller in size, and differ from western lowland gorillas both in the shape of their cranium and in their genetics⁶,⁷. They can be found in a small region of highland forests located along the borders of Cameroon and Nigeria. Cross River gorillas are vastly understudied, but data suggests they have diets similar to other western lowland gorillas⁹.

The eastern gorillas consist of two subspecies: mountain gorillas (*Gorilla beringei beringei*; Matschie, 1903), and eastern lowland gorillas, or Grauer's gorillas (*Gorilla beringei graueri*; Matschie, 1914)⁷. Eastern gorillas are characterised by long black hair, by narrow nostrils, and by the absence of a lip above the nasal septum⁶,⁷. The

<footnote n="1">¹Biomedicine Discovery Institute, Department of Anatomy and Developmental Biology, Monash University, Melbourne, VIC, Australia. ²Division of Mechanical Engineering and Design, London South Bank University, London, UK. ³Department of Archaeology and Anthropology, Institución Milá y Fontanals de Investigación en Humanidades, Barcelona, Spain. ⁴Culture and Socio-Ecological Dynamics, Department of Humanities, Universitat Pompeu Fabra Ramon Trias Fargas, Barcelona, Spain. ⁵School of Biological Sciences, Monash University, Melbourne, VIC, Australia. ⁶Geosciences, Museums Victoria, Melbourne, VIC, Australia. ⁷Department of Anthropology, Center for the Advanced Study of Human Paleobiology, The George Washington University, Washington, DC, USA. ⁸Department of Clinical Psychology and Psychobiology, University of Barcelona, Barcelona, Spain. ⁹Laboratory PALEVOPRIM, UMR CNRS 7262, University of Poitiers, Poitiers, France. ¹⁰Department of Paleoanthropology, Senckenberg Research Institute and Natural History Museum Frankfurt, Frankfurt, Germany. ¹²Department of Paleobiology and Environment, Institute of Ecology, Evolution, and Diversity, Goethe University, Frankfurt, Germany. ¹⁲Departamento de Biotecnología, Universidad de Alicante, 03690 Alicante, Spain. ¹³Instituto Universitario de Investigación en Arqueología y Patrimonio Histórico (INAPH), Universidad de Alicante, 03690 Alicante, Spain. ¹email: luca.fiorenza@monash.edu</footnote>
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Map of Sub-Saharan Africa showing the distribution of western lowland gorillas (*Gorilla gorilla gorilla*, in light blue), cross river gorillas (*Gorilla gorilla diehli*, in dark blue), Grauer's gorillas (*Gorilla beringei graueri*, in dark red), and mountain gorillas (*Gorilla beringei beringei*, in orange). Adapted from Xue et al.⁵.</figcaption></figure>

mountain gorillas have a large head, wide facial skeleton, and shorter and stockier forelimbs than the western gorillas⁶,⁷. They live in the high-altitude rainforests of Virunga Volcanoes, along the borders between Rwanda, Democratic Republic of Congo and Uganda, and in the Bwindi Impenetrable Forest between the southwest of Uganda and northwest Rwanda⁶. Grauer's gorillas have a narrower facial skeleton compared to mountain gorillas and a shorter pelage on the scalp. They can be found in eastern forests of the Democratic Republic of Congo⁶,⁸.

When rainfall is low and fruit is scarce, gorillas inhabiting lowlands rely on more herbaceous foods, including leaves, stems, pith and barks¹⁰. These foods fall into the general categories of terrestrial herbaceous vegetation (THV) and aquatic herbaceous vegetation (AHV). Less commonly, other foods, like the woody endocarps of the *Coula edulis* seeds, are incorporated in the diets of some populations of western lowland gorillas, indicating a broader dietary spectrum than previously thought¹¹. Both western lowland gorillas and Grauer's gorillas prefer relatively soft, ripe fruits¹⁰,¹². However, the percentage of THV consumed by western lowland gorillas and Grauer's gorillas differs, with a larger amount of fibrous food eaten by the eastern populations¹³. On the other hand, mountain gorillas mostly rely on THV foods, which are available year-round, while fruit, especially in the high-altitude forests of the Virunga Mountains, are relatively scarce throughout the year²,¹⁴,¹⁵.

Molar dental studies have found that mountain and Grauer's gorillas are characterised by molars with duller cusps compared to western lowland gorillas¹⁶. When duller cusps wear, they create longer compensatory crests, increasing the contact areas between the sharp blade-like crests on the molars and the food particles. This provides a functional advantage compared to molars with sharper cusps, as it increases the cutting surfaces for processing fibrous foods¹⁷. Another study utilised the Dirichlet normal energy (DNE; a dental topographic method) to analyse the diet of seven groups of great apes, finding a positive correlation between DNE and fibre content in apes among sympatric pairs¹⁸. DNE is a measure of surface curvature, meaning teeth with higher DNE values are curvier, or sharper, potentially making it easier to break down foods with a high fibre content¹⁹. The presence of character displacement in DNE suggests the relationship between tooth curvature and fibre content is affected not only by diet, but also other parameters, such as indirect interspecific competition over food resources.

Moreover, a study based on buccal microwear signals found geographic differences between various populations of Grauer's gorillas which were attributed to ecological conditions²⁰. Dental microwear texture analyses have compared western lowland gorillas with mountain gorillas without finding any clear differences between the two groups²¹. This is probably due to dietary overlap between mountain and lowland gorillas, which may be large enough that microwear differences are not observable. In addition to this, tooth microwear signatures can change very rapidly, yielding information only about an individual's diet in the weeks or even days before its death²². In turn, a recent study in chimpanzees revealed that dental microwear texture can provide information about short-term dietary changes, such as seasonality, but are less effective in providing long-term dietary signals²³. Elgart²⁴ correlated the toughness of food with the amount of dental wear in *Gorilla*, finding some inconclusive results in relation to their diet, potentially due to complications in gathering dietary mechanical property data²⁵. Galbany and colleagues²⁶ measured the percentage of dentine exposure (PDE) in permanent molars of mountain gorillas from the Volcanoes National Park (Rwanda), finding a significant relationship between the degree of tooth wear with the consumption of plant roots at the individual level.

However, we still do not comprehensively understand if ecological constraints and dietary differences observed in gorilla subspecies are reflected in tooth wear. The aim of this study is to explore if differences in diet in gorillas can be detected at the subspecies level by using molar microwear analyses in combination with various dental measurements.
<table><caption>**Table 1.** Summary statistics (median and SD) of relative wear facet areas divided by wear stages.</caption><thead><tr><th rowspan="2"></th><th rowspan="2">N</th><th colspan="2">Buccal P I</th><th colspan="2">Lingual P I</th><th colspan="2">PII</th></tr><tr><th>Median</th><th>SD</th><th>Median</th><th>SD</th><th>Median</th><th>SD</th></tr></thead><tbody><tr><td colspan="8"><strong>Wear 1</strong></td></tr><tr><td>G. b. beringei</td><td>5</td><td>0.04</td><td>0.05</td><td>0.39</td><td>0.07</td><td>0.51</td><td>0.08</td></tr><tr><td>G. b. graueri</td><td>3</td><td>0.11</td><td>0.06</td><td>0.33</td><td>0.05</td><td>0.49</td><td>0.03</td></tr><tr><td>G. g. gorilla</td><td>13</td><td>0.05</td><td>0.04</td><td>0.38</td><td>0.08</td><td>0.48</td><td>0.08</td></tr><tr><td colspan="8"><strong>Wear 2 and 3</strong></td></tr><tr><td>G. b. beringei</td><td>4</td><td>0.06</td><td>0.07</td><td>0.40</td><td>0.04</td><td>0.51</td><td>0.09</td></tr><tr><td>G. b. graueri</td><td>18</td><td>0.11</td><td>0.04</td><td>0.45</td><td>0.04</td><td>0.45</td><td>0.06</td></tr><tr><td>G. g. gorilla</td><td>54</td><td>0.09</td><td>0.06</td><td>0.44</td><td>0.06</td><td>0.47</td><td>0.05</td></tr><tr><td colspan="8"><strong>Wear 4</strong></td></tr><tr><td>G. b. beringei</td><td>1</td><td>0.03</td><td>-</td><td>0.30</td><td>-</td><td>0.67</td><td>-</td></tr><tr><td>G. b. graueri</td><td>4</td><td>0.07</td><td>0.06</td><td>0.53</td><td>0.05</td><td>0.42</td><td>0.03</td></tr><tr><td>G. g. gorilla</td><td>10</td><td>0.08</td><td>0.13</td><td>0.44</td><td>0.06</td><td>0.48</td><td>0.08</td></tr></tbody></table>

Dental macrowear, as opposed to dental microwear, is a cumulative process, which occurs during the individual's lifetime and thus reflects long-term dietary and environmental history<sup>27,28</sup>. Here we examine the macrowear patterns of second lower molars (M2s) in a large sample of western lowland gorilla, Grauer's gorillas and mountain gorillas by combining the Occlusal Fingerprint Analysis method (OFA)<sup>29</sup>, a well-established digital approach tracking changes in dental function by examining occlusal wear facets (planar areas with well-defined edged boundaries produced by the attritional and abrasional contact between upper and lower teeth), with other dental measurements, including the occlusal relief index (OR), PDE and percentage of enamel wear (PEW).

OFA has been used to reconstruct the diet and chewing behaviour of past human populations and in extinct and living non-human primates<sup>28,30-35</sup>. The proportions of crushing, grinding and shearing wear can be used to better discriminate between population with different diets, while tooth wear inclination provides information about the abrasiveness of the diet<sup>31,34,36</sup>. For example, Stuhlträger and colleagues<sup>35</sup> found that the inclination and size of wear facets in two chimpanzee populations were positively correlated with the abrasives in their diet, while a study of molar wear angles of phase II facets revealed significant differences among great ape species, with higher degree of inclination correlated with the consumption of tough and pliant foods such as leaves<sup>28</sup>. Moreover, OFA has also been used to identify unusual types of wear, revealing an unexpected case of bruxism (grinding the teeth and clenching the jaw) in a wild adult male gorilla characterised by the presence of a supernumerary maxillary premolar<sup>33</sup>.

OR has been widely used to reconstruct primate diets, and it is based on the assumption that folivorous species should be characterised by a high relief because it increases the efficacy in processing fibrous and tough foods such as leaves<sup>37</sup>. PDE is a method that can be employed to obtain information about the interaction between dental tissue with various ecological factors, including diet<sup>26,38</sup>, while PEW is a new measurement that aims to quantify the volume of wear in a different way.

Previous molar macrowear studies in Pleistocene and Holocene humans, colobus monkeys, and ungulates, have shown that the more fibrous the diet, the greater the proportion of buccal phase I facets<sup>27,30</sup>, while the consumption of softer foods in primates, such as fruit, is generally associated with the presence of larger lingual phase I facets<sup>27</sup>. Thus, we expect to find larger buccal shearing wear in the folivorous mountain gorillas, and larger lingual shearing wear in the more frugivorous Grauer's and western lowland gorillas. We also expect to find a more complex occlusal morphology (higher OR values) with sharper cusps (higher wear inclinations) in mountain gorillas enhancing comminution of fibrous foods such as leaves and stems. PDE in primates seems to be proportionally correlated with the amount of grit and dust inadvertently ingested with foods<sup>26,38</sup>. Because mountain gorillas consume a larger amount of THV foods, which are covered by soil and dust from the sediments, they should exhibit higher PDE values compared to Grauer's and western lowland gorillas.

We will also test, using PEW, if western lowland gorillas are characterised by the highest degree of wear as suggested by previous studies<sup>24</sup>. Finally, we also explore if there are any differences in molar macrowear between males and females of western lowland gorillas. Gorillas are characterised by a marked sexual dimorphism, with adult males often weighing twice as much as females<sup>39</sup>. This marked body size difference is reflected also in the nutritional requirements, with males spending more time in feeding than females<sup>40</sup>. Because prolonged chewing times have an effect on tooth wear, we expect to find a greater amount of wear in males than in females.

## Results

Overall, gorillas are characterised by small buccal phase I facets, and larger phase II and lingual phase I facets (Table 1). During the earlier wear stages, the macrowear patterns of the three gorilla groups look similar, with differences emerging only at the later wear stages (Fig. 2).

The macrowear pattern of mountain gorillas is dominated by large phase II facets (between 51 and 67%), followed by lingual (between 30 and 40%) and buccal phase I facets (between 3 and 6%). The macrowear pattern of Grauer's gorillas is rather homogenous compared to the other two groups and is characterised by large lingual phase I facets (between 33 and 53%), followed by phase II facets (between 42 and 49%), and by buccal
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Ternary diagrams showing the proportions (in %) of relative wear areas of buccal phase I facets, lingual phase I facets, and phase II facets, which are positioned in an equilateral triangle. Each base of the triangle represents a ratio of 0% while the vertices correspond to a percentage of 100%. Tip crush areas were included in phase II facets. *Gorilla beringei graueri* in green, *Gorilla gorilla gorilla* in orange, and *Gorilla beringei beringei* in blue. Specimens 1 = AMNH 115609, 2 = AMNH 200506, 3 = SMF 59159, 4 = ZMB 47526, and 5 = BM 64.12.1.5.</figcaption></figure>

<table><caption>**Table 2.** Summary statistics (median and SD) of wear facet inclinations divided by wear stages.</caption><thead><tr><th rowspan="2"></th><th rowspan="2">N</th><th colspan="2">Buccal P I</th><th colspan="2">Lingual P I</th><th colspan="2">Phase II</th><th colspan="2">Tip crush</th></tr><tr><th>Median</th><th>SD</th><th>Median</th><th>SD</th><th>Median</th><th>SD</th><th>Median</th><th>SD</th></tr></thead><tbody><tr><td colspan="10"><strong>Wear 1</strong></td></tr><tr><td>G. b. beringei</td><td>5</td><td>32.76</td><td>6.54</td><td>38.03</td><td>4.56</td><td>46.5</td><td>6.01</td><td>16.42</td><td>7.32</td></tr><tr><td>G. b. graueri</td><td>3</td><td>32.00</td><td>7.84</td><td>36.66</td><td>3.90</td><td>41.94</td><td>1.88</td><td>0.00</td><td>12.48</td></tr><tr><td>G. g. gorilla</td><td>13</td><td>29.43</td><td>6.02</td><td>40.55</td><td>4.40</td><td>44.75</td><td>4.70</td><td>0.00</td><td>9.54</td></tr><tr><td colspan="10"><strong>Wear 2</strong></td></tr><tr><td>G. b. beringei</td><td>2</td><td>24.23</td><td>0.35</td><td>31.49</td><td>1.44</td><td>41.98</td><td>1.040</td><td>11.87</td><td>16.78</td></tr><tr><td>G. b. graueri</td><td>9</td><td>35.01</td><td>2.55</td><td>36.74</td><td>5.75</td><td>38.97</td><td>5.08</td><td>11.36</td><td>8.17</td></tr><tr><td>G. g. gorilla</td><td>27</td><td>29.40</td><td>5.17</td><td>38.62</td><td>5.40</td><td>37.4</td><td>4.63</td><td>15.94</td><td>7.37</td></tr><tr><td colspan="10"><strong>Wear 3</strong></td></tr><tr><td>G. b. beringei</td><td>2</td><td>39.58</td><td>10.63</td><td>37.00</td><td>14.26</td><td>32.83</td><td>2.46</td><td>10.06</td><td>3.49</td></tr><tr><td>G. b. graueri</td><td>9</td><td>31.90</td><td>8.96</td><td>35.68</td><td>2.56</td><td>34.83</td><td>4.62</td><td>12.12</td><td>1.73</td></tr><tr><td>G. g. gorilla</td><td>27</td><td>32.73</td><td>7.44</td><td>37.30</td><td>3.72</td><td>32.23</td><td>5.53</td><td>13.13</td><td>5.99</td></tr><tr><td colspan="10"><strong>Wear 4</strong></td></tr><tr><td>G. b. beringei</td><td>1</td><td>36.81</td><td>-</td><td>30.64</td><td>-</td><td>33.40</td><td>-</td><td>7.75</td><td>-</td></tr><tr><td>G. b. graueri</td><td>4</td><td>28.40</td><td>15.49</td><td>37.57</td><td>2.64</td><td>28.70</td><td>1.59</td><td>11.46</td><td>5.50</td></tr><tr><td>G. g. gorilla</td><td>10</td><td>26.64</td><td>8.12</td><td>34.71</td><td>3.42</td><td>29.28</td><td>7.73</td><td>14.49</td><td>4.67</td></tr></tbody></table>

phase facets (between 7 and 11%). Finally, western lowland gorillas show the most variable molar macrowear pattern, with large lingual phase I facets (between 38 and 44%) and phase II facets (between 47 and 48%), and small buccal phase I facets (between 5 and 9%). This was confirmed by the Levene's test of equality of variance among the three gorilla groups (Table S1).

No statistically significant differences were found when comparing molars with wear stage 1 (Table S2 and S3).

Statistically significant differences were found when comparing lingual phase I facets of mountain gorillas with those of Grauer's gorillas (p = 0.040; Wear stage 2 and 3). No differences have been found between molars with wear stage 4. However, the only mountain gorilla with wear 4 included in this analysis looks very different compared to Grauer's and western lowland gorillas (Fig. 2).

Gorilla molars are generally characterised by steep facets, with average angles of phase I and phase facets comprised between 26 and 47 degrees (Table 2). Tip crush areas are relatively flat, with values comprised between 8 and 16 degrees. Phase II facets show much steeper angles if compared to phase I facets, especially in slightly worn molars. The opposite situation is found in more heavily worn molars with phase I facets characterised by higher values if compared to the planes of phase II facets. The mountain gorillas generally display the steepest wear facets, while western lowland and Grauer's gorillas tend to have more variable inclinations depending on the wear stage considered (Table S4 and S5).

We could not statistically test if mountain gorillas facet inclinations differ from those of Grauer's and western lowland gorillas, because we could not group together molars with wear stage 2 with molars with wear stage 3. We could only compare the wear facet angles of molars with wear stage 1, which did not produce any statistically
<table><caption>**Table 3.** Summary statistics (median and SD) of occlusal relief index (OR), percentage of dentine exposure (PDE) and percentage enamel wear (PEW).</caption><thead><tr><th rowspan="2">Groups</th><th rowspan="2">N</th><th colspan="2">OR</th><th colspan="2">PDE</th><th colspan="2">PEW</th></tr><tr><th>Median</th><th>SD</th><th>Median</th><th>SD</th><th>Median</th><th>SD</th></tr></thead><tbody><tr><td>G. b. beringei</td><td>10</td><td>1.895</td><td>0.133</td><td>0.245</td><td>3.156</td><td>20.215</td><td>8.540</td></tr><tr><td>G. b. graueri</td><td>25</td><td>1.990</td><td>0.190</td><td>1.660</td><td>3.550</td><td>24.290</td><td>6.600</td></tr><tr><td>G. g. gorilla</td><td>77</td><td>2.040</td><td>0.268</td><td>0.760</td><td>3.400</td><td>26.830</td><td>7.760</td></tr></tbody></table>

significant difference (Table S6 and S7). The only statistically significant difference has been found between the buccal phase I facets of Grauer's and western lowland gorillas ($p = 0.016$; Wear 2).

Mountain gorillas display lower OR values compared to those of western lowland gorillas (Table 3).

Overall, western lowland gorillas are also characterised by the highest degree of variation in OR, while the least variable are the mountain gorillas. Grauer's gorillas display intermediate OR values. No statistically significant differences have been found in OR values among the three gorilla groups examined in this study (Table S8).

In terms of PDE, Grauer's gorilla show the highest values, followed by western lowland gorillas and mountain gorillas (Table 3). No statistically significant differences have been found when comparing PDE values between the gorillas examined in this study (Table S8).

Western lowland gorillas are generally characterised by the highest PEW values, followed by Grauer's and mountain gorillas (Table 3). We found statistically significant differences in PEW values among the three gorilla groups ($p = 0.016$; Table S8). Finally, no statistically significant differences have been found between males and females of western lowland gorillas for all variables considered in this study (Tables S10–S13).

## Discussion

We have found significant differences in molar macrowear patterns between mountain gorillas and lowland gorillas, confirming our initial hypothesis. Grauer's and western lowland gorilla molar macrowear patterns are dominated by lingual phase I facets, which is probably associated with a larger intake of fruit in their diet. Previous studies of molar macrowear on various ungulate mammals and of colobus monkeys suggested that lingual phase I facets are linked to the consumption of softer food items such as fruit²⁷. However, we found that Grauer's gorillas display the largest lingual phase I facets, especially at the most advanced wear stages. We would have expected to find the largest lingual phase I facets in western lowland gorillas, considering their more frugivorous diet⁴¹. We observe this trend only when examining slight worn molars, while at intermediate levels of wear the proportion of lingual phase I facets in Grauer's and western lowland gorillas is similar. Larger lingual phase I facets are possibly correlated with greater transverse mandibular movements, which may also indicate the consumption of foods other than fruit such as roots, seeds and gums²⁷,³⁰. Mountain gorilla molars are characterised by the largest phase II facets which seem to be associated with a greater consumption of THV foods. Contrary to previous findings, where higher degrees of folivory was linked to larger buccal phase I facets²⁷, here we found that the consumption of herbaceous food is associated with an increase of the proportion of phase II facets and tip crush areas. Buccal phase I facets are very small in all gorillas examined in this study at any given wear stage. This seems to be very common in all great apes, which would indicate a distinctive masticatory behaviour compared to early and late hominins⁴². The functional significance of phase II facets is still highly debated⁴³. Phase II facets are used during both phases of mastication, including crushing (phase I) and grinding (phase II)⁴⁴. Grinding is a combination of forces that act parallel and perpendicular to the contact planes producing a more horizontal shear force⁴⁵. Consequently, the consumption of THV foods, which requires prolonged exposure to repetitive mastication⁴⁶, probably promote the formation of large phase II facets.

As expected, the western lowland gorillas show the highest degree of macrowear variation, suggesting a more flexible diet if compared to those of eastern gorillas. Western lowland gorillas occupy a large geographical area, ranging between the southern part of Cameroon to North Angola⁸. They also inhabit a variety of lowland forest types, including open and closed-canopy forests, and vast swamps in northern Congo⁶. Western lowland gorillas eat a wide variety of fruit, preferring ripe, succulent, sweet fruits that are low in proteins and fat⁴. However, fruit productivity is highly seasonal, and therefore the preferred food of western lowland gorillas is not always available. During these periods of fruit scarcity, western lowland gorillas consume larger quantity of THV, including pith, shoots, leaves and barks⁴¹. Moreover, some populations of western lowland gorillas regularly consume AHV, which are high in essential vitamins and minerals⁴⁸, while evidence of seasonal hard-object feeding has been documented in one population of western lowland gorillas from Loango National Park, Gabon¹¹.

Grauer's gorillas and the mountain gorillas show a less variable macrowear pattern suggesting the intake of a less diverse diet if compared to western lowland gorillas. Mountain gorillas occupy high altitude forests, which are characterised by a significantly lower diversity of plants and fewer fruits⁴⁶. Grauer's gorilla's diet is more variable than those of mountain gorillas. At higher altitude the diet mostly consists of THV foods, like mountain gorillas, while the populations from lower elevations incorporate more fruit in their diets³,⁶.

Mountain gorillas are generally characterised by steeper phase II facets, probably increasing the shearing capabilities, which is functionally advantageous in processing tough and pliant foods such as leaves²⁸. This is very different from the wear patterns observed in late *Homo*, where phase II facets are significantly flatter than buccal and lingual phase I facets³¹. Significant differences in inclination of buccal phase I facets between Grauer's and western lowland gorillas are probably linked to the consumption of THV foods. Although the buccal phase I
facets play a little role in gorilla chewing behaviour, their steeper planes in Grauer's gorillas increase their shearing capacity, more suitable for processing herbaceous vegetation. We observe the same pattern in mountain gorillas, but because of the small sample size we could not test if there were any statistically significant differences with western lowland gorillas.

While the percentage of dentine exposure is negligible in slightly worn molars, we noticed larger PDE values in more advanced wear stages in the mountain and Grauer's gorillas. This could be correlated with the increased amount of THV foods included in their diet. Feeding time in primates increases as toughness increases⁴¹. A previous study has shown that PDE in mountain gorillas from the Virunga National Park was positively correlated with the time spent consuming plant roots²⁶. This was also probably correlated with the amount of exogenous material accidentally introduced in their diet, which is more likely to cause increased tooth wear⁴². It has been observed that mountain gorillas pull the plants out from the soil to access the roots without properly eliminating all the gritty particles that cover the plants²⁶. We have also observed that the molars of western lowland gorillas are generally characterised by a higher occlusal relief index compared to those of mountain and Grauer's gorillas, partially confirming the results of Berthaume's study¹⁶ who found that western lowland gorillas were characterised by sharper cusps compared to mountain and Grauer's gorillas.

Mountain gorillas are characterised by the least amount of wear (PEW) compared to western lowland and Grauer's gorillas, confirming the results of previous studies²⁴,²⁶. This would indicate that mountain gorilla molars maintain a good chewing efficiency throughout the life of an individual despite the mechanically demanding diet²⁴. Similarly, the high PEW values of western lowland gorillas are in agreement with Elgart's study²⁴, where it was found that the western populations showed the highest quantities of wear.

Finally, dental topographic and macrowear analyses of western lowland gorilla molars did not reveal any difference between males and females, rejecting thus our initial hypothesis. Previous studies support these results as it was found that there were no differences in occlusal relief and in the amount of wear between males and females of western lowland gorillas²⁴,³⁷. Previous studies have also identified that western lowland gorilla males in Bai Hokü, Central African Republic, spend less time feeding and more time resting compared to females and immatures in their group⁴¹.

It is important to consider some limitations in our current study. First of all, the sample size of mountain gorillas is rather small, and therefore our interpretations need to be taken with caution. Unfortunately, most of the museum collections we examined included only a small number of mountain gorillas. Their teeth were often heavily worn and chipped, preventing thus their inclusion for the occlusal fingerprint analysis.

Statistically testing the variability across samples of unequal sizes can be problematic. The largest variation found in the western lowland gorillas could be simply driven by the larger sample size, and may not reflect the ecology of this taxon. To overcome this limitation, we have performed additional statistical tests, such as the Permutational Multivariate Analysis of Variance, which confirmed that western lowland gorillas are characterised by the most variable macrowear pattern. We also verified that the variance found in the groups analysed in this study is homogenous, which should have partially minimised this problem.

Because our sample mostly consists of museum specimens, the age of the individual was often unknown, and the geography was frequently uncertain. This additional information could have been very useful with the interpretation of the results. For example, a recent study of dental wear patterns in chimpanzee populations from precise geographic locations revealed previous information for the ecological interpretations of their results²³. In relation to this, more detailed information on how the sample has been acquired could have helped to better interpret the results. Because most of the sample analysed in this study has been wild-shot, and therefore did not live a full life, we should expect a less amount of wear compared to those individuals who have died of natural causes. Moreover, a couple of specimens, caught in the wild, actually died in zoos. A closer look at their macrowear patterns revealed some differences, especially for the specimens SMF 59159 and ZMB 47526, who died at the zoos of Duisburg and Berlin respectively⁴⁷ (Fig. 2). These specimens are characterised by a slightly larger percentage of lingual phase I facets, which may reflect a higher intake of fruit in their diet. The diet of great apes living in captivity has been traditionally based on cultivated fruits and vegetables, pellet concentrate and browse⁴⁸. It would be interesting to examine if macrowear patterns of apes grown in captivity differ from those caught in the wild considering they consume different diets. In addition, future studies could include a larger sample of mountain gorillas, and could also compare the analysis of permanent and deciduous dentition to examine if there are any dietary difference between adult and immature individuals.

# Materials and methods

**Ethical statement.** No live animals were handled in this study.

The sample examined consists of 112 gorilla specimens, including 10 *Gorilla beringei beringei*, 25 *Gorilla beringei graueri*, and 77 *Gorilla gorilla gorilla*. We selected only fully erupted M2s (either left or right, but no anti-meres) because they provide a good general overview of the development of masticatory function in primates⁴⁹. We included only molars with a slight to moderate degree of wear because in heavily worn teeth occlusal facets tend to coalesce often preventing a clear identification³⁰. We qualitatively evaluated the level of wear based on the amount of dentine exposure and cusp removal (wear stages 1–4)³⁰. When available we included from museum records⁴⁷ information about age (subadult and adult), sex, locality and acquisition (electronic supplementary material, Table S14). The sample mostly consists of specimens captured in the wild (wild-shot). Unfortunately, it was difficult to obtain additional information about the way the sample has been collected. We did not include any information about their acquisition for the specimens that either did not have any acquisition date, or that have been acquired after 1950, because it was difficult to establish if the sample has been collected opportunistically or hunted by local people for bushmeat. There are only a handful specimens that died in zoos, but that have
<figure><img src="image_5.png" /><figcaption>**Figure 3.** Second right lower molar of a western lowland gorilla adult male (90,194) showing the macrowear pattern. Facets are numbered according to the Maier and Schneck's labelling system⁵¹. Buccal phase I (facets 1, 2, 2.1, 3 and 4; in blue), lingual phase I (facets 5, 5.1, 6, 7, 8 and 8.1; in green) and phase II (facets 9, 10, 11, 12 and 13; in red). We also identified tip crush areas (in orange) and dentine exposures (in dark grey). P = protoconid, H = hypoconid, Hyp = hypoconulid.</figcaption></figure>

been collected in the wild. Some of these specimens, such as BM 64.12.1.5, may have not survived the journey when they arrived at their destination⁴⁷.

**Occlusal fingerprint analysis (OFA).** Three-dimensional (3D) digital models of teeth were post-processed using Polyworks V12 (InnovMetric Software), a 3D metrology software. OFA consists of four consecutive steps: (1) model orientation, (2) facet identification, (3) facet area, and (4) facet inclination²⁹. The polygonal models are oriented using a reference plane that is created along the cervical line of the tooth through the least square best-fit method. Successively, the reference plane is rotated to the *xy* plane obtained from the original coordinate system.

The identification of wear facets follows the numbering system originally created by Maier and Schneck⁵¹ and later modified by Kullmer and colleagues²⁹. Facets 1, 1.1, 2, 2.1, 3 and 4 develop along the buccal slopes of protoconid, hypoconid and hypoconulid, while facets 5, 6, 7 and 8 form along the buccal slopes of the metaconid and entoconid (Fig. 3). These facets are created during the initial phase of the rhythmic chewing cycle (phase I), and they can be further divided into buccal (facets 1–4) and lingual (facets 5–8) facets⁴⁴. Facets 9, 10, 11, 12 and 13 are formed during the second phase of the rhythmic chewing cycle (phase II), when the mandibular molars move out of occlusion. These facets develop along the lingual slopes of protoconid, hypoconid and hypoconulid.

During the rhythmic chewing cycle, the food bolus is initially processed by shearing, which is generated by forces parallel to the contact plane (phase I), followed by crushing between basins and cusps of molars, where the occlusal force is perpendicular to the contact plane⁴⁴. After the molars reach centric occlusion (phase II), the food bolus is processed by grinding, which is the resulting action of the combination of perpendicular and parallel forces to the contact plane.

We also identify tip crush areas, which are generally formed during puncture-crushing, when food is initially pulped by a series of masticatory cycle where tooth-to-tooth contacts do not occur²⁷. Finally, we identify three additional facets in gorilla molars: facet 5.1, facet 8.1 and facet 10.1. These facets are created by the presence of additional dental traits, which are less common in *Homo*. More specifically, facets 5.1 forms in mandibular molars around the area of cusp 7 (C7), and it occludes with maxillary molars in proximity of the lingual cingulum³². Facet 8.1 forms around cusp 6 (C6) and it occludes with the lingual cingulum of the maxillary molars in correspondence of the mesiolingual aspect of the protocone³³. Facet 10.1 has been previously described in Neanderthal molars, and it is created by the contact between the distolingual slope of the metaconule of the maxillary molar with the mesiolingual slope of the protoconid⁵². Facets 5.1 and 8.1 are grouped with lingual phase I facets, while facet 10.1 is considered a phase II facet.

Once the facets have been identified, we automatically calculate the surface area using the *area* function available in Polyworks V12 (InnovMetric Software). To facilitate the analysis of the occlusal contact areas we grouped phase II facets together with tip crush areas, and divided the phase I facets into buccal and lingual facets³⁰. However, tip crush areas and phase II facets have been kept separated for the analysis of facet inclination, considering they are characterised by significantly different angles³¹. Because larger teeth generally develop larger facets than smaller teeth, we eliminate the size factor by using relative areas only. Relative areas are calculated by dividing the absolute area of a facet with the total wear area (TWA), which is the sum of absolute area of all facets³⁰. The facet inclination is calculated by measuring the angle between the reference, or cervical plane, with the facet plane, which is created by selecting all triangles within its perimeter and by applying the best-fit plane function of Polyworks V12 (InnovMetric Software) (electronic supplementary material, figure S1).
**Occlusal relief index (OR).** The occlusal relief index (OR) was calculated by dividing the 3D area by the two-dimensional (2D) area of the occlusal surface⁵³. The occlusal plane, parallel to the cervical plane, was translated along the y axis until it reached the deepest point of the occlusal surface, known as the central fossa. Next, the digital model was sliced, with respect to the occlusal plane, and the 2D area was calculated. The 3D area was calculated by selecting all triangles of the polygonal models above the occlusal plane (electronic supplementary material, figure S2).

**Percentage of dentine exposure (PDE) and enamel wear (PEW).** To calculate the percent of dentine exposure (PDE) is obtained by dividing the sum of total dentinal areas for each tooth by the 3D occlusal area and then multiplied by 100²⁶. The percentage of enamel wear (PEW) is obtained by dividing the total wear area (TWA) with the 3D occlusal area and then multiplied by 100.

**Statistical analysis.** We employed summary statistical analyses (median and standard deviation, SD) for each variable considered in this study. Because we further divided our sample according to wear stages 1 to 4⁵⁰ the comparative groups became relatively small. As the presence of a small sample size prevents the assumption of a normal distribution, we used the nonparametric Mann-Whitney pairwise test, to test whether two univariate samples are taken from populations with equal medians⁵⁴. For the analysis of facet areas, we grouped together molars characterised by wear stage 2 and wear stage 3 in order to maximize the sample size. We then tested if we could combine these two groups by comparing molars wear stage 2 with molars wear stage 3 of western lowland gorillas without finding any statistically significant result (Table S15). However, for the inclination we kept the molars separated by their wear stages, because the level of wear strongly influences the values of these variables³¹. OR, PDE and PWE were not separated into different wear stage groups.

Difference across groups in the distribution of individual variables expressed as percentages was calculated after running arcsine transformation. This method was preferred over logit transformation because of the presence of zeros that would have otherwise generated negative infinite values. Because of the heavily unbalanced sample sizes, and of the potential violation of parametric assumptions, we only relied on non-parametric tests. Homoscedasticity across groups was ascertained through the Levene's Test of equality of variances using the robust Brown-Forsythe Levene-type procedure (i.e. the median)⁵⁵ and the combination of O'Brien's correction factor and Hines-Hines structural zero removal⁵⁶ contained in the package lawstat in R version 4.1⁵⁷,⁵⁸. After having controlled for homogeneity of variances due to unequal sample sizes across the sample, a Kruskal-Wallis test was used to test for the presence of significant differences in the distribution of wear areas within each wear stage and between the gorilla groups. Significant cases were further explored via a Dunn test with Benjamini-Hochberg (BH) correction to quantify pairwise relationships. The same procedure was used to test for differences in the distribution of OR, PDE, and PEW between geographic groups.

All variables expressing wear facet inclinations were used to calculate a pairwise Euclidean distance. Heteroscedasticity within the sample was then calculated on the distance matrix using the analysis of multivariate homogeneity of group dispersion (function betadisper in the package vegan in R, 999 permutations)⁵⁹, Tukey's Honest Significant Differences test, and permutation test for homogeneity of multivariate dispersion (function permutest in package vegan, *n* = 999). Difference between the gorilla groups within each wear stage class was then calculated through a Permutational Multivariate Analysis of Variance (Permanova; 999 permutations) using the function adonis in the package vegan in R.

Differences between male and female individuals of the western lowland gorilla were tested via two-tailed Mann-Whitney tests for independent sample design, while differences in wear facet area of western lowland gorillas across wear stage classes was ascertained through a pairwise Permanova using the package pairwise Adonis⁶⁰. The summary statistical analyses, ternary diagram and Man-Whitney test were performed using the software PAST v.3.22 (PALaeontological Statistics)⁶¹, while for all other statistical analyses we used the R software⁵⁸.

Received: 8 December 2021; Accepted: 10 May 2022

Published online: 02 June 2022

**References**

1. Fossey, D. & Harcourt, D. H. Feeding ecology of free-ranging mountain gorilla (Gorilla gorilla beringei). In *Primate ecology* (ed. Clutton-Brock, T. H.) 415–447 (Academy Press, New York, 1977).

2. Watts, D. P. Composition and variability of mountain gorilla diets in the central Virungas. *Am. J. Primatol.* **7**, 323–356 (1984).

3. Watts, D. Comparative socio-ecology of gorillas. In *Great Ape Societies* (eds. McGrew, W. C., Marchant, L. F. & Nishida, T.) 16–28 (Cambridge University Press, Cambridge, 1986).

4. Doran, D. M. & McNeilage, A. Gorilla ecology and behavior. *Evol. Anthrop.* **6**, 120–131 (1988).

5. Xue, Y. et al. Mountain gorilla genomes reveal the impact of long-term population decline and inbreeding. *Science* **348**, 242–245 (2015).

6. Mittermeier, R. A., Rylands, A. B. & Wilson, D. E. *Handbook of the mammals of the world*. Primates Vol. 3 (Lynx Edicions, 2013).

7. Groves, C. *Primate Taxonomy* (Smithsonian Institution Press, 2001).

8. Cooksey, K. E. & Morgan, D. B. Gorilla (Gorilla). In *The International encyclopedia of primatology*, Vol. 1 (ed. Fuentes, A.) 472–477 (Wiley, Hoboken, 2017).

9. McFarland, K. L. Ecology of cross river gorillas (*Gorilla gorilla diehli*) on Afi mountain, Cross River State, Nigeria. Ph.D. Dissertation. City University of New York, USA (2007).

10. Rogers, M. E., Maisels, F., Williamson, E. A., Fernandez, M. & Tutin, C. E. G. Gorilla diet in the Lopé Reserve, Gabon: a nutritional analysis. *Oecologia* **84**, 326–339 (1990).
11. van Casteren, A., Wright, E., Kupczik, K. & Robbins, M. M. Unexpected hard-object feeding in Western lowland gorillas. *Am. J. Phys. Anthropol.* **170**, 433–438 (2019).

12. Masi, S., Cipolletta, C. & Robbins, M. M. Western lowland gorillas (*Gorilla gorilla gorilla*) change their activity patterns in response to frugivory. *Am. J. Primatol.* **71**, 91–100 (2009).

13. Yamagiwa, J., Basabose, A. K., Kaleme, K. & Yumoto, T. Diet of grauer's gorillas in montane forest of Kahuzi, Democratic Republic of Congo. *Int. J. Primatol.* **26**, 1345–1373 (2005).

14. Grueter, C. C. et al. Long-term temporal and spatial dynamics of food availability for endangered mountain gorillas in Volcanoes National Park, Rwanda. *Am. J. Primatol.* **75**, 267–280 (2013).

15. Ostrofsky, K. R. & Robbins, M. M. Fruit-feeding and activity patterns of mountain gorillas (*Gorilla beringei beringei*) in Bwindi Impenetrable National Park, Uganda. *Am. J. Phys. Anthropol.* **173**, 3–20 (2020).

16. Berthaume, M. A. Tooth cusp sharpness as a dietary correlate in great apes. *Am. J. Phys. Anthropol.* **153**, 226–235 (2014).

17. King, S. J. et al. Dental senescence in a long-lived primate links infant survival to rainfall. *Proc. Natl. Acad. Sci. USA* **102**, 16579–16583 (2005).

18. Berthaume, M. A. & Schroer, K. Extant ape dental topography and its implications for reconstructing the emergence of early *Homo*. *J. Hum. Evol.* **112**, 15–29 (2017).

19. Sheine, W. S. & Kay, R. F. An analysis of chewed food particle size and its relationship to molar structure in the primates *Cheirogaleus medius* and *Galago senegalensis* and the insectivoran Tupaia glis. *Am. J. Phys. Anthropol.* **47**, 15–20 (1977).

20. Galbany, J., Estebaranz, F., Martinez, L. M. & Pérez-Pérez, A. Buccal dental microwear variability in extant African Hominoida: taxonomy versus ecology. *Primates* **50**, 221–230 (2009).

21. Scott, R. S., Teaford, M. F. & Ungar, P. S. Dental microwear texture and anthropoid diets. *Am. J. Phys. Anthropol.* **147**, 551–579 (2012).

22. Teaford, M. F. & Oyen, O. J. In vivo and in vitro turnover in dental microwear. *Am. J. Phys. Anthropol.* **80**, 447–460 (1989).

23. Stuhlträger, J. et al. Dental wear patterns reveal dietary ecology and season of death in a historical chimpanzee population. *PLoS ONE* **16**, e0251309 (2021).

24. Elgart, A. A. Dental wear, wear rate, and dental disease in the African apes. *Am. J. Primatol.* **72**, 481–491 (2010).

25. Berthaume, M. A. Food mechanical properties and dietary ecology. *Am. J. Phys. Anthropol.* **159**, 79–104 (2016).

26. Galbany, J. et al. Tooth wear and feeding ecology in mountain gorillas from Volcanoes National Park, Rwanda. *Am. J. Phys. Anthropol.* **159**, 457–465 (2016).

27. Janis, C. M. The correlation between diet and dental wear in herbivorous mammals, and its relationship to the determination of diets of extinct species, in *Evolutionary paleobiology of behavior and coevolution* (ed. Boucot, A. J.) 241–259 (Elsevier, Amsterdam, 1990).

28. Knight-Sadler, J. & Fiorenza, L. Tooth wear inclination in great ape molars. *Folia Primatol.* **88**, 223–236 (2017).

29. Kullmer, O. et al. Technical note: Occlusal fingerprint analysis: Quantification of tooth wear pattern. *Am. J. Phys. Anthropol.* **139**, 600–605 (2009).

30. Fiorenza, L. et al. Molar macrowear reveals Neanderthal eco-geographical dietary variation. *PLoS ONE* **6**, e14769 (2011).

31. Fiorenza, L., Benazzi, S., Oxilia, G. & Kullmer, O. Functional relationship between dental macrowear and diet in Late Pleistocene and recent modern human populations. *Int. J. Osteoarchaeol.* **28**, 153–161 (2018).

32. Fiorenza, L. et al. The functional role of the Carabelli trait in early and late hominins. *J. Hum. Evol.* **145**, 102816 (2020).

33. Fiorenza, L. & Kullmer, O. Occlusion in an adult male gorilla with a supernumerary maxillary premolar. *Int. J. Primatol.* **37**, 762–777 (2016).

34. Kullmer, O., Menz, U., & Fiorenza, L. Occlusal fingerprint analysis (OFA) reveal dental occlusal behaviour in primate teeth. In T. Martin & W. von Koenigswald (Eds.), *T. Martin, W von Koenigswald, K-H Südekum), Mammalian teeth: form and function*. (pp. 25–43). Munich, Germany: Dr. F. Pfeil (2020).

35. Stuhlträger, J. et al. Dental wear patterns reveal dietary ecology and season of death in a historical chimpanzee population. *PLoS ONE* **16**, e0251309 (2021).

36. Fiorenza, L., Benazzi, S., Estalrrich, A., & Kullmer, O. Diet and cultural diversity in Neanderthals and modern humans from dental macrowear analyses. In C. Schmidt & J. T. Watson (Eds.), *Dental wear in evolutionary and biocultural contexts* (pp. 39–72). London, UK: Academic Press (2020).

37. M'Kirera, F. & Ungar, P. S. Occlusal relief changes with molar wear in *Pan troglodytes troglodytes* and *Gorilla gorilla gorilla*. *Am. J. Primatol.* **60**, 31–41 (2003).

38. Galbany, J. et al. Age-related tooth wear differs between forest and savanna primates. *PLoS ONE* **9**, e94938 (2014).

39. Leigh, S. R. & Shea, B. T. Ontogeny and the evolution of adult body size dimorphism in apes. *Am. J. Primatol.* **36**, 37–60 (1995).

40. Watts, D. P. Environmental influences on mountain gorilla time budgets. *J. Am. Primatol.* **15**, 195–211 (1988).

41. Doran, D. M. et al. Western lowland gorilla diet and resource availability: New evidence, cross-site comparisons, and reflections on indirect sampling methods. *Am. J. Primatol.* **58**, 91–116 (2002).

42. Zanolli, C. et al. Evidence of increased hominid diversity in the Early and Middle Pleistocene of Indonesia. *Nat. Ecol. Evol.* **3**, 755–764 (2019).

43. Krueger, K. L., Scott, J. R., Kay, R. F. & Ungar, P. S. Dental microwear textures of “phase I” and “phase II” facets. *Am. J. Phys. Anthropol.* **137**, 485–490 (2008).

44. Kay, R. F. & Hiiemae, K. M. Jaw movement and tooth use in recent and fossil primates. *Am. J. Phys. Anthropol.* **40**, 227–256 (1974).

45. Wall, C. E., Vinyard, C. J., Johnson, K. R., Williams, S. H. & Hylander, W. L. Phase II jaw movements and masseter muscle activity during chewing in *Papio anubis*. *Am. J. Phys. Anthropol.* **129**, 215–224 (2006).

46. Glowacka, H. et al. Toughness of the Virunga mountain gorilla (*Gorilla beringei beringei*) diet across an altitudinal gradient. *Am. J. Primatol.* **79**, e22661 (2017).

47. Cooper, J. E. & Hull, G. *Gorilla pathology and health* (Academic Press, 2017).

48. Hammerton, R., Hunt, K. A. & Riley, L. M. An investigation into keeper opinions of great apes diet and abnormal behaviour. *J. Zoo Aquar. Res.* **7**, 170–178 (2019).

49. Kay, R. F. *Mastication, molar tooth structure and diet in primates*. Ph.D. thesis, Yale University, New Haven, CT (1973).

50. Smith, B. H. Patterns of molar wear in hunter-gatherers and agriculturalists. *Am. J. Phys. Anthropol.* **63**, 39–56 (1984).

51. Maier, W. & Schneck, G. Konstruktionsmorphologische Untersuchungen am Gebiß der hominoiden Primaten. *Z. Morphol. Anthropol.* **72**, 127–169 (1981).

52. Fiorenza, L., Benazzi, S., Tausch, J., Kullmer, O. & Schrenk, F. Identification reassessment of the isolated tooth Krapina D58 through Occlusal Fingerprint Analysis. *Am. J. Phys. Anthropol.* **143**, 306–312 (2010).

53. Kullmer, O., Huck, M., Engel, K., Schrenk, F. & Bromage, T. Hominid Tooth Pattern Database (HOTPAD) derived from optical 3D topometry. In *Three-dimensional imaging in paleoanthropology and prehistoric archaeology* (eds. Mafart, B. & Delingette, H.) 71–82 (Acts of the XIVth UISPP Congress, BAR Int. Ser.1049, 2002).

54. Hammer, Ø. & Harper, D. *Paleontological data analysis* (Blackwell Publishing, 2006).

55. Brown, M. B. & Forsythe, A. B. Robust tests for the equality of variances. *J. Am. Stat. Assoc.* **69**, 364–367 (1974).

56. Noguchi, K. & Gel, Y. R. Combination of Levene-type tests and a finite-intersection method for testing equality of variances against ordered alternatives. *J. Nonparam. Stat.* **22**, 897–913 (2010).
57. Gastwirth, J. L., et al. Lawstat: tools for biostatistics, public policy, and law. R package version 3.4. https://CRAN.R-project.org/package=lawstat (2020).

58. R Core Team. R: A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/. (2021)

59. Oksanen, J. et al. Vegan: Community Ecology Package. R package version 2.5-7. https://CRAN.R-project.org/package=vegan (2020).

60. Martinez Arbizu, P. PairwiseAdonis: pairwise multilevel comparison using adonis. R package version 0.4 (2017).

61. Hammer, Ø., Harper, D. A. T. & Ryan, P. D. PAST: Palaeontological statistics software package for education and data analysis. *Palaeontol. Electron.* **4**, 9 (2001).

## Acknowledgements

We would like to thank curators and institutions that gave access to the great ape specimens: Irina Ruf and Katrin Krohmann (Senckenberg Research Institute and Natural History Museum Frankfurt, Germany), Richard Sabin and Roberto Portela Miguez (Natural History Museum, London, UK), Jenkins and Eileen Westwig of the American Museum of Natural History (New York, New York), Yuzuru Hamada from the Primate Research Institute, Kyoto University (Inuyama, Japan), Christiane Funk and Steffen Bock (Museum für Naturkunde, Berlin, Germany) and Emmanuel Gilissen (Royal Museum for Central Africa, Tervuren, Belgium). We would also like to thank Jordan Knight-Sadler, David Rex Mitchell, Ada Klinkhamer, Zeinab Alimohamadi, Arnaud Mazurier (IC2MP, Poitiers, France) and Adélaïde Euriat (PALEVOPRIM, Poitiers, France) for generating the 3D digital surface models of teeth used in this study. This study was supported by the Australian Research Council (Grant Number: DP190100465), by the French National Research Agency (Grant Number: ANR-17-CE02-0010-01), by the Primate Research Institute Cooperative Research Program (Grant Number: 2016-B-91), and by the Spanish Ministerio de Ciencia e Innovación (Grant Number: PID2020-114517GB-I00).

## Author contributions

T.H. and L.F wrote the main manuscript text and prepared all the figures and tables. E.B. prepared the additional statistical analyses. M.A.B, O.K., F.G., V.L., J.G. and A.R. provided additional samples. All authors reviewed the manuscript.

## Competing interests

The authors declare no competing interests.

## Additional information

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-022-12488-8.

**Correspondence** and requests for materials should be addressed to L.F.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2022