Research article | Open | Published:
Global climate changes drive ecological specialization of mammal faunas: trends in rodent assemblages from the Iberian Plio-Pleistocene
BMC Evolutionary Biologyvolume 13, Article number: 94 (2013)
Several macroevolutionary hypotheses propose a synchrony between climatic changes and variations in the structure of faunal communities. Some of them focus on the importance of the species ecological specialization because of its effects on evolutionary processes and the resultant patterns. Particularly, Vrba’s turnover pulse hypothesis and resource-use hypothesis revolve around the importance of biome inhabitation. In order to test these hypotheses, we used the Biomic Specialization Index, which is based on the number of biomes occupied by each species, and evaluated the changes in the relative importance of generalist and specialist rodents in more than forty fossil sites from the Iberian Plio-Pleistocene.
Our results indicate that there was a decrease in the specialization degree of rodent faunas during the Pliocene due to the global cooling that triggered the onset of the glacial events of the Cenozoic (around 2.75 Ma). The subsequent faunal transition after this critical paleoenvironmental event was characterized by an increase of specialization related to the adaptation to the new environmental conditions, which was mainly associated with the Pleistocene radiation of Arvicolinae (voles).
The pattern of faunal turnover is correlated with the development of the modern glaciations in the Northern Hemisphere around 2.75 Ma, and represents a reorganization of the rodent communities, as suggested by the turnover pulse hypothesis. Our data also support the resource-use hypothesis, which presumes the role of the degree of specialization in resources specifically related to particular biomes as a driver of differential speciation and extinction rates. These results stress the intimate connection between ecological and evolutionary changes.
Almost since the first publication of the seminal contribution about the tendency of species to form varieties , there has been a debate between models that consider the competition between species as the key for evolutionary changes [2, 3] and the ones that regard factors external to the species, usually biogeography or climatic changes, as the main drivers responsible of biotic evolution [4–7]. Nevertheless, an increasing consensus has appeared during the last decades on the critical importance of changes in the physical environment, rather than biotic interactions themselves, for the evolution of organisms and ecosystems at large spatial and time scales [8–11].
The habitat theory proposed by Elisabeth S. Vrba [12, 13] is one of the best known of such evolutionary models. This theory comprises a set of hypotheses showing the influence of global climatic shifts and the subsequent environmental changes on the turnover of species assemblages as a result of the drifting of geographic distributions, lineage originations and extinctions.
As a core part of the habitat theory, Vrba [7, 14] developed the resource-use hypothesis, which stresses the relationship between the ecological specialization of species and the processes that regulate their evolution. The different conditions imposed by the ecological characteristics of species, and especially the biomes that they inhabit, have an effect on the macroevolutionary patterns that are observed through time and space [7, 14–17]. In this way, the resource-use hypothesis gives a great value to the differences between the species that inhabit one biome exclusively (biome specialists or stenobiomic species) and others that obtain their resources from more than one biome (biome generalists or eurybiomic species). In order to clarify the concepts used for this work, it is important to indicate that biome specialization is not necessarily related to specialization in other ecological traits of the species. For example, Rhynchomys soricoides (Murinae, Rodentia), which has been described as vermivore , exhibits a high grade of specialization in its dietary habits but, since it inhabits three different biomes (equatorial rainforest, tropical deciduous woodland and temperate evergreen forest) [19–21], is a biome generalist. On the other hand, Stochomys longicaudatus (Murinae, Rodentia) has an omnivore diet , but it is restricted to the equatorial rainforest of Central Africa .
According to the resource-use hypothesis, during episodes of climatic triggering of habitat change, specialist species are more prone to suffer limitation of their resources and, consequently, they are more susceptible to habitat fragmentation, vicariance and directional selection. Since environmental changes have stronger effects on biome specialists than on generalist species, which can find their resources in different biomes, the former are predicted to have higher speciation and extinction rates than the latter [14, 17]. Such phenomenon results in an increase of the specialist species versus generalists in the global fauna, which has been observed in mammalian assemblages from both Africa [14, 15] and South America  as well as in the ruminants at the global scale .
Importantly, because of the constraints on the evolutionary history of species imposed by changes in the physical environment, most of the evolutionary changes in biotic lineages (including speciation, extinction or dispersion) should occur synchronically with global climatic changes. Such an evolutionary scenario is developed in the turnover-pulse hypothesis [12, 24, 25]. According to this hypothesis, most speciations and extinctions across diverse groups of organisms are not randomly distributed in time, but show statistically significant concentrations near times of major physical change. While most of these turnover-pulses affect few lineages and/or restricted geographic areas, some of them are massive and of global extent .
The combination of all these issues, along with the hierarchical organization of ecological and evolutionary processes , raises three corollary predictions (Figure 1): (1) due to their large ecological preferences, biome generalists might constitute a predominant part of the set of species that survive during the moments of significant global change; (2) nevertheless, once the critical periods of global change have finished, the preponderance of generalist species will decrease and a set of new specialist species will develop through speciation of the surviving lineages as the environmental conditions are stabilized; (3) this would entail a progressive specialization of generalist clades through niche filling within the newly generated environments. On the other hand, due to ecological constraints, clades dominated by specialist species before the crisis are expected to severely decrease in importance, although their specialization degree may not change.
To test these predictions, here we focus on the rodent faunas from the Iberian Plio-Pleistocene. During this period, successive cooling pulses culminated with the establishment of continental northern-hemisphere glaciations and the modern ice age climate around 2.75 Ma [27–29]. This severe climatic event at the global scale is an ideal scenario for testing Vrba’s hypotheses. We used rodents as our study group because they are widespread, highly diverse and habitat-sensitive, which makes them one of the most environmental and climatically informative groups of mammals [30–34]. Besides, their Iberian fossil record is vast and extensively documented . Additionally, this group has been used to characterize intervals of great faunal change throughout the Cenozoic, usually associated with global climate fluctuations (e.g. [36, 37]). For these reasons they provide a suitable faunal set to test the predictions exposed above.
In this work we used the faunal lists of rodent communities from 44 fossil sites from the Iberian Plio-Pleistocene (Table 1) dated between 5.25 and 0.01 Ma . These fossil sites have been subject to intensive sampling during the last fifty years (see references in ). The minimum sample required to include a fossil site in our study was 100 first and second upper and lower molars, which is considered the minimum number necessary to render a representative sample of the original paleocommunity [30, 36, 39].
The specialization degree of each species was measured using the Biomic Specialization Index (BSI) developed by Hernández Fernández and Vrba . This index indicates the number of biomes inhabited by the species, following the climatic classification of Walter , which recognizes ten biomes. Therefore, BSI equals 1 for most specialized species whereas generalist species could exhibit a BSI as high as 10. The data on the biome residence for all the rodent species were obtained from , who derived biome residences from identifying their living ecological analogues as estimated by ecomorphological studies of the dentition [41, 42].
Finally, following Hernández Fernández and Vrba , we calculated the relative frequency of specialist and generalist species in each fossil site in terms of the average value of the BSI of the species found in the site. Taking into account that some taxa on our lists were not identified to species level, we decided to conduct the analysis by two different ways. First, for the undetermined taxa, the BSI value was calculated as the average of all the species that belong to their upper taxonomic level. Second, to avoid the potential noise in the data due to unidentified taxa, we only analysed the taxa that were determined at species level in each fossil site. Since both analyses yielded very similar results, here we only discuss those corresponding to the latter version of the analysis (all taxa at the species level).
Following our predictions derived from the interaction between the resource-use and the turnover-pulse hypotheses, the beginning of the modern glaciations in the Northern Hemisphere (around 2.75 Ma) should be coincident with the rising of mammal faunas dominated by generalists (high average BSI values). After this critical event we should find a progressive increase in the specialization degree of the rodent faunas from the Iberian Pleistocene (decrease of average BSI). In order to identify significant changes in the average BSI of fossil assemblages through time, we carried out Multivariate Adaptive Regression Splines (MARS; ) on the data set. MARS is a method that identifies hinge points in a time series that automatically minimize the residual sum of squares (RSS). In this way time shifts in BSI trends do not need to be fixed a priori. Progressively adding hinge points to the model increases the number of parameters. Since overly complex models may result in stochastic error (inflated variance), we need to find an equilibrium between fit and complexity. MARS does this automatically using the Akaike Information Criterion (AIC) scores of each model, which measures the goodness of the fit of a statistical model while penalizing the number of parameters (the complexity).
Additionally, we expected to find different evolutionary responses to this event in different rodent clades according to their degree of specialization. More generalist groups are predicted to survive and proliferate after the environmental crisis at 2.75 Ma, becoming progressively more specialized. Meanwhile clades dominated by specialist species before the crisis are expected to decrease in importance, although the specialization degree may not change. These predictions were tested by means of t-Student comparisons between species pools before and after the inflection point yielded by the MARS analysis within different taxonomic groups. Following Wilson & Reeder , we studied the specialization trend of Sciuridae, Gliridae, Castoridae, Arvicolinae, Cricetinae, Gebillinae, Murinae and Hystricidae separately. Subfamilies within Cricetidae (arvicolines and cricetines) and Muridae (gerbillines and murines)  were studied independently due to their importance in Iberian Plio-Pleistocene rodent faunas, and because they have been traditionally taken as independent families in paleontological studies  and their monophyly has been demonstrated by molecular studies [46, 47].
Through the time span considered here, our results show the existence of two opposed patterns with an inflection point coincident with the onset of the Pleistocene glaciations. This involves substantial changes in the average BSI value of the Iberian rodent faunas (Figure 2), supporting the prediction of a change in the pattern of ecological specialization in concert with the cooling pulse of the Plio-Pleistocene around 2.75 Ma.
We observed a considerable increase of the average BSI throughout the Pliocene, reaching maximum values close to the end of the period. This is to say; during the Pliocene there was a transition from faunas with a higher prevalence of specialist species to faunas where biome generalists were predominant. The best-fit model identified by MARS analysis (Table 2) includes a single hinge point in the BSI time-series that divides the time series in two sub-sets differing in their trend. The first sub-set of data ranges from 5.3 Ma to 2.9 Ma, depicting an increase in BSI values (towards less specialists; see Figure 2). The second sub-set ranges from 2.6 Ma to the present and shows a decrease of the BSI values (towards more specialists; see Figure 2). Adding this hinge point to the model resulted in a significant improvement of the fitness over the null model where specialization (BSI) remains linear through time (ΔAIC = 13.28, see Table 2). Increasing the number of hinge points did not improve the fit of the models (Table 2).
The t-student tests for different rodent groups indicated that the radiation of arvicolines (voles) is related to a significant change in biomic specialization, from few generalists to many specialist species, associated to the critical environmental change at 2.75 Ma (Figure 3). Other Pliocene generalist groups (beavers, gerbils and porcupines) appear to be scarce in Iberian faunas. More specialized groups (sciurids, glirids, cricetines and particularly murines) showed much lower number of species during the Pleistocene than in the Pliocene, although they did not show differences between their species BSI before and after the development of the first glaciations of the Pleistocene (Figure 3). These results corroborate the third prediction tested here.
The gradual decrease in the relative importance of specialist species during the Pliocene (increase in average BSI) is probably related to the progressive global cooling predominant during the Pliocene, which triggered modifications in the climate of southern Europe, from subtropical conditions with minor fluctuations of temperature to temperate conditions with noticeable annual thermal seasonality [49, 50].
Our results are consistent with a scenario where species dwelling in more than one biome (generalists) were more able to tolerate this change due to their ability to find resources in different environments. After the development of the first modern glacial events at the Northern Hemisphere the specialization degree of the faunas progressively increased (the value of average BSI decreases). This trend reflects a gradual “recovery” of the rodent communities after the ecological disturbance caused by the onset of the first glaciations. Finally, during the Holocene the ecological specialization (average BSI) of the faunas reached similar values to the ones found during the early Pliocene.
There were profound changes in the faunal composition of the rodent assemblages (Figure 4) stemming from the development of the temperate climates of the Plio-Pleistocene . During the Pliocene, Iberian rodent faunas were dominated by murines and cricetines, with a set of companion species within terrestrial squirrels and gerbils. In contrast, the dominant group during the Pleistocene was Arvicolinae (voles) . Most of the correlations between the percentage of species in each family and the average BSI value of each fossil site before and after the modern glaciations are not significant (Table 3). Nevertheless, it is noteworthy that Arvicolinae shows a significant relationship with different sign in each period. Before 2.75 Ma, the moderate increase in the number of arvicoline species is associated with higher average BSI (more generalists). Yet, after the climatic crisis, the flowering of arvicoline faunas is associated to a significant decrease in average BSI of the rodent communities. This signifies that the increase in arvicoline diversity likely resulted from an adaptive radiation involving highly stenobiomic taxa. Furthermore, the radiation of voles had a preponderant role in the reorganization of the rodent faunas from the Iberian Plio-Pleistocene. During the Pliocene, this group is mainly represented by generalist species in Promimomys, Dolomys and Mimomys, and an important part of the survivors during the climatic “deterioration” of the middle-late Pliocene belong to the latter genus. During the Pleistocene the large radiation of arvicolines is dominated by the evolution of specialist species in many different genera, especially in Microtus[53, 54]. Probably, this is related to the development of several new environments in the temperate latitudes , where resources were available for the exploitation by the generalist survivors of this family, which could start a process of progressive specialization along the Pleistocene. This is possibly one of the keys for the evolutionary success of this family (Figure 3), which now inhabits all the ecosystems of the Holarctic and is the dominant group of rodents in most of them.
On the contrary, the other rodent groups of the Iberian Plio-Pleistocene were not able to take advantage of the new environments that were associated with the glacial-interglacial cyclicity of the Pleistocene climate. Particularly, in the case of Murinae there was a substantial decrease in the relative importance of this group in the rodent faunas from the Iberian Peninsula (Figure 4). From communities with representatives of eight genera during the Pliocene , the Pleistocene faunas only preserved species of Apodemus, Castillomys, Micromys and Stephanomys, the latter only surviving until 2.0 Ma [57, 58]. Most murine species were adapted to forested and warm environments [34, 42], which disappeared from Iberian latitudes with the onset of the modern glaciations. Today this family is predominantly distributed in tropical areas of the Old World, with only a few genera and species in temperate regions . Murine demise in the Iberian Peninsula would be due not only to the decrease in temperature, but probably also to changes in precipitation amount and seasonality, associate to the development of Mediterranean climate . Although such species loss was conspicuous in the Iberian Pleistocene, the specialization degree of this group did not show significant differences after 2.75 Ma (Figure 3).
Cricetinae was the other group showing a clear impoverishment over the analyzed interval. Whereas cricetid communities during the Pliocene held the genera Apocricetus, Blancomys, Ruscinomys, Celadensia or Trilophomys, Pleistocene cricetid faunas were only represented by the immigrant Allocricetus bursae. This taxon is closely related to the few genera that occupy today the arid environments present in the temperate regions of Eurasia since the onset of the glaciations. Apparently, the global cooling at the Plio-Pleistocene boundary imposed strong ecological limitations to the species adapted to arid and open environments, such as the Iberian cricetids [42, 61]. Similarly, the terrestrial squirrels, characteristic of many Pliocene faunas from the Iberian Peninsula (Atlantoxerus), belong to Xerini, a group that today is restricted to the African savannas and semideserts, and disappeared completely from the Iberian record due to the cooling of the climate. They were replaced in the Iberian Pleistocene by Sciurini (Sciurus) and Marmotini (Marmota); these groups are predominant respectively in forested and open temperate environments. In the case of glirids, there were no large faunal changes before and after the beginning of the Pleistocene northern glaciations. They maintained the same genera, which are still living, and there are no significant differences in the BSI of their species over the period studied here (Figure 3).
Finally, three minor Pliocene generalist groups, gerbils, beavers and porcupines, were unable to proliferate in the Iberian Pleistocene. In the first case this is probably due to their specific specialization to arid climates mostly developed in other regions of Eurasia. The other two families show extremely low species number, which may have hampered their diversification.
It seems that the Pliocene witnessed a process of disassembly within the Iberian rodent communities. The impoverishment of the Iberian rodent communities throughout the Pliocene stemmed from a loss of diversity that affected most of the rodent families. The recovery of species diversity during the Pleistocene was linked to the radiation of stenobiomic arvicolines that resulted in the establishment of a new rodent fauna (Figure 4). This global radiation is evidenced within the Iberian fossil record by local speciation (e.g. development of the endemic lineage Iberomys within Microtus[62–64]) as well as by immigration of new species from other Eurasian regions, being both processes spurred by global environmental change. Overall, such ecological reorganization of assemblages appears to be triggered by global climatic changes and modulated by the differences in ecological specialization of the implied species, similarly to what was observed in earlier periods of faunal replacement in Spain [36, 65].
Our results offer support for some hypotheses included in Vrba’s habitat theory, which predicts a proportional decrease of specialist species associated with severe global climatic changes, and a later recovery of this kind of species associated to a complete faunal turnover. The pattern of faunal turnover is correlated with the development of the modern glaciations in the Northern Hemisphere around 2.75 Ma, which triggered a reorganization of the rodent communities as predicted by the turnover pulse hypothesis. In the same way our data support the resource-use hypothesis, which stresses the role of the degree of specialization in resources specifically related to particular biomes as a driver of differential speciation and extinction rates. In this way, ecological and evolutionary changes are intimately connected.
Finally, this work shows that the exceptional quality of the fossil record in the rodent assemblages from the Iberian Peninsula makes this group a good case for the study of the relevance of ecological characteristics of species in the development of macroevolutionary processes.
Darwin CR, Wallace AR: On the tendency of species to form varieties; and on the perpetuation of varieties and species by natural means of selection. J. Proc. Linn. Soc. 1858, 3: 46-50.
Darwin CR: On the Origin of the Species by Means of Natural Selection, or the Preservation of Favoured Races in the Struggle for Life. 1859, London: John Murray
Van Valen L: A new evolutionary law. Evolutionary Theory. 1973, 1: 1-30.
Gulick JT: On the variation of species as related to their geographical distribution, Illustrated by the Achatinellinae. Nature. 1872, 6: 222-224.
Gulick JT: On the diversity of evolution under one set of external conditions. J. Proc. Linn. Soc. 1873, 11: 496-505. 10.1111/j.1096-3642.1873.tb01670.x.
Wallace AR: Darwinism: An exposition of the Theory of Natural Selection, with some of its applications. 1912, London: Macmillan
Vrba ES: Evolution, species and fossils: how does life evolve?. S Afr J Sci. 1980, 76: 61-84.
Stenseth NC, Maynard Smith J: Coevolution in ecosystems: red queen evolution or stasis?. Evolution. 1984, 38: 870-880. 10.2307/2408397.
Barnosky AD: Distinguishing the effects of the red queen and court jester on Miocene mammal evolution in the Northern Rocky Mountains. J Vertebr Paleontol. 2001, 21: 172-185. 10.1671/0272-4634(2001)021[0172:DTEOTR]2.0.CO;2.
Badgley C, Barry JC, Morgan ME, Nelson SV, Behrensmeyer AK, Cerling TE, Pilbeam D: Ecological changes in Miocene mammalian record show impact of prolonged climatic forcing. PNAS. 2008, 105: 12145-12149. 10.1073/pnas.0805592105.
Benton J: The Red Queen and the Court Jester: species diversity and the role of biotic and abiotic factors through time. Science. 2009, 323: 728-732. 10.1126/science.1157719.
Vrba ES: Mammals as a key to evolutionary theory. J Mammal. 1992, 73: 1-28. 10.2307/1381862.
Vrba ES: On the connections between paleoclimate and evolution. Paleoclimate and Evolution with emphasis on Human Origins. Edited by: Vrba E, Denton GH, Partridge TC, Burkle LH. 1995, New Haven: Yale University Press, 24-45.
Vrba ES: Ecology in relation to speciation rates: some case histories of Miocene-Recent mammal clades. Evol Ecol. 1987, 1: 283-300. 10.1007/BF02071554.
Hernández Fernández M, Vrba ES: Macroevolutionary processes and biomic specialization: testing the resource-use hypothesis. Evol Ecol. 2005, 19: 199-219. 10.1007/s10682-004-8152-7.
Moreno Bofarull A, Arias Royo A, Hernández Fernández M, Ortiz-Jaureguizar E, Morales J: Influence of continental history on the ecological specialization and macroevolutionary processes in the mammalian assemblage of South America: differences between small and large mammals. BMC Evol Biol. 2008, 8: 97-10.1186/1471-2148-8-97.
Cantalapiedra JL, Hernández Fernández M, Morales J: Biomic specialization and speciation rates in ruminants (Cetartiodactyla, Mammalia): a test of the resource-use hypothesis at the global scale. PLoS One. 2011, 6 (12): 1-10. e28749
Rowe KC, Reno ML, Richmond DM, Adkins RM, Steppan SJ: Pliocene colonization and adaptive radiations in Australia and New Guinea (Sahul): Multilocus systematics of the old endemic rodents (Muroidea: Murinae). Mol Phylogenet Evol. 2008, 47: 84-101. 10.1016/j.ympev.2008.01.001.
Musser GG, Carleton MD: Superfamily Muroidea. Mammal Species of the World: A Taxonomic and Geographic Reference. Edited by: Wilson DE, Reeder DM. 2005, Baltimore: The Johns Hopkins University Press, 894-1531.
Heaney L, Balete DS: Rhynchomys isarogensis. IUCN 2012: IUCN Red List of Threatened Species Version 2012.2. 2012, http://www.iucnredlist.org,
Heaney L, Tabaranza B, Balete DS: Rhynchomys soricoides. IUCN 2012: IUCN Red List of Threatened Species Version 2012.2. 2012, http://www.iucnredlist.org,
Gomes Rodrigues H, Merceron G, Viriot L: Dental microwear patterns of extant and extinct Muridae (Rodentia, Mammalia): ecological implications. Naturwissenschaften. 2009, 96: 537-542. 10.1007/s00114-008-0501-x.
Boitani L: Stochomys longicaudatus. IUCN 2012: IUCN Red List of Threatened Species Version 2012.2. 2012, http://www.iucnredlist.org,
Vrba ES: Environment and evolution: alternative causes of the temporal distribution of evolutionary events. S Afr J Sci. 1985, 81: 229-236.
Vrba ES: Turnover-pulses, the Red Queen, and related topics. Am J Sci. 1993, 293: 418-452. 10.2475/ajs.293.A.418.
Vrba ES: What are the biotic hierarchies of integration and linkage?. Complex organismal functions: integration and evolution in vertebrates. Edited by: Wake DB, Roth G. 1989, New York: John Wiley & Sons, 379-401.
Zachos J, Pagani M, Sloan L, Thomas E, Billups K: Trends, rhythms, and aberrations in global climate 65 Ma to present. Science. 2001, 293: 686-693.
Miller KG, Kominz MA, Browning JV, Wright JD, Mountain GS, Katz ME, Sugarman PJ, Cramer BS, Christie-Blick N, Pekar SF: The Phanerozoic record of global sea-level change. Science. 2005, 310: 1293-1298. 10.1126/science.1116412.
Ruddimann WF: A Paleoclimatic Enigma?. Science. 2010, 328: 838-839. 10.1126/science.1188292.
van der Meulen AJ, Daams R: Evolution of Early-Middle Miocene rodent faunas in relation to long-term palaeoenvironmental changes. Palaeogeogr Palaeoclimatol Palaeoecol. 1992, 93: 227-253. 10.1016/0031-0182(92)90099-Q.
van Dam JA, Weltje GJ: Reconstruction of the Late Miocene climate of Spain using rodent palaeocommunity successions: an application of end-member modelling. Palaeogeogr Palaeoclimatol Palaeoecol. 1999, 151: 267-305. 10.1016/S0031-0182(99)00015-2.
Hernández Fernández M: Bioclimatic discriminant capacity of terrestrial mammal faunas. Global Ecol Biogeogr. 2001, 10: 189-204. 10.1046/j.1466-822x.2001.00218.x.
Hernández Fernández M: Rodent paleofaunas as indicators of climatic change in Europe during the last 125,000 years. Quaternary Res. 2006, 65: 308-323. 10.1016/j.yqres.2005.08.022.
Hernández Fernández M, Álvarez-Sierra MA, Peláez-Campomanes P: Bioclimatic analysis of rodent palaeofaunas reveals severe climatic changes in Southwestern Europe during the Plio-Pleistocene. Palaeogeogr Palaeoclimatol Palaeoecol. 2007, 251: 500-526. 10.1016/j.palaeo.2007.04.015.
Sesé C: Los roedores y lagomorfos del Neógeno de España. Estudios Geol. 2006, 62: 429-480.
Daams R, van der Meulen AJ, Peláez-Campomanes P, Álvarez-Sierra MA: Trends in rodent assemblages from the Aragonian (Early-Middle Miocene) of the Calatayud-Daroca Basin, Aragon, Spain. Hominoid evolution and climatic change in Europe Vol 1 The evolution of Neogene terrestrial ecosystems in Europe. Edited by: Agustí J, Rook L, Andrews P. 1999, Cambridge: Cambridge University Press, 127-139.
van Dam JA, Aziz HA, Álvarez-Sierra MA, Hilgen FJ, van den Hoek Ostende LW, Lourens LJ, Mein P, van der Meulen AJ, Peláez-Campomanes P: Long-period astronomical forcing of mammal turnover. Nature. 2006, 443: 687-691. 10.1038/nature05163.
Hernández Fernández M, Azanza B, Álvarez-Sierra MA: Iberian Plio-Pleistocene biochronology: micromammalian evidence for MNs and ELMAs calibration in southwestern Europe. J Quaternary Sci. 2004, 19: 605-616. 10.1002/jqs.862.
van de Weerd A, Daams R: Quantitative composition of rodent faunas in the Spanish Neogene and paleoecological implications. I & II. P K Ned Akad B Phys. 1978, 81: 448-473.
Walter H: Vegetationszonen und klima. 1970, Stuttgart: Eugen Ulmer
Daams R, van der Meulen AJ: Paleoecological interpretation of micromammal faunal successions in the Upper Oligocene and Miocene of north central Spain. Paléobiol Continent. 1984, 14: 241-257.
Hernández Fernández M, Peláez-Campomanes P: Ecomorphological characterization of Murinae and hypsodont “Cricetidae” (Rodentia) from the Iberian Plio-Pleistocene. Col Pa Vol. Extraordinario. 2003, 1: 237-251.
Hernández Fernández M, Vrba ES: Rapoport effect and biomic specialization in African mammals: revisiting the climatic variability hypothesis. J Biogeogr. 2005, 32: 903-918. 10.1111/j.1365-2699.2004.01188.x.
Friedman JH: Multivariate adaptive regression splines. Am Nat. 1991, 19: 1-67.
Wilson DE, Reeder DM: Mammal Species of the World. A Taxonomic and Geographic Reference. 2005, Baltimore: Johns Hopkins University Press
Steppan S, Adkins RM, Anderson J: Phylogeny and Divergence-Date Estimates of Rapid Radiations in Muroid Rodents Based on Multiple Nuclear Genes. Sys Biol. 2004, 53: 533-553. 10.1080/10635150490468701.
Fabre P-H, Hautier L, Dimitrov D, Douzery EJP: A glimpse on the pattern of rodent diversification: a phylogenetic approach. BMC Evolutionary Biology. 2012, 12: 88-10.1186/1471-2148-12-88.
Mein P: European Miocene Mammal Biochronology. The Miocene Land Mammals of Europe. Edited by: Heisig GRK. 1999, München: Verlag Dr. Friedrich Pfeil, 25-38.
Suc J-P, Zagwijn WH: Plio-Pleistocene correlations between the northwestern Mediterranean region and northwestern Europe according to recent biostratigraphic and palaeoclimatic data. Boreas. 1983, 12: 153-166.
Drivaliari A, Ticleanu N, Marinescu F, Marunteanu M, Suc J-P: A Pliocene climatic record at Ticleni (southwestern Romania). The Pliocene: Time of Change. Edited by: Wrenn JH, Suc J-P, Leroy SAG. 1999, Dallas: American Association of Stratigraphic Palynologists Foundation, 227-240.
Jaeger J-J: The evolution of biodiversity among the Southwest European Neogene rodent (Mammalia, Rodentia) communities: pattern and process of diversification and extinction. Palaeogeogr Palaeoclimatol Palaeoecol. 1994, 111: 305-336. 10.1016/0031-0182(94)90070-1.
Chaline J, Brunet-Lecomte P, Montuire S, Viriot L, Courant F: Anatomy of the arvicoline radiation (Rodentia): palaeogeographical, palaeoecological history and evolutionary data. Ann Zool Fenn. 1999, 36: 239-267.
Chaline J: Arvicolid data (Arvicolidae, Rodentia) and evolutionary concepts. Evolutionary Biology. 1987, 21: 237-310.
Martin RA, Siefker A, Marcolini F: Modelling the morphology and evolution of the linea sinuosa (crown-root junction) in arvicolid rodents; a test with Pliocene Ogmodontomys from Kansas, USA. Lethaia. 2009, 42: 155-166. 10.1111/j.1502-3931.2008.00139.x.
Potts R, Behrensmeyer AK: Late Cenozoic terrestrial ecosystems. Terrestrial Ecosystems trough Time. Edited by: Behrensmeyer AK, Damuth JD, DiMichele WA, Sues H-D, Potts R, Wing SL. 1992, Chicago: The University of Chicago Press, 419-541.
Mein P, Moissenet E, Adrover R: L’extension et l’âge des formations continentales Pliocènes du fossé de Teruel (Espagne). CR Acad Sci. 1983, 296: 1603-1610.
López-Martinez N, Michaux J, Hutterer R: The skull of Stephanomys and a review of Malpaisomys relationships (Rodentia: Muridae): Taxonomic incongruente in Murids. J Mamm Evol. 1988, 5: 185-215.
Renaud S, Michaux J, Schmidt D, Aguilar J, Mein P, Auffray J: Morphological evolution, ecological diversification and climate change in rodents. P Roy Soc B-Biol Sci. 2005, 272: 609-617. 10.1098/rspb.2004.2992.
Misonne X: African and Indo-Australian Muridae: Evolutionary Trends. Ann Mus Royal Afr Centr. 1969, 172: 1-219.
Fejfar O, Heinrich W-D, Kordos L, Maul LC: Microtoid cricetids and the early history of arvicolids (Mammalia, Rodentia). Palaeontol Electron. 2011, 14 (3 27A): 1-38.
García-Alix A, Minwer-Barakat R, Martín-Suárez E, Freudenthal M: Cricetidae and Gliridae (Rodentia, Mammalia) from the Miocene and Pliocene of southern Spain. Scripta Geologica. 2008, 136: 1-37.
Ayarzagüena J, López-Martínez N: Estudio filogenético y comparativo de Microtus cabrerae y Microtus brecciensis. Doñana Acta Vertebrata. 1976, 3: 181-204.
Garrido-García JA, Soriguer-Escofet RC: Cabrera’s Vole Microtus cabrerae Thomas, 1906 and the subgenus Iberomys during the Quaternary: Evolutionary implications and conservation. Geobios. 2012, 45: 437-444. 10.1016/j.geobios.2011.10.014.
Laplana C, Sevilla P: Documenting the biogeographic history of Microtus cabrerae through its fossil record. Mammal Rev. 2013, in press
van der Meulen A, Peláez-Campomanes P, Levin SA: Age structure, residents, and transients of Miocene rodent communities. Am Nat. 2005, 165: E108-E125. 10.1086/428683.
This work is dedicated to Pierre Mein (Université Claude Bernard Lyon 1); without his research on small mammal faunas and his contributions to the biochronology of European mammal faunas this kind of studies would not be possible. The Editor Davide Pisani (The Natural History Museum, London), Cesar Laplana (Museo Arqueológico Regional, Alcalá de Henares), Raef Minwer-Barakat (Institut Català de Paleontologia, Sabadell), Pablo Peláez-Campomanes (Museo Nacional de Ciencias Naturales, Madrid), Jan A. van Dam (Universiteit Utrecht) and several anonymous referees are acknowledged for their comments on preliminary versions of the manuscript. This is a contribution of the Team of Paleoclimatology, Macroecology and Macroevolution of Vertebrates (http://www.pmmv.com.es) of the Complutensian University of Madrid as a part of the Research Group UCM 910607 on Evolution of Cenozoic Mammals and Continental Paleoenvironments. We acknowledge financial support from the Spanish Ministry of Science and Innovation, through the projects CGL2006-01773/BTE and CGL2010-19116/BOS.
The authors declare that they have no competing interests.
ARGC JLC and MHF conceived, designed and performed the data analyses, wrote the manuscript, co-edited all drafts and prepared the final version of the manuscript. AM and AMB co-edited early drafts. AM and MHF conducted the gathering of data and developed the idea for the manuscript. MHF conceived, designed and coordinated the study, initiated the project, facilitated the gathering of contributors, refined the intellectual content, coordinated the authorship survey and is the guarantor for the integrity of the article as a whole. All authors read and approved the final manuscript.
Ana R Gómez Cano, Juan L Cantalapiedra and Manuel Hernández Fernández contributed equally to this work.