Pollen and Molecular Biomarkers from Sedimentary Archives in the Central Po...

[Full title: Pollen and Molecular Biomarkers from Sedimentary Archives in the Central Po Plain (N Italy): Assessing Their Potential to Deepen Changes in Natural and Agricultural Systems]

[Find the meaning and references behind the names: Natural, Pollen, Change, Plant, Doi, June, Aim, Basel, Local, Luigi, Evidence, Clo, Life, Roots, Yvette, Trace, Present, Emilia, Bronze, Cases, Central, Land, Cea, Rosa, Jeremy, Under, High, Age, Warm, Campi, Eleonora, Span, July, Open, Last, Sur, Culture, Assunta, France, Area, Record, Modena, Ledda, Grass, Season, Knowledge, April, Past, Study, Strong, Jacob, Millet, Reggio, Orme, Shown]

Citation: Florenzano, A.; Clò, E.; Jacob, J. Pollen and Molecular Biomarkers from Sedimentary Archives in the Central Po Plain (N Italy): Assessing Their Potential to Deepen Changes in Natural and Agricultural Systems Sustainability 2023 , 15 , 10408. https://doi.org/ 10.3390/su 151310408 Academic Editor: Luigi Ledda Received: 29 April 2023 Revised: 21 June 2023 Accepted: 29 June 2023 Published: 1 July 2023 Copyright: © 2023 by the authors Licensee MDPI, Basel, Switzerland This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https:// creativecommons.org/licenses/by/ 4.0/) sustainability Article Pollen and Molecular Biomarkers from Sedimentary Archives in the Central Po Plain (N Italy): Assessing Their Potential to Deepen Changes in Natural and Agricultural Systems Assunta Florenzano 1, * , Eleonora Cl ò 1 and J é r é my Jacob 2, * 1 Laboratory of Palynology and Paleobotany, Department of Life Sciences, University of Modena and Reggio Emilia, Via G. Campi 287, 41125 Modena, Italy; eleonora.clo@unimore.it 2 Laboratoire des Sciences du Climat et de l’Environnement (UMR CEA-CNRS-UVSQ 8212), Orme des Merisiers, CEA/Saclay B â t. 714, 91191 Gif-sur-Yvette CEDEX, France * Correspondence: assunta.florenzano@unimore.it (A.F.); jeremy.jacob@lsce.ipsl.fr (J.J.) Abstract: This paper proposes to improve the information provided by biological indicators from sedimentary archives by integrating biomolecular techniques and botanical skills. This study represents a first proposal for combining pollen and biomolecular markers to detect land use and improve knowledge of past environmental change drivers. The specific aim of the research is to verify the relationship between miliacin (a pentacyclic triterpene methyl ether, usually interpreted as a broomcorn millet biomarker) and Panicum pollen in three near-site stratigraphic sequences of the Terramara S. Rosa di Poviglio (Po Plain, N Italy). The three cores span the last ~15,000 years and potentially record the beginning of Panicum miliaceum cultivation attested in the area since at least the Bronze Age within the Terramare culture. Despite the fact that Panicum pollen grains were rare in the spectra and miliacin was barely detectable in most of the 31 samples selected for biomolecular analyses, their combined evidence testifies to the local presence of the plant Panicum pollen and sedimentary miliacin suggest the adoption of millet crops during the Recent Bronze Age by the Terramare culture, when climatic instability led to the diversification of cereal crops and the shift to drought-tolerant varieties Keywords: Panicum ; miliacin; pollen; millet crops; palaeoenvironment; archaeological site; pentacyclic triterpenes methyl ethers 1. Introduction Molecular biomarkers detected in sedimentary archives are being increasingly utilised as proxies for examining past environmental conditions and anthropic impacts (e.g., [ 1 – 5 ]). These organic compounds possess, in specific cases, high taxonomic specificity and are reputedly resistant to biodegradation [ 6 ], providing information on both the overall ecological conditions and the local presence of species. Tricyclic diterpenes are specific molecular biomarkers of vascular plants, mostly produced by gymnosperms, and pentacyclic triterpenes are mainly derived from angiosperms [ 7 , 8 ]. Among them, pentacyclic triterpene methyl ethers (PTMEs) have been shown to be specific mostly to Poaceae [ 9 ], whereas triterpenyl acetates are mostly produced by Asteraceae [ 10 ]. These plant-derived biomolecules could be particularly informative of the past local presence of a plant, even in the absence of any vegetal macroremains in the sedimentary archives Recent research is focused on the detection of molecular biomarkers specific to Panicum in sedimentary archives as traces of the past presence of millets in paleosols and sediments. Increased attention has been paid to the studies of miliacin, PTME specific to Panicum , and some other Poaceae. Miliacin is mainly concentrated in seeds, and it is also normally present in leaves and stems but absent from hulls and roots [ 11 ]. In sedimentary deposits, the detection of miliacin is commonly regarded as trace evidence of the past presence of the broomcorn millet [ 11 – 14 ]. Panicum miliaceum L. (Poaceae) is a C 4 warm-season annual grass Sustainability 2023 , 15 , 10408. https://doi.org/10.3390/su 151310408 https://www.mdpi.com/journal/sustainability

[Find the meaning and references behind the names: Foods, Transport, Path, Low, Europe, Hot, Poor, Sole, Topic, Show, Field, Time, China, Tool, Dry, Data, Hand, Situ, Few, Evolution, Due, Cover, Still, Arid, Self, Wind, General, Small, Ams, Need, Lake, Cal, Early, Mid, Short]

Sustainability 2023 , 15 , 10408 2 of 17 that is highly tolerant of poor soils, drought, and high temperatures. This small-grained cereal is used in foods and as a forage plant as it requires very little water and has a short growing cycle. Today, it is extensively cultivated in arid regions, which have poor soils and hot, dry weather [ 15 ]. The earliest finds of P. miliaceum were recorded in China in the form of charred grains from Early Neolithic sites dated to the early–mid 6 th millennium BC [ 16 ]. Archaeobotanical research indicated that millet was widely spread in Europe by the 2 nd millennium BC [ 17 – 19 ]. Uncertainties about the timing and diffusion path of its adoption by European societies were mostly resolved thanks to an extensive program of AMS dating of charred broomcorn millet grains from prehistoric sites in Europe [ 20 ]. The detection of P. miliaceum caryopses in archaeobotanical records, which is crucial for the reconstruction of the time and spread of its cultivation, could also be used to support the interpretation of miliacin as a plant-specific biomarker. In fact, some controversy remains over its interpretation as evidence of broomcorn millet since miliacin is produced by a variety of plants (53 species), mainly Poaceae (45) [ 21 ]. Recent biomolecular studies [ 9 , 11 , 13 ] have restricted the production of miliacin, in combination with other PTMEs, to a few plants belonging to Panicum and Pennisetum species. Nevertheless, if miliacin is found as the sole PTME and in high concentration, it is very likely that P. miliaceum is the source [ 11 ]. Nevertheless, the need for an unequivocal association with P. miliaceum alone remains Actually, P. miliaceum caryopses are hardly preserved in archaeological deposits due to biodegradative or taphonomic processes. Additionally, the identification of broomcorn millet in phytolith records [ 22 , 23 ] can support the attribution of miliacin as a biomarker specific to P. miliaceum , but this is still an underexplored research field in European countries (e.g., [ 24 ]). Furthermore, when present, both of these P. miliaceum macroand microremains show evidence of the use of this plant at a site but not necessarily its local cultivation When coupled with pollen, could miliacin be considered reliable evidence of broomcorn millet growing locally? Indeed, palynology is the most powerful tool to reconstruct vegetation cover and environmental changes, providing detailed information on vegetation composition and distribution and allowing for the identification of anthropogenic environments and cultural landscape evolution (e.g., [ 25 – 29 ]). Especially when related to archaeological contexts, pollen evidence could elucidate aspects of agricultural systems, land-use types, and crop cultivation (e.g., [ 30 – 32 ]). Panicum is mainly a self-pollinating plant and a low-pollen producer [ 33 ]. When found in pollen samples from sedimentary contexts, millet pollen grains confirm the local past presence of the plant Biomolecular investigations integrated with pollen studies are innovative and yet rarely applied methods for certifying the presence of broomcorn millet and its use in the agricultural activities of an archaeological site. In stratigraphic sequence samples, the comparison of miliacin concentrations and pollen data for Panicum could be particularly informative regarding this topic. This is conditioned by several factors that affect both the production, transport, and archiving of these two tracers from the same plant. On the one hand, in P. miliaceum , miliacin is very abundant only in seeds (0.5‰ weight; [ 11 ]), and, in general, its detection is dependent on the site. For example, detecting miliacin in lake sediments supposes that a sufficient (depending on the lake’s catchment area and surface area) number of seeds were produced in the catchment and that the mechanism of transportation of miliacin from the field to the lake sediment (more plausibly after ingestion and excretion by humans) was efficient [ 2 ]. On the other hand, millet pollen is transported by the wind a short distance from the source plant and attests to its presence in situ [ 34 ]. This paper aims to verify the relationship between miliacin and Panicum pollen content in the near-site cores of the Terramara S. Rosa di Poviglio (3500–3120 cal yr BP, Po Plain, N Italy), where sediment records continuously span at least the last ~15,000 years and provide detailed and high-resolution information on vegetation history and land-use changes in the area [ 35 ]. Pollen analyses at the Terramara Bronze Age sites (3500-1950 cal yr BP) in the Po Plain have indicated the presence of Panicum during increasing climate aridity [ 36 – 38 ], and the first biomolecular studies attest to traces of miliacin in contexts inside Terramara

[Find the meaning and references behind the names: Santa Rosa, Forest, Map, Wine, Day, Rivers, Wood, Santa, Peer, Long, Large, Woods, Areas, Piccolo, Century, Far, Mean, River, Lat, Semi, Milk, Grande, Canals, Location, Springs, Common]

Sustainability 2023 , 15 , 10408 3 of 17 settlements [ 39 ]. The integration of biomolecular techniques and pollen analyses may improve the information on broomcorn millet adoption by the Terramare culture 2. Materials and Methods 2.1. Study Area The Terramara Santa Rosa di Poviglio (lat. 44 ◦ 52 0 21 00 N; long. 10 ◦ 34 0 3 00 E) is located in the alluvial plain of the Po River, about 3 km southward from the present-day course of the river (Figure 1 ). The climate is semi-continental, with mean annual temperatures of 12–14 ◦ C and mean annual rainfall between 600 and 800 mm. This area was occupied by a swamp until the 15 th century [ 40 ], and even today, it is poorly drained. The territory is characterized by intensive industrial agriculture with the production of crops (mainly cereals for food, tomatoes for industrial use, and Medicago for fodder), orchards (Rosaceae– Prunoideae, including cherries, plums, apricots, and grapes for wine production), and animal breeding (especially pigs for meat and cattle for milk). Hygrophilous forests are spread near rivers, canals, and springs, while sparse woods consisting of mixed oakwood are the most common forest cover, as attested from past pollen diagrams in the region [ 41 – 44 ]. Almost 40 years of archaeological excavations of the site and a geophysical survey [ 45 ] have revealed that the archaeological structures of the Bronze Age settlement consist of an earthen rampart, wood fences, and two dwelling areas named “Villaggio Piccolo”— “Small Village” (PVG-VP) and “Villaggio Grande”— “Large Village” (PVG-VG; Figure 1 ) [ 37 ] Sustainability 2023 , 15 , x FOR PEER REVIEW 3 of 19 changes in the area [35]. Pollen analyses at the Terramara Bronze Age sites (3500-1950 cal yr BP) in the Po Plain have indicated the presence of Panicum during increasing climate aridity [36 – 38], and the first biomolecular studies attest to traces of miliacin in contexts inside Terramara settlements [39]. The integration of biomolecular techniques and pollen analyses may improve the information on broomcorn millet adoption by the Terramare culture. 2. Materials and Methods 2.1. Study Area The Terramara Santa Rosa di Poviglio (lat. 44 ° 5 2′ 21 ″ N; long. 10 ° 3 4′3″ E) is located in the alluvial plain of the Po River, about 3 km southward from the present-day course of the river (Figure 1). The climate is semi-continental, with mean annual temperatures of 12 – 14 °C and mean annual rainfall between 600 and 800 mm. This area was occupied by a swamp until the 15 th century [40], and even today, it is poorly drained. The territory is characterized by intensive industrial agriculture with the production of crops (mainly cereals for food, tomatoes for industrial use, and Medicago for fodder), orchards (Rosaceae – Prunoideae, including cherries, plums, apricots, and grapes for wine production), and animal breeding (especially pigs for meat and cattle for milk). Hygrophilous forests are spread near rivers, canals, and springs, while sparse woods consisting of mixed oakwood are the most common forest cover, as attested from past pollen diagrams in the region [41 – 44]. Almost 40 years of archaeological excavations of the site and a geophysical survey [45] have revealed that the archaeological structures of the Bronze Age settlement consist of an earthen rampart, wood fences, and two dwelling areas named “ Villaggio Piccolo ”— “ Small Village ” (PVG-VP) and “ Villaggio Grande ”— “ Large Village ” (PVG-VG; Figure 1) [37]. Figure 1. Location map of the Terramara Santa Rosa di Poviglio (PVG) site (VP = “ Villaggio Piccolo ”— “ Small Village ” ; VG = “ Villaggio Grande ”— “ Large Village ” ) in the central Po Plain, N Italy, and position of the three cores (N-S 3, C-S 1, and F-S 2) drilled in a SW-NE direction at ~150 m, ~320 m, and ~525 m, respectively, from the Terramara (GoogleEarth ™ ). 2.2. Sediment Samples Figure 1. Location map of the Terramara Santa Rosa di Poviglio (PVG) site (VP = “Villaggio Piccolo”— “Small Village”; VG = “Villaggio Grande”— “Large Village”) in the central Po Plain, N Italy, and position of the three cores (N-S 3, C-S 1, and F-S 2) drilled in a SW-NE direction at ~150 m, ~320 m, and ~525 m, respectively, from the Terramara (GoogleEarth™) 2.2. Sediment Samples During the excavation campaign carried out in 2018, three continuous near-site cores (Near-S 3, Central-S 1, and Far-S 2) were collected in a SW-NE transect at different distances (~150 m, ~320 m, and ~525 m) from the Terramara S. Rosa di Poviglio [ 35 ] (Figure 1 ). A total of 292 pollen samples were collected (N-S 3: 77 samples; C-S 1: 150 samples; F-S 2: 65 samples) at about 10–30 cm intervals, according to the stratigraphy and archaeological phases. These

[Find the meaning and references behind the names: Sand, Step, Range, Northern, Sheets, Csa, Size, Middle, Sample, Clay, Core]

Sustainability 2023 , 15 , 10408 4 of 17 cores span about ~15,000 years; the chronology relies on radiocarbon dates obtained from organic sediments collected along the stratigraphic sequence of the N-S 3 and C-S 1 cores and ceramic fragments (Figure 2 ). The organic layer of the N-S 3 core is archaeologically dated to the Middle Bronze Age (BM) and stratigraphically correlates to the occupation layers of the “Villaggio Piccolo” stratigraphy that has been archaeologically and radiometrically dated to a well-defined window of the Italian Middle Bronze Age [ 46 ]. The F-S 2 core has been stratigraphically correlated to the other two cores Sustainability 2023 , 15 , x FOR PEER REVIEW 4 of 19 During the excavation campaign carried out in 2018, three continuous near-site cores (Near-S 3, Central-S 1, and Far-S 2) were collected in a SW-NE transect at different distances (~150 m, ~320 m, and ~525 m) from the Terramara S. Rosa di Poviglio [35] (Figure 1). A total of 292 pollen samples were collected (N-S 3: 77 samples; C-S 1: 150 samples; F-S 2: 65 samples) at about 10 – 30 cm intervals, according to the stratigraphy and archaeological phases. These cores span about ~15,000 years; the chronology relies on radiocarbon dates obtained from organic sediments collected along the stratigraphic sequence of the N-S 3 and C-S 1 cores and ceramic fragments (Figure 2). The organic layer of the N-S 3 core is archaeologically dated to the Middle Bronze Age (BM) and stratigraphically correlates to the occupation layers of the “ Villaggio Piccolo ” stratigraphy that has been archaeologically and radiometrically dated to a well-defined window of the Italian Middle Bronze Age [46]. The F-S 2 core has been stratigraphically correlated to the other two cores. In 2020, a subsampling step was performed for biomarker analysis with a selection of 31 samples mainly related to Holocene chronology. Thick sediment discs ranging 1 – 2 cm in size were taken in correspondence with the 2018 pollen sampling, and then a subsampling of 1 cm 3 for each disc was taken and wrapped into aluminum sheets to avoid any contamination with plastic materials. Figure 2. Near-site cores (N-S 3, C-S 1, and F-S 2) of the Terramara S. Rosa di Poviglio: ( a ) stratigraphy and uncalibrated age BP marked by *; C = clay, Si = silt, Sa = sand, CSa = coarse sand; ( b ) AMS – 14 C ages for N-S 3 and C-S 1; the calibration into calendar years (2σ range) was calculated with R Clam using IntCal 20 Northern Hemisphere 14 C calibration curves [47]. 2.3. Pollen Treatment and Analysis Figure 2. Near-site cores (N-S 3, C-S 1, and F-S 2) of the Terramara S. Rosa di Poviglio: ( a ) stratigraphy and uncalibrated age BP marked by *; C = clay, Si = silt, Sa = sand, CSa = coarse sand; ( b ) AMS– 14 C ages for N-S 3 and C-S 1; the calibration into calendar years (2 σ range) was calculated with R Clam using IntCal 20 Northern Hemisphere 14 C calibration curves [ 47 ]. In 2020, a subsampling step was performed for biomarker analysis with a selection of 31 samples mainly related to Holocene chronology. Thick sediment discs ranging 1–2 cm in size were taken in correspondence with the 2018 pollen sampling, and then a subsampling of 1 cm 3 for each disc was taken and wrapped into aluminum sheets to avoid any contamination with plastic materials 2.3. Pollen Treatment and Analysis Pollen was extracted from 195 selected samples (N-S 3: 77 samples; C-S 1: 90 samples; F-S 2: 28 samples), according to the routine method in use in the Laboratory of Palynology and Palaeobotany of Modena-LPP [ 48 , 49 ]. About 2 g of sediment from each sample were treated through sieving and heavy-liquid flotation with sodium metatungstate hydrate.

[Find the meaning and references behind the names: Keys, Sum, Jelly, Plus, Wild, Avena, Gram, Zea, Mays, Factor, Pre]

Sustainability 2023 , 15 , 10408 5 of 17 Lycopodium spore tablets were added to calculate the concentration (pollen per gram = p/g) Pollen was identified at 1000 × magnification with the help of atlases/keys (e.g., [ 50 , 51 ]) and the LPP reference pollen collection Particular attention was paid to the identification of different types of Poaceae pollen grains (Figure 3 ), which are divided into morphological groups according to the literature [ 52 , 53 ], with the correction factor for glycerol jelly [ 25 ]. Pollen grains smaller than 37 µ m (only a few wild grasses are larger) with annulus diameters smaller than 8 µ m were included in the wild grass group [ 53 , 54 ]. According to Beug [ 52 ], cereal pollen is larger than 37 µ m. The Hordeum group includes pollen between 37 and 50 µ m in size with an annulus diameter between 8 and 12 µ m, while the Avena/Triticum group includes pollen larger than 50 µ m, with a mean annulus diameter larger than 12 µ m [ 25 , 52 ]. The presence of Zea mays (mean pollen size larger than 60 µ m) was identified according to Beug [ 55 ]. Morphological characterization of Panicum pollen can be difficult, as it includes pollen with medium to large sizes (32–62 µ m in diameter) and porus plus annulus between 8 and 13 µ m; as such, it overlaps, especially with the Hordeum group. However, it could be identified by the presence of insulate exine, which is absent in the Hordeum group [ 52 , 53 , 55 , 56 ]. Different species of Panicum and other genera (e.g., Setaria and Echinochloa ) could be included in the same pollen type. Therefore, the term Panicum s.l. was used to indicate this pollen type, which includes both cultivated millet ( P. miliaceum ) and wild grasses. In this paper, “wild cereals” refer to large Poaceae pollen (including the Avena/Tricticum group, Hordeum group, and Panicum s.l.) observed in pre-Holocene samples Sustainability 2023 , 15 , x FOR PEER REVIEW 5 of 19 Pollen was extracted from 195 selected samples (N-S 3: 77 samples; C-S 1: 90 samples; F-S 2: 28 samples), according to the routine method in use in the Laboratory of Palynology and Palaeobotany of Modena-LPP [48,49]. About 2 g of sediment from each sample were treated through sieving and heavy-liquid flotation with sodium metatungstate hydrate. Lycopodium spore tablets were added to calculate the concentration (pollen per gram = p/g). Pollen was identified at 1000× magnification with the help of atlases/keys (e.g., [50,51]) and the LPP reference pollen collection. Particular attention was paid to the identification of different types of Poaceae pollen grains (Figure 3), which are divided into morphological groups according to the literature [52,53], with the correction factor for glycerol jelly [25]. Pollen grains smaller than 37 μm (only a few wild grasses are larger) with annulus diameters smaller than 8 μm were included in the wild grass group [53,54]. According to Beug [52], cereal pollen is larger than 37 μm. The Hordeum group includes pollen between 37 and 50 μm in size with an annulus diameter between 8 and 12 μm, while the Avena/Triticum group includes pollen larger than 50 μm , with a mean annulus diameter larger than 12 μm [25, 52]. The presence of Zea mays (mean pollen size larger than 60 μm) was identified according to Beug [55]. Morphological characterization of Panicum pollen can be difficult, as it includes pollen with medium to large sizes (32 – 62 µ m in diameter) and porus plus annulus between 8 and 13 µ m; as such, it overlaps, especially with the Hordeum group. However, it could be identified by the presence of insulate exine, which is absent in the Hordeum group [52,53,55,56]. Different species of Panicum and other genera (e.g., Setaria and Echinochloa ) could be included in the same pollen type. Therefore, the term Panicums .l. was used to indicate this pollen type, which includes both cultivated millet ( P. miliaceum ) and wild grasses. In this paper, “ wild cereals ” refer to large Poaceae pollen (including the Avena/Tricticum group, Hordeum group, and Panicums .l.) observed in pre-Holocene samples. Figure 3. Poaceae pollen grains: ( a ) wild grass group (23 μ m); ( b ) Panicums .l. (43 μ m) and details of its verrucate exine; ( c ) Hordeum group (40 μ m); ( d ) Avena/Triticum group (51 μ m); ( e ) Zea mays (101 μ m). The scale is 10 µ m. Figure 3. Poaceae pollen grains: ( a ) wild grass group (23 µ m); ( b ) Panicum s.l. (43 µ m) and details of its verrucate exine; ( c ) Hordeum group (40 µ m); ( d ) Avena/Triticum group (51 µ m); ( e ) Zea mays (101 µ m). The scale is 10 µ m On average, about 350 pollen grains per sample were counted [ 35 ]. Percentages were calculated based on a pollen sum, including both identified and unidentified pollen grains. Only data on Poaceae (wild grass group, cereals— Avena/Triticum and Hordeum groups, and Panicum s.l.) were reported in this paper.

[Find the meaning and references behind the names: Mode, Carrier, Better, Gas, Standard, Final, Set, Prior, Hold, Api, Held, Mass, Nap, Energy, Artemisia, Non, Rather, Lower, Half, Min, Peak, Lesser]

Sustainability 2023 , 15 , 10408 6 of 17 Pollen diagrams, drawn with Tilia 3.0.1 [ 57 ], also show selected pollen curves as follows: Arboreal Pollen/Non-Arboreal Pollen (AP/NAP), hygrophilous herbs (Cyperaceae, Juncus type, Phragmites , Sparganium emersum type, and Typha latifolia type), and Anthropogenic Pollen Indicators—API group ( Artemisia , Centaurea , Cichorieae, Plantago , Trifolium type, Urtica , and cereals; [ 58 ]). 2.4. Biomolecular Analysis The lipid extraction and separation were based on Jacob et al. [ 9 ]. The samples were crushed to powder after being dried in an oven at 50 ◦ C for 48 h. About 2–3 g of sample powder were ultrasonically extracted with a mixture of dichloromethane and isopropanol (2:1 v : v ), which was repeated three times. After being combined, the extracts were separated into neutral and acidic fractions by solid-phase extraction using silica bonded with aminopropyl groups. The neutral compounds, including PTMEs and triterpenyl acetates, were eluted with dichloromethane and isopropanol (2:1 v : v ). Additionally, 5 α -cholestane was added as a quantitation standard prior to the analysis. The identification and quantification analyses were performed using a gas chromatography instrument (GC; Trace 1300, ThermoScientific, Waltham, MA, USA) coupled to a mass spectrometer (MS; ISQ QD, ThermoScientific). The gas chromatograph was fitted with a BPX 5 capillary column (60 m; 0.32 mm i.d.; 0.25 µ m film thickness). The GC operating conditions were as follows: the temperature was held at 80 ◦ C for 1 min, then increased from 40 to 120 ◦ C at 30 ◦ C min − 1 and 120 to 310 ◦ C at 3 ◦ C min − 1 , with a final isothermal hold at 310 ◦ C over 40 min The sample was injected splitless, with the injector temperature set at 280 ◦ C. Helium was the carrier gas at 1 mL/min. The mass spectrometer was operated in the electron ionization (EI) mode at 70 eV ionization energy and scanned from 50 to 650 Da. PTMEs, such as miliacin and triterpenyl acetates, were identified using authentic standards [ 10 , 11 ]. The concentrations were determined by comparing the peak areas of specific biomarkers on the GC trace with those of 5 α -cholestane and reported to the mass of the sample to provide a concentration value. Replicate injections of some samples lead to a reproducibility better than 3 ng/g 3. Results 3.1. Pollen of P. Miliaceum and Other Poaceae Poaceae pollen grains are ubiquitous in the three cores. The Poaceae wild grass group is dominant among herbs (N-S 3: 15.3%, C-S 1: 12.6%, and F-S 2: 12%, on average; NAP values: N-S 3: 54.7%, C-S 1: 56%, and F-S 2: 50.2%). Among the cereals, the Hordeum group is the most common in the pollen spectra, and it is present in at least half of the samples (N-S 3: 0.8%, C-S 1: 0.9%, and F-S 2: 0.6%) In the C-S 1 core, the almost continuous occurrence of high values of Panicum s.l. in the pre-Holocene samples (up to 3.3% at a depth of 846 cm) is probably not indicative of the presence of P. miliaceum but rather wild grasses producing large pollen (e.g., Echinochloa , Glyceria , and Setaria ), as no pre-Holocene broomcorn millet (or any other) cultivation is attested. However, in lower percentages, Panicum s.l. pollen is also present in the pre- Holocene levels of the N-S 3 (0.3% at 1013 cm) and F-S 2 cores (0.9% at 830 cm). Similarly, the Hordeum group (C-S 1: up to 4.6% at 825 cm) and, to a lesser extent, the Avena/Triticum group (N-S 3: up to 0.6% at 1165 cm) also show significant percentages in pre-Holocene levels, both in N-S 3 and C-S 1 Even for the Holocene levels, the C-S 1 core is the one with the highest presence of Panicum s.l. pollen, although always in lower percentages compared to the samples described above. The Hordeum group and, to a lesser extent, the Avena/Triticum group have been observed continuously during the Middle–Late Neolithic period (ca. 7050/6450- 5450 cal yr BP) and the pre-Bronze Age in all cores. The higher percentages of cereals characterize the Middle Bronze Age levels ( Hordeum group: 0.8% in N-S 3, 1.2% in C-S 1, and 1.4% in F-S 2; Avena/Triticum group: 0.4% in N-S 3, 0.1% in C-S 1, and 0.7% in F-S 2; mean values). The peak of the Hordeum group is reached in the N-S 3 core (2.2% at 128 cm), while

[Find the meaning and references behind the names: Bauer, Level, Reach, Date, Barley, Wide, Original, Back, Wheat, Ion, Ones]

Sustainability 2023 , 15 , 10408 7 of 17 the highest value of the Avena/Triticum group is attested in the F-S 2 core (1.4% at 281 cm) The traces of Panicum s.l. pollen in the C-S 1 core (0.3% at 177 cm—combined evidence of miliacin up to 17.2 ng/g) indicate an increased crop diversity (mainly wheat and barley until then). The increase in cereals is evident in the most recent levels of the cores, with especially high percentages for the Hordeum group (e.g., N-S 3: 5.6% at 99 cm and 4.8% at 45 cm; C-S 1: 4.9% at 81 cm and 6.2% at 61 cm; F-S 2: 1% at 73 cm) Zea mays pollen is present in the C-S 1 core (0.8% at 61 cm and 0.3% at 81 cm) 3.2. Molecular Biomarkers Apart from classical molecular biomarkers, such as n -alkanes and n -alcohols, the neutral fraction of the lipid extracts of Poviglio samples contained a wide variety of pentacyclic triterpenes, of which several original ones could be detected and identified (e.g., Figure 4 ). Miliacin (olean-18-en-3 β -ol methyl ether) and crusgallin (taraxer-14-en-3 β -ol methyl ether) were the only PTMEs detected. Pentacyclic triterpenyl acetates specific for Asteraceae were also found in some samples. These comprised taraxer-14-en-3 β -yl (taraxeryl), olean-12-en-3 β -yl ( β -amyrenyl), olean-13(18)-en-3 β -yl (germanicyl), bauer-8- en-3 β -yl (isobauerenyl), bauer-7-en-3 β -yl (bauerenyl), pichier-8-en-3 β -yl (isopichierenyl), and pichier-9(11)-en-3 β -yl (pichierenyl) acetates Sustainability 2023 , 15 , x FOR PEER REVIEW 7 of 19 been observed continuously during the Middle – Late Neolithic period (ca. 7050/6450-5450 cal yr BP) and the pre-Bronze Age in all cores. The higher percentages of cereals characterize the Middle Bronze Age levels ( Hordeum group: 0.8% in N-S 3, 1.2% in C-S 1, and 1.4% in F-S 2; Avena/Triticum group: 0.4% in N-S 3, 0.1% in C-S 1, and 0.7% in F-S 2; mean values). The peak of the Hordeum group is reached in the N-S 3 core (2.2% at 128 cm), while the highest value of the Avena/Triticum group is attested in the F-S 2 core (1.4% at 281 cm). The traces of Panicums .l. pollen in the C-S 1 core (0.3% at 177 cm — combined evidence of miliacin up to 17.2 ng/g) indicate an increased crop diversity (mainly wheat and barley until then). The increase in cereals is evident in the most recent levels of the cores, with especially high percentages for the Hordeum group (e.g., N-S 3: 5.6% at 99 cm and 4.8% at 45 cm; C-S 1: 4.9% at 81 cm and 6.2% at 61 cm; F-S 2: 1% at 73 cm). Zea mays pollen is present in the C-S 1 core (0.8% at 61 cm and 0.3% at 81 cm). 3.2. Molecular Biomarkers Apart from classical molecular biomarkers, such as n -alkanes and n -alcohols, the neutral fraction of the lipid extracts of Poviglio samples contained a wide variety of pentacyclic triterpenes, of which several original ones could be detected and identified (e.g., Figure 4). Miliacin (olean-18-en- 3β -ol methyl ether) and crusgallin (taraxer-14-en- 3β -ol methyl ether) were the only PTMEs detected. Pentacyclic triterpenyl acetates specific for Asteraceae were also found in some samples. These comprised taraxer-14-en- 3β -yl (taraxeryl), olean-12-en- 3β -yl ( β -amyrenyl), olean-13(18)-en- 3β -yl (germanicyl), bauer-8-en- 3β - yl (isobauerenyl), bauer-7-en- 3β -yl (bauerenyl), pichier-8-en- 3β -yl (isopichierenyl), and pichier-9(11)-en- 3β -yl (pichierenyl) acetates. Figure 4. Partial chromatograms of the neutral fraction of the lipid extract of sample M 04-S 3 (113 cm). ( a ) Total Ion Current chromatogram showing the peak of friedelin (0); ( b ) m/z 425 + 440 ionspecific chromatogram revealing the peaks of pentacyclic triterpene methyl ethers (1 — crusgallin; 2 — miliacin); ( c ) m/z 453 + 468 ion-specific chromatogram illustrating the peaks of pentacyclic triterpenyl acetates (3 — taraxeryl acetate; 4 —β -amyrenyl acetate; 5 — germanicyl acetate; 6 — isobauerenyl acetate; 7 — bauerenyl acetate; 8 — isopichierenyl acetate; 9 — pichierenyl acetate). Up until now, bauerenyl, pichierenyl, and isopichierenyl acetates were only described in Picris hieracioides (Asteraceae) [59 – 62]; thus, they are specific at a species level. Figure 4. Partial chromatograms of the neutral fraction of the lipid extract of sample M 04-S 3 (113 cm) ( a ) Total Ion Current chromatogram showing the peak of friedelin (0); ( b ) m/z 425 + 440 ion-specific chromatogram revealing the peaks of pentacyclic triterpene methyl ethers (1—crusgallin; 2—miliacin); ( c ) m/z 453 + 468 ion-specific chromatogram illustrating the peaks of pentacyclic triterpenyl acetates (3—taraxeryl acetate; 4— β -amyrenyl acetate; 5—germanicyl acetate; 6—isobauerenyl acetate; 7—bauerenyl acetate; 8—isopichierenyl acetate; 9—pichierenyl acetate) Up until now, bauerenyl, pichierenyl, and isopichierenyl acetates were only described in Picris hieracioides (Asteraceae) [ 59 – 62 ]; thus, they are specific at a species level. Apart from Lavrieux [ 63 ], this study constitutes their first detection in sediments. The samples also contained friedelin (friedelan-3-one), a pentacyclic triterpene that is found in many angiosperms but is not specific at the genus or even species level Miliacin concentrations are very low (<4 ng/g), below 250 cm, in the N-S 3 core (pre- Bronze Age), and then increase upward to reach 6 to 15 ng/g in Bronze Age sediments. The same pattern is observed in the C-S 1 core, with concentrations maximizing at 18 ng/g at 177 cm, i.e., for samples deposited during the Bronze Age (Figure 5 ). Crusgallin was only detected in one sample from N-S 3 and two samples from C-S 1, all of which date back to the Bronze Age. Only bauerenyl, pichierenyl, and isopichierenyl acetates were quantified.

[Find the meaning and references behind the names: New, Human, Sake, Key, Bot, Mark, Banks, Main, Marsh, Orange]

Sustainability 2023 , 15 , 10408 8 of 17 Significant concentrations (>5 ng/g) in bauerenyl, pichierenyl, and isopichierenyl acetates were only found in samples corresponding to the Bronze Age Sustainability 2023 , 15 , x FOR PEER REVIEW 8 of 19 Apart from Lavrieux [63], this study constitutes their first detection in sediments. The samples also contained friedelin (friedelan-3-one), a pentacyclic triterpene that is found in many angiosperms but is not specific at the genus or even species level. Miliacin concentrations are very low (<4 ng/g), below 250 cm, in the N-S 3 core (pre- Bronze Age), and then increase upward to reach 6 to 15 ng/g in Bronze Age sediments. The same pattern is observed in the C-S 1 core, with concentrations maximizing at 18 ng/g at 177 cm, i.e., for samples deposited during the Bronze Age (Figure 5). Crusgallin was only detected in one sample from N-S 3 and two samples from C-S 1, all of which date back to the Bronze Age. Only bauerenyl, pichierenyl, and isopichierenyl acetates were quantified. Significant concentrations (>5 ng/g) in bauerenyl, pichierenyl, and isopichierenyl acetates were only found in samples corresponding to the Bronze Age. Figure 5. Near-site cores (N-S 3, C-S 1, and F-S 2) of the Terramara S. Rosa di Poviglio: concentration of miliacin (ng/g) and Panicums .l. pollen percentages from the 31 samples studied for lipid biomarkers. The selected samples come from the Holocene levels; the orange bands mark the Bronze Age phase in each core. 4. Discussion The integrated analysis of pollen and molecular biomarkers from the near-site cores of the Terramara S. Rosa di Poviglio provided new data for improving knowledge on agriculture and human impact in the central Po Plain. Despite the low concentrations of miliacin along the cores, its specific presence in key phases of the local environmental transformation described by pollen records may contribute to detailing the reconstruction of the spread of Panicum miliaceum in N Italy. For the sake of clarity, in the following section, the pollen and biomarker samples are indicated with the depth of sampling. 4.1. The Dualistic Interpretation of the Open Environment Inferred from Poaceae Pollen Curves Pollen analyses from the N-S 3, C-S 1, and F-S 2 cores provide information on vegetation dynamics in the central Po Plain from a long-term perspective (from at least the last ~15,000 years; [35]). According to the stratigraphy, radiocarbon dates, and possible comparison with the main reference pollen records for the study area (e.g., [64 – 66]), the bot- Figure 5. Near-site cores (N-S 3, C-S 1, and F-S 2) of the Terramara S. Rosa di Poviglio: concentration of miliacin (ng/g) and Panicum s.l. pollen percentages from the 31 samples studied for lipid biomarkers The selected samples come from the Holocene levels; the orange bands mark the Bronze Age phase in each core 4. Discussion The integrated analysis of pollen and molecular biomarkers from the near-site cores of the Terramara S. Rosa di Poviglio provided new data for improving knowledge on agriculture and human impact in the central Po Plain. Despite the low concentrations of miliacin along the cores, its specific presence in key phases of the local environmental transformation described by pollen records may contribute to detailing the reconstruction of the spread of Panicum miliaceum in N Italy. For the sake of clarity, in the following section, the pollen and biomarker samples are indicated with the depth of sampling 4.1. The Dualistic Interpretation of the Open Environment Inferred from Poaceae Pollen Curves Pollen analyses from the N-S 3, C-S 1, and F-S 2 cores provide information on vegetation dynamics in the central Po Plain from a long-term perspective (from at least the last ~15,000 years; [ 35 ]). According to the stratigraphy, radiocarbon dates, and possible comparison with the main reference pollen records for the study area (e.g., [ 64 – 66 ]), the bottom of the cores intercepted the pre-Holocene phases. Pollen spectra from the pre-Holocene levels of the cores indicate the presence of an open vegetation dominated by Poaceae (Figure 6 ). According to the presence of plants in wet environments, the occurrence of increasing values of Panicum s.l. at the bottom of the cores is probably not indicative of the presence of Panicum itself but rather wild grasses producing large pollen (e.g., Echinochloa , Glyceria , and Setaria ), which are typical of wet and ruderal environments, such as banks, riverbeds, and swamps, and warm temperatures [ 67 ]. Marsh and alluvial plain deposits from almost the whole Late Pleniglacial/Lateglacial have also been identified at other sites in the lower central Po Plain [ 68 ].

Sustainability 2023 , 15 , 10408 9 of 17 1 Figure 6. Near-site cores (N-S 3, C-S 1, and F-S 2) of the Terramara S. Rosa di Poviglio: percentage pollen diagram of selected sums for palaeoenvironmental reconstruction (AP/NAP, hygrophilous herbs, API group [ 58 ], wild and cultivated Poaceae). Enhanced curves × 10.

[Find the meaning and references behind the names: Empire, Rapid, Own, Villa, Western, Fields, Lakes, Loss, Six, Part, Meadows, Roman, Romagna, End, Worth, Cremaschi]

Sustainability 2023 , 15 , 10408 10 of 17 Pollen records indicate that vegetation changes in the early Holocene were mainly related to climate oscillations, and fairly open riparian communities (dominated by hygrophilous herbs: Cyperaceae, Sparganium emersum type, and Typha latifolia type) characterized the study area [ 69 ]. The open landscape is also indicated in the upper part of the cores and is associated with the increasing pressure of human activities in the area, as evidenced by the rapid and significant increase in Anthropogenic Pollen Indicators (API group [ 58 ] in Figure 6 ). Palynological data show vegetation changes linked to the interplay between climatic and human forces acting on the environment. The openness of the landscape strongly increased toward the Bronze Age, as indicated by an expansion of the space dedicated to cereal crops [ 37 ]. Cereals were mainly barley (represented by the Hordeum group), but wheat ( Avena/Triticum group) was also cultivated in the surrounding area (Figure 6 ). Pollen curves suggest that the distribution of barley fields was almost ubiquitous, while wheat is documented mainly in the N-S 3 and F-S 2 cores. In the C-S 1 core, the presence of Panicum s.l. could be an indication that crops were probably more diversified According to Cremaschi et al. [ 37 ], broomcorn millet spread during the Bronze Age in the Po Plain, and it was probably introduced by the population of the Terramara S. Rosa di Poviglio during drier phases (Bronze Age Warm Period, or BAWP; [ 70 ]) or sown at the end of the season in the same fields where other cereal crops had already been harvested Widespread arid conditions and poor nutrient supply are also indicated by the decrease in cyanobacteria colonies and riparian vegetation [ 37 , 69 ] and by the increase in Cichorieae and other Asteraceae [ 37 ], suggesting an expansion of dry meadows and pasturelands in the area. The biomolecular analyses also indicated the presence of triterpenyl acetates specific for Asteraceae- Picris hieracioides (such as bauerenyl and pichierenyl acetates) in the Bronze Age samples During the Roman Period, the territory was densely occupied, and the plain was organized into a rural system known as “centuriation” (from about the 2 nd century BC until the 5 th century AD; [ 71 ]). The high percentages of cereals (up to 6.4% at 99 cm in N-S 3) agree with this new human occupation. It is worth noting that a Roman villa rustica with its own agricultural productive area was located in the same territory occupied by the Bronze Age settlement [ 72 , 73 ]. In the Middle Ages (6 th–10 th centuries AD), a palustrine environment spread in the territory following the “Roman empire crisis” [ 40 , 73 ], with a large diffusion of hygrophilous herbs, including plants belonging to Paniceae [ 35 ]. During the Renaissance, a rapid loss of wet environments occurred after land reclamation (15 th–16 th centuries AD). The agrarian landscape began to characterize the Po Plain, dominated by cereal crops in the lowland. The recovery of Zea mays pollen in the upper part of the core (Figure 5 ) matches historical sources documenting the introduction of this cultivation in the Emilia-Romagna Region between the end of the 16 th century and the first half of the 17 th century [ 74 ]. 4.2. The Combined Evidence of Panicum Pollen and Sedimentary Miliacin during the Holocene Panicum pollen was identified in the three Poviglio cores, but only in 3 of the 31 samples analyzed for lipid biomarkers (N-S 3: 224 cm—pre-Bronze Age; C-S 1: 350 cm—pre-Bronze Age, 117 cm—Bronze Age; Figure 5 ). Panicum s.l. pollen grains were identified in the Middle–Late Neolithic samples (C-S 1 and F-S 2; Figure 6 ). Miliacin was detected in low abundance in at least three Neolithic samples in the N-S 3 core (Figure 6 ). In a recent synthesis of miliacin sedimentary records from six lakes in the European Western Alps by Jacob et al. [ 14 ], a low concentration of miliacin before the supposed introduction of millet in the Alps was considered an indication of the presence of another miliacin producer, such as wild Panicum sp. or Pennisetum sp. At that time, broomcorn millet was already a “Neolithic package of founder crops” in Europe (e.g., [ 75 – 78 ]). Based on Panicum macroremains and biomolecular markers, the challenging topic of reconstructing the timing and spread of millet cultivation in the study area requires further investigation.

[Find the meaning and references behind the names: Dal Corso, Quite, Fill, Iron, Dung, Ash, Minor, Base, Feed, Trend, Diet, Verona, Agro, Corso, Quality, Dal]

Sustainability 2023 , 15 , 10408 11 of 17 Miliacin concentrations increased slightly in the pre-Bronze Age samples, indicating the presence of miliacin producers. According to our hypothesis, the presence of broomcorn millet is suggested when the combined evidence of pollen and miliacin occurs (N-S 3: 224 cm; C-S 1: 350 cm; Figure 5 ). Panicum pollen grains are almost absent in the Bronze Age levels Panicum s.l. has been found only in the C-S 1 core (177 cm), whereas miliacin was documented in six samples of the three cores (Figure 5 ). As the occurrence of Panicum pollen grains was quite rare in the sediment samples (the plant is a low-pollen producer, as are the other cereal species), when it is recorded, it testifies to the occurrence of the plant. Archaeobotanical evidence reveals that millet was widespread in Europe in the 2 nd millennium BC, but there is uncertainty about the timing and diffusion path of its adoption by European societies [ 20 ]. Current research indicates that millet cultivation was spread for alimentary purposes by the Terramare civilization [ 20 , 79 , 80 ]). Archaeobotanical remains from Terramare sites (e.g., Montale [ 36 ]; Fondo Paviani [ 39 ]) indicate the diffusion of millets from the Early Bronze Age with an increasing trend of spread up to the Recent Bronze Age. During this last phase, a climatic crisis, including the depletion of water, affected the agro-sylvo-pastoral system of Terramare, mainly based on irrigated agriculture [ 37 ]. The increasing arid conditions toward the Recent Bronze Age are reflected by the diversification of cereal crops, shifting to drought-tolerant varieties as millets adopted by the Terramare people [ 37 , 81 ]. The combined evidence of Panicum pollen and miliacin in the Bronze Age level of C-S 1 (177 cm) suggests that broomcorn millet could have been cultivated at some distance (~300 m, i.e., C-S 1 position; Figure 1 ) from the settlement. Interestingly, in the same sample, the Avena/Triticum and Hordeum pollen groups decreased. Nevertheless, here as in other Terramare sites, intercropping practices are attested (e.g., [ 34 , 36 , 37 ]. In the N-S 3 core, the continuous presence of miliacin in the last phases of the Bronze Age can be considered an indicator of the local presence of Panicum sp Miliacin concentrations in the three Poviglio cores are always lower than those found in sediment samples from the Terramara of Fondo Paviani (Verona, N-E Italy; [ 39 ]). Many factors can control miliacin concentrations, such as dilution by the mineral matrix and transport of millet from cultivation or consumption sites [ 14 ]. Even at Fondo Paviani, the lowest miliacin yields came from a minerogenic unit constituting the base of the peripheral earthen rampart. The highest miliacin yields in the study by Dal Corso et al. [ 39 ] came from ash lenses of anthropogenic origin in the ditch fill. The presence of phytoliths derived from the chaff of pooid and panicoid cereals, grass leaves and stems, and panicoid inflorescences have suggested that the ashes were derived from the burning of dung from domestic animals [ 24 , 82 ]. As for the more recent chronological phases, the ubiquitous presence of miliacin seems to suggest a growth of Panicum in the area. However, the presence of this plant is not reflected in the pollen record. Broomcorn millet became a common culture during the Iron Age in northern and central Italy (ca. 3100-2400 cal yr BP; [ 83 ]), while in Roman times, the consumption of millet was reserved for animal feed and it was deemed suitable only for the poor as it was considered a low-quality cereal [ 84 ]. In general, archaeobotanical evidence suggests a minor presence of millet compared to other cereal crops at Roman sites [ 85 ]. A recent synthesis of Roman food from about one hundred Roman-period archaeological sites (from the 3 rd century BC to the 6 th century AD) located in northern Italy [ 86 ] highlights the presence of broomcorn millet grains in inhabited places and infrastructures and, in lesser amounts, in cemeteries and cult places. Carpological remains indicate the importance of “minor” cereals, especially broomcorn millet, during the Early Middle Ages in northern Italy, motivated by environmental and cultural factors that modified habits concerning food production and diet [ 84 , 87 ]. 4.3. Wild-Cereal Hints from Large-Size Cereal Pollen (Pre-Holocene Phases) Pollen spectra from the three cores, together with those from previous palynological research carried out in stratigraphic records from the Terramara S. Rosa di Poviglio, revealed

[Find the meaning and references behind the names: Less, Mercuri, Woodland, Macro, Might, Seem, Ability, Major, Meadow, Lack, Desert, Andrieu, Proto, Sahara, Role, See, Able]

Sustainability 2023 , 15 , 10408 12 of 17 the presence of wild cereal pollen grains even in phases prior to the settlement. This feature is particularly evident in the pre-Holocene levels (Figure 6 ), when common meadow grasses covered most of the lawns. Poaceae exploit a wide variety of environments as pioneer taxa; as such, Poaceae spread in woodland clearings and swamps in the alluvial plain The same environmental conditions allowed for the spread of cereals, mainly represented by the Hordeum group and Panicum s.l., produced by wild and not domesticated grass species. However, the morphology of wild cereal pollen cannot be distinguished from that of domesticated species [ 53 ]. In particular, the pollen size of wild Paniceae (e.g., pollen of Setaria and Echinochloa ) is mainly medium to large (e.g., [ 56 ]), partially overlapping with cereal pollen’s diameter measurements Wild cereals (e.g., most of the wild progenitors of cultivated cereals) are characterized by key adaptive features that allow them to thrive in disturbed soils and easily adapt to past environmental changes and conditions of environmental instability [ 56 ]. These plants are often described as pioneering, ruderal, weedy, and generalist and, thus, are able to colonize a wide range of habitats, from swamps to open grasslands to semi-desert habitats. The same behavior that allowed these plants to survive in a human-disturbed environment, described by Mercuri et al. [ 88 ], could also be suggested for Paniceae and other wild cereals in the central Po Plain during the pre-Holocene phases. At that time, unstable environments subject to alluvial events were present in the lowlands [ 35 ]. Pollen records from the plain identify a large open area dominated mostly by shrub-herbaceous communities related to alluvial processes that led to less-arid soil conditions and the temporary formation of wetlands and wet meadows (e.g., [ 68 ]). Arid grasslands and ephemeral pioneer meadows were probably present, defining the xerophilous herbaceous component. Open spaces and large meadows may have favored the frequentation of the area by grazing wild animals, as also suggested by the presence of coprophilous fungi in the palynological record [ 35 ]. However, the low quantities of coprophilous fungal spores [ 35 ] cannot support the hypothesis on the genetic modifications of the Poaceae taxa and the appearance of “proto-cereals” as a consequence of continuous overgrazing, as described in a different context by Andrieu-Ponel et al. [ 89 ]. In our record, the presence of large Poaceae pollen should probably be linked to climate change in the central Po Plain during the pre-Holocene [ 42 ], when the atmospheric drought may have favored the appearance of polyploidy in some species of Poaceae [ 89 – 92 ]. A recent study from the Early–Middle Holocene site of Takarkori in central Sahara [ 93 ] suggests that the occurrence of pollen anomalies, mainly in Paniceae, has derived from multiple genomes of some wild species as well as from plants’ ability to respond to environmental stresses. The Poviglio pollen record suggests that climate could have played a major role in the spread of Poaceae with large pollen in pre-domestication phases. The presence of such proto-cereals and evidence of polyploid grasses seem to indicate that a genetic modification of Poaceae might be related to adaptive responses to environmental changes that appeared long before the emergence of agriculture 5. Conclusions The results obtained from this research provide a new contribution to the reconstruction of the history of the broomcorn millet crop in northern Italy by combining the molecular tracer of Panicum miliaceum and its pollen evidence. Archaeocarpological records often lack grains of broomcorn millet due to preservation issues; hence, further approaches to investigate its introduction and adoption are necessary. Whereas on their own, Panicum pollen and miliacin as proxies are ambiguous, and the co-occurrence of this pollen type and miliacin in the same deposits can be taken as evidence of the presence and local cultivation of Panicum miliaceum . This is what we see in the sediment samples from the near-site cores (N-S 3, C-S 1, and F-S 2) of the Terramara S. Rosa di Poviglio. Our study demonstrates the potential of considering both types of materials preserved in the sediments and the importance of such an approach at sites where other (microor macro-) botanical evidence is poor or lacking altogether. In the cores of the Terramara S. Rosa di Poviglio, miliacin

[Find the meaning and references behind the names: Marta, Carry, Resources, Milan, Andrea, Board, Carl, Heron, Terra, Zanon, Marco, Desio, British, Prin, Net, Read, Dell, Kiel, Mic, Milano, Novel, Fondi, State, Della, Constant, Christian, Thank, Case, Author, Padua, Mauro, Principal]

Sustainability 2023 , 15 , 10408 13 of 17 shows the highest concentration in the samples attributed to the Bronze Age, confirming that the Terramare population cultivated millet fields. When recorded, the joint occurrence of Panicum pollen and sedimentary miliacin provides evidence of the local presence of the plant and suggests the adoption of millet crops, especially during the Recent Bronze Age, when the diversification of cereal crops shifting to drought-tolerant varieties has been documented [ 37 , 81 ]. In addition to the information obtained by the combination of pollen and biomolecular proxies, the detailed study of Poaceae curves provides new information on the presence of open areas related to different environmental conditions in the central Po Plain from a long-term perspective. Pollen spectra from the three cores indicate the presence of wild cereal pollen grains in the pre-Holocene levels. The adaptation of wild cereals to past environmental changes and environmental instability is attributed to their key adaptive features [ 88 ]. In this context, wild Poaceae with large pollen could result from adaptive responses (e.g., polyploidy) to climate changes, and this might be linked to phases of natural environmental changes occurring during the pre-Holocene phases in the central Po Plain The novel approach based on the combined evidence of Panicum pollen and miliacin indicated unequivocally the presence of broomcorn millet in sedimentary archives. In the current state of research, the precise chronology of the introduction and diffusion of millet crops in several large regions of Europe has been addressed [ 20 ], but some uncertainties remain about their local cultivation. Tracing millet through the two proxies could provide new, essential information on the timing of the earliest introduction of broomcorn millet in northern Italy Author Contributions: Conceptualization, A.F.; methodology, A.F., E.C. and J.J.; formal analysis, E.C and J.J.; investigation, A.F. and E.C.; resources, A.F.; data curation, A.F., E.C. and J.J.; writing—original draft preparation, A.F., E.C. and J.J.; writing—review and editing, A.F., E.C. and J.J.; visualization, A.F.; supervision, A.F. and J.J.; funding acquisition, A.F. All authors have read and agreed to the published version of the manuscript Funding: This research was funded by FAR 2020—Departmental Projects of the Dept. of Life Sciences, University of Modena and Reggio Emilia-UNIMORE (PI: A.F.). The cores were included in the SUCCESSO-TERRA Project (PRIN-20158 KBLNB—Human societies, climate/environment changes and resource exploitation/sustainability in the Po Plain in the mid-Holocene: the case study of the Terramare culture; Principal Investigators M. Cremaschi—A. Zerboni, Universit à degli Studi di Milano; https://www.successoterra.net (accessed on 1 June 2023) Institutional Review Board Statement: Not applicable Informed Consent Statement: Not applicable Data Availability Statement: The data presented in this study are available on request Acknowledgments: We thank Marta Dal Corso (University of Padua), Marco Zanon (Christian- Albrechts-Universität zu Kiel), and Carl Heron (British Museum) for their support in the early stages of the research. We are grateful to Mauro Cremaschi for the helpful comments on the stratigraphy and chronology of the cores. We are also indebted to Andrea Zerboni, who contributed to the interpretation of pollen data within the depositional context of the site and provided technical support for the sampling procedures. We thank the Municipality of Poviglio for its constant support. Funds for the excavations came from the University of Milan (Fondi Speciali di Ateneo per le Ricerche Archeologiche dell’Universit à degli Studi di Milano). All of the authorizations to carry out the excavation and sampling were released by the SABAP-BO (MiC) to the Dipartimento di Scienze della Terra “A. Desio” of the Universit à degli Studi di Milano. We thank the anonymous reviewers, whose constructive remarks helped improve the clarity of the manuscript Conflicts of Interest: The authors declare no conflict of interest.

[Find the meaning and references behind the names: Van Mourik, Van Loon, De Goede, Zhang, Liu, Dalla, Leonardi, Modern, Catena, Doyen, Wiley, Albuquerque, Brazil, Siddique, Press, Rolland, Heiss, Longa, Eds, Parisot, Arch, Veg, Billaud, Simonneau, Masuda, Hunt, Hatcher, Francis, Wang, Int, Sons, Vander, Sci, Revel, Schouten, Chapron, John, Chem, Desmet, Mes, Harding, Nicosia, Mani, Carbon, Staff, Evol, Arnaud, Dubois, Bosi, Loon, French, Front, London, Future, Jansen, Sifeddine, Bossard, Rajasekaran, Rep, Lallier, Histoire, Farooq, Cox, Ganesan, Track, Abe, Spengler, Smart, Miller, Messager, Fen, Linden, Oxford, Peng, Stika, Amsterdam, Jones, Iversen, Bourget, Schlesinger, Dans]

Sustainability 2023 , 15 , 10408 14 of 17 References 1 Castañeda, I.S.; Schouten, S. A review of molecular organic proxies for examining modern and ancient lacustrine environments Quat. Sci. Rev 2011 , 30 , 2581–2891. [ CrossRef ] 2 Dubois, N.; Jacob, J. Molecular Biomarkers of Anthropic Impacts in Natural Archives: A Review Front. Ecol. Evol 2016 , 4 , 92 [ CrossRef ] 3 Jansen, B.; Hooghiemstra, H.; de Goede, S.P.C.; van Mourik, J.M. Chapter 5—Biomarker analysis of soil archives. In Developments in Quaternary Sciences ; van Mourik, J.M., Van Der Meer, J.J.M., Eds.; Elsevier: Amsterdam, The Netherlands, 2019; Volume 18, pp 163–222 4 Zhang, Y.; Peng, T.; Yu, J.; Zhengwen, L. Natural and anthropogenically driven changes in the n-alkanols of lake sediments and implications for their use in paleoenvironmental studies of lakes Catena 2021 , 207 , 105591. [ CrossRef ] 5 Jansen, B.; van Loon, E.E.; Hooghiemstra, H.; Verstraten, J.B. Improved reconstruction of palaeo-environments through unravelling of preserved vegetation biomarker patterns Palaeogeogr. Palaeoclimatol. Palaeoecol 2010 , 285 , 119–130. [ CrossRef ] 6 Brocks, J.; Summons, R. Sedimentary hydrocarbons, biomarkers for early life. In Treatise on Geochemistry ; Schlesinger, W.H., Ed.; Elsevier: Oxford, UK, 2003; Volume 8, pp. 63–115 7 Simoneit, B.R.T.; Grimalt, J.O.; Wang, T.G.; Cox, R.E.; Hatcher, P.G.; Nissenbaum, A. Cyclic terpenoids of contemporary resinous plant detritus and of fossil woods, ambers and coals Org. Geochem 1985 , 10 , 877–889. [ CrossRef ] 8 Abe, I.; Rohmer, M.; Prestwich, G.D. Enzymatic cyclization of squalene and oxidosqualene to sterols and triterpens Chem. Rev 1993 , 93 , 2189–2208. [ CrossRef ] 9 Jacob, J.; Disnar, J.-R.; Boussafir, M.; Spadano Albuquerque, A.L.; Sifeddine, A.; Turcq, B. Pentacyclic triterpene methyl ethers in recent lacustrine sediments (Lagoado Caç ó , Brazil) Org. Geochem 2005 , 36 , 449–461. [ CrossRef ] 10 Lavrieux, M.; Jacob, J.; Le Milbeau, C.; Zocatelli, R.; Masuda, K.; Br é heret, J.G.; Disnar, J.R. Occurrence of triterpenyl acetates in soil and their potential as chemotaxonomical markers of Asteraceae Org. Geochem 2011 , 42 , 1315–1323. [ CrossRef ] 11 Bossard, N.; Jacob, J.; Le Milbeau, C.; Sauze, J.; Terwilliger, V.; Poissonnier, B.; Verg è s, E. Distribution of miliacin (olean-18-en-3 β -ol methyl ether) and related compounds in broomcorn millet ( Panicum miliaceum ) and other reputed sources: Implications for the use of sedimentary miliacin as a tracer of millet Org. Geochem 2013 , 63 , 48–55. [ CrossRef ] 12 Jacob, J.; Disnar, J.-R.; Arnaud, F.; Chapron, E.; Debret, M.; Lallier-Verg è s, E.; Desmet, M.; Revel-Rolland, M. Millet cultivation history in the French Alps as evidenced by a sedimentary molecule J. Archaeol. Sci 2008 , 35 , 814–820. [ CrossRef ] 13 Jacob, J.; Disnar, J.-R.; Bardoux, G. Carbon isotope evidence for sedimentary miliacin as a tracer of Panicum miliaceum (broomcorn millet) in the sediments of Lake le Bourget (French Alps) Org. Geochem 2008 , 39 , 1077–1080. [ CrossRef ] 14 Jacob, J.; Simonneau, A.; Bajard, M.; Messager, E.; Billaud, Y.; Le Milbeau, C.; Poulenard, J.; Giguet-Covex, C.; Rapuc, W.; Doyen, E.; et al. Diffusion et dynamique de la culture du millet dans les Alpes ouest-europ é ennes sur la base des concentrations en miliacine dans les s é diments lacustres. In Actes des 41 è mes Rencontres Internationales d’Arch é ologie et d’Histoire—Biodiversit és , Environnements et Soci é t és Depuis la Pr é histoire: Nouveaux Marqueurs et Approches Int é gr é es ; Nicoud, E., Balasse, M., Desclaux, E., Th é ry-Parisot, I., Eds.; É ditions APDCA: Antibes, France, 2021; pp. 123–134 15 Rajasekaran, R.; Francis, N.; Mani, V.; Ganesan, J. Proso millet ( Panicum miliaceum L.). In Neglected and Underutilized Crops: Future Smart Food ; Farooq, M., Siddique, K.H.M., Eds.; Academic Press: London, UK, 2023; pp. 247–278. [ CrossRef ] 16 Miller, N.F.; Spengler, R.N.; Frachetti, M. Millet cultivation across Eurasia: Origins, spread, and the influence of seasonal climate Holocene 2016 , 26 , 1566–1575. [ CrossRef ] 17 Hunt, H.V.; Vander Linden, M.; Liu, X.; Motuzaite-Matuzeviciute, G.; College, S.; Jones, M.K. Millet across Eurasia: Chronology and context of early records of the genera Panicum and Setaria from archaeological sites in the Old World Veg. Hist. Archaeobotany 2008 , 17 , 5–18. [ CrossRef ] [ PubMed ] 18 Motuzaite-Matuzeviciute, G.; Staff, R.A.; Hunt, H.V.; Liu, X.; Jones, M.K. The early chronology of broomcorn millet ( Panicum miliaceum ) in Europe Antiquity 2013 , 338 , 1073–1085. [ CrossRef ] 19 Stika, H.-P.; Heiss, A.G. Plant cultivation in the bronze age. In The Oxford Handbook of the European Bronze Age ; Fokkens, H., Harding, A., Eds.; Oxford University Press: Oxford, UK, 2013; pp. 369–384 20 Filipovi´c, D.; Meadows, J.; Dal Corso, M.; Kirleis, W.; Alsleben, A.; Akeret, Ö.; Bittmann, F.; Bosi, G.; Ciut ă , B.; Dreslerov á , D.; et al New AMS 14 C dates track the arrival and spread of broomcorn millet cultivation and agricultural change in prehistoric Europe Sci. Rep 2020 , 10 , 13698. [ CrossRef ] 21 Ohmoto, T.; Ikuse, M.; Natori, S. Triterpenoids of the Gramineae Phytochemistry 1970 , 6 , 2137–2148. [ CrossRef ] 22 Lu, H.; Zhang, J.; Wu, N.; Liu, K.-B.; Xu, D.; Li, Q. Phytoliths Analysis for the Discrimination of Foxtail Millet ( Setaria italica ) and Common Millet ( Panicum miliaceum ) PLoS ONE 2009 , 4 , e 4448. [ CrossRef ] [ PubMed ] 23 Out, W.A.; Madella, M. Morphometric distinction between bilobate phytoliths from Panicum miliaceum and Setaria italica leaves Archaeol. Anthropol. Sci 2016 , 8 , 505–521. [ CrossRef ] 24 Dal Corso, M.; Nicosia, C.; Balista, C.; Cupit ò , M.; Dalla Longa, E.; Leonardi, G.; Kirleis, W. Bronze Age crop processing evidence in the phytolith assemblages from the ditch and fen around Fondo Paviani, northern Italy Veg. Hist. Archaeobotany 2017 , 26 , 5–24 [ CrossRef ] 25 Faegri, K.; Iversen, J Textbook of Pollen Analysis ; John Wiley & Sons: Chichester, UK, 1989.

[Find the meaning and references behind the names: Fruits, Creation, Gran, Richard, Torri, Brea, Alexandria, Accorsi, Valley, Ammann, Nel, Rinaldi, Mazzanti, Marine, Bandini, Giardini, Masi, Urban, Finn, Marchesini, Janick, Mele, Brandolini, Development, Italiano, Lago, Lippi, Germany, Vescovi, Giudici, Morel, Limousin, Siala, Vives, Bassi, Montecchi, Vignola, Flood, Gardens, Vigliotti, Olmi, Pizzi, Bonnier, Valsecchi, Trevisan, Genesis, Else, Progress, Taylor, Mariotti, Plateau, Miras, Firenze, Bari, Cardarelli, Gelabert, Labate, Socio, March, Massamba, Grandi]

Sustainability 2023 , 15 , 10408 15 of 17 26 Miras, Y.; Guenet, P.; Richard, H. Holocene vegetation, landscape, and reconstruction of human activity from prehistory to the roman period based on new pollen data performed in “the plateau de Millevaches” (Limousin, Massif Central, France) Quaternaire 2011 , 22 , 147–164. [ CrossRef ] 27 Mercuri, A.M. Genesis and evolution of the cultural landscape in central Mediterranean: The ‘where, when and how’ through the palynological approach Landsc. Ecol 2014 , 29 , 1799–1810. [ CrossRef ] 28 Hjelle, K.; Overland, A.; Gran, M.M.; Romundset, A.; Ystgaard, I. Two thousand years of Landscape—Human interactions at a coastal peninsula in Norway revealed through pollen analysis, shoreline reconstruction, and radiocarbon dates from archaeological sites Front. Ecol. Evol 2022 , 10 , 911780. [ CrossRef ] 29 Vignola, C.; Hättestrand, M.; Bonnier, A.; Finn é , M.; Izdebski, A.; Katrantsiotis, C.; Kouli, K.; Liakopoulos, G.C.; Norström, E.; Papadaki, M.; et al. Mid-late Holocene vegetation history of the Argive Plain (Peloponnese, Greece) as inferred from a pollen record from ancient Lake Lerna PLoS ONE 2022 , 17 , e 0271548. [ CrossRef ] 30 Sadori, L.; Mercuri, A.M.; Mariotti Lippi, M. Reconstructing past cultural landscape and human impact using pollen and plant macroremains Plant Biosyst 2010 , 144 , 940–951. [ CrossRef ] 31 Mercuri, A.M.; Florenzano, A.; Burjachs, F.; Giardini, M.; Kouli, K.; Masi, A.; Picornell-Gelabert, L.; Revelles, J.; Servera-Vives, G.; Sadori, L.; et al. From influence to impact: The multifunctional land-use in Mediterranean prehistory emerging from palynology of Archaeological sites (8.2–2.8 ka BP) Holocene 2019 , 29 , 830–846. [ CrossRef ] 32 Revelles, J. The Role of Palynology in Archaeoecological Research: Reconstructing Human-Environment Interactions during Neolithic in the Western Mediterranean Appl. Sci 2021 , 11 , 4073. [ CrossRef ] 33 Baltensperger, D.D. Foxtail and proso millet. In Progress in New Crops ; Janick, J., Ed.; ASHS Press: Alexandria, VA, USA, 1996; pp. 182–190 34 Mercuri, A.M.; Accorsi, C.A.; Bandini Mazzanti, M.; Bosi, G.; Trevisan Grandi, G.; Cardarelli, A.; Labate, D.; Marchesini, M.; Olmi, L.; Torri, P. Cereal fields from the middle-recent Bronze Age, as found in the Terramara di Montale, in the Po Plain (Emilia-Romagna, Northern Italy), based on pollen, seeds/fruits and microcharchoals. In The Archaeology of Crop Fields and Gardens ; Morel, J.-P., Tresserras, J., Matalama, J.C., Eds.; Edipuglia: Bari, Italy, 2006; pp. 251–270 35 Cl ò , E. Climate Change and Human Impact in a Long-Term Perspective: Palynology of Central Po Plain Deposits Since the Lateglacial. Ph.D. Thesis, Universit à degli Studi di Modena e Reggio Emilia, Modena, Italy, 2022 36 Mercuri, A.M.; Accorsi, C.A.; Bandini Mazzanti, M.; Bosi, G.; Cardarelli, A.; Labate, D.; Marchesini, M.; Trevisan Grandi, G. Economy and environment of Bronze Age settlements e Terramaras e in the Po Plain (Northern Italy): First results of the archaeobotanical research at the Terramara di Montale Veg. Hist. Archaeobotany 2006 , 16 , 43–60. [ CrossRef ] 37 Cremaschi, M.; Mercuri, A.M.; Torri, P.; Florenzano, A.; Pizzi, C.; Marchesini, M.; Zerboni, A. Climate change versus land management in the Po Plain (Northern Italy) during the Bronze Age: New insights from the VP/VG sequence of the Terramara Santa Rosa di Poviglio Quat. Sci. Rev 2016 , 136 , 153–172. [ CrossRef ] 38 Dal Corso, M.; Marchesini, M.; Leonardi, G.; Kirleis, W. Environmental changes and human impact during Bronze Age in northern Italy: On-site palynological investigation at Fondo Paviani (Legnago, Verona). In Collapse or Continuity? Environment and Development of Bronze Age Human Landscapes. Proceedings of the International Workshop “Socio-Environmental Dynamics over the Last 12,000 Years: The Creation of Landscapes II (14–18 March 2011)” in Kiel ; Kneisel, J., Kirleis, W., Dal Corso, M., Taylor, N., Tiedtke, V., Eds.; Kiel Archaeology: Kiel, Germany, 2012; Volume 1, pp. 71–84 39 Dal Corso, M.; Zanon, M.; Heron, C.; Rottoli, M.; Cupit ò , M.; Dalla Longa, E.; Kirleis, W. Tracing millet through biomarker analysis in archaeological sites in alluvial plains: The first miliacin data from the northern Italian Bronze Age. In Millet and What Else? The Wider Context of the Adoption of Millet Cultivation in Europe ; Kirleis, W., Dal Corso, M., Filipovi´c, D., Eds.; Sidestone Press: Leiden, The Netherlands, 2022; pp. 231–252 40 Brandolini, F.; Cremaschi, M. The Impact of Late Holocene Flood Management on the Central Po Plain (Northern Italy) Sustainability 2018 , 10 , 3968. [ CrossRef ] 41 Valsecchi, V.; Tinner, W.; Finsinger, W.; Ammann, B. Human impact during the Bronze Age in the vegetation at Lago Lucone (northern Italy) Veg. Hist. Archaeobotany 2006 , 15 , 99–113. [ CrossRef ] 42 Vescovi, E.; Kaltenrieder, P.; Tinner, W. Late-glacial and Holocene vegetation history of Pavullo nel Frignano (Northern Apennines, Italy) Rev. Palaeobot. Palynol 2010 , 160 , 32–45. [ CrossRef ] 43 Bosi, G.; Mercuri, A.M.; Bandini Mazzanti, M.; Florenzano, A.; Montecchi, M.C.; Torri, P.; Labate, D.; Rinaldi, R. The evolution of Roman urban environments through the archaeobotanical remains in Modena e Northern Italy J. Archaeol. Sci 2015 , 53 , 19–31 [ CrossRef ] 44 Mercuri, A.M.; Bandini Mazzanti, M.; Torri, P.; Vigliotti, L.; Bosi, G.; Florenzano, A.; Olmi, L.; Massamba N’siala, I. A marine/terrestrial integration for mid-late Holocene vegetation history and the development of the cultural landscape in the Po Valley as a result of human impact and climate change Veg. Hist. Archaeobotany 2012 , 21 , 353–372. [ CrossRef ] 45 Mele, M.; Cremaschi, M.; Giudici, M.; Lozej, A.; Pizzi, C.; Bassi, A. The terramare and the surrounding hydraulic structures: A geophysical survey of the Santa Rosa site at Poviglio (Bronze Age; Northern Italy) J. Archaeol. Sci 2013 , 40 , 4648–4662. [ CrossRef ] 46 Bernab ò Brea, M.; Cremaschi, M Il Villaggio Piccolo Della Terramara di S. Rosa a Poviglio. Scavi 1987/1992 ; Istituto Italiano Preistoria e Protostoria: Firenze, Italy, 2004.

[Find the meaning and references behind the names: Du Nord, Di Giuseppe, Di Lernia, De Amicis, Mutti, Files, Trees, Mol, Bologna, Pharm, Andersen, Joly, Inoue, Italia, Bayliss, Carta, Fournier, Lipp, Zone, Ales, Comune, Penny, Nord, Grain, Pini, Bard, Weiss, Ann, Denmark, Jakupovic, West, Hunter, Bianchini, Gli, Adda, Centro, Moore, Marchetti, December, Piazza, Cave, Giuseppe, Bottazzi, Friedrich, Asia, Nelle, Lernia, Afrique, Grimm, Motte, Alp, Garibaldi, Hahn, Pfeil, Falsone, Flora, Cazzola, Marseille, Rawa, Alpi, Java, Collins, Dam, Tours, Blackwell, Sibilia, Vittori, Edwards, Reimer, Sec, Sua, Hopf, Pederzoli, Routledge, Webb, Munich, Lett, Suzuki, Fischer, Trombino, Lac, Puy, Pignatti, Barreau, Ishida, Beni, Serio, Pessina, Mais, Magny, Sol, Vezzoli, Gustav, Guida, Cremonini, Rome, Nella, Iris, Giglio, Ramsey, Lin, Mariani, Austin, Zanchi, Vianello, Cheng, Parma, Ravello, Florence]

Sustainability 2023 , 15 , 10408 16 of 17 47 Reimer, P.J.; Austin, W.E.N.; Bard, E.; Bayliss, A.; Blackwell, P.G.; Bronk Ramsey, C.; Butzin, M.; Cheng, H.; Edwards, R.L.; Friedrich, M.; et al. The IntCal 20 Northern Hemisphere radiocarbon age calibration curve (0–55 kcal BP) Radiocarbon 2020 , 62 , 725–757. [ CrossRef ] 48 van der Kaars, S.; Penny, D.; Tibby, J.; Fluin, J.; Dam, R.; Suparan, P. Late Quaternary palaeoecology, palynology and palaeolimnology of a tropical lowland swamp: Rawa Danau, West Java, Indonesia Palaeogeogr. Palaeoclimatol. Palaeoecol 2001 , 171 , 185–212 [ CrossRef ] 49 Florenzano, A.; Mercuri, A.M.; Pederzoli, A.; Torri, P.; Bosi, G.; Olmi, L.; Rinaldi, R.; Bandini Mazzanti, M. The significance of intestinal parasite remains in pollen samples from Medieval pits in the Piazza Garibaldi of Parma, Emilia Romagna, Northern Italy Geoarchaeology 2012 , 27 , 34–47. [ CrossRef ] 50 Moore, P.D.; Webb, J.A.; Collins, M.E Pollen Analysis , 2 nd ed.; Blackwell Scientific Publications: Oxford, UK, 1991 51 Reille, M Pollen et Spores d’Europe et d’Afrique du Nord ; Laboratoire de Botanique Historique et Palinologie: Marseille, France, 1992 52 Beug, H.J Leitfaden der Pollenbestimmung für Mitteleuropa und Angrezende Gebiete ; Gustav Fischer: Stuttgart, Germany, 1961 53 Andersen, S.T Identification of Wild Grasses and Cereal Pollen ; Danmarks Geologiske Undersøgelse: Copenhagen, Denmark, 1979 54 Joly, C.; Barill é , L.; Barreau, M.; Mancheron, A.; Visset, L. Grain and annulus diameter as criteria for distinguishing pollen grains of cereals from wild grasses Rev. Palaeobot. Palynol 2007 , 146 , 221–233. [ CrossRef ] 55 Beug, H.J Leitfaden der Pollenbestimmung für Mitteleuropa und Angrenzende Gebiete ; Dr. Friedrich Pfeil: Munich, Germany, 2015 56 Mercuri, A.M. Palynological analysis of the Early Holocene sequence. In The Uan Afuda Cave Hunter-gatherer Societies of Central Sahara, Arid Zone Archaeology, Monographs 1 ; Di Lernia, S., Ed.; All’Insegna del Giglio: Florence, Italy, 1999; pp. 149–181 57 Grimm, E.C TILIA and TGView ; Illinois State Museum: Springfield, IL, USA, 2004 58 Mercuri, A.M.; Bandini Mazzanti, M.; Florenzano, A.; Montecchi, M.C.; Rattighieri, E.; Torri, P. Anthropogenic Pollen Indicators (API) from archaeological sites as local evidence of human-induced environments in the Italian peninsula Ann. Bot. (Rome) 2013 , 3 , 143–153. [ CrossRef ] 59 Shiojima, K.; Masuda, K.; Lin, T.; Suzuki, H.; Ageta, H.; Inoue, M.; Ishida, T. Composite constituents: Three gammacer-16-ene derivatives, novel triterpenoids isolated from roots of Picris hieracioides L. subsp japonica. Tetrahedron Lett 1989 , 30 , 4977–4980 [ CrossRef ] 60 Shiojima, K.; Masuda, K.; Ooishi, Y.; Suzuki, H.; Ageta, H. Composite constituents: New migrated gammacerane triterpenoids from roots of Picris hieracioides L. subsp japonica. Tetrahedron Lett 1989 , 30 , 6873–6874. [ CrossRef ] 61 Shiojima, K.; Masuda, K.; Suzuki, H.; Lin, T.; Ooishi, Y.; Ageta, H. Composite constituents: Forty-two triterpenoids including eight novel compounds isolated from Picris hieracioides subsp japonica. Chem. Pharm. Bull 1995 , 43 , 1634–1639. [ CrossRef ] 62 Jakupovic, J.; Tsichrintzis, F.; Zdero, C.; Hahn, A.; Lipp, W.; Smalla, B.; Berendsohn, W.G. Natural Substances in the Compositae: The Bohlmann Files. Available online: http://www.bgbm.org/BioDivInf/projects/bohlmannfiles/ (accessed on 4 December 2011) 63 Lavrieux, M. Biomarqueurs Mol é culaires d’Occupation des Sols, du Sol au S é diment: Exemple du Bassin-Versant et du Lac d’Aydat (Puy-de-D ô me). Ph.D. Thesis, Universit é de Tours, Tours, France, 2011 64 Ravazzi, C.; Peresani, M.; Pini, R.; Vescovi, E. Il Tardoglaciale nelle alpi italiane e in pianura padana. Evoluzione stratigrafica, storia della vegetazione e del popolamento antropico Alp. Mediterr. Quat 2007 , 20 , 163–184 65 Ravazzi, C.; Deaddis, M.; De Amicis, M.; Marchetti, M.; Vezzoli, G.; Zanchi, A. The last 40 ka evolution of the Central Po Plain between the Adda and Serio rivers G é omorphologie 2012 , 18 , 131–154. [ CrossRef ] 66 Vittori Antisari, L.; Bianchini, G.; Cremonini, S.; Di Giuseppe, D.; Falsone, G.; Marchesini, M.; Marvelli, S.; Vianello, G. Multidisciplinary study of a Lateglacial-Holocene sedimentary sequence near Bologna (Italy): Insights on natural and anthropogenic impacts on the landscape dynamics J. Soils Sediments 2016 , 16 , 645–662. [ CrossRef ] 67 Pignatti, S Flora d’Italia ; Edagricole: Bologna, Italy, 1982 68 Ravazzi, C.; Donegana, M.; Vescovi, E.; Arpenti, E.; Caccianiga, M.; Kaltenrieder, P.; Londeix, L.; Marabini, S.; Mariani, S.; Pini, R.; et al. A new Lateglacial site with Picea abies in the Northern Apennine foothills: A population failing the model of glacial refugia trees Veg. Hist. Archaeobotany 2006 , 15 , 357–371. [ CrossRef ] 69 Cl ò , E.; Florenzano, A. Heterocysts of Rivularia Type for Interpreting a Palaeoenvironmental Context of the Late Quaternary in Northern Italy Sustainability 2022 , 14 , 15332. [ CrossRef ] 70 Magny, M. Holocene climatic variability as reflected by mid-European lake-level fluctuations, and its probable impact on prehistoric human settlements Quat. Int 2004 , 113 , 65–79. [ CrossRef ] 71 Bottazzi, G.; Bronzoni, L.; Mutti, A Carta Archeologica del Comune di Poviglio 1986–1989 ; IRIS: Poviglio, Italy, 1995 72 Mutti, A. Gli scavi nella terramara Santa Rosa a Fodico di Poviglio. In Guida All’esposizione ; Fondazione P. Manodori: Reggio Emilia, Italy, 2007 73 Brandolini, F.; Trombino, L.; Sibilia, E.; Cremaschi, M. Micromorphology and site formation processes in the Castrum Popilii Medieval Motte (N Italy) J. Archaeol. Sci. Rep 2018 , 20 , 18–32. [ CrossRef ] 74 Cazzola, F. L’introduzione del mais in Italia e la sua utilizzazione alimentare (sec. XVI–XVIII). In La Pr é paration Alimentaire des C é r é ales ; Fournier, S., Sigaut, F., Eds.; Centro Universitario Europeo per i Beni Culturali: Ravello, Italy, 1992; pp. 109–127 75 Zohary, D.; Hopf, M.; Weiss, E Domestication of Plants in the Old World , 4 th ed.; Clarendon Press: Oxford, UK, 2012 76 Rottoli, M.; Pessina, A. Neolithic agriculture in Italy: An update of archaeobotanical data with particular emphasis on northern settlements. In The Origins and Spread of Domestic Plants in Southwest Asia and Europe ; Colledge, S., Conolly, J., Eds.; Routledge: London, UK, 2016; pp. 141–154.

[Find the meaning and references behind the names: Ashman, Mcdonald, Canci, Oliver, Fornaciari, Fuller, Cambridge, Force, Magri, Rudolf, Dei, North, Chavarria, Fox, Pozzi, Band, Delle, Bonn, Kariyawasam, Martin, Michaud, Valenti, Tafuri, Lightfoot, Zappa, Idea, Kluska, Mitchell, Rey, Alexander, Costa, Proserpio, Severi, Dias, Sulla, Barbaro, Craig, Ideas, Arc, Cell, Gallinaro, Honour, Rochette, Bastida, Olds, Nat, Lehman, Carro, Murphy, Vanin, Castiglioni, Furia, Atlas, Manzaneda, Property, Mayda, Raskin]

Sustainability 2023 , 15 , 10408 17 of 17 77 Kirleis, W Atlas of Neolithic Plant Remains from Northern Central Europe ; Barkhuis Publishing: Groningen, The Netherlands, 2019 78 Liu, X.Y.; Motuzaite-Matuzeviciute, G.; Hunt, H.V. From a fertile idea to a fertile arc: The origins of broomcorn millet 15 years on. In Far from the Hearth: Essays in Honour of Martin K. Jones ; Lightfoot, E., Liu, X.Y., Fuller, D.Q., Eds.; McDonald Institute of Archaeology: Cambridge, UK, 2019; pp. 155–164 79 Tafuri, M.A.; Craig, O.E.; Canci, A. Stable isotope evidence for the consumption of millet and other plants in Bronze Age Italy Assoc. Phys. Anthropol 2009 , 139 , 146–153. [ CrossRef ] 80 Chiaffarelli, G.; Indipolage Kariyawasam, T.; Proserpio, B.; Rottoli, M. Indipolage Kariyawasam, T.; Proserpio, B.; Rottoli, M. I macroresti vegetali nelle strutture al margine del Villaggio Grande della Terramara Santa Rosa di poviglio. In Terramara Santa Rosa di Poviglio—Le Strutture Idrauliche al Margine del Villaggio Grande. Scavi 1998–2011 ; Cremaschi, M., Pizzi, C., Eds.; Istituto Italiano di Preistoria e Protostoria: Florence, Italy, 2021; Origines, 35; pp. 389–401 81 Florenzano, A.; Cl ò , E.; Zappa, J.; Montecchi, M.C.; Furia, E.; Torri, P.; Mercuri, A.M. Paesaggio vegetale sulla base delle analisi del riempimento dei pozzi al margine del Villaggio Grande della Terramara Santa Rosa di Poviglio. In Terramara Santa Rosa di Poviglio—Le Strutture Idrauliche al Margine del Villaggio Grande. Scavi 1998–2011 ; Cremaschi, M., Pizzi, C., Eds.; Istituto Italiano di Preistoria e Protostoria: Florence, Italy, 2021; Origines, 35; pp. 373–387 82 Dal Corso, M. Environmental History and Development of the Human Landscape in a Northeastern Italian Lowland during the Bronze Age: A Multidisciplinary Case-Study. In Universitätsforschungen zur Prehistorischen Archäologie—Band 312, HDL 12 ; Dr Rudolf Habelt GmbH: Bonn, Germany, 2018 83 Fundurulic, A.; Valenti, I.; Celant, A.; Barbaro, B.; Costa, M.; Manhita, A.; Severi, E.; Barrocas Dias, C.; Magri, D. Millets and Cereal Meals from the Early Iron Age Underwater Settlement of “Gran Carro” (Bolsena Lake, Central Italy) Sustainability 2022 , 14 , 3941. [ CrossRef ] 84 Castiglioni, E.; Rottoli, M. Broomcorn millet, foxtail millet and sorghum in north Italian Early Medieval sites PCA 2013 , 3 , 131–144 85 Murphy, C. Finding millet in the Roman world Archaeol. Anthropol. Sci 2016 , 8 , 65–78. [ CrossRef ] 86 Bosi, G.; Castiglioni, E.; Rinaldi, R.; Mazzanti, M.; Marchesini, M.; Rottoli, M. Archaeobotanical evidence of food plants in Northern Italy during the Roman period Veg. Hist. Archaeobotany 2020 , 29 , 681–697. [ CrossRef ] 87 Ganzarolli, G.; Alexander, M.; Chavarria, A.; Oliver, A.; Craig, E. Direct evidence from lipid residue analysis for the routine consumption of millet in Early Medieval Italy J. Archaeol. Sci 2018 , 96 , 124–130. [ CrossRef ] 88 Mercuri, A.M.; Fornaciari, R.; Gallinaro, M.; Vanin, S.; di Lernia, S. Plant behaviour from human imprints and the cultivation of wild cereals in Holocene Sahara Nat. Plants 2018 , 4 , 71–81. [ CrossRef ] 89 Andrieu-Ponel, V.; Rochette, P.; Demory, F.; Alçiçek, H.; Boulbes, N.; Bourl è s, D.; Helvacı, C.; Lebatard, A.-E.; Mayda, S.; Michaud, H.; et al. Continuous presence of proto-cereals in Anatolia since 2.3 Ma, and their possible co-evolution with large herbivores and hominins Sci. Rep 2021 , 11 , 8914. [ CrossRef ] 90 Ejsmond, M.J.; Wro ´nska-Pilarek, D.; Ejsmond, A.; Dragosz-Kluska, D.; Karpi ´nska-Kołaczek, M.; Kołaczek, P.; Kozłowski, J. Does climate affect pollen morphology? Optimal size and shape of pollen grains under various desiccation intensity Ecosphere 2011 , 2 , 117. [ CrossRef ] 91 Manzaneda, A.J.; Rey, P.J.; Bastida, J.M.; Weiss-Lehman, C.; Raskin, E.; Mitchell-Olds, T. Environmental aridity is associated with cytotype segregation and polyploidy occurrence in Brachypodium distachyon (Poaceae) New Phytol 2012 , 193 , 797–805. [ CrossRef ] 92 Fox, D.T.; Soltis, D.E.; Soltis, P.S.; Ashman, T.-L.; Van de Peer, Y. Polyploidy: A biological force from cells to ecosystems Trends Cell Biol 2020 , 30 , 688–694. [ CrossRef ] 93 Mercuri, A.M.; Cl ò , E.; Florenzano, A. Multiporate Pollen of Poaceae as Bioindicator of Environmental Stress: First Archaeobotanical Evidence from the Early–Middle Holocene Site of Takarkori in the Central Sahara Quaternary 2022 , 5 , 41. [ CrossRef ] Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.