Abstract Pleistocene mammalian communities display unique features which differ from present-day faunas. The paleocommunities were characterized by the extraordinarily large body size of herbivores and predators and by their unique structure consisting of species now inhabiting geographically and ecologically distinct natural zones. These features were probably the result of the unique environmental conditions of ice age ecosystems. To analyze the ecological structure of Last Glacial and Recent mammal communities we classified the species into biome and trophic-size categories, using Principal Component analysis. We found a marked similarity in ecological structure between Recent eastern Altai-Sayan mammalian assemblages and comparable Pleistocene faunas. The composition of Last Glacial and Recent eastern Altai-Sayan assemblages were characterized by the occurrence of large herbivore and predator species associated with steppe, desert and alpine biomes. These three modern biomes harbor most of the surviving Pleistocene mammals. None of the analyzed Palearctic Last Glacial faunas showed affinity to the temperate forest, taiga, or tundra biome. The Eastern part of the Altai-Sayan region could be considered a refugium of the Last Glacial-like mammalian assemblages. Glacial fauna seems to persist up to present in those areas where the forest belt does not separate alpine vegetation from the steppes and deserts.

Citation: Pavelková Řičánková V, Robovský J, Riegert J (2014) Ecological Structure of Recent and Last Glacial Mammalian Faunas in Northern Eurasia: The Case of Altai-Sayan Refugium. PLoS ONE 9(1): e85056. https://doi.org/10.1371/journal.pone.0085056 Editor: Michael Hofreiter, University of York, United Kingdom Received: April 22, 2013; Accepted: November 25, 2013; Published: January 13, 2014 Copyright: © 2014 Pavelková Řičánková et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. Funding: The study was supported by Czech Ministry of Education (MSM 6007665801) http://www.msmt.cz/index.php?lang=2 and Czech Science Foundation (# P504/11/0454) http://www.gacr.cz/en/. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing interests: The authors have declared that no competing interests exist.

Introduction The unique structure of Pleistocene mammalian communities has drawn the attention of scientists for many decades [1]–[7]. The extraordinary large body size of some Pleistocene mammals (e.g., mammoths, giant deer, or cave bear) and composition of the ice age communities have no analogies in the present-day faunas [8]–[9]. The Pleistocene communities consisted of species which now inhabit geographically and ecologically distinct natural zones (tundra, forest, steppe, savanna) [5], [8]. Arctic reindeer and musk-ox thus lived in sympatry with e.g. steppe horse and bison and/or with forest elk and roe deer [4], [10]. The reason for such a peculiar structure of Last Glacial (corresponding with Weichselian Glaciation) mammalian assemblages is probably associable with the unique environmental conditions of ice age ecosystems [8], [10], [11]. The non-taxonomic, ecological structure of a mammalian community (∼ its guild composition according to e.g. Simberloff & Dayan [12]) is determined mainly by environmental factors such as climate, type of biome or vegetation heterogeneity [13]–[15]. Each type of ecosystem (e.g. woodland or arid ones) is characterized by a specific trophic-size structure of its mammalian community [5], [13], [16]. Mammalian communities from areas with similar regional climates tend to converge to similar community structures [15]. Historical factors also play an important role in forming community structure, but they operate on a different, probably much longer timescale than environmental factors [14]–[15]. The non-analogue Last Glacial communities evolved in a cold and dry continental climate, which supported highly heterogeneous vegetation and landscape structure, usually described as tundra-steppe or mammoth steppe [4], [8], [10], [17]. The Pleistocene tundra-steppe ecosystem was quite heterogeneous locally but displayed a relatively high degree of homogeneity on the continental scale. This ecosystem covered wide areas of the northern part of the globe, thrived for approximately 100,000 years without major changes, and then suddenly went extinct about 12,000 years ago [4], [10]. Surprisingly, environmental conditions similar to the Last Glacial period have been found in the Central Eurasian Altai-Sayan mountains [18]. This climatic analogue has recently been supported by biological data. Recent findings of the paleo-biome reconstruction [19]–[20] and pollen-analytical research [21]–[23] suggest that present-day Altai-Sayan landscapes could be considered the closest modern analogy to the Last Glacial environments. The area is currently inhabited by mollusc assemblages that were characteristic of full-glacial environments across large areas in Eurasia but went extinct in the regions that experienced considerable climatic change, namely in Europe [24]. Simulated paleovegetation maps based on paleoclimatic models and plant functional types have also suggested considerable stability in Central Eurasia over the last 40,000 years [25], [26]. Detailed analysis of the Altai Late Pleistocene assemblages of small mammals revealed that no significant changes occurred between the cold phase of the Pleistocene and the Holocene [27]. The environment of this region can thus be considered as conservative and stable. In this study, we compared the ecological structure of Recent Altai-Sayan mammalian assemblages to the ecological structure of Last Glacial fauna of Altai-Sayan and several adjacent regions as well as to Recent mammalian communities from various natural zones of northern Eurasia. In order to examine the most important structural characteristics of Pleistocene assemblages, we assigned individual species according to their biome and the trophic-size categories. Given the analysis of Willis et al. [28] and Rodrigues et al.[7], [29] showing that glacial vegetation consisted of steppe, tundra and forest and mammalian fauna was characterized by large herbivores and predators, we would expect that glacial communities are characterized (and differ from Recent communities) by the co-occurrence of steppe, tundra, and forest species [28], and by high a proportion of large herbivores and predators [29]. We hypothesize that, given the environmental stability of central Eurasia [18], Altai-Sayan Recent assemblages will be more similar to the glacial communities than to any other Recent community [22]–[23].

Discussion We found a marked similarity in ecological structure between Recent Altai-Sayan mammalian assemblages and Last Glacial paleocommunities. The ecological structure of mammalian communities confirmed the possible persistence of the Last Glacial-like fauna in the present-day Altai-Sayan region, as the Recent assemblages of eastern Altai-Sayan (i.e. Ukok-Sailiugem, Khar Us Nuur, and Uvs Nuur) were more similar to various glacial localities than to the communities of any other Recent areas. Glacial communities have close modern analogues in the three eastern Altai-Sayan areas where e.g. reindeer and saiga antelope still live in sympatry [38]. Our results are congruent with other evidence supporting the persistence of Pleistocene biota in the Altai-Sayan region. Rodríguez [7] reported a similarity between the ecological structure of present-day mammalian communities from Central Eurasia (Uvs Nuur and Great Gobi) and Iberian Pleistocene communities. The Altai Mountains are an important refugium for full-glacial snail faunas, as recently documented by Horsák et al. [24] and Hoffmann et al. [39]. Kuneš et al. [22] and Pelánková & Chytrý [23] demonstrated a close similarity between glacial pollen samples from central Europe and modern surface-pollen spectra from the Altai-Sayan region. Fossil pollen spectra from the Altai and adjacent regions indicate little difference between modern biomes of this region and those reconstructed for the Last Glacial Maximum [20]. Similar conditions possibly occur in climatically stable areas of North American cold deserts. Fossil evidence from the Great Basin indicate that Pleistocene plant assemblages are comparable with the modern ones [40]. According to vegetation studies, three major biomes occurred widely in the Pleistocene: steppe, tundra, and taiga [28], [41]–[42]. In contrast, the composition of the Last Glacial (and Recent eastern Altai-Sayan) mammalian faunas was characterized by the co-occurrence of steppe, desert and alpine species. These three modern biomes harbor most of the surviving Last Glacial mammals. The importance of the desert biome was probably more pronounced in the examined region of southern Russian plain and central Asia in comparison to the northern and western part of Eurasia. Glacial-like communities still persist in areas where the forest belt does not separate alpine vegetation from the steppes and (semi)deserts. High aridity in eastern Altai-Sayan restricts forests to isolated patches with higher soil moisture [43]–[44]. Holocene fragmentation of the alpine grasslands to isolated patches surrounded by forests could have lead to the extinction of some large mammals typical of the mammoth steppe. Sharp separation of originally intermixed faunal elements into distinct biomes seems to be the major pattern of the Last Glacial/Holocene transition [8], [45]. Most previous attempts to find a modern analogue of the mammoth steppe have been focused on regions of the Arctic tundra, i.e., Yakutia, NE Russia, and Alaska [8], [9], [46], [47]. However, present-day faunas of the Arctic regions have diverged greatly in their ecological structure from Last Glacial ones, showing more affinity to the taiga biome. The divergent ecological structure of modern Yakutian fauna is determined mainly by Holocene immigrant species, characterized by a high proportion of arboreal herbivores and predators of invertebrates associated with forest and taiga biomes (Fig. 4). The modern arctic tundra is characterized by low productivity and relatively high homogeneity [4], [25], in contrast to the alpine biome where tundra vegetation occurs in association with forests and steppes, e.g. [48]–[49]. During the Last Glacial, tundra vegetation was confined to places with higher precipitation or lower evapotranspiration, e.g. in the mountains or at non-glaciated higher latitudes, often in the close vicinity of steppe [20]. Taiga and temperate forest species occurred in all of the examined Last Glacial assemblages; however their percentage was very low, with the exception of the Last Glacial assemblage of the Carpathian Mountains. This assemblage holds a special position as the site resembles Recent rather than any other Last Glacial community. The Carpathian Mountains of eastern Europe possibly represented a glacial refugium of forest vegetation and forest-dwelling animals, as suggested by Sommer & Nadachowski [50], Jankovská and Pokorný [21], Markova et al. [51], Willis et al. [28], and Willis & van Andel [52]. The main difference in trophic-size structure between present-day and fossil assemblages is the higher richness of large mammals and proportionally lower richness of small mammals in the latter. The high proportion of large herbivores observed in the Last Glacial as well as in some present-day communities is generally typical of areas with low tree cover [15] and could be further supported by the year-round availability of high quality food in the glacial steppe [10], [53]. Fossil assemblages are probably biased against small mammal species owing to fossil record incompleteness [54], [7], [29], [55]–[56]. However, the low proportion of small mammals cannot be considered completely artifactual because the present-day assemblages from the eastern Altai-Sayan region show very similar composition to Last Glacial communities. The low species richness of insectivores and aquatic predators in Last Glacial-like Altai-Sayan assemblages (and probably in other Last Glacial assemblages as well) could be due to the dry and cold climate associated with permafrost which strongly limits insectivores’ food sources [57, 48 58]. These results were confirmed using two different ecological classifications and using two independent statistical methods (NMDS, PCA). Moreover, our fossil datasets are based on large areas and a long time period to avoid taphonomic biases in species occurrence. Therefore, we suggest that potential bias in the Pleistocene data subset cannot significantly affect our results. In general, Pleistocene assemblages were characterized by the occurrence of large herbivore and predator species associated with steppe, desert and alpine biomes. Mammalian paleocommunities classified according to biome type are relatively homogeneous, confirming the view of the mammoth steppe as a single Last Glacial biome [10]. In contrast to biome classification, the trophic-size structure of mammalian paleocommunities shows a degree of heterogeneity, comparable to the Recent localities. The trophic-size structure of communities is probably less influenced by historical factors [16]. Historical processes are considered to be the main factor promoting differences between communities from similar environments [15], [29]. Our results open new research possibilities for many aspects of Quaternary paleoecology. The Altai-Sayan region offers a possibility to study factors shaping the structure of so called non-analogue communities and explore vegetation and faunal changes associated with Pleistocene/Holocene transition. Research of soil nutrient availability and cycling in glacial-like localities can provide an insight into the carrying capacity of ice age ecosystems supporting numerous large herbivore species. Modeling the impact of climate changes on the glacial-like landscape may elucidate the process of biome diversification in the Holocene. Our result can be confirmed by thorough paleontological research of the as yet unexplored eastern Altai-Sayan region as well as by phylogeographical analyses of typical glacial species (e.g. steppe lemming or pika) including Altai-Sayan populations.

Acknowledgments We thank David Storch, Martin Konvička, Daniel Frynta, Ivan Horáček, Vojtěch Novotný, Michal Horsák, Jan Zrzavý, Milan Chytrý, John Stewart and anonymous referee who helped greatly with useful comments. We thank Ingrid Steenbergen for language corrections and Martin Hais for maps.

Author Contributions Conceived and designed the experiments: VPŘ JRi. Analyzed the data: JRi. Wrote the paper: VPŘ Jro JRi. Performed the literature search and data extraction: JRo VPŘ.