# Possible Interbreeding in Late Italian Neanderthals?
## New Data from the Mezzena Jaw (Monti Lessini, Verona, Italy)

Silvana Consemi, Aurélien Mounier, Paolo Giunti, Martina Lari, David Caramelli, Laura Longo

► To cite this version:

Silvana Consemi, Aurélien Mounier, Paolo Giunti, Martina Lari, David Caramelli, et al.. Possible Interbreeding in Late Italian Neanderthals? New Data from the Mezzena Jaw (Monti Lessini, Verona, Italy). PLoS ONE, 2013, 8 (3), pp.e59781. 10.1371/journal.pone.0059781. halshs-00849140

HAL Id: halshs-00849140

https://shs.hal.science/halshs-00849140v1

Submitted on 19 Sep 2018

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 CC BY 4.0 - Attribution - International License
# Possible Interbreeding in Late Italian Neanderthals? New Data from the Mezzena Jaw (Monti Lessini, Verona, Italy)

Silvana Condemi¹*, Aurélien Mounier¹,², Paolo Giunti³, Martina Lari⁴, David Caramelli⁴, Laura Longo⁵*

1 UMR 7268 CNRS/Aix-Marseille Université/EFS ADES - Anthropologie bioculturelle, Droit, Ethique et Santé Faculté de Médecine - Secteur Nord Aix-Marseille Université, Marseille, France, 2 The Leverhulme Centre for Human Evolutionary Studies Biological, Anthropology Division, Department of Archaeology and Anthropology, University of Cambridge, Cambridge, United Kingdom, 3 Istituto Italiano di Preistoria e Protostoria, Firenze, Italy, 4 Università di Firenze, Dipartimento di Biologia Evoluzionistica, Laboratorio di Antropologia, Unità di Antropologia Molecolare/Paleogenetica, Firenze, Italy, 5 Musei Civici Fiorentini, Firenze, Italy

## Abstract

In this article we examine the mandible of Riparo Mezzena a Middle Paleolithic rockshelter in the Monti Lessini (NE Italy, Verona) found in 1957 in association with Charentian Mousterian lithic assemblages. Mitochondrial DNA analysis performed on this jaw and on other cranial fragments found at the same stratigraphic level has led to the identification of the only genetically typed Neanderthal of the Italian peninsula and has confirmed through direct dating that it belongs to a late Neanderthal. Our aim here is to re-evaluate the taxonomic affinities of the Mezzena mandible in a wide comparative framework using both comparative morphology and geometric morphometrics. The comparative sample includes mid-Pleistocene fossils, Neanderthals and anatomically modern humans. This study of the Mezzena jaw shows that the chin region is similar to that of other late Neanderthals which display a much more modern morphology with an incipient mental trigone (e.g. Spy 1, La Ferrassie, Saint-Césaire). In our view, this change in morphology among late Neanderthals supports the hypothesis of anatomical change of late Neanderthals and the hypothesis of a certain degree of interbreeding with AMHs that, as the dating shows, was already present in the European territory. Our observations on the chin of the Mezzena mandible lead us to support a non abrupt phylogenetic transition for this period in Europe.

**Citation:** Condemi S, Mounier A, Giunti P, Lari M, Caramelli D, et al. (2013) Possible Interbreeding in Late Italian Neanderthals? New Data from the Mezzena Jaw (Monti Lessini, Verona, Italy). PLoS ONE 8(3): e59781. doi:10.1371/journal.pone.0059781

**Editor:** David Frayer, University of Kansas, United States of America

**Received** October 30, 2012; **Accepted** February 18, 2013; **Published** March 27, 2013

**Copyright:** © 2013 Condemi et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

**Funding:** The authors have no support or funding to report.

**Competing Interests:** Please note that in accord with PLOS ONE's guidelines, the fact that one of the authors of this article, David Caramelli is a PLOS ONE Editorial Board member does not alter the authors' adherence to all the PLOS ONE policies on sharing data and materials.

* E-mail: silvana.condemi@univ-amu.fr (SC); laura.longo@comune.fi.it (LL)

## Introduction

One of the most frequently debated questions in paleoanthro-pology concerns the arrival in Europe of our species, *Homo sapiens*, anatomically modern humans (AMHs), and the fate of the humans who lived in this territory prior to their arrival, the Neanderthals.

For many decades there were two different responses to this question: according to one point of view the Neanderthals did not really disappear, but were incorporated into the new incoming modern humans. Authors who support this hypothesis have argued that there was a change in morphology of late Neander-thals [1–3] and have interpreted certain anatomic features observed among early AMHs in Europe as a result of a continuity with Neanderthals [4,5]. A certain degree of continuity has also been supported by archaeologists who have identified much more complex forms of behavior among Neanderthals than was previously acknowledged [6,7]. European Neanderthals have been considered to be not only the producers of Mousterian assemblages but also the makers of the later so-called "transitional assemblages" (Chatelperronian, Uluzzian, Bohunician, Szeletian) [8,9], either by internal modification [7,9,10] or through acculturation by AMHs [11,12].

An opposing model has claimed that there is great discontinuity between Neanderthals and modern humans [13,14] and relates the demise of the Neanderthals to the territorial expansion of AMHs from Africa through the Near East.

The scenarios which have generally been accepted argue that this territorial expansion occurred during a period of great climatic change [15,16]. According to this hypothesis, the expansion of AMHs, identified primarily through their association with Aurignacian assemblages [13,17], pushed the Neanderthals associated with Mousterian assemblages toward southern Europe and, in particular, toward the Iberian and Italian peninsulas in the Mediterranean area [11,18]. This view was reinforced by genetic data which have shown that there is no contribution of Neanderthals to the mitochondrial DNA of *H. sapiens* [19,20].

During recent years, data collected in Europe that seemed to support this view have been questioned. First, Neanderthal nuclear DNA shows a low level of interbreeding (4%) with sapiens [21]. Furthermore, *H. sapiens* is now associated with local (Uluzzian) so-called "transitional assemblages" at Grotta del Cavallo in the southern Italian peninsula while the human remains were previously thought to be Neanderthals [9,22]. The presence of AMHs in Grotta del Cavallo has been demonstrated based on the morphological pattern of the enamel on human deciduous teeth, and the age of Uluzzian artifacts associated with the teeth (levels of unit E) has been re-analysed [23]. The new dating shows that the AMHs reached the southern Italian peninsula at around 45–43 ka BP, which is at least 7000 years earlier than was previously supposed. This study indicated the difficulty of advancing a general explanation [11] valid for all of Europe, since the
replacement of Neanderthals by AMHs on the Italian peninsula took place earlier [23,24] and was probably different than that which occurred in Iberia [25–27].

In this article we examine the morphology of the Mezzena mandible (Figure 1) found in 1957 [28]. We argue that the mandibular morphology of late Italian Neanderthals, in particular the chin, may help us better understand the transition between the two human groups. The study of the human remains of Middle Paleolithic Riparo Mezzena, a rockshelter in the Monti Lessini (Venetian region -NE Italy) associated with Charentian Moustérian lithic assemblages [29,30] has led to the identification of the only genetically typed Neanderthal of the Italian peninsula (cf. [31–33] and this study) and has confirmed through dating that it belongs to a late Neanderthal (i.e. $34.5 \pm 655$ ka) [30]. Our aim is to re-evaluate the taxonomic affinities of the Mezzena jaw in a wide comparative framework using both comparative morphology and geometric morphometrics analyses. The comparative sample includes mid-Pleistocene fossils, Neanderthals and anatomically modern humans (cf. [34], Tables 1 and S1). This study of the Mezzena mandible shows that the chin region is similar to that of other late Neanderthals which display a much more modern morphology with an incipient mental trigone (e.g. Spy 1, Saint Césaire). In our view, this change in morphology among late Neanderthals reopens the debate on the "more modern like" morphology of late Neanderthals and can lend support to the hypothesis of a certain degree of continuity with AMHs or a possible interbreeding with them.

## Results

### Comparative Morphology

The mandible of Mezzena (IG VR 203334) is incomplete (Figure 1). All of the fractures are old, the two vertical branches and the left side of the mandibular corpus (from P1) are broken. However, the symphyseal region is complete. On its right side the body of the mandible is conserved up to the level of the second molar. No teeth are present in the mandible: most of them were lost *ante-mortem* (Figure 1D). Destruction and pathological remodeling of internal bone in the vicinity of the right premolar and molar teeth was revealed through x-ray and computerized-tomodensitometric examinations. This lesion has been interpreted as subsequent to an infection due to bacterial invasion developed from the right premolar [35].

The Mezzena corpus is somewhat robust. However, this robustness cannot be evaluated with great precision since the alveolar rim is damaged throughout its length, in other words from the level of the first left premolar until the second right molar (Figure 1D). In spite of this damage we have evaluated the robusticity of the mandible. The remaining height is 25 mm at level of the symphysis, 27.4 mm at the level of the second molar (M2) and, where the alveolar rim is less damaged, at the level of the mental foramen, the height is 34 mm. (Table 1). Thus, even if not completely accurate due to damage to the alveolar rim, the height of Mezzena at the mental foramina and at the M2 is situated within the range of the variation recognized for European Neanderthals (mental foramen: $v = 36–24$ mm, N = 18; second molar: $v = 28–33,5$ mm, N = 15). The thickness of the symphysis measures 14 mm, whereas it is 13 mm at the level of the mental foramen and 16.5 mm at the level of the M2. Thus, the thickness of the mandibular corpus of Mezzena is also close to those of European Neanderthals (mental foramen: $v = 14–16.2$ mm; N = 15; second molar: $v = 12.7–19$ mm; N = 15). The index of robustness (i.e. RI) at the level of the mental foramen is 38.23. Due to damage of the alveolar rim in Mezzena, this index is situated

slightly below the lower range of variation of European Neanderthals ($v = 60.4–39.37$ mm; N = 18; cf. Table 1 and [36]). It is of particular interest to note that of all the European Neanderthal jaws, the Mezzena mandible index of robustness is closest to that of St Césaire (RI = 39.37) and Guattari III (RI = 40) and not far from Spy (RI = 42.4) and Arcy sur Cure (RI = 42.6) (Table 1). The index of robustness at the level of M2 is 60. It is situated in the upper range of variation of classical European Neanderthals ($v = 73.9–44$; N = 18; cf. Table 1). It should be noted that variation of index of robustness among modern humans is very wide and can include the Neanderthal variation [36].

The external face of the right side of the mandibular corpus is present up to the alveolus of the second molar (Figure 1C). The principal mental foramen is small (with an opening measuring 3.7 mm) and positioned under the second molar. A smaller secondary mental foramen is found below the second premolar and the first molar. This feature is important since it is generally acknowledged that a mental foramen positioned under the first molar or between the first molar and the second premolars is a typically Neanderthal character. Indeed, in 25% of the European Neanderthals the mental foramina are situated between the second premolar and the first molar, in 65% they are located below the first molar and, in the 10%, below the second premolar [36–38]. The position of this feature in the Mezzena mandible suggests that it cannot be excluded from the Neanderthals and that the Mezzena jaw was moderately elongated similar to Guattari III, Saint-Césaire, Zafarraya, and Palomas 59 [36,39,40]. Additionally, the mental foramina of Mezzena are situated half way up the body of the bone which is usually considered to be an archaic feature [36].

On the lateral surface of the mandibular body of Mezzena (Figure 1C) there is a slight swelling, the *prominentia lateralis*, situated below the location of the second molar and the bone fracture. Near the base of the jaw there is a marginal anterior tubercle (*tuberculum marginalis anterior*). Its size is relatively small. As for the position of the *prominentia lateralis*, the presence of this tubercle is considered to be a diagnostic feature of Neanderthals [36–38].

The interior face of the mandibular body of Mezzena (Figure 1B) displays a well-defined oblique internal line or *linea mylohyoidea*, as on Neanderthals. This line is situated, as on Neanderthals, in a lower position than in modern humans. Toward the front and the upper part there is a small *fossa sublingualis*. Due to the state of preservation of the mandible it is not possible to follow this line throughout its length. Thus it is only possible to note toward the back the presence of a *fossa subalveolaris posterior*. The clear relief of the oblique internal line shows that the mandible of Mezzena had powerful mylo-hyoid muscles. In general, this region, as with other parts of the mandible, resembles the morphology found on Neanderthals.

The symphysal region of the Mezzena mandible (i.e. the region delimited by the mesial rim of the canines, Figure 1A) is relatively well-preserved and provides important information. However, the fractures of the alveolar rim do not make it possible to calculate the angles in order to evaluate the inclination of the symphysis. The bone surface displays a swelling in the region of the *trigonum mentale* which is composed of a very small *tuber symphyseo* and of two *tubercula lateralia*, the latter of which are distinctly separated from the inferior margin of the mandible. This surface does not display an *incurvatio mandibulae*. The *incisura submentalis* is present and, in the basal region, it forms a slight concavity in the shape of an arc with a maximum height of 5 mm. In lateral view the symphyseal region does not appear to be concave as among modern humans, nor convex, as among ancient European fossils, but vertical with a slight swelling. Both morphologies of the lateral profile of the
<figure><img src="image_2.png" /><figcaption>**Figure 1.** The mandible from Mezzena. Frontal view : A, internal view : B, lateral view: C, superior view : D), inferior view : E.
doi:10.1371/journal.pone.0059781.g001</figcaption></figure>

symphysis and of the *incisura* are similar to that found among classic Neanderthals, such as Guattari III or Regourdou from France. But this morphology is present in particular among late Neanderthals (e.g. St Césaire, Spy 1, La Ferrassie 1, Las Palomas 59 and Vindija [5,36]) and to a lesser extent among Neanderthals of the Near East (especially Tabūn II and Amud 1). All these fossils have an incipient *mentum osseum*. On the internal face of the symphysis (Figure 1B), the alveolar rim is severely damaged but the mental spines (*spinae mentales*) can be noticed below the fracture. They are clearly separated as on the Neanderthal La Ferrassie 1. Above these two spines a *foramen spinosum* is clearly visible. Under the upper mental spine there is a very slight half-moon shaped notch comprising the *fossa genioglossa*.

The inferior margin of the mandible (Figure 1D) is very thick and presents visible digastric muscles imprints. They form two digastric fossae which are well delimited and distinct. They are large and ellipsed shaped. At the point of junction of the two fossae there is a marked *crista intergastrica* in front of which is the *trigonum basale* of Toldt. Here too, this region bears a similarity to the Neanderthals.

## Geometric Morphometrics (Shape Analysis)

The M Box test results (M = 207.445, F = 1.150, ddl1 = 110, ddl2 = 3741.902, p = 0.139, Table S3) indicates that the covariance matrices are homogenous, and therefore a linear Discriminant Function Analysis (DFA) is appropriate.

The first discriminant function (F1) of the DFA accounts for 81.0% of the total variance of the discrimination; it separates the three pre-defined groups: *H. sapiens*, *H. neanderthalensis* and mid-Pleistocene specimens which have been previously attributed to *H. heidelbergensis* [35]. This discrimination is supported by a significant Wilks' lambda value (Wilk's λ = 0.097, chi-square = 80.504, df = 20, p<0.0001, Table S4). Most of the intra-group shape variation is represented along the second function (F2:19.0% of variance), which, coherently Wilk's lambda value is less significant (Wilk's λ = 0.504, chi-square = 23.616, df = 9, p = 0.005, Table S4). Nevertheless, these results suggest that the variables can be used to distinguish between the groups. Results of the validation procedure (i.e. cross validation, see Method section) indicate that 78.6% of cross validated grouped specimens were correctly re-classified (i.e. 80% of the modern Humans, 73.3% of the Neanderthals and 83.3% of *H. heidelbergensis*) compare to 95.2%
<table><caption>**Table 1.** Specimens of the study.</caption><thead><tr><th>Specimens</th><th>Hgt mf</th><th>Thc mf</th><th>RI mf</th><th>Hgt M2</th><th>Thc M2</th><th>RI M2</th><th>Labels in Figure 2</th></tr></thead><tbody><tr><td colspan="8">Middle Pleistocene specimens</td></tr><tr><td colspan="8">H. heidelbergensis</td></tr><tr><td>Tighenif 1</td><td>35.5</td><td>18.1</td><td>51</td><td>35.4</td><td>22.2</td><td>62.6</td><td>T1</td></tr><tr><td>Tighenif 2</td><td>31</td><td>15</td><td>48.4</td><td>33</td><td>20</td><td>60.6</td><td>T2</td></tr><tr><td>Tighenif 3</td><td>32.5</td><td>19</td><td>58.5</td><td>-</td><td>22.2</td><td>-</td><td>T3</td></tr><tr><td>Mauer</td><td>33</td><td>19.4</td><td>58.8</td><td>34</td><td>20.5</td><td>60.2</td><td>Ma</td></tr><tr><td>Montmaurin</td><td>28.8</td><td>15</td><td>52.1</td><td>28</td><td>159</td><td>56.9</td><td>Mt</td></tr><tr><td>Arago II</td><td>33.1</td><td>16</td><td>53.1</td><td>30.5</td><td>16.5</td><td>54</td><td>ArII</td></tr><tr><td>Arago XIII</td><td>31.2</td><td>22</td><td>68.7</td><td>29</td><td>23</td><td>79.3</td><td>ArXIII</td></tr><tr><td>Sima de los Huesos (ATB- 1)</td><td>30.2</td><td>16.9</td><td>55.9</td><td>28.9</td><td>17.7</td><td>61.2</td><td>ATB-1</td></tr><tr><td>Sima de los Huesos (AT-75)</td><td>29.1</td><td>15.6</td><td>53.6</td><td>28.7</td><td>16.7</td><td>58.1</td><td>*</td></tr><tr><td>Sima de los Huesos (AT- 250)</td><td>31</td><td>16.9</td><td>54.5</td><td>29.5</td><td>17.9</td><td>60.6</td><td>*</td></tr><tr><td>Sima de los Huesos (AT-300)</td><td>34.3</td><td>17.1</td><td>49.8</td><td>30.1</td><td>18.0</td><td>59.8</td><td>*</td></tr><tr><td>Sima de los Huesos (AT-605)</td><td>37.1</td><td>16.7</td><td>45.0</td><td>32.9</td><td>18.4</td><td>55.9</td><td>AT-605</td></tr><tr><td>Sima de los Huesos (AT-607)</td><td>27.1</td><td>15.2</td><td>56.0</td><td>27.2</td><td>17.1</td><td>62.8</td><td>AT-607</td></tr><tr><td>KNM-BK 67</td><td>32.5</td><td>14.9</td><td>45.9</td><td>32.3</td><td>20.4</td><td>63.1</td><td>BK67</td></tr><tr><td>Ehringsdorf F</td><td>25.5</td><td>16.5</td><td>63.5</td><td>27</td><td>16</td><td>59.3</td><td>EhF</td></tr><tr><td colspan="8">Late Pleistocene</td></tr><tr><td colspan="8">H. neanderthalensis</td></tr><tr><td>Bañolas</td><td>28.9</td><td>16</td><td>55.5</td><td>33</td><td>19</td><td>57.5</td><td>Ba</td></tr><tr><td>La Naulette</td><td>25</td><td>14</td><td>56</td><td>23</td><td>17</td><td>73.9</td><td>Na</td></tr><tr><td>Malarnaud</td><td>24</td><td>14</td><td>60.4</td><td>22</td><td>15</td><td>68.2</td><td>*</td></tr><tr><td>La Chaise BD1</td><td>33</td><td>15.1</td><td>45.7</td><td>30.5</td><td>15</td><td>50</td><td>*</td></tr><tr><td>Krapina G</td><td>30</td><td>15</td><td>50</td><td>28</td><td>14.5</td><td>51.8</td><td>*</td></tr><tr><td>Krapina H</td><td>35</td><td>15</td><td>42.8</td><td>33.5</td><td>15</td><td>44.8</td><td>*</td></tr><tr><td>Krapina J</td><td>33.5</td><td>16</td><td>47.7</td><td>33.2</td><td>16.2</td><td>48.8</td><td>KJ</td></tr><tr><td>Krapina D</td><td>27</td><td>13</td><td>44.4</td><td>-</td><td>-</td><td>-</td><td>*</td></tr><tr><td>Regourdou</td><td>33</td><td>16</td><td>48.5</td><td>32</td><td>15</td><td>46.8</td><td>Reg</td></tr><tr><td>La Quina H5</td><td>34</td><td>15</td><td>44.1</td><td>34</td><td>16</td><td>47</td><td>*</td></tr><tr><td>La Quina H9</td><td>37</td><td>16</td><td>43.2</td><td>-</td><td>-</td><td>-</td><td>QH5</td></tr><tr><td>Arcy-sur-Cure</td><td>38</td><td>16.2</td><td>42.6</td><td>33</td><td>19</td><td>57.6</td><td>*</td></tr><tr><td>La Ferrassie 1</td><td>33</td><td>15</td><td>45.4</td><td>32</td><td>14</td><td>43.7</td><td>LF1</td></tr><tr><td>Spy I</td><td>33</td><td>14</td><td>42.4</td><td>31</td><td>16</td><td>51.6</td><td>Spy1</td></tr><tr><td>Zafarraya</td><td>33.3</td><td>16</td><td>47.9</td><td>32.3</td><td>14</td><td>43.2</td><td>Zaf</td></tr><tr><td>Saint Césaire</td><td>32</td><td>12.6</td><td>39.37</td><td>28.8</td><td>12.7</td><td>44.09</td><td>StC</td></tr><tr><td>Guattari II</td><td>36</td><td>17</td><td>47.2</td><td>-</td><td>-</td><td>-</td><td>Gull</td></tr><tr><td>Guattari III</td><td>35</td><td>14</td><td>40</td><td>35</td><td>16</td><td>45.7</td><td>GullII</td></tr><tr><td>Mezzena</td><td>(34)</td><td>13</td><td>(38.3)</td><td>27.5</td><td>16.5</td><td>60</td><td>Mezzena</td></tr><tr><td>Tabün C1</td><td>27.5</td><td>15</td><td>54.5</td><td>26.2</td><td>16.2</td><td>58</td><td>TC1</td></tr><tr><td>Tabün II</td><td>42.5</td><td>16.4</td><td>55.5</td><td>38.5</td><td>18</td><td>58</td><td>TII</td></tr><tr><td>Amud 1</td><td>34</td><td>15</td><td>44.1</td><td>34</td><td>16.5</td><td>48.5</td><td>*</td></tr><tr><td colspan="8">H. sapiens</td></tr><tr><td>Skhül V</td><td>(36)</td><td>13.2</td><td>36.7</td><td>(34.5)</td><td>13</td><td>37.7</td><td>SV</td></tr><tr><td>Qafzeh 9</td><td>35</td><td>16.6</td><td>47.4</td><td>29</td><td>17</td><td>58.6</td><td>Q9</td></tr><tr><td>Ohalo II</td><td>30.33</td><td>12</td><td>39.7</td><td>28.54</td><td>14.8</td><td>51.8</td><td>OII</td></tr><tr><td>Cro-Magnon I</td><td>28</td><td>13.5</td><td>48.2</td><td>26.3</td><td>13.3</td><td>50.6</td><td>CMI</td></tr><tr><td>Abri Pataud 1</td><td>27.9</td><td>11.7</td><td>41.9</td><td>28.1</td><td>17.2</td><td>61.2</td><td>AP1</td></tr></tbody></table>

Specimens used in the study. * indicates specimens that were not included in the geometrics morphometric analysis due to the state of preservation of the fossils or to their inaccessibility. Bold types indicate when original fossil was examined. Robustness index (RI) of the mandibular corpus are calculated from the thicknesses (Thc) and the heights (Hgt) at the mental foramen level (mf) and at the M2 level. All measurements have been taken in accord with standard procedures defined by [64] and are derived from [36], with the exception of Zafarraya [39] and the Mezzena mandible (italic, present study).

doi:10.1371/journal.pone.0059781.t001
of correct classification during the original DFA procedure (Table S7 and S8). Again, this value suggests that the groups can be distinguished by the DFA.

The linear regressions (PC1: $R^2=0.190, p=0.004$; PC2: $R^2=0.007, p=0.591$; PC3: $R^2=0.005, p=0.653$; PC4: $R^2=0.203, p=0.002$; PC5: $R^2=0.016, p=0.425$; PC6: $R^2=0.006, p=0.630$; PC7: $R^2=0.013, p=0.468$; PC8: $R^2=0.005, p=0.657$; PC9: $R^2=0.047, p=0.164$; PC10: $R^2=0.013, p=0.459$) indicate that centroid size does not significantly impact specimens' shape (see, Table S6). Thus, differences in shape between specimens are not due to allometry.

The F1 is responsible for most of the dispersion of the cloud of points (Figure 2A). Neanderthals and AMHs groups overlap at the centre of the chart. Late *H. neanderthalensis* specimens (i.e. Saint-Césaire and Spy 1), Near-Eastern Neanderthals (i.e. Tabūn II and Amud 1), as well as the more classic Neanderthal specimens La Ferrassie 1 and Guattari III, tend to be positioned at the left margin of the Neanderthals cloud of points in the overlapping area with AMHs *H. sapiens* fossil specimens are positioned well within of the recent human cloud of points and are not similar in shape to Neanderthals to the exception of Abri Pataud 1. On the contrary, it is two 20th century specimens (i.e. Java 1 and Nigeria 2) which share more similarities in shape with Neanderthals. The mid-Pleistocene specimens are quite homogeneous on F1 where they are segregated from Neanderthals. The cloud of points nevertheless overlaps slightly with Neanderthals. This is mainly due to the position of the Arago XIII and Tighenif 3 mandibles. Most of the dispersion of the cloud of points is observed on F2, and the African and European specimens show strong similarities in their shape, especially the Mauer and KNM-BK 67 specimens.

We can elaborate on the distributions of mandibular shape among hominins if we look at Figure 2B. Extreme shapes for modern humans show the presence of a well-developed chin (*tuber symphyseos*, pogonion, landmark #6), the absence of a *planum*

*alveolare* (planum alveolare, #9) and a relatively gracile mandibular corpus. The Neanderthals are characterized by a more vertical symphyseal profile with no well-developed *tuber symphyseos* (pogonion, #6), a slightly thicker mandibular corpus at the mental foramen level, which is also positioned slightly more posteriorly (foramina mentale, #4). Mid-Pleistocene fossils show a strongly receding symphyseal profile with no chin development (pogonion, #6). The antero-posterior thickness of the symphysis at the level of the dental arch is important (infradentale, #7 and infadentale posterior, #8), and the *planum alveolare* is well-developed (planum alveolare, #9). A wide *incisura submentalis* is present under the lower rim of the symphysis (tuberculus marginalis superius, #3 and gnathion, #5) which is absent in AMHs and weakly developed in Neanderthals.

The position on the scatter plot of our specimen of interest, Mezzena, has been calculated *a posteriori*. Unsurprisingly, the Mezzena mandible does not present any particular affinities with mid-Pleistocene specimens. It is most similar to AMHs being positioned within the *H. sapiens* cloud of points and the DFA classifies the specimen with modern humans (Table S7). Especially its shape is similar to that of Ohalo II and to a lesser extent to the recent modern human specimen China5. However, it should be noted that its position also indicates affinities with some Neanderthal specimens: the late Neanderthal Spy 1 and Saint-Césaire, the Near-East specimens Tabūn II and Amud 1, and to a lesser extent the classic Neanderthals La Ferrassie 1 and Guattari III (Figure 2).

## Genetic Analysis

About 100 mg of bone powder were removed by drilling the bone of the Mezzena mandibular corpus; DNA was extracted in a DNA laboratory exclusively dedicated to ancient DNA work. We performed three different extractions and two different PCR for each extracts. After purification of these PCR products the

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Discriminant Function Analysis based on landmarks data (A) and associated cranial shapes (B). A) Crosses indicate centroids of each a priori sample. Triangles = modern humans (1–10: Holocene specimens with indication of sex when known, labeled: fossil specimens); circles = Neanderthals; squares = *H. heidelbergensis* sample; diamond = Mezzena included a posteriori in the analysis. Labels descriptions are provided in Table 2. B) The configuration of landmarks is indicated by circles superimposed on views of the Mezzena mandible (full, visible landmarks; empty, landmarks non visible in the current view); shapes in *norma lateralis*(upper left), *norma verticalis*(lower left) and *norma facialis*(right) are portrayed for the extremities of each axis (full lines, shape change; dashed lines, consensus). Modern humans, Neanderthals and middle Pleistocene specimens are discriminated on Function 1. Function 2 mostly shows intra-group morphological variation. The architectural shape of Mezzena is closer to modern human specimens, particularly to Abri Pataud 1. Mezzena's shape shows also strong similarities with Neanderthals, especially late Neanderthals (e.g. La Ferrassie 1, Spy 1).</figcaption></figure>

doi:10.1371/journal.pone.0059781.g002
<table><caption>**Table 2.** Landmarks used in for the geometric morphometrics analysis.</caption><thead><tr><th>n°</th><th>Name</th><th>Description</th><th>type</th></tr></thead><tbody><tr><td>1</td><td>Basal</td><td>point vertically positioned under the centre of the <i>proeminentia lateralis</i></td><td>III</td></tr><tr><td>2</td><td><i>proeminentia lateralis</i></td><td>centre of the <i>proeminentia lateralis</i></td><td>II</td></tr><tr><td>3</td><td><i>tuberculus marginalis superius</i></td><td>most latero-inferior point of the tubercle</td><td>II</td></tr><tr><td>4</td><td><i>foramina mentale</i></td><td>point at the centre of the <i>foramen mentale</i></td><td>I</td></tr><tr><td>5</td><td>Gnathion</td><td>most inferior midline point on the mandible</td><td>II</td></tr><tr><td>6</td><td>Pogonion</td><td>most anterior midline point on the chin</td><td>II</td></tr><tr><td>7</td><td>Infradentale</td><td>most anterior midline point between the alveoli of the incisors</td><td>II</td></tr><tr><td>8</td><td>infradentale posterior</td><td>most posterior midline point between the alveoli of the incisors</td><td>II</td></tr><tr><td>9</td><td><i>planum alveolare</i></td><td>most posterior midline point of the <i>planum alveolare</i></td><td>II</td></tr><tr><td>10</td><td>Geni</td><td>most posterior midline point of the mental spines</td><td>II</td></tr></tbody></table>

Number, name, description and type for each landmark.

doi:10.1371/journal.pone.0059781.t002

sequence of the hypervariable region I of mtDNA was divided into seven amplicons but due to high DNA degradation we obtain only the fragment NL 16230 NH16262, in which all sequences were endogenous (Figure S1). All the substitutions observed in the MLS mtDNA jaw (determined between positions 16230 nt to 16262 nt) have been consistently reproduced in different amplifications and in three different extracts.

The comparison with other 30 partial mitochondrial Neanderthal sequences are currently available: Feldhofer 1 and Feldhofer 2 from Germany; Mezmaiskaya from the Russian Caucasus; Okladnikov from Siberia, Russia; Vindija 75, Vindija 77 and Vindija 80 from Croatia; Engis 2 and Scladina from Belgium; La Chapelle-aux-Saints and Rochers-de-Villeneuve from France; El Sidrón 441, El Sidron 1252, El Sidron R011, El Sidron 331c, El Sidron 1327 h, El Sidron 753, El Sidron 1161, El Sidron 763a, El Sidron 566, El Sidron 500, El Sidron 1634, El Sidron 763b, El Sidron 634 and Valdegoba and Cova del Gegant from Spain; Teshik Tash from Uzbekistan plus other five complete mtDNA (in three of them the HVR-I region were previously typed Feldhofer 1 and 2, Vindija, 33.16 previous called V80) and two new samples, El Sidrón 1253 from Spain, and another one from Croatia (Vindija 33.25) (for more details and references see, Table S10).

The new Mezzena mitochondrial sample of the mandible presents a classic Neanderthal motif (16234 T, 16244 A, 16256 A, 16258 G) with the diagnostic transversion 16256 C/A (see, [41] and Table S10). The nucleotides at these sites are very unlikely to reflect contamination, because they were consistently observed in amplicons also showing mutations typical of Neanderthals and not of modern humans. Moreover these substitutions were previously observed in 5 other Neanderthals (Feldhofer 1, Vindija 75, El Sidrón 441, Vindija 80 (33.16)) and in the previous Mezzena (MLS 1) cranial fragments examined [31,32].

## Discussion

The genetic analysis of the small fragment of hypervariable region 1 of the mtDNA with the well-known diagnostic Neanderthal substitutions (determined between positions 16230 nt to 16262 nt) presents a classic Neanderthal motif (16234 T, 16244 A, 16256 A, 16258 G) with the diagnostic transversion 16256 C/A which classifies the Mezzena mandible as a Neanderthal. These results are further supported by the comparative morphology analysis which shows similarities between the Mezzena mandible and Neanderthals. The symphysis of

the Mezzena mandible is very close to the European classic Neanderthals (i.e. Regourdou and Guattari III) and late Neanderthals (Spy, La Ferrassie 1, Saint-Césaire, Vindija and Las Palomas) showing an incipient mental trigone and more vertical symphyses than earlier Neanderthals [5,36,40].

However, the pattern is not quite clear, especially considering the results of the DFA analysis where Mezzena mandibular shape is more similar to that of modern humans and is classified as *H. sapiens* (Figure 2, Table S7). The DFA based on geometric morphometrics and Procrustes analysis distinguishes between the three pre-defined groups: mid-Pleistocene fossils (most specimens being attributed to *H. heidelbergensis*, [34]), Neanderthals (*H. neanderthalensis*), and AMHs (*H. sapiens*). The first axis (81.0% of variance) presents a taxonomic-based distribution of the specimens, with *H. sapiens*, *H. neanderthalensis*, and *H. heidelbergensis* being separated along this function. The second axis (19.0% of variance) accounts for most of the intra-group variation.

Nevertheless, overlaps occur between the three groups of specimens. This might be partly due to the chosen landmarks which failed to describe the full mandibular shape. They are concentrated on the symphysis and on the most anterior part of the mandibular corpus. Most of them are type II or III landmarks (see, [42]) which cannot be fully considered as “discrete anatomical loci that can be recognized as the same loci in all specimens” ([43]:23). This can explain the relatively high Wilk’s lambda value (see, Table S7). However, the preservation state of the Mezzena mandible does not allow a better description of the full mandibular shape.

This methodological problem does not rule out the fact that the DFA results are taxonomically coherent (i.e. they provide a classification which succeeds in attributing most of the specimens of our sample to their correct alleged species, Figure 2 and Table S7) and support the existence of different species among our data: *H. sapiens* and *H. neanderthalensis* of course, and to a lesser extent *H. heidelbergensis*. Additionally, studies of morphological and metrical variations in AMHs and Neanderthals already show some level of overlap (e.g. [5,44–47]) although the two populations are overall morphologically distinct (e.g. [44,45,47–49]). Finally, we must keep in mind that most of the Neanderthals diagnostic features have been identified on the mandibular corpus and especially on the ramus (e.g. retromolar space, truncated gonial angle, medial position of intersection between mandibular notch and condyle and deeply excavated pterygoid fossa, see for instance, [36,38,50]). Thus, the absence of a ramus and of part of the mandibular corpus
on the Mezzena specimen might have artificially enhanced its
resemblance to AMHS in the geometric morphometrics analysis.

In this light, we can interpret the position of the Mezzena
mandible which stands within the modern human shape space,
while presenting strong shape similarities with some Neanderthal
specimens. Such a conflicting taxonomical position is not
surprising, considering the geological age of the mandible [30].
Indeed, numerous late Neanderthals such as Spy 1, Saint Césaire
and the Near-East mandibles Amud 1 and Tabun II possess hints
of a chin (i.e. tuber symphyseo) though not a true modern human
morphology [37,51]. Late Neanderthals lived in area where
AMHs might have been already present [2,23,52], while the
Levantine fossils are displaying a less derived Neanderthal
morphology [35,36].

Therefore, in our view, this change in morphology of the
mandibular chin among the fossils of Mezzena and other late
Neanderthals could have been the result of a small degree of
interbreeding with AMHs. We must nevertheless keep in mind
that this data cannot exclude the possibility that the estimated
genetic admixture between *Homo sapiens* and *Homo neanderthalensis*
might be due to a sub-structure of an ancient African ancestor of
archaic human and present-day human populations [53], or a
more complex model recently published [54].

Thus, morphological and geomorphological analyses of the
mandible of Riparo Mezzena strongly support the hypothesis of
change in morphology on this genetic typed late Mousterian
Neanderthal. This study confirms that simple models of abrupt
behavioral and phylogenetic transition for this period in Europe
should be abandoned, at least in certain geographical areas
[5,44,55]. In Italy while AMHs with Uluzzian assemblages
reached the south of the peninsula at Grotta del Cavallo at
around 45-43 ka BP [23] and settled in Northeastern Italy close to
the Mezzena rockshelter at the site of Fumane with proto-
Aurignacian industries at 41/40 ka cal BP [56], in Riparo
Mezzena [30], morphologically and genetically identifiable
Neanderthals still making Mousterian industries had not yet
disappeared.

**Materials and Methods**

Materials

The fossil sample was selected in order to encompass a large
part of morphological variation of the middle and late Pleistocene
fossil record. 48 fossils from Africa, Asia, and Europe were studied
(middle Pleistocene specimens among which most of the individ-
uals were attributed to *H. heidelbergensis* [34]:18, *H. neanderthalensis*:
24 and *H. sapiens*: 5) (Table 1). Additionally, 10 modern humans
from Africa, Europe and Asia (4 Neolithic, 6 extant modern
humans among which 3 males, 3 females and 4 non-sexed
specimens) were included in the geometric morphometrics analysis
in order to: 1- provide a sufficient sample of modern humans
spread out over a span of time similar to that of the Neanderthals
(i.e. 130,000 years); 2- take into account the margin of error in
dating the fossil sample; and 3- test the reliability of the character
data set and the statistical method used in the study. The reduce
number of modern humans used in the DFA is due to the
obligation to respect the hypothesis of equality of the co-variant
matrices of the three groups. The three groups must be of roughly
equal size (i.e. *H. heidelbergensis* N = 12, *H. neanderthalensis* N = 15, *H.
sapiens* N = 15) to be able to interpret the results of the DFA.

Methods

**Geometric morphometrics shape analysis.** (see, [42]) is
based on 10 landmarks (Figure 2, Table 2) chosen to best describe

the mandible morphology while taking into consideration the state
of preservation of the fossils and especially of the Mezzena
specimen. The method follows the protocol described in [57]. We
ran a Generalized Procrustes Analysis, a Principal Component
Analysis (PCA) based on the procrustes residuals and a Discrim-
inant Function Analysis (DFA) to discriminate three pre-defined
groups (*H. sapiens*, *H. neanderthalensis* and *H. heidelbergensis* sensu
[34]). The number of variables must be lower than 12 (smallest
group number of specimens), thus this analysis uses the first 10
Principle Components (PC) which represent 90.02% of the total
variance (Table S2, S4 and S5). The discrimination between these
groups is used as a “pattern” to study the Mezzena mandible
which is introduced *a posteriori* in the analysis. The Wilks’ lambda
statistics (see, Table S4), used to validate the discrimination,
necessitates covariance matrices equality of each group which can
be tested using a Box’s M test (Table S4). A cross validated
classification was then ran. It successively classifies all cases but one
to develop a more reliable discriminant function and then
categorizes the case that was left out (Table S7 and S8).
Additionally, we tested the impact of size on specimen shape
modifications in order to identify a possible allometric trend in our
data. We used linear regression, which was calculated for each PCs
involved in the computation of the discriminant functions when
compared to centroid size (Table S6). We used Morphologika 2
v2.5 [58] (APG, ACP, linear regression) and SPSS v11.5 ©SPSS
Inc. 1989-2002 (linear Discriminant Function Analysis).

**Genetic analysis, experimental procedures.** We performed three different extractions and two different PCRs for each extracts. All of the most stringent protocols and procedures for the analysis of ancient DNA have been followed [59-61]. Extraction was performed as described in [62], with UNG treatment [63] in order to minimize postmortem damage. mtDNA sequences were generated by using 60 cycles of PCR and 5 ul of extract. The strong inhibitory effect of the extract required a final 1 to 100 dilution to obtain positive amplifications. Different primer pairs were used, some of them designed to match Neanderthal-specific substitutions (Figure S1, Table S9). The PCR products were cloned using the TOPO TA Cloning kit (Invitrogen), according to the manufacturer's instructions. Screening of white recombinant colonies was accomplished by PCR. The colonies were transferred into a 30-µl reaction mix (67 mM Tris HCl [pH 8.8], 2 mM MgCl2, 1 µM of each primer, 0.125 µM of each dNTP, and 0.75 U of Taq polymerase) containing M13 forward and reverse universal primers. After 5 min at 92°C, 30 cycles of PCR (30 s at 90°C, 1 min at 50°C, 1 min at 72°C) were performed and clones with an insert of the expected size were identified by agarose-gel electrophoresis. After purification of these PCR products with Microcon PCR devices (Amicon), a volume of 1.5 µl was cycle-sequenced, according to the BigDye Terminator kit (Applied Biosystems) supplier's instructions. The sequence was determined using an Applied BioSystems 3100 DNA sequencer. The hypervariable region I of mtDNA was divided in seven amplicons (L 15995- H16132; L16022- H16095; L 16106- H16282; NL 16223 -H16385; NL16230- NH16262; L16299-H16400; 3; L16311- H16402) but due to high DNA degradation we obtain only the fragment NL 16230 NH16262, in which all sequences were endogenous (Figure S1). All the substitutions observed in the MLS mtDNA jaw (determined between positions 16230 nt to 16262 nt) have been consistently reproduced in different amplifications and in three different extracts.
## Supporting Information

**Figure S1** Sequences of the 31 clones from which the consensus sequence was determined in the MLS Neanderthal jaw sample. The first line reports the human reference sequence (CRS) with the numbering of the nucleotide positions. Second line reports the sequences of primers used. Nucleotides identical to the reference sequence are indicated by dots. Clones are identified by an abbreviation and three numbers: the first number indicates the extraction; the second number indicates the PCR, the third number indicates the amplicon.
(DOCX)

**Table S1** Holocene modern humans included in the geometrics morphometric analysis.
(DOCX)

**Table S2** Main Principal Components from the procrustes shape analysis. Eigenvalues. percentage of variance and percentage of cumulated variance for each principal component.
(DOC)

**Table S3** Discriminant Function Analysis: Box's M results on the covariance matrices of the three predefined groups. p>0.5, the hypothesis of equality of the covariance matrices is accepted. The covariance matrices of the three groups are considered to be equal.
(DOC)

**Table S4** Discriminant Function Analysis: quality of the discrimination. The Wilks' lambda results validate the discrimination for function 1 at p<0.0001. Function 2 is less discriminating (i.e. Wilks' lambda = 0.504, p = 0.005).
(DOC)

**Table S5** Discriminant Function Analysis: Principal Component contribution to each discriminant function and coefficient for each function. CP2 and CP1 contribute the most to the first discriminant function while CP3 and 7 contribute the most to the second discriminant function (in bold).
(DOC)

## References

1. Smith FH, Falsetti AB, Donnelly SM (1989) Modern human origins. Yearbook of Physical Anthropology 32: 35-68.

2. Frayer D (1992) Evolution at the European edge: Neanderthal and Upper Paleolithic relationships. Préhistoire Européenne 2: 9-69.

3. Wolpoff MH (1999) The systematics of Homo. Science 284: 1774-1775.

4. Trinkaus E, Moldovan O, Milota S, Bilgar A, Sarcina L, et al. (2003) An early modern human from the Pestera cu Oase, Romania. Proceedings of the National Academy of Science of the USA 100: 11231-11236.

5. Janković I, Karavanić I, Ahern JCM, Brajkovic D, Mauch Lenardić J, et al. (2011) Archaeological, Paleontological and Genomic Perspectives on Late European Neandertals at Vindija Cave, Croatia. In: Condeki S, Weniger GC, editors. Continuity and Discontinuity in the Peopling of Europe: One Hundred Fifty Years of Neanderthal Study, 299 Vertebrate Paleobiology and Paleoanthropology: Springer. 299-313.

6. d'Errico F, Zilhão J, Julien M, Baffier D, Pelegrín J (1998) Neanderthal Acculturation in Western Europe? A Critical Review of the Evidence and Its Interpretation. Current Anthropology 39: S1-S44.

7. Zilhão J (2012) Personal Ornaments and Symbolism Among the Neanderthals. Developments in Quaternary Science 16: 35-49.

8. Lévéque F, Vandermeersch B (1980) Découverte des restes humains dans un niveau castelperronien à Saint-Césaire (Charente Maritime). Comptes Rendus de l'Académie des Sciences, Paris 291, Série D: 187-189.

9. Churchill SE, Smith FH (2000) Makers of the early Aurignacian of Europe. Am J Phys Anthropol Suppl 31: 61-115.

10. Zilhão J (2006) Neandertals and Moderns Mixed, and It Matters. Evolutionary Anthropology 15: 183-195.

11. Mellars P (2005) The Impossible Coincidence. A Single-Species Model for the Origins of Modern Human Behavior in Europe. Evolutionary Anthropology 14: 12-27.

12. Hublin J-J, Talamo S, Julien M, David F, Connet N, et al. (2012) Radiocarbon dates from the Grotte du Renne and Saint-Césaire support a Neandertal origin

**Table S6** Linear regression results for the ten first principal components when compared to centroide size. The R² values indicate that the linear regression is not a good approximation of the data with a maximum of 19.0% for PC1 and 20.3% PC4 of the data explained by the linear regression. Additionally, Fisher's tests are not significant
(DOC)

**Table S7** Classification of the specimens from the original DFA and from the cross validation procedure. Specimens presented with their original assigned group (i.e. #1: *H. sapiens*, #2: *H. neanderthalensis* and #3: *H. heidelbergensis*) and their predicted group. Discriminant function scores are indicated for each specimen. Incorrectly classified specimens are signalled with **. Note the attribution of the Mezzena mandible to the *H. sapiens* group (#1).
(DOC)

**Table S8** Summary of the original classification and of the cross validated classification. 95.2% of original cases are correctly classified and 78.6% of cross validated groups cases are correctly classified. Groups numbers: #1: *H. sapiens*, #2: *H. neanderthalensis* and #3: *H. heidelbergensis*.
(DOC)

**Table S9** Primers sequences used in this study. NL, NH: primers designed to match with Neanderthal-specific substitutions.
(DOC)

**Table S10** Mitochondrial DNA sequences. MtDNA sequences showing Neanderthal diagnostic positions in HVRI (hypervariable region I).
(DOC)

## Author Contributions

Edited, discussed, and commented on the manuscript: SC AM PG ML DC LL. Conceived and designed the experiments: SC AM LL DC. Performed the experiments: SC AM ML DC. Analyzed the data: SC AM PG DC LL. Wrote the paper: SC AM PG ML DC LL.

for the Châtelperronian. Proceedings of the National Academy of Science of the USA 109: 18743-18748.

13. Mellars P, Stringer C (1989) The Human Revolution Behavioural and Biological Perspectives on the Origins of Modern Humans. Princeton: Princeton University Press. 800 p.

14. Bräuer G, Rimbach KW (1990) Late archaic and modern Homo sapiens from Europe, Africa, and Southwest Asia: Craniometric comparisons and phylogenetic implications. Journal of Human Evolution 19: 789-807.

15. Van Andel TH, Davies W, editors (2003) Neanderthals and Modern Humans in the European Landscape of the Last Glaciation: Archaeological Results of the Stage 3 Project. Cambridge: McDonald Institute for Archaeological Research. 278 p.

16. Stewart JR, Stringer CB (2012) Human evolution out of Africa: the role of refugia and climate change. Science 335: 1317-1321.

17. Davies SWG (2001) A very model of a modern human industry: New perspectives on the origins and spread of the Aurignacian in Europe. Proceedings of the Prehistoric Society 67: 195-217.

18. Finlayson C, Carrion JS (2007) Rapid ecological turnover and its impact on Neanderthal and other human populations. Trends Ecol Evol 22: 213-222.

19. Serre D, Paabo S (2004) Evidence for gradients of human genetic diversity within and among continents. Genome Res 14: 1679-1685.

20. Green RE, Krause J, Ptak SE, Briggs AW, Ronan MT, et al. (2006) Analysis of one million base pairs of Neanderthal DNA. Nature 444: 330-336.

21. Green RE, Krause J, Briggs AW, Marici T, Stenzel U, et al. (2010) A Draft Sequence of the Neandertal Genome. Science 328: 710-722.

22. Palma di Cesnola A (1993) Il paleolitico superiore in Italia. Firenze: Garlatti e Razai Editori.

23. Benazzi S, Douka K, Fornai C, Bauer CC, Kullmer O, et al. (2011) Early dispersal of modern humans in Europe and implications for Neanderthal behaviour. Nature 479: 525-528.
24. Douka K, Grimaldi S, Boschian G, del Lucchese A, Higham TF (2012) A new chronostratigraphic framework for the Upper Palaeolithic of Riparo Mochi (Italy). *J Hum Evol* 62: 286-299.

25. Camps M, Higham T (2012) Chronology of the Middle to Upper Palaeolithic transition at Abric Romani, Catalunya, *J Hum Evol* 62: 89-103.

26. Maroto J, Vaquero M, Arrizabalaga Á, Baena J, Baquedano E, et al. (2012) Current issues in late Middle Palaeolithic chronology: New assessments from Northern Iberia. Quaternary International 247: 15-25.

27. Walker MJ, López-Martínez MV, Ortega-Rodrigáñez J, Haber-Uriarte M, López-Jiménez A, et al. (2012) The excavation of buried articulated Neanderthal skeletons at Sima de las Palomas (Murcia, SE Spain). Quaternary International 259: 7-21.

28. Corrain C (1968) Resti scheletrici umani del "Riparo Mezzena". Memorie del Museo di Storia Naturale di Verona, Sez Scienze dell'Uomo 16: 97-101.

29. Giunti P, Longo L (2011) Production system of the mousterian lithic industry of layer 3 from Riparo Mezzena (Verona, northern Italy). In: Longo L, editor. Integrated methodological approaches to the study of lithic industries Human Evolution Special Issue. 131-144.

30. Longo L, Boaretto E, Caramelli D, Giunti P, Lari M, et al. (2012) Did Neandertals and anatomically modern humans coexist in northern Italy during the late MIS 3? Quaternary International 259: 102-112.

31. Caramelli D, Lalueza-Fox C, Condemi S, Longo L, Milani L, et al. (2006) A highly divergent mtDNA sequence in a Neandertal individual from Italy. Curr Biol 16: R630-632.

32. Lalueza-Fox C, Rompler H, Caramelli D, Staubert C, Catalano G, et al. (2007) A melanocortin 1 receptor allele suggests varying pigmentation among Neanderthals. Science 318: 1453-1455.

33. Lari M, Rizzi E, Milani L, Corti G, Balsamo C, et al. (2010) The microcephalin ancestral allele in a Neanderthal individual. PLoS ONE 5: e10648.

34. Mounier A, Marchal F, Condemi S (2009) *Is *Homo heidelbergensis* a distinct species? New insight on the Mauer mandible. Journal of Human Evolution 56: 219-246.

35. Condemi S, Tardivo D, Foti B, Ricci S, Giunti P, et al. (2012) A case of an osteolytic lesion on an Italian Neanderthal jaw. Comptes Rendus Palevol 11: 79-83.

36. Condemi S (2001) Les Néanderthaliens de La Chaise. Paris: Comité des Travaux Historiques et Scientifiques. 178 p.

37. Stefan VH, Trinkaus E (1998) Discrete trait and dental morphometric affinities of the Tabun 2 mandible. Journal of Human Evolution 34: 443-468.

38. Rosas A (2001) Occurrence of Neanderthal features in mandibles from the Atapuerca-SH site. American Journal of Physical Anthropology 114: 74-91.

39. Barroso Ruiz C (2001) Le site moustérien de la grotte du Boquete de Zafarraya, province de Malaga en Andalousie, Espagne. Paris: Museum National d'Histoire Naturelle. 395 p.

40. Walker MJ, Lombardi AV, Zapata J, Trinkaus E (2010) Neandertal mandibles from the Sima de las Palomas del Cabezo Gordo, Murcia, southeastern Spain. Am J Phys Anthropol 142: 261-272.

41. Briggs AW, Good JM, Green RÉ, Krause J, Maricic T, et al. (2009) Targeted retrieval and analysis of five Neandertal mtDNA genomes. Science 325: 318-321.

42. Bookstein FL (1991) Morphometric Tools for Landmark Data: Geometry and Biology. Cambridge: Cambridge University Press. 456 p.

43. Zelditch ML, Swiderski DL, Sheets SD, Fink WL (2004) Geometric Morphometrics for Biologists: A Primer. New York and London: Elsevier Academic Press. 443 p.

44. Smith FH, Jankovic I, Karavanić I (2005) The assimilation model, modern human origins in Europe, and the extinction of Neandertals. Quaternary International 137: 7-19.

45. Trinkaus E (2006) Modern Human versus Neandertal Evolutionary Distinctiveness. Current Anthropology 47: 597-620.

46. Weaver TD (2009) The meaning of Neandertal skeletal morphology. Proceedings of the National Academy of Sciences 106: 16028-16033.

47. Trinkaus E (2011) Late Pleistocene adult mortality patterns and modern human establishment. Proc Natl Acad Sci U S A 108: 1267-1271.

48. Holliday TW (2006) Neanderthals and modern humans: an example of a mammalian syngameon? In: Harvati K, Harrison T, editors. Neanderthals Revisited. New York: Springer. 336.

49. Schwartz JH, Tattersall I (2006) Morphology, Variability, and Systematics: Lesson from Krapina. Periodicum Biologorum 108: 389-401.

50. Creed-Miles M, Rosas A, Kruszynski R (1996) Issues in the identification of Neandertal derivative traits at early post-natal stages. Journal of Human Evolution 30: 147-153.

51. Schwartz GT, Tattersall I (2000) The human chin revisited: what is it and who has it? Journal of Human Evolution 38: 367-409.

52. Higham T, Jacobi R, Julien M, David F, Basell L, et al. (2010) Chronology of the Grotte du Renne (France) and implications for the context of ornaments and human remains within the Chatelperronian. Proc Natl Acad Sci U S A 107: 20234-20239.

53. Eriksson A, Manica A (2012) Effect of ancient population structure on the degree of polymorphism shared between modern human populations and ancient hominins. Proceedings of the National Academy of Science of the USA 109: 13956-13960.

54. Sánchez-Quinto F, Botigué LR, Civit S, Arenas C, Ávila-Arcos MC, et al. (2012) North African Populations Carry the Signature of Admixture with Neandertals. PLoS ONE 7: e47765.

55. Trinkaus E (2011) Late Neandertals and Early Modern Humans in Europe, Population Dynamics and Paleobiology. In: Condemi S, Weniger GC, editors. Continuity and Discontinuity in the Peopling of Europe: One Hundred Fifty Years of Neanderthal Study, 299 Vertebrate Paleobiology and Paleoanthropology: Springer. 315-329.

56. Higham T, Brock F, Peresani M, Broglio A, Wood R, et al. (2009) Problems with radiocarbon dating the Middle to Upper Palaeolithic transition in Italy. Quaternary Science Reviews 28: 1257-1267.

57. Mounier A, Condemi S, Manzi G (2011) The Stem Species of Our Species: A Place for the Archaic Human Cranium from Ceprano, Italy. PLoS ONE 6: e18821.

58. O'Higgins P, Jones N (2006) Tools for statistical shape analysis. York: Hull York Medical School.

59. Cooper A, Poinar HN (2000) Ancient DNA: Do It Right or Not at All. Science 289: 1139.

60. Pääbo S, Poinar H, Serre D, Jaenicke-Després V, Hebler J, et al. (2004) GENETIC ANALYSES FROM ANCIENT DNA. Annual Review of Genetics 38: 645-679.

61. Willerslev E, Cooper A (2005) Ancient DNA. Proc Biol Sci 272: 3-16.

62. Caramelli D, Lalueza-Fox C, Vernesi C, Lari M, Casoli A, et al. (2003) Evidence for a genetic discontinuity between Neandertals and 24,000-year-old anatomically modern Europeans. Proceedings of the National Academy of Science of the USA 100: 6593-6597.

63. Sampietro ML, Caramelli D, Lao O, Calafell F, Comas D, et al. (2005) The genetics of the pre-Roman Iberian Peninsula: a mtDNA study of ancient Iberians. Ann Hum Genet 69: 535-548.

64. Martin R, Saller K (1957) Lehrbuch der Anthropologie in Systematischer Darstellung mit Besonderer Berücksichtigung der Anthropologischen Methoden; Verlag GF, editor. Stuttgart: Jena. 2999 p.