Dietary Adaptations of Early and Middle Pleistocene Equids From the Anagni Basin (Frosinone, Central Italy)

The intermontane Anagni Basin (Frosinone, central Italy) is an important region for Italian biochronology and paleoecology due to the presence of two rich fossil assemblages dated to the Early (Coste San Giacomo) and Middle Pleistocene (Fontana Ranuccio). These sites have yielded a vast collection of large fossil mammals with a well-documented presence of fossil equids in both localities (represented mostly by isolated teeth). Coste San Giacomo is dated to around 2.1 Ma, thereby having recorded the effects of the onset of the Quaternary glacial cycles, which led to a gradual trend toward colder and more arid conditions in the Northern Hemisphere. The fossil equids of this site belong to the ﬁrst group of grazing stenonid equids of the genus Equus that spread to the Italian Peninsula during the so called “Elephant- Equus” event, which marked the appearance of new large mammals living in herds in open and arid environments. The site of Fontana Ranuccio is dated to around 400 ka, close to the MIS 12–11 succession (the “Mid-Brunhes event”), which marked the end of the Middle Pleistocene Transition. The fossil horses from Fontana Ranuccio represent one of the oldest caballoid (or “true horses”) populations of the Italian Peninsula. The Anagni Basin, thus, provides important data to investigate paleoecological adaptations of these groups of equids in response to two critical environmental and climatic shifts of the Pleistocene. We explore their niche occupation by examining long-term dental wear patterns and tooth enamel carbon and oxygen stable isotopic composition. Both taxa appear to have exhibited a narrow dietary niche, displaying a clear abrasive (highly specialized) grass-rich diet. In particular, caballoid equids from Fontana Ranuccio show a more abrasion-dominated mesowear signature. Stenonid equids from Coste San Giacomo exploited broader and more diverse landscapes during the Early Pleistocene, whereas caballoid horses from Fontana Ranuccio appeared to have limited their dietary adaptations to a stricter grazing behavior in more closed environments.


INTRODUCTION
The intermontane Anagni Basin, situated in the central part of the Italian Peninsula, developed between the Late Pliocene and the early Middle Pleistocene (Carrara et al., 1995;Galadini and Messina, 2004). Two rich vertebrate fossil sites, the Early Pleistocene locality of Coste San Giacomo (CSG) and the Middle Pleistocene site of Fontana Ranuccio (FR), were discovered in the 1970s, following different surveys carried out by the Italian Institute of Human Paleontology (IsIPU) (Figure 1). CSG has been dated to around 2.1 Ma based on magnetostratigraphy, pollen, and its record of micro mammals recovered from alluvial deposits (Bellucci et al., 2014). It is one of the few Italian localities that record the terrestrial habitat conditions after the onset of the Pleistocene glaciations that marked the Plio-Pleistocene transition at 2.6 Ma Raymo, 2005, 2007). FR fossiliferous layer located within a succession of pyroclastic deposits has been dated to around 0.4 Ma (Pereira et al., 2018), close to the climatic transition of the "Mid-Brunhes Event"  with the definitive association of the glacial cycles ruled by a 100 kyr periodicity and the consolidation of the Early-Middle Pleistocene Transition (Head and Gibbard, 2015).
Importantly, both localities have yielded abundant fossil mammal material, among which equid remains are wellrepresented (Segre and Ascenzi, 1984;Bellucci et al., 2012Bellucci et al., , 2014Strani et al., 2018a). Fossil human teeth were also unearthed in FR (Rubini et al., 2014). The equid remains from CSG were attributed to Equus stenonis by Biddittu et al. (1979), but a recent study (Palombo et al., 2017) points to stronger affinities with the middle-size Equus senezensis and proposes a new provisional identification of the CSG material as E. senezensis aff. E. sen. stehlini. This taxon is among the first groups of stenonid equids (monodactyl horses with a "V" shaped linguaflexid between the metaconid and metasylid of the lower cheek teeth ;Forsten, 1988) to reach the Italian Peninsula during the so-called "Elephant-Equus" event with the appearance of new large mammals living in herds in drier habitats (Azzaroli, 1977;Lindsay et al., 1980). Equids' remains from FR have been tentatively classified as Equus cf. E. mosbachensis (Segre and Ascenzi, 1984), belonging to one of the first caballoid (or "true horses, " monodactyl equids with a "U" shaped linguaflexid ;Forsten, 1988) populations recorded on the Italian Peninsula after the disappearance of the early Middle Pleistocene stenonid forms (Caloi, 1997;Alberdi and Palombo, 2013).
Both modern and fossil Equus taxa are usually considered highly specialized grazers that consistently occupied similar habitats through the Pleistocene and Holocene (i.e., open grasslands), despite it has been observed that extant species can also include a noticeable amount of browse in their diet (Roeder, 1999) and that some fossil Equus populations displayed less abrasive feeding behaviors (Kaiser and Franz-Odendaal, 2004). The material from the Anagni basin represents an exceptional case to investigate the habitat occupation and resource use of two groups of fossil equids and test the hypothesis that both stenonid and caballoid horses exploited exclusively open habitats and fed mostly on grasses in Central Italian Peninsula. In this study we compare dental remains of E. senezensis aff. E. sen. stehlini from CSG and Equus cf. E. mosbachensis from FR in occlusal surface mesowear, a proxy for dietary abrasiveness in herbivorous mammals (Fortelius and Solounias, 2000), and in stable carbon and oxygen isotope abundances, proxies for habitat openness (Van der Merwe and Medina, 1991) and temperature (Dansgaard, 1964), to shed light on the niche occupation and resource exploitation mechanisms of stenonid and caballoid equids on central Italian Peninsula.
Our data will allow us to obtain comprehensive information on the paleoecological adaptations of the Equus genus in the region following two major climatic events, and to understand if this group could be more flexible in terms of feeding behavior and habitat occupation during the Early and Middle Pleistocene.

Dental Mesowear
Mesowear reflects the cumulative effects at a macroscopical level of the items ingested by herbivores (foods and exogenous particles such as grit or dust) on tooth morphology that are produced in a long period of time compared to the lifespan of the animal (Fortelius and Solounias, 2000). Dental mesowear is a direct signal of a species' diet that varies between browsers, which feed mostly on soft plant resources (e.g., leaves, twigs, green stems), and grazers that clip abrasive vegetation (e.g., grasses) at or near ground level.
Attrition (tooth-to-tooth contact) and abrasion (tooth-tofood contact) are the main factors that influence the occlusal morphology of the teeth, determining the cusps sharpness and the height of the occlusal relief. In browsing animals attrition mainly determines the tooth wear leading to sharper cusps and a higher occlusal relief. In species that feed on more abrasive plant resources the higher level of tooth-to-food contact produces blunter cusps and a lower occlusal relief (Fortelius and Solounias, 2000). Traditional mesowear is limited to upper second molars (M2), but it can be extended to other tooth positions (Kaiser and Solounias, 2003;DeMiguel et al., 2010DeMiguel et al., , 2012Strani et al., 2015Strani et al., , 2018a. The method has been developed to examine the diet-related dental wear patterns in ungulates and has been successfully employed to investigate dietary adaptations of different groups of both modern and fossil artiodactyls and perissodactyls (Fortelius and Solounias, 2000;Schulz et al., 2007;DeMiguel et al., 2008DeMiguel et al., , 2011DeMiguel et al., , 2018Bernor et al., 2014;DeMiguel, 2016;Marín-Leyva et al., 2016;Maniakas and Kostopoulos, 2017). Dental mesowear also proved to be an optimal tool to study equid paleoecological adaptations, allowing to obtain information about the dietary change and evolution of this group from the early forms of the Eocene to the horses of the Pleistocene (Kaiser and Franz-Odendaal, 2004;Mihlbachler et al., 2011;Bernor et al., 2017).

Carbon Stable Isotopes
Carbon isotopic (δ 13 C) abundance in herbivores reflects the presence of C3 and C4 vegetation in lower and mid-latitude habitats as well as a degree of canopy closure vs. openness in wooded high-latitude habitats. C3 plants, or woody vegetation and grasses that have a cool growing season, have δ 13 C values ranging between −34 and −23‰, with a suggested mean of −27‰ (Calvin and Benson, 1948;Hatch et al., 1967;Smith and Epstein, 1971;Koch, 1998). The canopy effect can result in up to 5‰ differences in δ 13 C in the same species at one site in pure C3 forest environments between the herbivores feeding in a heavily enclosed wooded environment or on plants from different tree elevation, and the herbivores feeding in more open parkland or grassland environments with higher evapotranspiration (Van der Merwe and Medina, 1991;Bocherens et al., 1999;Feranec and MacFadden, 2006;Drucker et al., 2008;Hofman-Kaminska et al., 2018).
The δ 13 C values of the carbonate fraction of bioapatite in large mammalian herbivores are higher than or enriched by 14‰ compared to their diet: δ 13 C bioapatite = δ 13 C diet + 14‰ (Cerling et al., 1999). In warm temperate and cold boreal environments δ 13 C values higher than −13‰ in herbivore enamel indicates open landscapes, such as grasslands and steppes, in which the plants have δ 13 C values higher than −27‰ in diet or ecosystem δ 13 C values, whereas tooth enamel δ 13 C values lower than −13‰ indicate woodlands and denser forests (Bocherens, 2014).

Oxygen Stable Isotopes
Oxygen isotopic (δ 18 O) abundance of the herbivores' enamel apatite reflects the source and amount of body-water ingested, either free drinking water or obtained from food and plants. The δ 18 O in the enamel apatite is a proxy for δ 18 O values of meteoric water, which varies due to the source of precipitation and effects of latitude, altitude, continentality, temperature, and evaporation (Dansgaard, 1964;Longinelli, 1984;Luz et al., 1984;Bryant and Froelich, 1995;Bryant et al., 1996;Kohn, 1996;Sponheimer and Lee-Thorp, 1999;Levin et al., 2006). Crudely, δ 18 O values of meteoric water are higher in warmer or dryer environments (Dansgaard, 1964). Thus, as a general rule, herbivores feeding in open steppe habitat are expected to have more enriched δ 18 O values than species from forested cooler and more humid habitat (Sponheimer and Lee-Thorp, 1999). Obligate or water-dependent drinkers usually have δ 18 O values lower than non-obligate drinkers, being dependent on water availability and seasonality of rainfall. Seasonal variation in rainwater in mid-and high latitude environments and species migration can be significant, producing intra-tooth variations of 3-4‰, or more (e.g., in high-crowned ungulates, horses and bison, Hoppe et al., 2004;Feranec et al., 2009). Oxygen isotopic composition can be affected not only by daily ecology and migration but also by life history and physiology of mammals. Oxygen isotopic compositions are higher in nursing animals during lactation and weaning periods because of milk consumption, which is related to the body water in the mother, and more enriched in 18 O due to the preferential loss of 16 O through sweat, urine, and expired water vapor (Bryant and Froelich, 1995;Kohn, 1996).

MATERIALS AND METHODS
The material studied consists of upper cheek teeth belonging to E. senezensis aff. E. sen. stehlini from CSG and to Equus cf. E. mosbachensis from FR. All the examined fossil material is housed at the Italian Institute of Human Paleontology (IsIPU) (Anagni, Italy). The comparative material was obtained from literature (see in text).

Dental Mesowear
Past studies on the dental wear patterns of E. senezensis aff. E. sen. stehlini and Equus cf. E. mosbachensis from the Anagni Basin examined both upper and lower cheek teeth (Strani et al., 2015(Strani et al., , 2018a. In this study we examined only upper teeth [P4-M3, following Kaiser and Solounias (2003)] in their occlusal relief (high or low) and cusp shape (sharp, rounded, or blunt) of the apex of the paracone or metacone and scored for a better comparison with available data of both modern and fossil equids because mesowear analysis is often restricted to the upper dentition. Occlusal relief and cusp shape scores were also converted to a single mesowear score (MWS) following the "mesowear ruler" developed for scoring dental mesowear on fossil equids by Mihlbachler et al. (2011). The method is based on seven cusp types (numbered from 0 to 6), ranging in shape from high and sharp (stage 0) to completely blunt with no relief (stage 6). Additionally, a "stage 7" is given to teeth with a convex cusp apex. A total of 30 specimens were scored using this method (E. senezensis aff. E. sen. stehlini N = 20; Equus cf. E. mosbachensis N = 10). Previously unavailable upper teeth of E. senezensis aff. E. sen. stehlini (N = 9) have been added to the original sample analyzed in Strani et al. (2015Strani et al. ( , 2018b. Results were compared with mesowear data of four extant equids: plains zebras (Equus quagga) (data from Solounias, 2000 andSemprebon, 2010), Grévy's zebras (Equus grevyi) (data from Fortelius and Solounias, 2000), onagers (Equus hemionus), and African free ranging donkeys (Equus asinus) (data from Schulz et al., 2007). To distinguish fossil and modern equid taxa based on mesowear variables (percentages of high relief, rounded, and blunt cusps) a hierarchical cluster analysis was performed using the Ward's method and the Euclidean distance. E. senezensis aff. E. sen. stehlini dental mesowear data were compared to those of E. stenonis from the Early Pleistocene locality of Olivola (∼1.8 Ma, north central Italy) (Strani et al., 2018c) and of the Britain sites of Mundesley and East Runton (∼1.8 Ma) (Rivals and Lister, 2016). Equus cf. E. mosbachensis dental wear patterns were also compared to those of the fossil Middle Pleistocene caballoid equid population from the central European localities (Steinheim and Heppenloch,MIS 11, Germany) (data from Rivals and Ziegler, 2018), from different Britain sites (Clacton, Hoxne and Swanscombe; MIS 11) (Rivals and Lister, 2016) and from Caune de l'Arago (438 ± 31 ka, France) (data from Rivals et al., 2015a).
Discriminant analysis was performed to examine the resolution of mesowear variables applied to the fossil. The percentages of high relief, rounded and blunt cusps were used as independent variables and two dietary (conservative and radical) classifications of modern ungulate taxa as grouping variables (Fortelius and Solounias, 2000). Equus quagga from Rumuruti (Rivals and Semprebon, 2010), Equus hemionus and Equus asinus (Schulz et al., 2007) were added to the modern taxa with known diets and treated as grazers in both conservative and radical classifications. All analyses were performed using IBM SPSS Statistics 24.

Stable Carbon and Oxygen Isotope Analyses and Preparation of Teeth Enamel Samples
We obtained stable carbon (δ 13 C) and oxygen isotopic (δ 18 O) abundances from the upper cheek teeth (upper first or second molars). In horses the isotopic composition in M1 or M2 would reflect nutrients from mother's milk and should be slightly more depleted in 13 C and more enriched in 18 O compared to isotopic Number of specimens (N); percentage of specimens with high (%High) and low (%Low) occlusal relief; percentage of specimens with sharp (%Sharp); rounded (%Rounded), and blunt (%rBlunt); mesowear score (MWS).
Frontiers in Ecology and Evolution | www.frontiersin.org composition of subsequently erupting P2, P3, M3, and P4 that would better reflect solid diet of an adult animal (Hoppe et al., 2004). Here we made comparisons with M1 or M2 from German sites and took into account this difference when comparing samples taken from premolars. We used 3 teeth samples of Equus cf. E. mosbachensis from FR and two E. senezensis aff. E. sen. stehlini from CSG. One FR horse sample was from PM4 and exhibited enrichment by 1‰ in 13 C and more enriched in 18 O compared to two FR M2. In German sites only E. hydruntinus had permanent teeth compared and the obtained values were within variation ranges of 1 and 2 ‰ in 13 C and 18 O, respectively. Enamel samples were either mechanically separated and ground to a finer condition in an agate mortar or drilled with a diamond-impregnated rotary tool. The uppermost surface of the tooth was removed before sampling usually by drilling with a dremel © borer. Samples of approximately 20 mg were analyzed. Fossil enamel powders were pre-treated following the method described by Bocherens et al. (1996) but using the Continuous Flow (CF) for collecting carbon dioxide Bocherens et al., 2009). Our experience with fossil material from limestone cave has shown that the pre-treatment used here is strong enough to remove exogenous carbonates but we are aware of possible slight changes of the oxygen isotopic composition of tooth enamel carbonate.
Pretreated (Sample /Enamel) carbonate was reacted with 100% H 3 PO 4 for 4 h at 70 • C using a MultiFlow-Geo interfaced with the Elementar IsoPrime 100 IRMS. Final isotopic ratios are reported per mil (‰) calibrated with international standards (IAEA-603 δ 13 C: 2.46/δ 18 O: −2.37 and NBS-18 δ 13 C: −5.014/δ 18 O: −23.2), as well as three in-house standards. IonOS software (Version 3.2) by Elementar was used to carry out multi-point standard isotope calibration by generating a trend line (y = mx+c) that maps measured vs. expected isotopic results, which is then used to calibrate sample results. The measurement uncertainty was monitored using three in-house standards. The overall analytical precision is higher than 0.1‰ for carbon and better than 0.2‰ for oxygen isotopic values.
Stable isotopic results are expressed as the following standard δ-notation: X = [(Rsample/Rstandard)-1] * 1,000, where X is referred to δ 13 C and δ 18 O values and R is equivalent to 13 C/ 12 C or 18 O/ 16 O, respectively. The recorded delta values follow the international reference standards, a "Vienna PeeDee Belemnite" (VPDB) for the carbon and oxygen. Additionally, δ 18 O values relative to "Vienna Standard Mean Ocean Water" (VSMOW) are given. In addition, the carbonate content (CaCO 3 %) was determined using the ratio between amount of CO 2 released by the reaction, as detected from the peak intensity for mass 44 and the weight of pure carbonate used as a standard, with an analytical error of 0.3%, based on multiple analysis of reference enamel samples.
To calculate values for temperature and drinking water from oxygen isotope values, and convert oxygen carbonate and phosphate fractions, we used and equation from Lécuyer et al. (2010) δ 18 O c = 1.035 δ 18 O P + 8.33 The estimates of the δ 18 O values of drinking water consumed by species were calculated, first, by converting the oxygen isotope composition of structural carbonate to its equivalent in phosphate, using equations that are robust for many different species: δ 18 O P = 0.71 × δ 18 O W +22.60 (Delgado Huertas et al., 1995).
Then the oxygen isotope composition of phosphate was converted to that of drinking water of different groups of species, according to formulas below and reformulated based on original data from Skrzypek et al. (2011) based on current relationship between air and water. We calculated temperature ( • Celcius) using the formula for continental Europe: We used both parametric and non-parametric tests with α = 0.05. To identify the sources of significant differences we used ANOVA and among unequal variances we used post hoc Tamhane tests.
The results were compared with the late middle and late Pleistocene horse data from Steinheim, Mauer, Bockstein and Vogelherd (Pushkina et al., 2014).

Dental Mesowear
E. senezensis aff. E. sen. stehlini upper teeth from CSG showed mostly low occlusal relief (75%) and a higher percentage of rounded cusps (45%) over sharp (35%) and blunt cusps (20%) ( Table 1). Low occlusal relief and high occurrence of rounded cusps E. senezensis aff. E. sen. stehlini from CSG point to a generally abrasive diet. Equus cf. E. mosbachensis teeth also displayed low occlusal relief and mostly rounded cusps (50%) ( Table 1), a pattern pointing to a dominance of abrasion over attrition.

Stable Carbon and Oxygen Isotope Analyses
The δ 13 C values of Equus cf. E. mosbachensis from the Middle Pleistocene FR ranged from −13.93 to −12.63‰ (Table 2 The only peculiar difference between the Italian and German sites is demonstrated in the temperature calculations, indicating much higher temperatures for both CSG and FR (6.2-12 • C) than for the German sites (−0-6.8 • C) ( Table 2).

DISCUSSION
Mesowear and stable isotopic analyses suggest that horses in both sites inhabited rather open and dry environments with abrasive C 3 vegetation as a dietary source for the majority. As obligate drinkers, however, horses could have used different water sources in more humid and closed, or more open and arid environments, which could have been at slightly different elevation and subjected to seasonal variation in water availability. Both E. senezensis aff. E. sen. stehlini from CSG and Equus cf. E. mosbachensis from FR display dental mesowear patterns highly consistent with a grazing feeding behavior, according to the results reported in Strani et al. (2015Strani et al. ( , 2018a, where both upper and lower cheek teeth were analyzed. Differences, however, are observed in terms of niche occupation in the Early Pleistocene stenonid equid, displaying a less abrasion-dominated diet in a more open landscape of steppe, shrubland, or light forest compared to the Middle Pleistocene caballoid species that fed strictly on abrasive sources in a more closed habitat of woodland or shrubland. Temperature-wise these localities did not appear to have differed. These findings allow us to provide information on the paleoecology of the two groups of fossil equids of the Anagni Basin, hereby shedding light on their ecological response to the major climatic events that marked the Early and Middle Pleistocene: the onset of the Pleistocene glaciation, which led to a gradual trend toward cooler climatic conditions and subsequent reduction of sub-tropical vegetation in Eurasia (Fortelius et al., 2006;Combourieu-Nebout et al., 2015), and strengthening of the Early-Middle Pleistocene Transition, which led to an increase of interglacial temperatures (Head and Gibbard, 2015 (Pushkina et al., 2014). Dental mesowear patterns are compatible with a grazing behavior similarly to E. stenonis populations, e.g., from Sésklo (MN17 biozone, Early Pleistocene, Greek Peninsula) (Rivals and Athanassiou, 2008) and from Olivola (∼1.8 Ma, Italian Peninsula) (Strani et al., 2018b). E. stenonis is a common species of both Italian and Greek Peninsulas associated with fossil ungulates showing generally mixed or grazing dietary adaptations in open and dry environments (Rivals and Athanassiou, 2008;Strani et al., 2018b). In contrast, E. senezensis aff. E. sen. stehlini from CSG is associated with artiodactyls showing a wider range of feeding behaviors in a more heterogeneous landscape comprised of wetlands, forests and grasslands (Strani et al., 2015(Strani et al., , 2018c, while displaying a less abrasion-dominated dental mesowear pattern than E. stenonis from Olivola (Figure 4). Similar conditions are observed in East Runton, where E. bressanus displayed FIGURE 6 | δ 13 C (A) and δ 18 O (B) values of the studied taxa compared to fossil equids of central Europe (data from Pushkina et al., 2014). The red line represents the boundary between δ 13 C values that indicate open (δ 13 C > −13‰) or closed landscapes (δ 13 C < −13‰) in warm temperate and cold boreal environments (Bocherens, 2014).
long-term grazing (according to molar mesowear) and shortterm browsing diet (according to tooth microwear) and is grouped with other ungulates with a variety of feeding behaviors (Rivals and Lister, 2016).
The stenonid Equus altidens from the late Early Pleistocene Southern European site of Vallparadís Estació (layer EVT7) also displayed a certain dietary plasticity with a short-term mixed feeding behavior (Strani et al., 2019). E. senezensis aff. E. sen. stehlini is also one of the best represented ungulates from CSG, whereas fewer individuals have been recorded in Olivola (Strani et al., 2018b). This suggests that this equid, while feeding mostly on grasses, may have also accessed more diverse food items leading to a greater intake of softer plant resources compared to the similarly aged E. stenonis from Sésklo and to the later stenonid equids from Olivola. This scenario is consistent with the onset of a gradual trend toward cooler climatic conditions after the initiation at the Plio-Pleistocene boundary (2.6 Ma) of the Quaternary glaciations regulated by a 41 kyr periodicity (Lisiecki and Raymo, 2005), which led to the augmentation of open habitats and the disappearance of sub-tropical vegetation in the European continent (Bertini, 2003;Bertini, 2006, Fortelius et al., 2006Bertini, 2010Bertini, , 2013Kahlke et al., 2011;Combourieu-Nebout et al., 2015). This process probably affected later Northern Europe as the paleoenvironmental conditions of East Runton appear to be more heterogeneous compared to the coeval locality of Olivola.
The more abrasion-dominated mesowear patterns of Equus cf. E. mosbachensis point to a stricter grazing diet, and are comparable to the ones observed in modern zebra populations, which feed almost exclusively on grasses (Bodenstein et al., 2000;Fortelius and Solounias, 2000;Rivals and Semprebon, 2010;Kartzinel et al., 2015). The same trend is observed also in the Middle Pleistocene caballoid equids from other European localities, which also display abrasion-dominated mesowear signatures (Rivals et al., , 2015aRivals and Lister, 2016;Rivals and Ziegler, 2018). A similar trend from less abrasion-dominated to strictly grazing diet is observed even in Early and Middle Pleistocene equid populations of Britain (Rivals and Lister, 2016).
Isotope analysis points instead to a preference for slightly more closed environments in Equus cf. E. mosbachensis from FR similarly to E. mosbachensis from Mauer and Middle and Late Pleistocene equids from Steinheim (Pushkina et al., 2014) in comparison to E. senezensis aff. E. sen. stehlini. Long-term grazing behavior and short-term less abrasive diets are observed in E. ferus from the MIS 11 German localities of Hoppenloch and Steinheim (Rivals and Ziegler, 2018) and from the MIS 9 site Schöningen (Rivals et al., 2015b), suggesting that also in central Europe Middle Pleistocene horses appeared to have successfully exploited less open habitats.
FR Equus cf. E. mosbachensis is associated with ungulates displaying a wide range of feeding behavior with a noticeable abundance of browsing and mixed feedings cervids (Strani et al., 2018a,d). A similar condition is observed in Steinheim, where herbivores with diverse diets in humid and more closed habitats are recorded (Pushkina et al., 2014;Rivals and Ziegler, 2018). Isotopic analysis from the ungulates of Mauer also suggest warm and densely wooded landscape for this locality (Pushkina et al., 2014). On the contrary, fossil horses from the layer G (MIS 12) of Caune de l'Arago Cave (France) and from Hoppenloch are associated with a lower range of dietary traits pointing to drier environmental settings Rivals and Ziegler, 2018). This suggests that Equus cf. E. mosbachensis from the Anagni Basin occupied an ecological niche similar to the fossil horses from Steinheim in comparable environmental conditions, as grazers that feed also in more closed habitats. The higher MWS recorded in the FR population (4.8) compared to that of the Steinheim (3.8) (Figure 5), suggests that the fossil horses of central Italian Peninsula had a stricter grazing behavior compared to central European populations. FR archeopaleontological layer with its fossil assemblage is dated to around "Mid-Brunhes event" (Biddittu et al., 1979;Muttoni et al., 2009;Pereira et al., 2018), which is associated with an increase in the amplitude of the interglacials and glacial cycles (Lang and Wolff, 2011;McClymont et al., 2013), and with one of the warmest and longest interglacial periods of the last 400 ka, registered during MIS 11 (Weirauch et al., 2008). According to palynological data from the MIS 11 Boiano Section (central Italian Peninsula) these warmer conditions led to the expansion of closed canopy forests with high summer precipitations (Combourieu-Nebout et al., 2015). The paleoecological information for the fossil horse of FR, suggests that in this scenario the Middle Pleistocene caballoid equids from central Italian Peninsula adapted to the reduction of open grasslands by exploiting also more closed habitats in heterogeneous landscapes. A possible competition with other large gregarious ungulates, such as Bos primigenius and Stephanorhinus sp., could have also played a key role in the niche occupation of this equid. The higher temperature gradient recorded by the isotopic analysis of the equid dental material compared to those observed in Steinheim, also can imply warmer conditions that characterized the Anagni Basin during this interval. This is concordant with the hypothesis of the Italian Peninsula acting as a southern refugium for mammal taxa and human populations during colder periods and extreme glacials of the Pleistocene (Manzi et al., 2011;Manzi, 2016).

CONCLUSIONS
This work offers important updated information on the paleoecological adaptations of Pleistocene equids of the Italian Peninsula. Overall Early Pleistocene stenonid horses grazed in open landscapes but demonstrated less abrasion-dominated diet compared to the later forms, suggesting the existence of more heterogeneous environmental conditions during the early Early Pleistocene before the Gelasian/Calabrian (middle/late Villafranchian) transition. In the same region, Middle Pleistocene caballoid populations displayed more abrasion-dominated longterm feeding behavior, having occupied less open habitats (at the edge between grasslands and woodlands) than stenonid forms of the Early Pleistocene, both in possibly warmer climatic conditions compared to French and German localities of the same age. Temperature calculations, however, should be considered with caution, since the sample size is quite small and based only on horses, and oxygen stable isotope variation apart from higher temperature can be associated with other factors, e.g., water sources from different elevation and subjected to seasonal variation in water availability.
Obtained data on the niche occupation of the stenonid Equus of the Anagni Basin support the spread of open landscapes in the Italian Peninsula following the initiation of the Pleistocene glacial cycles, suggesting, however, that these equids could have successfully adapted to exploit also more heterogeneous resources in the long-term timescale. Data from the caballoid Equus support the expansion of more closed environments in the region after the MIS 12-11 transition, though horses seemed to feed almost exclusively on abrasive vegetation possibly to avoid competition with other large ungulates that exploited softer food resources.

AUTHOR CONTRIBUTIONS
All authors have provided substantial contribution to the article. FS conceived the ideas and designed the approach. FS collected and elaborated the mesowear data. DP collected and performed the isotopic analysis with the help of HB, FS, and DP led the writing of the manuscript with the support of HB and DD. RS and LB provided helpful comments and inputs for the manuscript. All authors have approved the manuscript and agree with its submission to the journal.

ACKNOWLEDGMENTS
This research has been supported by Avvio alla Ricerca 2017 (Ref. AR11715C7CE78D4F to FS) and the Spanish Agencia Estatal de Investigación (Ministerio de Economía, Industria y Competitividad)-European Regional Development Fund of the European Union (CGL2016-76431-P, AEI/FEDER EU). We thank the Italian Institute of Human Paleontology for preserving the collection and allowing the study of the material object of this investigation. We are also very grateful to Gina Semprebon and Mikael Fortelius for providing the raw mesowear data of modern equids (Equus quagga and Equus grevyi) and other extant ungulates We also acknowledge to Raymond L. Bernor for inviting us to contribute to the Examining Evolutionary Trends in Equus and its Close Relatives from Five Continents Research Topic. Finally, we thank the two reviewers for their helpful comments and suggestions.