Exploring late Paleolithic and Mesolithic diet in the Eastern Alpine region of Italy through multiple proxies

Abstract Objectives The analysis of prehistoric human dietary habits is key for understanding the effects of paleoenvironmental changes on the evolution of cultural and social human behaviors. In this study, we compare results from zooarchaeological, stable isotope and dental calculus analyses as well as lower second molar macrowear patterns to gain a broader understanding of the diet of three individuals who lived between the end of the Late Pleistocene and the Early Holocene (ca., 17–8 ky cal BP) in the Eastern Alpine region of Italy. Materials and methods We analyze individuals buried at the sites of Riparo Tagliente (Verona), Riparo Villabruna, and Mondeval de Sora (Belluno). The three burials provide a unique dataset for diachronically exploring the influence of climatic changes on human subsistence strategies. Results Isotopic results indicate that all individuals likely relied on both terrestrial and freshwater animal proteins. Even though dental calculus analysis was, in part, hindered by the amount of mineral deposit available on the teeth, tooth macrowear study suggests that the dietary habits of the individuals included plant foods. Moreover, differences in macrowear patterns of lower second molars have been documented between Neanderthals and modern humans in the present sample, due to a prevalence of Buccal wear among the former as opposed to higher values of Lingual wear in modern human teeth. Discussion Isotopic analyses have emphasized the contribution of animal proteins in the diet of the three foragers from the Eastern Alpine region. The possible intake of carbohydrate‐rich plant foods, suggested by the retrieval of plant remains in dental calculus, is supported by the signal of macrowear analysis. Moreover, the latter method indicates that the distribution of macrowear in lower second molars (M2s) allows us to discriminate between Neanderthals and modern humans within the present reference sample. Overall, our results show these three prehistoric hunter‐gatherers were well adapted to the environment in which they lived exploiting many natural resources.

a unique dataset for diachronically exploring the influence of climatic changes on human subsistence strategies.
Results: Isotopic results indicate that all individuals likely relied on both terrestrial and freshwater animal proteins.Even though dental calculus analysis was, in part, hindered by the amount of mineral deposit available on the teeth, tooth macrowear study suggests that the dietary habits of the individuals included plant foods.Moreover, differences in macrowear patterns of lower second molars have been documented between Neanderthals and modern humans in the present sample, due to a prevalence of Buccal wear among the former as opposed to higher values of Lingual wear in modern human teeth.
Discussion: Isotopic analyses have emphasized the contribution of animal proteins in the diet of the three foragers from the Eastern Alpine region.The possible intake of carbohydrate-rich plant foods, suggested by the retrieval of plant remains in dental calculus, is supported by the signal of macrowear analysis.Moreover, the latter method indicates that the distribution of macrowear in lower second molars (M 2 s) allows us to discriminate between Neanderthals and modern humans within the present reference sample.Overall, our results show these three prehistoric huntergatherers were well adapted to the environment in which they lived exploiting many natural resources.
Investigating various strands of evidence can be particularly relevant in the case of individuals who lived in periods of climatic changes, which can trigger specific patterns of human-environment interaction in distinctive regional settings.Between the end of the Last Glacial Maximum (LGM-~16.5ka cal yr BP; Lambeck, Rouby, Purcell, Sun, & Sambridge, 2014) and the early Holocene (after 11.7 ka cal yr BP), climatic temperature oscillations set in motion significant changes in human culture, mobility, adaptive strategies, and population structure across continents.
After the LGM, the climate amelioration triggered a general withdrawal of the main glaciers in the Eastern Alpine region.The retreat of glaciers in the region brought to the repopulation of mountainous areas by animals and human groups (Bertola et al., 2007).The abundant archeo-osteological evidence and the long-lasting tradition of prehistoric studies in this area make the Eastern Alpine region significant for exploring the influence of climatic changes on subsistence strategies.
These sites span the period from ~16.5 to 8 ka cal BP which falls in the Late Upper Paleolithic and the end of the Late Mesolithic (Figure 1) and are crucial for understanding modalities of human-environment interaction after the LGM in the analyzed region.
To explore prehistoric dietary habits, we first provide the paleoenvironental background of the area between the end of the Late Pleistocene and the early Holocene.Then we combine data obtained through different analytical methods, generally performed independently, such as (a) analysis of faunal assemblages recovered in layers that are coeval to the burials; (b) stable isotope analysis; (c) dental macrowear analysis.Also, we analyzed dental calculus, limited by a small amount of mineralized biofilm preserved on the three individuals, for its potential to inform about dietary strategies.Last, we discuss dietary strategies of the three foragers who lived in the Eastern Alpine region against paleoenvironmental data already available for the region of interest.et al., 2008;Monegato, Scardia, Hajdas, Rizzini, & Piccin, 2017) and sea-level dropped of ~120 or 149 m (Pellegrini et al., 2015;Maselli et al., 2014).The pre-Alpine area underwent glacial and peri-glacial conditions with profound modifications of the ecosystems, the evolution of the physiographic outline of the valley floors (Rossato, 2013), extinction of the less adapted large herbivore and carnivore species (Sala, 1990;Terlato et al., 2018), setting the grounds for the development of the recent ecosystems.
The subsequent Late Glacial period (i.e., the interval between the end of the LGM and the beginning of the Holocene) was punctuated by a series of millennial-scale climate oscillations and dynamic glacial fluctuations (Seguinot et al., 2018).
After the major glaciers retreated, the evolution of the valley floor was controlled by slope-damming processes.Larch, scots pine, and mountain pine forests were limited to the foreland and the external belt of Pre-Alps (Ravazzi et al., 2014;Vescovi et al., 2007).During this period, groups of foragers moved on the southern slope of the Pre-Alps.
The oldest known archeological site on the Southern slope of the Alps to be re-occupied by human groups after the LGM was Riparo Tagliente (Fontana et al., 2008, Fontana, Cilli, et al., 2009;Fontana et al., 2018;Naudinot, Tomasso, Tozzi, & Peresani, 2014) spanning the first part of the Late Glacial until the beginning of the Bølling-Allerød interstadial (Fontana, Cilli, et al., 2009;Fontana et al., 2018).During the first part of the Late Glacial, an individual was buried in the area protected by the overhang of the shelter (Level 16-13; Tagliente 1).A pit opened in historical times partially destroyed the burial and removed the entire skeleton above the pelvis (Table 1).A hemi-mandible (Tagliente 2) was found within disturbed deposits just outside the rock-shelter.These deposits probably derived from the partially destroyed Late Upper Paleolithic sequence of the inner area of the shelter, something also suggested by the association of the mandible with Epigravettian archeological materials (Corrain, 1966).
The climatic warming at the onset of the Bølling-Allerød (B-A) interstadial, dated ca.14.7 ka cal BP, promoted the development of alpine habitats to higher altitudes.Such favorable climatic conditions promoted the increase in the number of Epigravettian sites and gradual colonization of high altitudes (i.e., above 1,000 m a.s.l.).The geographic distribution of the settlements is attested between the valley-bottom areas and the middle mountain range mostly in the vicinity of peat bogs or landmarks, also in the inner Alps and the Dolomitic region (Avigliano, Di Anastasio, Improta, Peresani, & Ravazzi, 2000;Bertola et al., 2007;Broglio & Impronta, 1994-1995;Montoya, Duches, Fontana, Peresani, & Visentin, 2018) and karst plateaux up to an altitude of 1,600 m a.s.l.(Bertola et al., 2007;Broglio, 2001;Broglio & Impronta, 1994-1995;Fiore & Tagliacozzo, 2005;Montoya et al., 2018;Peresani, Bertola, De Stefani, & Di Anastasio, 1999).Within this context, seasonal settlement systems based on the interplay of sites with complementary functionality developed.Such behavior is well attested at Riparo Tagliente as well as Riparo Villabruna (Aimar et al., 1992;Fontana, Cilli, et al., 2009;Fontana et al., 2018;Naudinot et al., 2014).The latter site, also object of investigation in this article, is located at 500 m a.s.l. in the Venetian Dolomites at the confluence of the Cismòn Valley with Rosna creek and characterized by a group of rock shelters situated in the proximity of the valley bottom (Broglio & Villabruna, 1991).Two meters of deposit have excavated at the site, and an Epigravettian burial was brought to light at the bottom of the stratigraphic sequence.Grave goods, represented by isolated blades, a bone point, a backed knife, a flint blade, a flint core, a stone retoucher and a ball of probable resin and wax, were positioned probably in a bag, on the left forearm near the pelvis (Broglio & Villabruna, 1991).Natural stones and large painted stones collected on the stream bed close to the shelter were used to cover the body after the internment.Radiocarbon dates obtained on a charcoal from the grave sediment (Layer 17A; R-2023:12,040 ± 150 yr BP; Aimar et al., 1992), and the human cranium (KIA-27004:12.140 ± 70;Vercellotti et al., 2008) place the burial at 14,190-13,770 B.P.calibrated using IntCal13 in the OxCal 3.0 program (Ramsey et al., 2009;Reimer et al., 2013; Table 1).
Faunal remains indicate that during the summer and autumn, Epigravettian groups used to occupy higher altitudes to exploit the variety of biological resources (Fiore & Tagliacozzo, 2005;Fiore & Tagliacozzo, 2006;Phoca-Cosmetatou, 2005a, 2005b, 2009).Overall, red deer and ibex were the most frequent herbivores hunted although chamois, roe deer, and wild boar were also common in some faunal assemblages, while aurochs and horse were less targeted.Specialized hunting of marmots was also practiced (Romandini, Peresani, Gurioli, & Sala, 2012).
The successive Younger Dryas (YD) cold event, which started at 12.7 ka cal BP, had strong effects on stands of thermophilus trees formerly expanded in the forelands of the Southern Alps.Treelines shifted 200-300 m downward, whereas steppe and alpine grasslands with xerophytic shrubs expanded.Unfortunately, the available archeological data are not detailed enough to advance an alternative explanation for supposed continuity in settlement dynamics during the YD (Mussi & Peresani, 2011;Naudinot et al., 2014).The YD oscillation preceded a sharp increase of temperature identifying the final transition to the Holocene.At the beginning of these interglacial conditions, within a few centuries, the timberline reached an altitude of about 2,100 m a.s.l.(Oeggl & Wahlmüller, 1994).Afterward, a rapid cooling and moisture increase occurred at ca., 8.2 ka cal BP, which is clearly defined in the Greenland ice core δ 18 O records, favoring the expansion of spruce (Picea sp.) and silver fir (Abies sp.) in the subalpine belt (Pini, 2002).
During the early Holocene, the whole of north-eastern Italy was covered by forests.In the Alps, the belt of territory between 1,700 and 2,100 m asl was occupied by Larix-Pinus cembra forests, and at around 10,500 BP, the tree-line reached 2,250 m (Drescher-Schneider, 2009;Wick, 1994).This territory was occupied by Mesolithic hunter-gatherers attested from the present coastal areas of the Northern Adriatic to the inner Alpine highlands.Intense colonization of the central-eastern Italian Alpine region is documented with evidence of occupation both in the mountains and on valley-bottoms.We can thus infer a social organization based on small size groups with mobility based on seasonal displacements from the plain/valley-bottoms to highland territories through the main drainage systems Highland sites were situated at recurrent locations (in the proximity of small lakes and passes, along panoramic crests, either open-air or under the overhang of erratic boulders) and altitudes, mostly between 1,900 and 2,300 m a.s.l. and occupied during the warmer season of the year (Broglio, 1980;Broglio & Lanzinger, 1990;Fontana & Visentin, 2016;Visentin et al., 2016) Sauveterrian, for instance, hunting of ungulates was particularly relevant besides other small prey such as pikes and cyprinids, sweet water mollusks, and marsh turtles in the valley bottoms (Wierer and Boscato, 2006;Bazzanella. Betti, & Wierer, 2007).
One important site for understanding the Late Mesolithic occupation of the region is Mondeval de Sora.This site is located at 2,150 m a.s.l.under the overhang of a large erratic boulder on a terrace in the high valley of the Cordevole River (Belluno Dolomites).
Traces of human occupation were unearthed under two sides of the boulder (Sectors I and III; Alciati et al., 1994;Fontana & Guerreschi, 2003;Fontana, Guerreschi et al., 2009;Fontana, Govoni, et al., 2009;Valletta, Fontana, Bertola, & Guerreschi, 2016).In Sector I, located along the south-western side, a complex stratigraphic series was explored over a surface of about 60 m 2 and yielded layers dated between the Early Mesolithic and Medieval age.Mesolithic layers were preserved only in the southern portion of the site.The Early Mesolithic (Sauveterrian) sequence documents some dwelling structures associated to two main layers containing a rich lithic assemblage along with faunal and charcoals remains (Alciati et al., 1994;Berto, Luzi, Guerreschi, Fontana, & Valletta, 2016;Colombo et al., 2016;Fontana, Govoni, et al., 2009;Fontana & Guerreschi, 2003;Fontana, Guerreschi, et al., 2009;Fontana & Vullo, 2000;Thun Hohenstein et al., 2016;Valletta et al., 2016).The Late Mesolithic (Castelnovian) layers, although disturbed by the later occupation of the site, yielded a well-preserved burial of a young individual.The skeleton lied in a supine position within a pit naturally delimited by two dolomite boulders with the lower part of the body covered with stones (Fontana, 2006;Fontana et al., 2016;Gerhardinger & Guerreschi, 2006).The individual was accompanied by a rich repertoire of burial goods consisting of 60 items of various typologies laid around the body, and on the right side of the skeleton, a small patch of red ochre was recovered.One direct AMS radiocarbon date on the skeleton yielded a result of 7,425 ± 55 years BP (OxA-7468) (8,380-8,070 cal BP;Fontana, Guerreschi, Bertola, Briois, & Ziggiotti, 2016) calibrated using IntCal13 in the OxCal 3.0 program (Ramsey et al., 2009;Reimer et al., 2013; Table 1).

| Riparo Tagliente
Two different human remains have been discovered: Tagliente 1 is a burial preserved from the pelvis to the feet and belonging to a ca.20-29-year-old young adult male (Bartolomei et al., 1974;Corrain, 1977), Tagliente 2, consists of a hemi-mandible with teeth in sockets (Figure S1).
Considering the uncertainty of the stratigraphic position of Tagliente 2, we believe these two human remains likely belonged to different individuals.In this study, we focus our attention on Tagliente 2 to implement morphometric and morphological information previously published (Corrain, 1966).The teeth have been virtually extracted (Figure S1) and analyzed (Figure S2) in a digital model generated from μCT data.Deposits of dental calculus on the teeth of the hemi-mandible of Tagliente 2 were removed from the original specimen for conservative cleaning procedures.A minimal quantity of the mineralized biofilm left on the first molar was sampled (ca., 0.4 mg).

| Riparo Villabruna
The complete skeleton discovered at Riparo Villabruna belongs to a ca.25-year-old adult male, 167 cm tall.Overall, the human remains are in an excellent state of preservation (Vercellotti et al., 2008).However, the distal portions of both lower limbs are incomplete due to damage that occurred during the work that led to the discovery of the site.Both the cranium and the mandible are complete, although the lower left first incisor has been lost postmortem, and the lower-left second premolar as well as the lower right second molar were not preserved.Dental calculus from the Villabruna skeleton was removed for conservative cleaning procedures.Only a minimal quantity of dental calculus was left and sampled for the analysis (ca., 0.5 mg).

| Mondeval de Sora
The buried individual from Mondeval de Sora belongs to a ca.40year-old male, 165 cm tall.The remains are incredibly well preserved, lacking only the most fragile components of the facial part and some distal bones of the left foot.The preservation of the skeleton allowed various studies to be performed, including stable isotope analysis (Gazzoni, 2011).Dental calculus was removed for conservative cleaning procedures, and the very small sample left was collected for micro-debris analysis (ca., 0.8 mg).

| Comparative sample for the macrowear analysis
The sample used in this study is conservative and consists of second mandibular molars (M 2 ) of hunter-gatherers.The specimens listed in  In this work, already available data on the analyzed human and faunal remains (e.g., the isotopic analysis and zooarcheological assemblages for Tagliente 1 and Riparo Villabruna) are discussed against the results obtained through the application of different methods (macro-wear, isotope, and dental calculus analyses) on Tagliente 2, Villabruna and Mondeval individuals.

| Data acquisition
As far as Villabruna specimen is concerned, the high resolution of the polygonal models of the M 2 , obtained by μCT images (Oxilia et al., 2015), allowed us to avoid the use of smoothing of the meshes preserving the original surface from the segmented μCT data.
The hemi-mandible from Riparo Tagliente structures were visualized by high-resolution μCT at the Multidisciplinary Laboratory of the Abdus Salam International Centre of Theoretical Physics (Trieste, Italy), using a system specifically designed for the study of paleontological and archeological materials (Tuniz et al., 2013).The μCT acquisition of the complete specimen was carried out by using a sealed X-ray source (Hamamatsu L8121-03) at a voltage of 110 kV, a current of 90 μA and with a focal spot size of 5 μm.The X-ray beam was filtered by a 1 mm-thick aluminum absorber.A set of 1,800 projections of the sample were recorded over a total scan angle of 360 by a flat panel detector (Hamamatsu C7942SK-25) with an exposure time/projection of 1.5 s.The resulting μCT slices were reconstructed using the commercial software DigiXCT (DIGISENS SAS) in 32-bit format and obtaining an isotropic voxel size of 41.33 μm.
To avoid any damage, we preferred to acquire macrowear data from a dental cast of the Mondeval de Sora individual.Original dentitions and high-resolution dental replicas were molded using a lowviscosity polyvinylsiloxane impression material (President Light-body; Coltene, Switzerland).The replicas were then produced using special gypsum (EverestH Rock, KaVo), possessing nonreflective properties and optimized for light surface scanning (Fiorenza, Benazzi, & Kullmer, 2009).Three-dimensional digital models were generated using a high-resolution (up to 50 μm) Artec Eva-S Spider 3D blue light technology scanner.Collection and alignment of the scan-data point clouds were carried out using Artec studio 12 Professional software.
Even though we used different approaches (μCT and surface scan of dental casts) to acquire the original sample, virtual models do not show significant differences in terms of reliability if compared to the original morphology.The 3D virtual models were further processed using Geomagic Design X (3D System Software) in order to remove errors and degenerate/duplicate triangles.

| Faunal assemblage analyses
Taxonomic and skeletal identifications were based on the reference collections stored at the Department of Humanities of the University of Ferrara and the Department of Anatomy, Pharmacology and Forensic Medicine of the University of Torino.Taphonomical analysis of faunal assemblages was performed by following protocols previously published (Aimar et al., 1992;Fontana, Cilli, et al., 2009;Thun Hohenstein et al., 2016; see Supporting information).As far as Mondeval de Sora and part of Villabruna assemblages are concerned, we only provided preliminary new data (Table 2).

| Stable isotopes of Mondeval de Sora and Tagliente 2
Stable isotopes of collagen from Tagliente 2 (from the root of the lower left first molar) were measured using a Thermo Finnigan Flash EA coupled to a Delta V isotope ratio mass spectrometer at the Max Planck Institute for Evolutionary Anthropology, Leipzig (MPI-EVA) (see Supporting information and Table 3).

| Dental calculus analysis
Most of the dental calculus deposits on the teeth of the three skeletons were removed before the time of our analysis for conservation/ cleaning purposes.As a consequence, very little mineralized deposit was left and sampled from the Mondeval de Sora individual (0,7 mg), Villabruna and Tagliente 2 (ca., 0,5 mg and 0,4 mg, respectively).Subgingival dental calculus was sampled with sterile disposable blades and saved in sterile centrifuge tubes.Powder-free gloves were worn at all times during the sampling (see Supporting information).Samples were transported to the DANTE-Diet and Ancient Technology laboratory (Sapienza University) for analyzing masticatory, dietary and environmental microfossils entrapped in their matrix.
Ancient samples were processed in a dedicated clean space and strict procedures were followed to avoid any modern contamination of the sample and the working area (see Supporting information).
Once in the laboratory, calculus samples were put into 1.5 ml sterile tubes, filled with ultrapure water, and agitated using a "vortex" machine to shake any soil particles adhering to their surfaces.The calculus was then removed, washed in ultrapure water and placed again into a sterile centrifuge tube.The sequence was repeated until water washed off clean.Calculus samples were then examined using a stereomicroscope.Any soil still adhering was removed using a fine acupuncture needle dampened with 0.6 N HCl.All the liquids produced during the decontamination were put aside and later analyzed as a control check.The "clean" calculus samples were then inserted in tubes and covered with HCl 6 N. Powder-free gloves and masks were always used during the procedure.Once dissolved, the calculus was pipetted out and mounted on a glass slide with a drop of 50:50 mixture of glycerol and water.Slides were observed using a Zeiss Imager cross-polarized transmitted light as well as a NIKON ECLIPSE cross-polarized microscopes, both equipped with DIC prism and examined at magnifications ranging from ×200 to ×630.Each type of microfossil recovered in dental calculi was described, recorded and photographed.Starch granule identification was achieved using methodologies and criteria known and accepted in the field of modern and ancient starch granules research (Henry et al., 2009;Henry, Brooks, & Piperno, 2011;Yang & Perry, 2013).In addition to this, plant micro-remains entrapped in dental calculus were interpreted using a modern reference collection of microdebris extracted from more than 200 species of plants native to the Mediterranean region stored at the DANTE laboratory and the Department of Archeology at the University of York.Our reference collection also included animal dietary debris (e.g., muscle, skin, fur particles, fish meat, and scales), inorganic nondietary debris that can be trapped in dental calculus during crafting activities (e.g., plant fibers used for cordage, wood particles, etc.), and minerals (e.g., ochre).Our interpretation of bird feather fragments is based on both published literature (Dove, 1998;Dove & Agreda, 2007;Dove & Koch, 2011;Harwood, 2011) and an extensive reference collection including aquatic bird species from the Anatidae, Strigidae, and Accipitridae families.
Feather fragments were primarily interpreted based on the node shape and width, node density per mm, and pigmentation.
Occlusal wear facet areas (polished enamel areas with well-defined margins) of M2s were analyzed using the occlusal fingerprint analysis (OFA) method, a digital approach used to reconstruct the major mandibular movements occurring during the rhythmic chewing cycle (Kullmer et al., 2009).Wear facets were manually mapped on each digital surface model, and labeled according to the numbering system created by Maier and Schneck (1981) and later modified by Kullmer et al. (2009).
Wear facets can be produced by occlusal tooth to tooth contact during masticatory processes.The vertical movement (puncturecrushing) is the first step of the masticatory cycle in mammals (Crompton & Hiiemae, 1970;Hiiemäe & Crompton, 1971;Kay & Hiiemae, 1974) during which the dental wear produces rough surfaces blunting the cusp tips, generated by the contact of the tooth surface with food and any other extrinsic materials (i.e., grit or sand; Kay, 1977;Stone, 1948).Once food has been partitioned, a rhythmic chewing phase (or power stroke) follows where the jaw movements are guided by the occlusal relief of the upper and lower teeth (Crompton & Hiiemae, 1970;Kay & Hiiemae, 1974) insalivating the food particles into a bolus.The power stroke is divided into two phases (Hiiemae & Kay, 1972;Kay & Hiiemae, 1974)

| Faunal assemblage analyses
The large mammal assemblage from the Late Epigravettian series of Riparo Tagliente shows differences in composition between the lower (17-13) and the upper layers (12-5).While the former layers reflect a dominance of species that lived in open environments, the latter ones document a sensible increase of taxa adapted to milder climates, which attest for the diffusion of woodland areas in the region (Fontana, Cilli, et al., 2009).The composition of the small mammal assemblage confirms such a pattern (Berto, Luzi, Montanari Canini, T A B L E 3 Bone samples from Riparo Tagliente (Gazzoni, 2011), Riparo Villabruna (Vercellotti et al., 2008), and Mondeval de Sora (Gazzoni, 2011) 2) (Fontana et al., 2018;Gazzoni et al., 2013).Carnivores, lagomorphs and rodents are also well represented, especially Marmota marmota (around 30%).As far as freshwater fish is concerned, bones were also present and currently under investigation.Butchery marks consisting of cut marks and scraping marks from skinning, defleshing, disarticulation, and tendon removal on 10% of ungulates (ibex, roe deer, red deer, and chamois) as well as percussion notches for marrow extraction are documented.Several marmot bones exhibit cut marks.Burned fragments are also frequent (Fontana, Cilli, et al., 2009;Rocci Ris, 2006).
The faunal assemblage from the Late Epigravettian occupation layers of Riparo Villabruna shows a prevalence of ungulates mainly represented by Capra ibex, Cervus elaphus, and Rupicapra rupicapra (Aimar et al., 1992; Table 2).The lower layers (17-10) are dominated by Caprinae, suggesting cold and dry climatic conditions.In the upper layers (9-3), the increase of red deer remains suggests a more temperate and wet climate.
As far as the layers of the burial (16 and 17) are concerned, the species mainly hunted were ungulates such as C. ibex (56,5%), R. rupicapra (22.8%), and C. elaphus (6.5%;Table 2).Butchery marks are present on 21% of the total analyzed bones with traces of skinning, defleshing, disarticulation, and removal of tendons.Only six burned fragments have been identified.
The faunal assemblage from layers coeval to the Castelnovian burial of Mondeval de Sora is very small, partially mixed due to postdepositional events.A prevalence of undeterminable fragmentary bones, not useful for paleoenvironmental and paleoeconomic analyses, is also documented.Nonetheless, considering Mesolithic occupation layers in their whole C. elaphus is the most represented taxon (Table 2).
Bone and antler remains were also found among the grave goods.

| Stable isotope analysis
As far as Tagliente 2 is concerned, the value of δ 13 C is −19.5‰ while the δ 15 N is equal to 11.5‰ (Figure 2 and Table 3).Compared with the fauna from the same site, Tagliente 2 δ 13 C values are 0.3 and 0.1‰,
Faunal data from Riparo Tagliente, Riparo Villabruna and Mondeval de Sora are also shown (see Table 3).In the inset, intrasite differences in δ 15 N between human and (mean) fauna data are summarized food held in their diet.Our results should be considered for their potential to inform about the awareness of certain plant species but not as the clear evidence of their regular use by the analyzed foragers.
In region and the Balkan peninsula, was conducted for previous and current analyses of dental calculus (Cristiani et al., 2016;Cristiani et al., 2018;Nava et al., in press).In our reference collection, the dimensions of Type A starch granule of Aegilops can vary from 10 μm in some species (e.g., in A. ovata and A. neglecta) to more than 30 μm in other (e.g., A. triuncialis, A. peregrina, A. geniculate, A. caudalis, etc.).
Based on our reference collection, and excluding any modern contamination in the analyzed sample or in the laboratory, we are confident the starch granules from different Aegilops species from our reference collection (see Figure S3).
Together with the morphotype above, several small irregular starch granules, other plant structures were identified in the calculus (Figure 3, 15-18) as well as animal remains.In particular, numerous barbule fragments, consistent with birds of the Anatidae family based on the shape, features, and distance of their nodes, were also found (Figure 3,(19)(20).

| Dental wear pattern
Relative wear areas of the three individuals are characterized by high values of Lingual Phase I facets (Table S3).Mondeval (BPI: 0.27) exhibits the lowest value of Buccal Phase I wear when compared to Tagliente 2 (BPI: 0.33) and Villabruna (BPI: 0.31), as opposed to the trend observed for Phase II facets (Figure 4a).
In order to understand the segregation which emerged among our individuals, we compared the three specimens with a reference sample (Table S1) created ex novo, comprising both Neanderthal and modern human M 2 s in wear stage 2/3 (Smith, 1984).Reference individuals are grouped based on Ecogeographic zones, MIS-groups, geography, and species.
The comparison of Phase II proportions separated into ecogeographic zones provided no significant differences, both considering the interspecific reference sample and the intraspecific modern human specimens only (Table S4).In contrast, Buccal Phase I exhibits significant differences, although after Bonferroni correction, no significant pairwise result was obtained from posthoc Mann-Whitney tests (Table S5).These results are confirmed by the lack of population structure as measured via AMOVA (Φ ST = 0.07, p-value = .054).
No significant differences were found after correction between MIS-based groups for any of the investigated macrowear variables (Table S6).The result once again is confirmed by the lack of population structure (Φ ST = 0.09, p-value = .09).
The impact of geography on the distribution of dental macrowear in the reference sample was also investigated (Table S7).No significant differences were found through univariate analyses, while geography explains a considerable amount of variability when all variables are considered at once (Φ ST = 0.14, p-value = .033).
The segregation between Neanderthals and modern humans emerging from the ternary diagram (Figure 4c) is confirmed by AMOVA (Φ ST = 0.28, p-value <.001; Table S8).Univariate analyses support significant differences between Neanderthal and modern humans for all masticatory phases (Figure 5).Considering the abovementioned results, we focused on modern human specimens alone for the following analyses.
Lingual Phase I of modern humans exhibits no significant difference between individuals adapted to Mediterranean evergreen (MED) and Steppe-Coniferous forest (SCF; Mann-Whitney W = 2, p-value = .07;Figure 6a), although effect size shows that adaptation to SCF has a more substantial impact on Buccal Phase I while adaptation to MED impacts more Lingual Phase I wear.On the other hand, a significant difference between MIS5 and MIS2 emerged for Buccal Phase I wear (W = 7.47, p-value = .02)as suggested by statistical analysis (effect size-Figure 6 B,C).

| DISCUSSION
By integrating different data from various sources, we achieved an indepth reconstruction of the subsistence and dietary preferences of three foragers who lived between the end of the Late Pleistocene and the Early Holocene in the Eastern Alpine region.
Overall, the isotopic analysis and the study of faunal remains indicate a significant consumption of terrestrial and, to a less extent, of aquatic animal proteins; however, the small number of faunal samples limits our interpretation.
The individual discovered at Riparo Tagliente (Tagliente 2) exhibits traces of a strong reliance on animal proteins while butchering marks on the faunal remains (Fontana et al., 2018;Fontana, Cilli, et al., 2009) point to ungulates as one of the likely sources of protein for this individual.Besides, the δ 15 N level measured in human remains is well above a trophic step shift (i.e., >3-5‰; Hedges & Reynard, 2007) compared to local herbivores and omnivores, suggesting an additional consumption of high trophic level animals (e.g., foxes) and/or freshwater resources.The composition of the faunal assemblage supports both hypotheses (see Bartolomei et al., 1982;Fontana, Cilli, et al., 2009;Gazzoni, 2011).
The δ 13 C value of −19.5‰ found at Riparo Tagliente (Tagliente 2) seems consistent with a typical terrestrial environment and freshwater fauna.As the δ 13 C ratio is generally predictive of a marine (higher δ 13 C) versus freshwater (lower δ 13 C) diet, the exploitation of freshwater fish seems the more conservative hypothesis for this individual, considering that Riparo Tagliente was located ca. 360 km away from the Adriatic Sea at the time of the Late Paleolithic occupation (Gazzoni et al., 2013).Accordingly, the consumption of freshwater resources has also been suggested for other Epigravettian human remains (e.g., Arene Candide, Liguria; Grotta del Romito, Calabria; and San Teodoro, Sicily).
In a wider perspective, the isotopic signature of Tagliente 2 is close to San Teodoro 7 (δ 15 N = 11.5 and δ 13 C = −19.1‰;Mannino, Di Salvo, et al. 2011) and in line with the high protein diet suggested for other Upper Paleolithic sites across Italy (Figure 2).
Interestingly, the other Epigravettian human from Riparo Tagliente (Tagliente 1) presents both higher δ 15 N (13‰) and δ 13 C (−18.4‰), indicating an even stronger reliance on aquatic resources (Gazzoni, 2011;Gazzoni et al., 2013;Lugli et al., 2019).Both these values may indicate that marine proteins cannot be entirely excluded from the diet of the individuals of Riparo Tagliente (Gazzoni, 2011;Gazzoni et al., 2013).Human displacements during the Late Glacial in Europe, as reflected in the archeological record, may cover comparable distances, even longer in Western and Central Europe (Roebroeks, 2001).Flint provenance also offers a unique opportunity Epigravettian (Bertola, 2012;Bertola, Fontana, & Visentin, 2018;Peresani, 2019).However, exclusive reliance on terrestrial animals cannot be entirely ruled out for these individuals.For example, the consumption of lactating juvenile animals, cooked or putrid meat may increase the δ 15 N value above the canonical trophic level shift (see, for example, Jaouen et al., 2019).
Tagliente 1 and Tagliente 2 show a δ 13 C and a δ 15 N differences of 1.1 and 1.5, respectively.Such differences might result from the different types of samples analyzed, thus suggesting a slightly different diet between juvenile and adult individuals.In particular, stable isotopes from Tagliente 2 individual were measured on the root of the Ternary diagrams showing the proportions (in %) of relative wear surface of Buccal Phase I areas, Lingual Phase I areas, and Phase II areas (Kay & Hiiemae, 1974), 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%.Our three individuals (Tagliente, Villabruna, and Mondeval) were identified in the ternary plots as black points.as the main cause of the isotopic variation (see Fahy, Deter, Pitfield, Miszkiewicz, & Mahoney, 2017;Olsen et al., 2014).
The δ 15 N and δ 13 C analysis performed on the individual of Villabruna shows clear evidence of the consumption of terrestrial proteins, even though the intake of alternative proteins (e.g., aquatic resources) cannot be excluded.δ 15 N and δ 13 C values are too low to fit in the marine protein range (Richards & Hedges, 1999).The low Human isotopic values of the Late Mesolithic individual from Mondeval de Sora suggest a strong reliance on terrestrial animal proteins again, with some potentially limited intake of fish and freshwater resources (see below).On average, the isotopic composition of herbivores from the site is −20.3 ± 1.2‰ (1σ) for carbon and 2.6 ± 1‰ (1σ) for nitrogen (n = 5;Gazzoni, 2011).The fact that the human collagen δ 15 N is well above a typical trophic step shift when compared to the herbivore mean (Δ human-herbivore = 6.5‰) indicates the likely additional intake of aquatic proteins for this individual (Gazzoni, 2011;Gazzoni et al., under review).However, considering a δ 13 C of −19.9‰ for the human collagen, such additional intake might be attributed to freshwater resources rather than marine resources.
Freshwater resources tend to have a wider δ 13 C range of values, overlapping the values commonly observed in terrestrial ecosystems F I G U R E 5 Univariate analyses result in significant differences between Neanderthal (N) and modern humans (S) for all masticatory phases (Phase II, Buccal Phase I and Lingual Phase I) (Fuller, Müldner, Van Neer, Ervynck, & Richards, 2012;Richards & Hedges, 1999).Future sulfur isotope analyses may help to untangle this issue, aiming to understand whether fish consumption occurred at Mondeval the Sora.
Overall, the isotopic data obtained for the three foragers reveal a diet primarily based on terrestrial animal proteins with contribution from freshwater and possibly marine resources.A similar pattern was already suggested for other Late Glacial and Mesolithic humans from Italy (Mannino & Richards, 2018).Overall, dietary data available for Mediterranean foragers (including Spain, Croatia and Corsica) indicate a heavy reliance on terrestrial proteins (Mannino & Richards, 2018).
The use of aquatic resources is often patchy and mostly observed in groups along the Atlantic (see Richards & Hedges, 1999) or during later phases of the Mesolithic (see Cristiani et al., 2018;Mannino et al., 2015).No specific dietary shift is detected isotopically between late Upper Paleolithic and early Mesolithic.Yet, a long-term trend (from ca., 30 to 8 ka) of both δ 15 N and δ 13 C values is evident in human collagen data (Figure 7).Such a pattern may have resulted shows isotope values compatible with a predominant terrestrial diet.
However, the presence of a harpoon among the grave goods and the δ 15 N value slightly above one-trophic step shift may suggest a role of freshwater resources in the diet (Hedges & Reynard, 2007).In the calculus of the Mondeval individual, feather barbules belonging to birds consistent with the Anatidae family were also recovered, supporting the possible consumption of such aquatic resource (as indicated by stable isotopes).This being said, the possibility that the feather fragments could have been ingested after the extra-masticatory use of the teeth cannot be ruled out (e.g., using teeth as a tool while hafting weapons, and/or in the production of ornaments).Higher exploitation of aquatic proteins as a seasonal resource occurred during the Late Glacial phases of climatic instability (Mannino & Richards, 2018).
Dental calculus and dental macrowear analyses provide a more balanced picture of the diet of Eastern Alpine foragers, underlining a possible contribution of plant species as food.In particular, starch granules belonging to grass grains known to be of dietary importance were recovered in the analyzed dental calculi, hence providing the first, direct evidence that local foragers consumed vegetal resources during their life.Interestingly, a higher familiarity with starchy foods is documented for Mondeval de Sora forager compared to the other individuals (Figure 3).However, we cannot exclude that these results are related to the different amounts of mineralized deposit analyzed in this work.Leoni, 1993;Castelletti & Leoni, 1984;Leoni, Castelletti, & Castiglioni, 2002;Mussi et al., 2000).In Sicily, wild legumes (Lathyruys sp., Pisum sp.), acorns (Quercus sp.), arboreal fruits (Arbutus unedo), and wild grapes (Vitis silvestris) have been found in the Mesolithic Phase I facets in groups of individuals who relied on mixed food resources, including a large percentage of plants (Fiorenza, Benazzi, Tausch, et al., 2011).A more recent study, however, has shown a mirrored occlusal relationship in the macrowear pattern of upper and lower molars in an Australian Aboriginal population from Yuendumu (Oxilia et al., 2018).(Borgognini Tarli & Repetto, 1985).A similar pattern is also visible in the ternary diagram (Figure 4), where Villabruna, Riparo Tagliente, and Mondeval de Sora individuals plot within the range of variability attributed to Neanderthals and modern humans from the MED group (and partly near the SCF group).Previous analyses of these groups suggested the exploitation of various food sources (Fiorenza, Benazzi, Tausch, et al., 2011), including large proportions of meat, as documented in Australian Aboriginal populations that consumed large quantities of animal proteins (O'Dea et al., 1991).When only modern humans are observed, univariate exploratory results suggest that all three individuals exhibit high values of Lingual wear consistent with SCF and out of the range of individuals adapted to MED ecozones (Figure 4a).
In particular, the Villabruna sample seems to be closer to steppe coniferous forest sample groups, as suggested by the distribution pattern of macrowear (Figure 4a), and from faunistic and palaeoenviromental data.In accordance with previous radiocarbon dating (Alciati, 1992), Villabruna falls in the fossil group that lived during MIS 2 (Figure 4b), while Neanderthals are clustered on the opposite side of the ternary plot (Figure 4c).
The ecogeographic zone that roughly defines the environment in which Tagliente 2 individual lived is attributable to steppe coniferous F I G U R E 7 Carbon and nitrogen stable isotope values (see Table S2) are plotted over time for Upper Paleolithic and Mesolithic huntergatherers from Italy.Oxygen isotope data from NGRIP (Rasmussen et al., 2014) are also reported to show the climatic variability forest (Figure 4a), in agreement with the relevant zooarchaeological findings.Similar to the individual from Villabruna, the forager from Tagliente 2 falls closer to humans of the MIS 2 stage (Figure 4b).
The Mondeval de Sora individual presents similar Lingual macrowear values (Figure 4a).This result cannot be supported or refuted on the basis of the available zooarcheological evidence.As far as MIS-based grouping is concerned, the Mondeval individual falls closer to humans of MIS 2 stage in the ternary diagram.However, its coordinates are closer to those of individuals adapted to warmer climates (Figure 4b).This is consistent with the direct radiocarbon dating performed on the individual (8,377-8,067 cal BP) (Fontana et al., 2016) as well as with paleoenvironmental data available for the site (Fontana, Govoni, et al., 2009).
The reference sample analyzed in this study provides additional information suggesting the possibility of discriminating between Neanderthals and modern humans by observing the distribution of macrowear in M 2 s.Considering, however, that the present reference sample is conservative and too small to generate inference at a population level, future studies should be focused on a wider sample to confirm this signal.
To sum up, through the application of a multi-proxy approach, we successfully unraveled additional details on the dietary strategies of

1. 1 |
Paleoenvironmental and archeological background The Eastern Alpine region of Italy is one of the best-suited areas to investigate the environmental and cultural influence on the subsistence and dietary preferences of hunter-gatherer groups of the last deglaciation.This region is located north of the Po valley, and the Friulian-Venetian plains in a wide pre-Alpine belt of mountains peaking over 2000 m, karst plateau elevated to 1,000-1,200 m, dissected by gorges, river valleys, and alpine lakes.During the LGM, European Alps were extensively covered by an ice-dome reaching a maximum ice extent around ca., 25 ka cal yr BP, when massive Piedmont glaciers advanced onto the Alpine foreland (Ivy-Ochs . The site, located at an altitude of 250 m a.s.l.along the corridor of the Valpantena valley giving access to the Lessini Plateau.A 4.5 m deep stratigraphic sequence, formed by two main deposits separated by a river escarpment, documents lower strata referred to MIS 4-3 occupation with Mousterian and Aurignacian assemblages and upper ones dated to the Late Glacial, attesting a Late Epigravettian occupation.According to radiocarbon dates (calibrated using IntCal09), which range from 17,219-16,687 cal BP (LTL4441A; layer 13a alpha) to 14,572-13,430 cal BP (R-371; layers 10-8), the Epigravettian series of Riparo Tagliente is one of the most complete in Northern Italy F I G U R E 1 Paleogeographic context of the three individuals buried at the sites of Riparo Tagliente (Verona), Riparo Villabruna and Modeval de Sora (Belluno).Snapshots showing Alpine glaciers extension modeled during three different time intervals.Frames modified from Seguinot et al., 2018 (Modeling last glacial cycle ice dynamics in the Alps: https://doi.org/10.5446/35164,fixed frame version) Middle/Upper Paleolithic (i.e., Neanderthals and anatomically modern humans) and Epipaleolithic Natufian individuals.Each individual was assigned to an ecogeographic group based on literature (Fiorenza, 2015), to geographic areas based on the region of discovery, to climatic classes based on the relevant Marine Isotope Stage (MIS), and to the respective species (Homo neanderthalensis and Homo sapiens).
: the first (Phase I, divided in Lingual and Buccal) happens when opposite molar cusps tend to slide past each other moving to centric occlusion (maximum intercuspation).The second phase (Phase II) is an anterior-medial movement, where the lower molars move out of occlusion.The chewing cycle terminates with the opening of the jaw.The relative surface area of the wear facets attributed to Buccal Phase I, Lingual Phase I, and Phase II of the occlusal power stroke was computed by summing the absolute areas (in mm 2 ) belonging to the same phase and dividing this sum by the total occlusal wear area.The resulting relative values (proportions) were visualized in a ternary plot, describing the relative proportions that sum to 1 or 100%.The presence of significant differences in the distribution of each macrowear variable was investigated among groups defined by species, ecogeographic zones, geographic areas, and MIS-based groups by applying Kruskal-Wallis tests (α = .05)and posthoc, two-tailed Mann-Whitney U tests for independent sample design (α = .025).In the case of species, one-tailed Mann-Whitney tests (α = .05)were run based on observation of variable distributions.To avoid problems of multiple testing, all Mann-Whitney p-values were treated with Bonferroni correction.We measured effect size to assess the magnitude of the difference between ecogeographic and MIS-based groups of modern humans using Cohen's d for unpaired samples with pooled standard deviation and Hedge's correction for small sample size via the Cohen's d function of the package effsize in R(Torchiano, 2018) The presence of population structure ascribable to the different grouping criteria was assessed by applying Analysis of Molecular Variance (AMOVA;Excoffier, Smouse, & Quattro, 1992) to the pairwise Euclidean distance computed on all reference individuals, in order to use all three macrowear variables at once.
the individual of Riparo Tagliente (Tagliente 2), six irregular polyhedral starch granules (ca.10 μm in length) with central hilum and central fissures were found (Figure 3, 1-4).While the interpretation of such irregular grains remains unspecific, our experimental reference collection and published data show that similar granules of around 10 μm in size are found in the Paniceae/Andropogoneae tribe of the Poaceae grass family.Tagliente 2 displayed very few plant micro remains in the calculus sample, including one small oval starch with central hilum and a few particles of charcoal.Interestingly, plenty of red mineral residues were found entombed in the calculus, the microscopic features of which are consistent with experimental particles of ochre when observed in transmitted polarized light (Figure 35-8).Plant micro-remains consistent with starch granules, wood, and other plant fibers as well as micro-particles of charcoal were well represented in the calculus matrix of the Mondeval de Sora individual.Within the starch granules, large granules (ca., 20 μm), oval to subround in 2D shape, lenticular in 3D shape, with a central hilum and visible lamellae were the most abundant.In the Mondeval sample, such morphotypes of granules occurred as single grains (Figure3,10-14) and, in one case, as in bimodal distribution, that is, with large circular granules with a clear central hilum and a higher density of deep lamellae toward the middle portion (also known as Type-A) appearing together with small almost spherical grains with a central hilum (≤10 μm, known as Type-B; Figure3,9).In our reference collection, such features appear more consistently in starch granules of the genus Aegilops within the Triticeae tribe of the Poaceae grass family.A total of 21 starch granules assigned to morphotype 1 were found in the sample of calculus from Mondeval de Sora.Few type-A granules also appeared damaged and having lost their original birefringence, possibly due to enzymatic damage or mechanical processing of the starch.Extensive work on 13 species of the genus Aegilops (Aegilops triuncialis,A.comosa, A. crassa, A. ovata, A. biuncialis,   A. cilidrica, A. speltoides, A. neglecta, A. columnaris, A. tawskii,   A. peregrina, A. caudata, A. geniculata), native of the Eastern Alpine

F
I G U R E 3 Starch and micro-debris recovered in dental calculus from Riparo Tagliente, Riparo Villabruna and Mondeval de Sora individuals.(1-4) Small irregular starch granules with central hilum and central fissures from Tagliente; (5-8) mineral debris from Villabruna; micro-debris from Mondeval de Sora: (9) large sub-round granule occurring in bimodal distribution; (10-14) large sub-round granules; (15) egg-shape aggregate of starch granules; (16) wood particle showing tracheids and bordered pits.The bar is 10 μm; (17) charcoal remains embedded in dental calculus; (18-19) plant structure enclosed in calculus; (20-21) feather particles with evident nodes.The bar is 20 μm.The bar is 10 μm (1-8) and 20 μm (9-21) in stating that the dimensions of the archeological starch granules (ca., 20 μm) recovered in the Mondeval de Sora dental calculus (Figure 3, 9-14) fall within the range of various species of the Aegilops genus.In order to provide a solid basis to our claim about the variability of Aegilops starch granules dimensions, we decided to add a figure with for tracing Early and Early-Late Epigravettian movements across the Po river valley (Po Plain).Remarkable evidence is provided by the reddish-brownish chert contained in the Upper Cretaceous to Middle Eocene limestone-marl formation Scaglia Rossa exploited in the Central Apennines and introduced as finished equipment in the Early Epigravettian sites of the Monti Berici and Istria during the LGM and northward to Riparo Tagliente during the early phase of the Late (a) Ecozones-based group; (b) MIS-based group; (c) species-based group (Red: Homo neanderthalensis, blue: Homo sapiens) lower left first molar, the mineralization of which ends at the age of ~0 years.These isotope values thus represent a broad period of the individual's life, starting from late childhood to the adult age, due to the presence of secondary dentine, which grows during the whole life.Differences in isotopic values between Tagliente 1 and 2 might be explained in three ways: (1) Tagliente 1 and 2 belong to two individuals with slightly diverse dietary inputs or provenance; (2) Ultrafiltrated collagen (Tagliente 2) versus nonultrafiltrated collagen (Tagliente 1) yielded different stable isotope results (see Jørkov, Heinemeier, & Lynnerup, 2007); (3) Tagliente 1 and 2 belong to the same individual, with bone turnover between mandible and rib (if any) 15 N value recorded in the deer bone collagen (1.6‰) from the same context could relate to climatic variations (e.g., physiological and/or dietary herbivore adaption to cold environments; differential volatilization rate in soil compounds or 15 N depletion in soil regenerated following deglaciation; see Richards & Hedges, 2003 for a detailed discussion).Accordingly, δ 15 N values of herbivore fluctuated from 40 ka to the present with the lowest peaks observed during colder phases.
from (a) global and/or local climatic variations, affecting plants and animal physiology, or behavior (see, for example, Bocherens, Drucker, & Madelaine, 2014; Iacumin, Nikolaev, & Ramigni, 2000 and Hedges, Stevens, & Richards, 2004); (b) a lack of isotopic data for some specific periods (e.g., only three Gravettian individuals were studied from Italian contexts); (c) a dietary shift connected to climatic adaption.Concerning the latter hypothesis, the only three Gravettian hunter-gatherers for whom isotopic data are available display both high δ 15 N and δ 13 C values, indicating a more substantial reliance on aquatic (maybe marine) resources.On the contrary, the lowest isotopic values are registered during the Holocene, possibly suggesting that Mesolithic individuals consumed primarily terrestrial foods, with a marginal role of aquatic foodstuff.Accordingly, our Mesolithic individual (Mondeval de Sora)

F
I G U R E 6 (a) Box plot of Lingual Phase I among Homo sapiens considering Mediterranean Evergreen (MED) and Steppe-Coniferous Forest (SCF) ecozones.(b) Effect size between MED and SCF ecozones for each masticatory Phase In the Homo sapiens reference sample.(c) Effect size between MIS2 and MIS5e-d for each masticatory Phase In the same sample layers of Grotta dell'Uzzo (Costantini, 1989).Interestingly, direct evidence of Late Epigravettian and Mesolithic foragers' familiarity with species of grass grains of the genus Aegilops and the Tribe Paniceae comes from the dental calculi of individuals buried at Grotta Continenza (Nava et al., submitted).Our results further support the idea that Paleolithic and Mesolithic foragers consumed carbohydrate-rich grass grains despite their low recovery rate in preagrarian contexts, in the absence of specific modalities of conservation (i.e., carbonization) and systematic flotation activity (see, e.g., Divišová & Šída, 2015 for an overview of the main plant remains recovered in Mesolithic Europe).Together with dental calculus, patterns of dental macrowear from molars support a possible reliance on relatively abrasive plant foods, as signaled by high values of Lingual Phase I(Fiorenza, 2015;Lucas, 2004) wear.A previous analysis of upper molars of recent hunter-gatherers revealed a wear pattern dominated by Buccal Phase I facets in meat-eater human groups, and a pattern with large Lingual According to this comparative information, the macrowear patterns of the two Italian Upper Paleolithic and one Late Mesolithic individuals suggest a diet mainly based on meat with only a moderate intake of plant foods.Interestingly, at Grotta dell'Uzzo and Grotta della Molara, both in Sicily, dental microwear analysis carried out on Mesolithic individuals confirms plant intake as a component of the local diet three hunter-gatherers, who lived in the Eastern Alpine region between the late Upper Paleolithic and the end of the Mesolithic.Epigravettian individuals (Tagliente and Villabruna) practiced specialized hunting and exploited wider variety of biological resources linked to the local biodiversity, and including terrestrial as well as freshwater resources.This conclusion is also supported by the dominance of species that lived in open environments and adapted to milder conditions in the associated faunal assemblages.In accordance with the still scanty paleobotanical data available for the region, plant foods were also part of the Epigravettian and Mesolithic diet, as demonstrated by starch granules identified from dental calculus and documented by Lingual wear dominance from macrowear.The intense colonization of the central-eastern Italian Alpine region, which occurred during the final transition and the early Holocene, was promoted by climatic oscillation, which produced forests increasing up to 2.250 m asl.More balanced exploitation of available resources (animal proteins and plant foods) was reflected by the dietary spectrum of the individual discovered at Mondeval de Sora.Human isotopic values suggest a strong reliance on terrestrial animal and freshwater resources supplemented by plant foods, the latter documented by macrowear and possibly by dental calculus analyses.5 | CONCLUSION Data discussed in this article support the idea that foragers who lived in the Eastern Alpine region during the Late Paleolithic and Mesolithic were familiar with animal, freshwater as well as plant resources, hence suggesting a quite versatile human adaptation to the local environment.The isotope data primarily emphasized the contribution of animal proteins in the prehistoric diet.However, the combination of macrowear data, dental calculus results and paleobotanical information, by stressing the potential familiarity of forager diet with plantfoods, underlines the significance of multiple proxies for obtaining a more exhaustive picture of prehistoric foragers' diet.Besides, macrowear study generated a considerable amount of information, including the potential differences between Neanderthals and modern humans in patterns distribution on the occlusal surfaces on M 2 s.ACKNOWLEDGMENTS This study was supported by the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No 639286-HIDDEN FOODS awarded to Prof. Emanuela Cristiani; grant agreement No 724046-SUCCESS awarded to Prof. Stefano Benazzi-www.erc-success.eu-).Studies on the Villabruna site and the burial are also supported by the Feltre and Sovramonte Municipalities.Isotopic studies on the Mondeval de Sora human remains were supported by Associazione Amici del Museo of Selva di Cadore.We thank Veneto Archeological Superintendency for the access to the human remains, in particular SABAP of Verona, Rovigo and Vicenza (Riparo Tagliente) and SABAP of Belluno, Padova and Treviso (Riparo Villabruna e Mondeval de Sora).We are grateful to Dr. Francesca Rossi and Dr Nicoletta Martinelli for giving access to the human remains from Riparo Tagliente and the support she provided during the sampling phase.We are indebted to Lysann Rädisch and Sven Steinbrenner of the Department of Human Evolution at the MPI-EVA for technical assistance.
Radiocarbon dates of the burials discovered at the Riparo Tagliente, Riparo Villabruna and Mondeval de Sora.