Multiporate Pollen of Poaceae as Bioindicator of Environmental Stress: First Archaeobotanical Evidence from the Early–Middle Holocene Site of Takarkori in the Central Sahara

: This paper reports on the most ancient unusual morphological trait of the apertures of Poaceae pollen found in archaeological layers. In Poaceae, high levels of hybridization, polyploidy, apomixis, and multiporate pollen are often related. Multiple genomes in polyploids are critical for the adaptation of plant species to stresses and could be revealed by anomalies in pollen development. Therefore, the paleoenvironmental research can gain great beneﬁts from identifying polyploids in past contexts by observing anomalous pollen morphology during pollen counts. The occurrence of multiporate pollen in Poaceae has also been related to special features of the ecology of the species showing this anomaly, as well as to climatic and environmental stresses experienced by Poaceae living in a given region. Multiporate and bi- or tri-porate instead of monoporate pollen grains have been observed in samples taken from Takarkori rockshelter, an archaeological site in southwestern Libya (central Sahara) that has been occupied between ~10,200 and ~4650 cal BP. Multiporate pollen was found in organic sands and coprolites of ovicaprines. On the basis of archaeobotanical research, this work aims to investigate whether the presence of supernumerary pores in Poaceae pollen may be an effect of both climatic/hydrological changes and continued anthropogenic pressure on the wild grasses living in the region. The presence of multiporate pollen reveals that Poaceae that lived in central Sahara tackled several kinds of stress during the early and middle Holocene. The Takarkori pollen record suggests that climate change could have played a major role in the early Holocene, while human pressure became stronger during the middle Holocene. The change in environmental conditions determined adaptive responses of polyploid grasses even in the form of multiporate pollen.


Introduction
Plant reproductive biology and archaeology have found unprecedented interconnection through the discovery of anomalous pollen in archaeological contexts. Multiporate grass pollen was recorded with unusual frequency at the Holocene archaeological site of Takarkori, in central Sahara. The fact that this pollen belongs to the cereal family-Poaceae-suggests the availability of resources in an area long frequented by people who depended on it. This area has simultaneously undergone extensive and sometimes rapid environmental transformations over the last few thousand years.

Developmental Conditions of Multiporate Pollen
The Poaceae family is classically endowed with ana-monoporate pollen, i.e., a single "pore" aperture develops distally in each monad during microsporogenesis. Although this is the most common and widespread condition of the plant cycle, multiple pores have been reported in the pollen of several grass species. Poaceae include wild, cultivated, what are cause and effect, the presence of multiporate pollen can be used for preliminary identification of apomixis [3,9].
To sum up, Poaceae multiporate pollen is an effect of reproductive cycle anomalies. It is often related to high levels of hybridization, polyploidy and amphiploidy, and apomixis [12]. The occurrence of this pollen anomaly in Poaceae has been connected to special properties of the plant species, as well as to the ability to reply to environmental stresses experienced by individuals living in a given region.

Aim of the Paper and the Archaeobotanical Evidence
This paper reports on multiporate pollen observed in an archaeological context with the aim of testing whether pollen with this anomaly could be a new type of environmental marker.
The Takarkori rockshelter, located in the Tadrart Acacus mountains in southwestern Libya (central Sahara), has been the subject of interdisciplinary research that has shed new light on the complexity of the relationships between humans and plants. Interdisciplinary studies have yielded a wealth of information about climate, vegetation, and people that lived in the region (e.g., [19,20]). Different cultures, with different degrees of social complexity, inhabited this desert area while challenging environmental transformations for millennia. The preservation of organic matter in this sheltered site is exceptional and allowed the preservation of molecular residues, lipids, and ancient DNA. They revealed, together with archaeology and archaeobotany, hitherto unknown details of plants consumption and uses [21,22]. The cp-DNA of some wild cereals preserved by dehydration was extracted, amplified, and sequenced, confirming the morphological identification of the African minor cereals Echinochloa colona (L.) Link., Panicum laetum Kunth., and Sorghum bicolor (L.) Moench [23]. Further evidence of the extraordinary state of preservation of the plant remains was in the findings of weavings, likely used for the construction of baskets to collect or store grains [24] and of wood with the colour still preserved from the pigments used for rock art [25]. Well preserved pollen grains were extracted by sediments accumulated into the site [20]. Among them, the multiporate pollen of Poaceae is an interesting anomaly that has never been reported in the Holocene Sahara. Our hypothesis is that this particular pollen is linked to a combination of environmental/anthropogenic stresses that occurred at the site and in the region at the time of occupation.

Materials and Methods
The site of Takarkori (N 24.53714, E 10.51969; 1100 m asl) has been excavated and studied in the framework of the Archaeological Mission in the Sahara (former Joint Italian-Libyan Archaeological Mission in Acacus and Messak) of the Department of Antiquities, Tripoli, and Sapienza University of Rome. In this area, the archaeological research identified a first exploitation by hunter-gatherer-fishers, which was followed by a long occupation of pastoral communities. These were herders of cattle, and of sheep and goats.
Archaeological chronology is supported by ceramics, lithics and 43 radiocarbon dates on heaths, floors, organic sand and burials (see Table 3 in [20]), besides 30 radiocarbon dates on seed accumulations (see Table 1 in [23]). The chronological phases of occupation of the rockshelter, based on radiocarbon dates, archaeological materials, geoarchaeology and stratigraphy, can be resumed in the four main phases reported below, from the oldest to the most recent (e.g., [26]): − Late Acacus (LA) = hunter-gatherer-fishers,~10,170-~8180 cal BP − Early Pastoral (EP) = pastoralists,~8300-~6890 cal BP − Middle Pastoral (MP) = pastoralists,~7160-~5610 cal BP − Late Pastoral (LP) = pastoralists,~5700-~4650 cal BP The excavation campaigns of Takarkori lasted from 2003 to 2006. The rockshelter showed a particularly well-preserved sedimentary sequence and was excavated in four sectors over 143 sqm of total surface ( Figure 1). The stratigraphy of the site is reported in the interdisciplinary paper by Cremaschi et al. [20] who described its deposit as made by sediments rich in organic matter, charcoal and dried plant parts, human and faunal remains. Archaeobotanical analyses including pollen and macroremains have been carried out on sediment samples collected from every stratigraphical unit unearthed during the excavation.

−
Middle Pastoral (MP) = pastoralists, ~7160-~5610 cal BP − Late Pastoral (LP) = pastoralists, ~5700-~4650 cal BP The excavation campaigns of Takarkori lasted from 2003 to 2006. The rockshelter showed a particularly well-preserved sedimentary sequence and was excavated in four sectors over 143 sqm of total surface (Figure 1). The stratigraphy of the site is reported in the interdisciplinary paper by Cremaschi et al. [20] who described its deposit as made by sediments rich in organic matter, charcoal and dried plant parts, human and faunal remains. Archaeobotanical analyses including pollen and macroremains have been carried out on sediment samples collected from every stratigraphical unit unearthed during the excavation. Pollen analyses were carried out on sediments and coprolites taken from all the sectors [20,27,28]. The analyses on Poaceae pollen reported here derive from new observations of the slides already studied for palaeoenvironmental and ethnographic purposes. These detailed morphological analyses refer to 30 samples taken from the wall of TK-NS, a 4 × 2 m trench. Seven Units, including loose, organic or laminated sand, form the stratigraphy (Table 1). Moreover, 55 coprolites from the Main sector were also reconsidered: 45 specimens belong to Late Acacus layers, and 10 specimens belong to Middle Pastoral layers. The list of samples, with descriptions of the stratigraphical units, cultural phases, and chronology, is reported in Table 1.
Pollen extraction from sediments was made with a method including sieving and heavy liquid separation [20], while coprolites were boiled in 10% NaOH solution [28]. All treatments included acetolysis. Lycopodium tablets were added to calculate pollen concentration. Residues were included in glycerol to mount permanent slides. For each sample, one or two whole slides were observed, counting the monoporate and multiporate pollen grains of Poaceae. In previous studies [20], pollen spectra have been calculated on a pollen sum that included all pollen. Large pollen grains are the grass records with a maximum diameter ≥40 µm that is often considered as "Cerealia type". Nomenclature largely follows the African Pollen Database, and Chenopodiaceae is preferred to Amaranthaceae because it was used in previous papers.
Data were analyzed by Pearson's correlation coefficients and Student's t-test (p < 0.05), performed to examine the relationship between the presence of multiporate Poaceae pollen and environmental/anthropogenic stresses (variables considered: >40 µm Poaceae Pollen analyses were carried out on sediments and coprolites taken from all the sectors [20,27,28]. The analyses on Poaceae pollen reported here derive from new observations of the slides already studied for palaeoenvironmental and ethnographic purposes. These detailed morphological analyses refer to 30 samples taken from the wall of TK-NS, a 4 × 2 m trench. Seven Units, including loose, organic or laminated sand, form the stratigraphy (Table 1). Moreover, 55 coprolites from the Main sector were also reconsidered: 45 specimens belong to Late Acacus layers, and 10 specimens belong to Middle Pastoral layers. The list of samples, with descriptions of the stratigraphical units, cultural phases, and chronology, is reported in Table 1. Pollen extraction from sediments was made with a method including sieving and heavy liquid separation [20], while coprolites were boiled in 10% NaOH solution [28]. All treatments included acetolysis. Lycopodium tablets were added to calculate pollen concentration. Residues were included in glycerol to mount permanent slides. For each sample, one or two whole slides were observed, counting the monoporate and multiporate pollen grains of Poaceae. In previous studies [20], pollen spectra have been calculated on a pollen sum that included all pollen. Large pollen grains are the grass records with a maximum diameter ≥40 µm that is often considered as "Cerealia type". Nomenclature largely follows the African Pollen Database, and Chenopodiaceae is preferred to Amaranthaceae because it was used in previous papers.
Data were analyzed by Pearson's correlation coefficients and Student's t-test (p < 0.05), performed to examine the relationship between the presence of multiporate Poaceae pollen and environmental/anthropogenic stresses (variables considered: >40 µm Poaceae pollen-Asteraceae, Chenopodiaceae, Cyperaceae, Typha, all other Poaceae). Pollen data [20] together with multiporate Poaceae pollen were ordinated by Principal Component Analysis (using XLSTAT) to obtain further information on differences between samples in terms of plant composition.
In TK-NS, multiporate pollen was in 53% of the samples representing 0.1% of the total grasses (Table 2). Diporate and triporate were the most common with pore sizes which varied significantly from each other. In general, the exine and size of the different records were quite variable suggesting intraspecific variability and the possible presence of several species sharing the multiporate pollen anomaly ( Figure 2). Table 2. Takarkori rockshelter, Libya: Multiporate pollen observed in the TK-NS sequence (last column) and the number of Poaceae observed; the table also reports on the percentage of Poaceae in pollen spectra, their concentration (p/g), and the presence of large pollen of grasses (X = > 40 µm are > 1% in the spectrum). Data elaboration shows that there is not linear relationship between multiporate pollen and many other records [>40 µm Poaceae pollen (Pearson correlation coefficient r = −0.11), Typha (r = −0.19), Cyperaceae (r = 0.04) and Chenopodiaceae (r = 0.06)], while a negative correlation resulted between the multiporate and all other Poaceae (r = −0.38), and a positive correlation exists with Asteraceae (r = 0.44). Although from the Student's t-test the difference is not statistically significant, the value of r suggests that multiporate pollen tends to decrease when all other grasses increase, while tends to increase when shrubs and herbs of the daisy family, common in the arid environments of the desert [29], increase. Quaternary 2022, 5, x FOR PEER REVIEW 7 of 16 Although from the Student's ttest the difference is not statistically significant, the value of r suggests that multiporate pollen tends to decrease when all other grasses increase, while tends to increase when shrubs and herbs of the daisy family, common in the arid environments of the desert [29], increase.

Discussion
Since the occurrence of multiporate pollen in normally monoporate species is related to hybridization causing changes in reproductive biology [3,9], we argue that most Takarkori's Poaceae were polyploids probably also equipped with apomixis. Actually, the tribes Andropogoneae and Paniceae, largely in the form of spikelets and caryopses in the archaeobotanical record of this rockshelter [23] and abundant in other sites of the area (e.g., [30,31]), are known to include polyploids and also some obligately apomictic species (e.g., [32,33]).

Discussion
Since the occurrence of multiporate pollen in normally monoporate species is related to hybridization causing changes in reproductive biology [3,9], we argue that most Takarkori's Poaceae were polyploids probably also equipped with apomixis. Actually, the tribes Andropogoneae and Paniceae, largely in the form of spikelets and caryopses in the archaeobotanical record of this rockshelter [23] and abundant in other sites of the area (e.g., [30,31]), are known to include polyploids and also some obligately apomictic species (e.g., [32,33]).
In Takarkori, anomalous pollen was common in the samples. Multiporate pollen was found in the pollen sequence and in coprolites dated to~9500-5700 cal BP [20,23,26]. They seem to have been distributed in layers of different ages, mainly concentrated in the Late Acacus and late Middle Pastoral (see Table 2). The presence of multiporate pollen reveals that Poaceae that lived in central Sahara have tackled several environmental stresses, under climate or anthropogenic change pressures, during the early and middle Holocene (Table 3). Table 3. Takarkori rockshelter, Libya: Pollen samples from TK-NS pollen sequence with main results on pollen zones [20], and plant macroremains (spot and mix accumulations: [23]; data on palaeoenvironment are from the palynological study of Wadi Teshuinat [31].

Sample
No.

Pollen Zones
Plant Accumulations of Caryopses and Spikelets

Polyploidy as Effect and Adaptation to Environmental or Anthropogenic Stresses
Thanks to the conservation of exine in fossil pollen, polyploids can be identified in past contexts, giving paleoenvironmental research the possibility of obtaining hitherto unknown information on plant genetics of some archaeological records. Multiple genomes, which can be revealed by anomalies in pollen development, are in fact critical for the adaptation of plant species to external stresses.
To the best of our knowledge, only two cases of the discovery of these anomalies in archaeological deposits have been described in the literature. Radaeski et al. [4], studying fossil pollen from the Late Holocene Àguas Claras site in Brazil, advanced the hypothesis that multiporate pollen in monoporate populations had some ecological significance. They observed diporate pollen of Poaceae in the late Holocene (~3420 cal BP); then, this pollen increased from~1340 cal BP, when a vegetation change with expansion of forest occurred at the expense of grassland. The presence of diporate pollen was interpreted as the evidence that grassland was scattered even where forest vegetation dominated pollen spectra; in that context, their increase was linked to the stress determined by an increase in atmospheric humidity. In contrast, the presence of pollen from large Poaceae with many germination pores has been correlated with anthropogenic stress by Fall [34] who pointed out the appearance of multiporate pollen following the arrival of Polynesian settlers on the island of Tonga. In the swamps of this island, Poaceae percentage was low (1-2%) before the arrival of settlers (~2850 cal BP); then, grass pollen increases abruptly (up to 85%) with the arrival of settlers who introduced the practice of controlled burning of the swamps, as indicated by the presence of microcharcoals. Concurrently with this increase, multiporate pollen appeared, and their large size suggested that these were plants cultivated and so, in this case, the appearance of abnormal pollen would be related to selective hybridization and agriculture.
In the Takarkori rockshelter context, the Principal Component Analysis proved effective in interpreting the data as the components obtained have a clear ecological-environmental connotation (Figure 3). The first PC explains 10.39% (F1) of the total variance, the second PC explains 9.72% (F2). The first PC has low loadings with similar value for many pollen taxa and describes a high biodiversity in a context characterized by plants with a synanthropic value (e.g., Tamarix, Artemisia, and Poaceae > 40 µm), separating pollen spectra in "more anthropic" and "less anthropic". The second PC has high loadings for pollen taxa typical of arid or semi-arid environments (e.g., Cichorieae), while pollen taxa related to wet conditions have low loadings (e.g., Typha and Lemna). The second component separates pollen spectra along a humidity gradient, from "prevalently wet" to "prevalently warm/dry" conditions. The position of multiporate pollen in the Quadrant I suggests that drier climate conditions might have triggered an adaptive reply (the anomalous pollen) in Poaceae of Takarkori, regardless of human's presence. Many samples have relatively low scores for the second PC. Overall, the arrangement of samples in the PCA space seems to be influenced more by anthropogenic action and by seasonal wet phases, as is visible during the transition between the Late Acacus and the Early Pastoral phase (mainly Quadrant II and III). Later, in the Middle Pastoral phase, heterogeneity is observed among the samples with different ecological conditions (Quadrant I and IV).   [20] together with multiporate Poaceae pollen. The graph shows the sample points from dataset projected on the max variance plane given by PCA. The first axis separates pollen spectra along an anthropic pressure gradient (from "more anthropic" to "less anthropic"), and the second one along a climate gradient from prevalently wet to prevalently warm/dry conditions.

Multiporate Poaceae and the Climate Stress
In the TK-NS sequence, the highest values of multiporate pollen have been observed as concomitant with high values of Asteraceae (Figure 4: bottom), that are dated approximately from the early Holocene to the transition to the middle Holocene (~9500-8400 cal BP: point a in Figure 4) and in the second part of the middle Holocene (~5700 cal BP: point b in Figure 4).  [20] together with multiporate Poaceae pollen. The graph shows the sample points from dataset projected on the max variance plane given by PCA. The first axis separates pollen spectra along an anthropic pressure gradient (from "more anthropic" to "less anthropic"), and the second one along a climate gradient from prevalently wet to prevalently warm/dry conditions.
In the archaeobotanical record of Takarkori probably the concomitant stresses of climate instability and human selection have acted on suitable grasses that responded with polyploidy and multiporate pollen. Below the two determinants are discussed to emphasize their role in the development of anomalous pollen in the plant populations living in the Sahara at early and middle Holocene.

Multiporate Poaceae and the Climate Stress
In the TK-NS sequence, the highest values of multiporate pollen have been observed as concomitant with high values of Asteraceae (Figure 4: bottom), that are dated approximately from the early Holocene to the transition to the middle Holocene (~9500-8400 cal BP: point a in Figure 4) and in the second part of the middle Holocene (~5700 cal BP: point b in Figure 4).   (Table 3) and radiocarbon dates along the sequence (Table 1) are reported. The presence of Typha, marking local availability of fresh-water places, is reported in the bottom diagram. Points "a" and "b" mark the occurrence of multiporate pollen in samples characterised by high values of Asteraceae, including many shrubby and herb species adapted to steppic and dry savannah in the Sahara.
From geoarchaeological, archaeobotanical and zooarchaeological data, several climatic phases have been identified, with alternating cool and warm, arid and humid oscillations in the region (Table 3). Between ~9500 and 8100 cal BP, most of the calcareous tufa deposits reflecting a prolonged period of humid climate due to the northward expansion of the SW African monsoon system, were observed in Tadrart Acacus [35]. This massif, like other Saharan mountains especially sensitive to climate change, hosted plant migrations as a refugium and was covered by grasslands [31] (and see [36] for Tibesti and Hoggar). During this African Humid Period (AHP), tropical plant taxa migrated some 500 km northwards compared to their present-day distribution [37], and Mediterranean taxa de-  (Table 3) and radiocarbon dates along the sequence (Table 1) are reported. The presence of Typha, marking local availability of fresh-water places, is reported in the bottom diagram. Points "a" and "b" mark the occurrence of multiporate pollen in samples characterised by high values of Asteraceae, including many shrubby and herb species adapted to steppic and dry savannah in the Sahara.
From geoarchaeological, archaeobotanical and zooarchaeological data, several climatic phases have been identified, with alternating cool and warm, arid and humid oscillations in the region (Table 3). Between~9500 and 8100 cal BP, most of the calcareous tufa deposits reflecting a prolonged period of humid climate due to the northward expansion of the SW African monsoon system, were observed in Tadrart Acacus [35]. This massif, like other Saharan mountains especially sensitive to climate change, hosted plant migrations as a refugium and was covered by grasslands [31] (and see [36] for Tibesti and Hoggar). During this African Humid Period (AHP), tropical plant taxa migrated some 500 km northwards compared to their present-day distribution [37], and Mediterranean taxa descended to the central Sahara mountains as shown by remnants in the Tadrart Acacus and in Tassili N'Agger [31,38]. Similar patterns have been recognized in the analysis of faunal, especially aquatic, remains [39].
During the early to middle Holocene, the transition between the Late Acacus and the Early Pastoral phase corresponds to the transition from a cool and dry to a warm oscillation moving from dry to wet conditions locally. Contrary to the situation described by Radaeski et al. [4], the transition to a drier climate might have triggered an adaptive reply in the wild cereals of Takarkori. It would seem, therefore, that it was not the presence of a humid or arid climate but the change in environmental conditions that determined adaptive responses of polyploid grasses. Then, in the TK-NS sequence, relatively significant multiporate pollen were found during periods of environmental instability.
In the Pastoral phases, the AHP gave way to a progressive drying of the climate. In the Central Sahara, this phase is known to have begun from~6000 cal BP and led to an expansion of grassland and reduction of forested and tree savanna vegetation. At the end of the AHP, the Sahara experienced a major environmental crisis with the southward displacement of its desert margin [40], rapid decrease of precipitation at 5800-4800 cal BP [41], decline of wetlands [42,43], and related loss in biodiversity [37,44]. Pollen data from Takarkori shows that large stands of grasses moved from grasslands to grow closer to the contracted wet environments: their genotypes were largely equipped with high levels of genetic variability and plasticity for climate-responsive traits [23]. This alone may have been able to give them the ability to survive dramatic environmental transformations by enduring water shortages, at least in the early stages of transformation processes.

Multiporate Poaceae and the Anthropogenic Stress
Wild grains that have been accumulated in the Takarkori rockshelter include genera whose multiporate pollen has already been documented in the literature, such as Panicum [3] and Echinochloa [4]. These grasses may have pollen with large size, between 40 and >60 µm, and may often be included in the Cerealia type [31]. In fact, wild cereals are often marked by large pollen as well as the domesticated species. In the TK-NS sequence, large pollen of Poaceae was found in almost all samples ( Table 2) but there is not a linear correlation between multiporate and large pollen of Poaceae in the sequence.
In order to answer the question of whether human action may have been a stress agent for the production of multiporate pollen in wild polyploids, it may be useful to consider data from coprolites and from plant macroremains. The highest amount of multiporate Poaceae pollen in coprolites was found in Layer 275 (LA2,~9500-8600 cal BP): the morphometric analysis of the coprolites in this sample demonstrated that young Barbary sheep, together with very few adults, were kept in a stone-based enclosure [27]. Pollen analysis proved that the enclosure was likely made of reeds and other wetland plants like Phragmites and Typha, reinforced with stones at the base. Grasses were a relatively abundant plant resource easily available. The high accumulation of stems and other plant remains in the infilling of the enclosure suggests that grass fodder was largely transported to the penning area. To do this, humans repeatedly selected grasses. They made a pressure on wild cereals, thus triggering the production of multiporate pollen. The combination of pollen and archaeological evidence suggested to Rotunno et al. [27] that the fodder collected to feed the young animals was selected from high stands of weeds and wild cereals, mostly growing close to wet environments, known to the gatherers and repeatedly visited to collect wild cereals for food [23].
The abundance of wild cereal caryopses and spikelets at the site also documents an undoubted predilection and selection ability of certain nutrient-rich wild cereal species by Takarkori's inhabitants. Human exploitation of wild grasses, their accumulation and cultivation can be regarded as forms of anthropogenic stress. Accumulations of wild grains used as food and, likely, as fodder and litter, related to the Late Acacus and Middle Pastoral periods were common in the area. In Takarkori, the concentration of seeds was visible as accumulations of different species during the Late Acacus and Pastoral phases.
In the four millennia that wild grains have been exploited by the people of the Central Sahara, there is evidence that systematic harvesting and processing of these grains evolved more than once into cultivation. As it seems that wild species have never acquired all the characteristics proper to domestication, a cultivation without domestication was suggested by Mercuri et al. [23]. Nevertheless, anthropogenic pressure on the grains themselves caused noticeable changes in the morphology of the seeds of Panicum laetum, which became rounder, and Echinochloa colona, which became more elongated [23]. Jungle rice-E. colona was probably brought to the Takarkori area by pastoralists. The harvested plants were also used as fodder for animals, a kind of pressure that reasonably increased during the Pastoral period [20,23]. At the beginning of the Middle Pastoral period, the plant coexisted with a very similar species, liverseed grass-Urochloa panicoides; later, an increasing expansion of jungle rice at the expense of liverseed grass was favoured by domesticated sheep and goats overgrazing that facilitated Echinochloa, a more competitive, opportunistic, and invasive species, and by human selection [23]. Interestingly, intermediate forms and hybrids are currently known for Urochloa and Brachiaria [45], which are so similar that ambiguous morphologies of caryopses are common in the archaeobotanical records of the Tadrart Acacus area (e.g., [30,46]). The presence of multiporate pollen of Poaceae in the Middle Pastoral period could therefore also be explained in the context of continuous and increasing anthropic pressure by pastoralists, an action that took place during a long phase of climatic instability.

Multiporate Pollen as Palaeoecological Marker
Takarkori's multiporate pollen of Poaceae is an important pollen marker that for the first time could offer the possibility of obtaining palaeoecological information from an anomaly rather than from the normal morphology of the identified pollen.
Recently, a similar consideration was made on gymnosperms found in a 320 million years old deposit. There, a 3% limit of anomalous pollen was taken as an indicator of increased environmental stress. Based on field observations, one sac or three sacs malformation rate < 3% is normal for conifers under low stress; however, the exposure to UV-B radiation induced trees to produce 3-sac pollen. In the paleo-deposit, pollen with 3-sacs instead of 2-sacs was common, but a spike in 3-sacs pollen was interpreted as an indicator of ozone weakening events [47]. As in the case study of Takarkori, this research showed that certain types of malformations may be useful to indicate what kind of stress those plants were experiencing. It is likely that these anomalies have brought, and still can bring, benefits to the species that are capable of them. Variability in pollen tetrad morphology and the size of Classopollis, another extinct gymnosperm pollen type, revealed evidence for unreduced 2n pollen: the polyploid hybrids, which were more stress tolerant, may have reduced the ancient extinction risk of these plants during the environmental crisis at the Triassic-Jurassic transition, 200 million years ago [48].
In general, doubling of the genome may be accompanied by diversification in gene expression and epigenetic changes, which are important for controlling phenotypic plasticity and promoting adaptation to environmental stress conditions, such as rapid climate change [4]. Accordingly, the benefits of polyploidy could promote the adaptation of species to changing environments and their continuous transformation [49]. In the particular case of multiporate instead of monoporate pollen, studies on the mechanical model of pollen grain shape and exine show that the pores deform significantly more than the exine as the pollen swells. Therefore 'the presence of the pores relieves the stress on the exine and reduces it below the values which it would attain if the pores were not present' [50]. This can be of strategic importance to maintain pollen viability during rapid environmental changes or continuous stresses.

Conclusions
This article shows for the first time the use of anomalous pollen as an environmental marker (bioindicator) applied to Holocene archaeological contexts. The coexistence of plant species with phenotypic plasticity and some kind of stress is indisputable in the Takarkori chrono-cultural-sequence. Many grass species were involved, as evidenced on the one hand by the abundance of archaeobotanical records, especially Paniceae and Andropogoneae, and on the other hand by the morphological variety of the multiporate pollen grains.
However, it is not easy to determine whether climate change or anthropogenic pressures were the cause of the stress. This could have been in agreement with phases of natural environmental changes that involved major transformations in wadi and lake water availability and vegetation cover in the Sahara during the Holocene. Stress could also have been caused by the continuous selection of wild cereals harvested for food or fodder. The Takarkori pollen record suggests that climate change could have played a major role in the early Holocene. The occurrence of multiporate pollen and the increase of shrubby vegetation during dry periods can be considered primarily as consequences of the climatic changes during that period. During the periods of grass reduction and dryness expansion, polyploids were more common and humans did not hesitate to collect and transport to the site these plants that ensured consistent and predictable characteristics. The continuative anthropogenic pressure on Poaceae in the region brought about by their exploitation, cultivation, and the spread of pastoralism, maintained polyploid populations assuring food and fodder throughout environmental crises. Human pressure became relatively stronger even on polyploids during the late middle Holocene, eventually resulting in overexploitation that irreversibly reduced these resources [31].
Polyploids and amphiploids, which are one of the determinants of multiporate pollen formation, largely reproduce by apomixis. This mode of reproduction naturally occurs in about 35 families of angiosperms, and today is chosen in agriculture because it ensures the perpetuation of clonal genotypes through self-reproduction of seeds (agamospermy) and without the loss of vigor of the polyploid hybrid [51]. In the present period of profound and rapid climate change, affecting both wild and cultivated plant species, with significant repercussions on all human populations, the Poaceae family continues to have the same enormous attractiveness for humans, considering that it includes the grains that form the basis of human nutrition in most cultures. Takarkori has once again made it possible to gain a deeper insight into the history of African wild grains. In its archaeobotanical record, it could be possible to observe the development of multiporate pollen under climatic and human stress conditions over millennia. Understanding the significance of Poaceae multiporate pollen grains may allow a better knowledge of how plants have responded and are responding to exogenous stresses and what the consequences are for the ecology of some of them, including wild cereals.
Fornaciari who studied cereal macroremains from Takarkori, and Paola Torri and Michela Boscaini who noticed the multiporate Poaceae pollen in the coprolites from the site. We also thank Linda Olmi who studied the TK-NS sequence within her PhD research. Thank you to Juan Manuel Rubiales and two other anonymous referees for their constructive reviews.

Conflicts of Interest:
The authors declare no conflict of interest.