|Home | About | Journals | Submit | Contact Us | Français|
Although the recent historical period is usually treated as a temporal base-line for understanding patterns of mammal extinction, mammalian biodiversity loss has also taken place throughout the Late Quaternary. We explore the spatial, taxonomic and phylogenetic patterns of 241 mammal species extinctions known to have occurred during the Holocene up to the present day. To assess whether our understanding of mammalian threat processes has been affected by excluding these taxa, we incorporate extinct species data into analyses of the impact of body mass on extinction risk. We find that Holocene extinctions have been phylogenetically and spatially concentrated in specific taxa and geographical regions, which are often not congruent with those disproportionately at risk today. Large-bodied mammals have also been more extinction-prone in most geographical regions across the Holocene. Our data support the extinction filter hypothesis, whereby regional faunas from which susceptible species have already become extinct now appear less threatened; they may also suggest that different processes are responsible for driving past and present extinctions. We also find overall incompleteness and inter-regional biases in extinction data from the recent fossil record. Although direct use of fossil data in future projections of extinction risk is therefore not straightforward, insights into extinction processes from the Holocene record are still useful in understanding mammalian threat.
Analysis of global patterns of threat across major taxonomic groups has become an important tool for understanding the ecological patterns, dynamics and magnitude of the current-day biodiversity crisis, and constitutes a key step for prioritizing management actions and allocation of finite resources in systematic conservation planning . Comprehensive assessment of the conservation status of the 5487 wild mammal species recognized as extant since AD 1500 has shown that 25 per cent of all species for which adequate data are available are currently threatened with extinction, and that the distribution of these threatened species is geographically biased . A range of other studies have used large-scale datasets of life-history parameters and phylogenetic ‘supertrees’ for threatened and non-threatened mammals [3,4] to identify intrinsic and extrinsic correlates of extinction risk in different mammal groups [5–10], and to predict future changes in the geographical distribution of threatened species .
These studies all treat the recent 500-year historical period, an interval which broadly corresponds with the concerted dispersal of European explorers, traders and colonists around the globe, as the temporal base-line for understanding mammalian threat and extinction. However, although human impacts on global mammal diversity have been documented relatively well for this recent interval [12,13], mammalian biodiversity loss has also occurred throughout much of the Late Quaternary in association with older historical and prehistoric human migration, population expansion and environmental modification. The best-studied series of prehistoric Late Quaternary extinctions is the loss of at least 97 genera of megafaunal vertebrates (more than 44 kg; sensu ), mostly mammals but also some birds and reptiles, across the world's continents and some island systems during the Late Pleistocene Epoch without subsequent ecological replacement . A further extensive series of mammalian species extinctions, population extirpations and range contractions continued across insular and continental regions into the Holocene Epoch, the current geological interval, through recent prehistory and into the historical period .
The staggered timing of Late Pleistocene mammal extinctions and their close temporal association with first human colonization of different continental regions is increasingly well understood for North America [17,18] and Australia [19–22], and direct or indirect human involvement in these events is now widely accepted by most palaeontologists. Meaningful biological and life-history data are also available for these extinct regional mammal faunas [23–25]. However, the Late Pleistocene was also characterized by large-scale climatic and associated vegetational shifts during the transition from glacial to interglacial conditions, and the nature and extent of human involvement even in these relatively well-studied regional extinctions remains the subject of considerable debate, with non-anthropogenic environmental factors often advocated as primary extinction drivers on the basis of climatic and vegetation models and archaeological data [26–32]. On a regional level, the relationship between prehistoric human activity, climate change and Late Pleistocene mammal extinction in northern Eurasia is apparently more complex than previously thought [33,34]. Late Pleistocene mammalian extinction chronologies and faunal turnover are more poorly resolved for South America, Africa, the insular Mediterranean, eastern and southeast Asia and Australasia, constraining our ability to understand the extent of human involvement in these regional extinction events [16,35,36].
In contrast, mammal extinctions during the Holocene occurred in an interval of modest or minimal climatic variation, and under broadly ‘modern’ climatic and environmental boundary conditions . Human involvement in these events is much less controversial, and only an extremely small proportion of Holocene species or population losses (e.g. disappearance of the mammoth population on St Paul Island; ) can even questionably be interpreted as non-anthropogenic events . Although high-resolution geographical data on Holocene range change and population-level extirpation are still only known for a few well-studied mammal species (e.g. ), a global dataset of all mammal species currently known or suspected to have died out during the Holocene is now available . However, the relatively well-preserved recent fossil record is still substantially incomplete, and although new insights may be gained through the use of novel methodological approaches [40–42], constraining extinction chronologies and identifying causative extinction drivers remains a major challenge .
Despite these potential concerns, the fossil and archaeological records are increasingly recognized as important sources of data for evidence-based conservation, as they can provide a unique temporal perspective on past patterns of ecosystem diversity and function and how global systems respond to changing climatic and anthropogenic drivers [43–46]. However, past human-caused extinctions have primarily been studied by palaeontologists and archaeologists rather than in the context of current mammal extinction risk. Although some studies have investigated biological correlates of extinction-proneness in extinct Late Quaternary mammals, these have only incorporated data for well-studied regional components of the Late Pleistocene megafauna [47–49]. Existing studies that have investigated patterns and correlates of current extinction risk in mammals have typically excluded the great majority of the most susceptible species that have already become extinct as a result of human activities (e.g. [7,8]). Even from such a dataset of extant mammal species, it has been recently suggested that an ‘extinction filter’  has acted throughout the Holocene: large mammals are primarily threatened today in tropical areas, where agricultural impacts in the past have been relatively low, implying that large mammals in temperate regions have already been reduced to small but currently stable populations or driven extinct . Patterns of reduced modern-day species loss in areas that have the longest records of human occupation, and for which large-scale undocumented past extinctions can be assumed, have also been suggested for other taxonomic groups . The possible existence of such an extinction filter implies that the findings of studies based exclusively on extant species may be incomplete, and highlights the need for a more inclusive assessment of human impacts on global mammal faunas over time.
Here, we explore the spatial, taxonomic and phylogenetic patterns of all known Holocene mammal extinctions, to assess whether our understanding of mammalian threat processes has been affected by excluding this category of taxa. Throughout this study, the Holocene is defined as the interval from approximately 11 500 years ago up to the present day (i.e. also including the recent historical period). In order to investigate whether data from the past can provide useful new insights for understanding current mammal extinctions, we incorporate extinct Holocene mammal species into analyses resembling the current comparative literature on extinction risk correlates in extant mammals, adding 241 extinct species to the latest global studies. We also assess the proposed extinction filter effect by investigating whether these Holocene extinctions have been selective with regard to large body mass. Both analyses assume that the fossil record correctly reflects past extinction events; to assess the validity of this assumption, the last part of this paper investigates the extent to which the incompleteness of the recent fossil record constrains our ability to identify extinction drivers and use past extinction data for informing present-day conservation efforts. We therefore use our expanded dataset to investigate whether past and present mammal extinction risks display phylogenetically, ecologically and geographically congruent patterns.
Our dataset of all mammal species known or suspected to have become extinct during the Holocene is based on Turvey , and updated following recent taxonomic revisions and new species descriptions [53–56]. This provided a total of 241 mammal species which went extinct after the Pleistocene–Holocene transition (see the electronic supplementary material, table S1, for a full dataset). We were able to incorporate 210 of these species into an existing supertree [3,51], either using published taxon-specific phylogenies or on the basis of taxonomy  if phylogenies were unavailable, creating polytomies if necessary (see the electronic supplementary material, table S1). We did not run new dating analyses as for the last taxonomic update  because sequence data do not yet exist for most extinct species. We were mostly interested in topological placement of species to assess the phylogenetic signal of extinctions and to control for phylogeny in comparative analyses; we therefore used known divergence dates where available in our references (see the electronic supplementary material, table S1), and otherwise placed new nodes temporally halfway between existing ones. The new supertree contained 5210 species, 210 of which were extinct, and was 51 per cent resolved. To assess the impact of our ad hoc dating, we also ran all phylogenetic analyses on the basis of topology alone, with branch lengths set to equal .
Body mass estimates were obtained for 220 of the extinct species using four different methods. (i) A small number of body masses for extinct mammals were already available in the PanTHERIA database . (ii) Body masses for species with extant or extinct congeners of known body mass were calculated as the genus-level mean, largely using the PanTHERIA database. Body masses for species with no extant congeners were either (iii) collated from published estimates that were calculated using predictive regression equations based on skeletal measurement parameters, or (iv) were newly calculated for this study using published regression equations for different taxonomic groups [59–62] and published or newly measured skeletal morphometric data (see the electronic supplementary material, table S1).
Geographical distributions for extinct species were only possible to record at the country level, because a large proportion of these species are poorly known and only recorded from a very limited number of point localities. This precluded accurate reconstruction of their past geographical ranges, so range size was unavailable for extinct species and not used in our analyses. Our country-level occurrences represented only Holocene records; we excluded additional Pleistocene records that are available for many species because climate-led environmental changes between these two time intervals are known to have driven major geographical range shifts in many taxa [33,63,64]. Species from islands that consist of more than one modern geopolitical country (e.g. Hispaniola, New Guinea) were interpreted as having occurred in all of the countries on the island, because in most cases it remains difficult to test whether proposed regions of intra-island endemism in extinct mammals (e.g. ) represent genuine biogeographical patterns or artefacts of limited fossil sampling.
Country occurrence and IUCN (International Union for Conservation of Nature) Red List status for extant mammals were taken from the Global Mammal Assessment . We modified their species range maps according to the taxonomy of , and excluded historical, introduced and questionable ranges. We also excluded marine species (cetaceans, sirenians, pinnipeds, and the polar bear, sea otter and marine otter) from our spatial analyses. Our final dataset contained 5564 species, including 128 marine species; country occurrence was known for 5380 non-marine species. Of these species, 236 were extinct, two were extinct in the wild, and 768 were ranked as Data Deficient (DD) . Body mass data for 3801 extant mammal species were taken from the PanTHERIA database . A recent study showed that this dataset is biased in that more large-bodied species have known body mass values ; to avoid this bias, we used their method to impute body mass data for an additional 1143 extant species as the value of their closest relative on the basis of the near-complete mammalian phylogeny. These imputed data were only used in the randomization when testing for size selectivity of extinction, not in the extinction risk correlate analyses.
Last-occurrence dates for extinct mammals and dates for regional first human and European arrival were largely taken from Refs. [16,68,69] (see the electronic supplementary material, table S2). A new Holocene radiometric last-occurrence date for the extinct Hispaniolan sloth Neocnus dousman was calculated for this study as 7141 ± 35 yr BP full size from newly collected fossil material (Oxford Radiocarbon Accelerator Unit, OxA-19 725).
We used R v. 2.10.1 for all analyses . We initially explored patterns of Holocene extinctions and current extinction risk by mapping the number and proportion of extinct and threatened species within countries. Proportions of extinct species were defined as: (number of extinct species)/(number of extant + number of extinct species); proportions of threatened species as: (number of species in any of the three threatened categories (Vulnerable, Endangered, or Critically Endangered; ))/(number of extant species − number of species classified as DD by IUCN). Throughout our analyses, we classified the two mammal species listed as Extinct in the Wild (Elaphurus davidianus and Oryx dammah; ) as being extinct. Congruence of spatial patterns was assessed with Pearson correlation coefficients; we did not perform statistical tests for their significance as these would be biased by spatial autocorrelation, but the estimates of the coefficients themselves are thought to be unbiased .
To investigate taxonomic patterns in past extinctions and present risk, we calculated the number and proportion of extinct and threatened species for each mammalian order, and tested for equal proportions with Pearson's χ2. We also tested each order's proportion against its expectation given a binomial distribution, with the probability of extinction or risk set as the proportion of all mammals extinct or threatened. These analyses were run on a dataset of 4454 extant and 241 extinct species (including marine but excluding DD species). Phylogenetic signal in extinction and risk was assessed using the D statistic, which tests whether extinct or threatened species are significantly more clustered than expected under phylogenetically random extinction or threat (D < 1) . Additionally, D tests whether the observed phylogenetic signal is as strong as expected if the tested binary trait (e.g. extinct versus extant) was generated by ranking species according to values of a continuous trait evolved under a Brownian motion model. For example, if the largest species were always threatened, and body mass evolved according to a Brownian model, D for threat would not be different from zero (its expectation under a Brownian threshold model). Sample size for phylogenetic signal was the 5210 species in the supertree, 210 of which were extinct in the Holocene (again including marine species).
To assess the ‘extinction filter’ effect of extinct species on modelling body mass as a correlate of extinction risk, we fitted models within countries including and excluding extinct species. The response variable was the IUCN Red List categories converted to a numerical scale (Least Concern = 0, Near Threatened = 1, Vulnerable = 2, Endangered = 3, Critically Endangered = 4, Extinct in the Wild or Extinct = 5) ; log10-transformed body mass was the explanatory variable. The dataset therefore contained the 3410 non-marine species which were in the phylogeny, not ranked as DD, and for which non-imputed body mass estimates were available; 217 of these were extinct. Previous studies have excluded species ranked as threatened under criteria other than A (population decline; ) to avoid circularity arising from inclusion of geographical range size as an explanatory variable for species that have been ranked as threatened due to small range size [5,7,8]. This circularity did not arise for our analysis, as we were unable to use geographical range size as a predictor; therefore, we present analyses for all species.
Models were fitted only for countries with more than 10 species with all data. We used phylogenetic comparative analysis with the CAIC (Comparative Analysis using Independent Contrasts) package for R  to control for the fact that related species will be more similar in evolved traits such as body size, but also in extinction risk rating given its strong phylogenetic signal [58,74]. Assuming our branch lengths for extinct species are acceptable estimates, we ran phylogenetic generalized least-squares (PGLS) regression, which controls for variable amounts of phylogenetic covariation within the country datasets . In order to assess effects of branch lengths, we also performed the analysis with phylogenetically independent contrasts setting all branch lengths to equal .
We simulated random extinction within countries to investigate whether extinct species were unusually large. The number of extinct species within each country was drawn at random 1000 times from the country's pool of species with body mass values, including those extinct during the Holocene. We then assessed whether the observed median body mass of extinct species within each country was as expected under this random-extinction scenario within countries with a one-tailed significance test. The dataset for this procedure included the extant species for which body mass was imputed (see above), so the global dataset contained 5115 species, 217 of these extinct.
Finally, we examined the resolution of the recent fossil record and the robustness of our understanding of mammalian extinction chronologies and extinction drivers by analysing the temporal pattern of last-occurrence dates in relation to the earliest known dates of prehistoric human arrival on each island. This was conducted for the diverse extinct endemic mammal faunas from Madagascar, Cuba and Hispaniola, the three large islands that have experienced a disproportionately high number of known Holocene mammal extinctions .
Most recorded Holocene mammal species extinctions have occurred on islands, whereas a high proportion of currently threatened species are also found on continents (figure 1). Australia, Indonesia, Madagascar and the larger Caribbean islands have the highest overall number of extinct species (figure 1a), while only the Caribbean and some Pacific islands stand out for their proportion of extinct species (figure 1b). Most threatened species are concentrated in Indonesia, other regions of Asia (India, China, peninsular southeast Asia), and the Neotropics (Mexico, Brazil; figure 1c). Globally, the proportions of threatened species within countries are much higher than that of documented Holocene species extinctions, especially for peninsular southeast Asia, Indonesia, Madagascar, and some other tropical islands (figure 1d). Congruence between the numbers of extinct and threatened species across countries was relatively weak (Pearson correlation coefficient = 0.37).
Across our whole mammal dataset, 4.3 per cent (241/5564) of known species are extinct and 24.3 per cent (1083/4454) are threatened with extinction (excluding DD species), but proportions differ between orders (table 1). Past extinction and current threat are not distributed randomly across mammalian orders (proportion of extinct species: χ2 = 183.95, d.f. = 29, p < 0.001; proportion of threatened species: χ2 = 169.47, d.f. = 27, p < 0.001). Assuming the binomial distribution, artiodactyls, bandicoots and bilbies, bibymalagasians, proboscideans and xenarthrans have significantly more extinct species than expected, and bats have significantly fewer (table 1). Conversely, significantly more species of artiodactyls, perissodactyls, primates and sirenians are threatened, while bats, opossums and rodents have fewer threatened species than expected.
Past extinctions on the mammal supertree show significant phylogenetic signal (figure 2), with extinct species being significantly phylogenetically clumped whether branch lengths are used (D = 0.28, p < 0.001) or not (D = 0.35, p < 0.001). Using branch lengths, D is not quite significantly different from the expectation under the Brownian motion threshold model (p = 0.06), but when branch lengths are discarded the difference becomes significant (p < 0.01).
Overall, these results confirm the hypothesis that looking exclusively at current extinction risk patterns may distort our understanding of mammalian extinction: past extinctions throughout the Holocene have affected some areas and taxa disproportionately, and these are often not the areas and taxa that are disproportionately at risk today. The low spatial congruence of within-country past extinction and current extinction risk provides some support for the extinction filter hypothesis, whereby the faunas of countries in which many susceptible species have already become extinct now appear less threatened [50,51]. Similarly, some higher-order taxa containing high proportions of extinct species (e.g. bandicoots and bilbies, xenarthrans) now show average levels of current extinction risk.
However, these results could also be interpreted as indicating that different extinction processes are responsible for driving past versus present species losses. This alternative hypothesis may also be probable given the unprecedented and escalating level of modern-day anthropogenic habitat destruction . In birds, the effects of introduced species as the primary global extinction driver have recently been superseded in importance by habitat loss . Although prehistoric human-driven habitat modification is considered likely to be responsible for some Holocene mammal species losses in regions such as Madagascar , many past mammalian island extinctions have been linked to direct overexploitation by early colonists or the effects of invasive species, rather than to direct habitat clearance [16,80].
Different implications of our analyses provide some support for both of these hypotheses. Firstly, the taxonomic and phylogenetic clustering shown by both past extinctions and current extinction risk suggests that conserved biological traits make some taxa or clades consistently more susceptible to extinction than others . Although some higher-order taxa differ in their past and current risk, as described above, this overall common pattern may indicate similar extinction processes operating over time for certain mammal groups. For example, bats are consistently less threatened than expected from the mammalian average, while artiodactyls are consistently more threatened (see ). An alternative explanation for this observation is consistent knowledge bias, as we may know more about risk status in artiodactyls than bats; in fact, DD status among extant mammals shows a phylogenetic signal, indicating taxonomic bias in our knowledge of whether mammal species are threatened or not . Additionally, the fact that closely related species tend to live in the same geographical regions and therefore experience similar threats could inflate the global taxonomic and phylogenetic signal .
Secondly, while taxonomic clustering of Holocene extinctions has previously been demonstrated , as have taxonomic and phylogenetic clustering of currently threatened species [72,81,82,85], our analyses enable a direct comparison to be made between past and present extinction risk. Phylogenetic signal as measured by D is much stronger for extinct species (this study) than for all threatened ones (D = 0.64) . Our results using phylogenetic branch lengths show that phylogenetic clumping of extinctions is indistinguishable from that caused by extinguishing the top cohort of species for an evolved trait such as body mass; this provides further support for the extinction filter hypothesis.
Thirdly, the amount of phylogenetic signal in extinction risk has previously been shown to differ depending on the underlying threat process . Extant species threatened by overexploitation show a very strong phylogenetic signal, presumably due to strong body size selectivity, while species threatened by habitat loss show a much lower phylogenetic signal, as the impact of this threat process is less dependent on evolutionarily conserved biological traits. Species threatened by introduced species show an intermediate level of phylogenetic signal. Comparing our D values for extinct species to D values for extant species  shows that phylogenetic signal in extinct species is not as strong as for species threatened by overexploitation (D = 0.08), and is slightly stronger than for those threatened by introduced species (D = 0.46), but it is much stronger than for extant species threatened by habitat loss (D = 0.60). This comparison may suggest that Holocene mammal extinctions up to date have been largely driven by a combination of overexploitation and invasive species rather than habitat loss.
Including extinct species, body mass is significantly correlated with extinction risk in 79 mostly tropical countries, out of the 174 countries for which we have models (figure 3a). Large body mass is almost exclusively associated with high extinction risk (negative slope in only 17 countries); only Australia shows a significant negative relationship between body mass and risk. Results excluding extinct species differ slightly, particularly in areas with many extinct species (figure 3b), with slope estimates for several Caribbean countries and Madagascar decreasing by large amounts. Excluding extinct species makes the relationship between large size and high risk non-significant in 22 countries, most notably in Russia, southern Europe and the Middle East, for which the smaller numbers of documented terrestrial Holocene species extinctions have only involved megafaunal mammals (e.g. Bison priscus, Bos primigenius, Equus hydruntinus, Mammuthus primigenius, Megaloceros giganteus). Interestingly, in Hispaniola, Brazil, Saudi Arabia and South Africa, extinction risk excluding extinct species is significantly related to large body mass, but not when including the extinct species, which were mainly either mesoherbivores (e.g. Antidorcas bondi, Gazella saudiya) or comprised both small-bodied and large-bodied sympatric taxa (e.g. nesophontid island-shrews, heteropsomyine and capromyid rodents, and megalonychid sloths on Hispaniola).
Results using phylogenetically independent contrasts on the phylogeny with branch lengths set to equal give very similar results to PGLS (see the electronic supplementary material, figure S1), and the spatial congruence across slope estimates between the two methods is high (Pearson correlation coefficient = 0.81 including extinct species, 0.74 excluding extinct species). Only a few countries differ in whether their slope estimates are significant or not, and these tend to have low sample sizes. This high congruence supports the assumption that exact dating of phylogenetic splits for extinct species is not an important factor for our analyses, so we only discuss PGLS results below.
In most countries, Holocene extinctions have been selective with respect to body size (figure 4): median body size of extinct species is larger than expected under random extinction nearly everywhere. Notable exceptions are the USA, Australia, and some Mediterranean and Latin American countries; most of these include both continental regions and also offshore islands (e.g. Hawaii and the California Channel Islands; Corsica and Sardinia), with the islands having been the focus of disproportionately high numbers of extinctions of comparatively small-bodied insular species (typically rodents). On average, extinct species were largest in continental Europe and Asia, but even the smaller-bodied extinct species in insular regions such as the Caribbean, Indonesia and New Guinea constituted a biased sample in terms of body size for these countries.
Overall, our results support previous studies in that large body size is an important correlate of high extinction risk , and that its effect varies spatially [8,51]. Large extant species are currently at risk almost exclusively in the tropics, but including extinct species in our analyses provides further evidence that there may have been an extinction filter acting throughout the Holocene. This filter seems to have removed large species during the Holocene from many geographical regions, such as southern Europe, Russia, the Caribbean, Madagascar and Indonesia. Additional evidence for a Holocene body-size extinction filter is provided by our simulations of random extinctions within countries, which reveal size selectivity of past extinction across much of the globe. Whereas the increased vulnerability of large mammals to extinction during the Late Pleistocene has been well established [14,49], and the early loss of large-bodied species has been proposed for specific regional mammal faunas during the Holocene (e.g. the Caribbean; [16,69]), increased extinction-proneness of large mammals has not previously been demonstrated for global mammal faunas across the Holocene.
On the other hand, we found a consistently negative relationship between body mass and extinction risk in Australia. It has been suggested that Australian mammals of intermediate body mass (within a proposed ‘critical weight range’ of 35–5500 g) have an elevated risk of extinction due to region-specific interactions between mammal ecology and human threat factors [86,87]. Since we did not fit nonlinear terms, the many more available body mass values for intermediate and large species may drive an apparent negative relationship in Australia, while the true relationship would be hump-shaped. Furthermore, even our inclusion of extinct Holocene species reveals no relationship between body size and extinction risk in the USA, Mexico, several South American countries, and most of central and eastern Asia, which could also be taken as evidence against a general extinction filter hypothesis. However, the non-significant simulation results for countries comprising both continental regions and offshore islands may reflect a mixture of different past extinction processes, notably the higher impacts of invasive species on small island mammals versus anthropogenic overexploitation of larger mammals on continents . It is also probable that older extinction filters have significantly affected our analyses for continental regions such as the Americas, Australia and much of Asia, where many large susceptible mammals were already removed during the Late Pleistocene through glacial cycling or earlier human agency. It would be interesting to test this hypothesis by incorporating country-level data on Late Pleistocene mammal extinctions into similar analyses once more information becomes available on these faunas.
Although we are dependent upon the recent fossil record as the primary source of information on past mammalian species losses, extinction data derived from the fossil record for all three of the major island systems that have experienced particularly high levels of Holocene extinctions (Madagascar, Cuba and Hispaniola) are still demonstrably unresolved. All three islands have incomplete extinction chronologies for their endemic mammal faunas: known radiometric last-occurrence dates (either directly or indirectly dated) are only available for 70 per cent of the extinct species known from the recent fossil record on Madagascar (18/26), 45 per cent on Cuba (9/20) and 42 per cent on Hispaniola (10/24). In addition, a higher proportion of these last-occurrence dates fall completely within the period of human occupancy for Madagascar (16/18) and Hispaniola (9/10) compared with Cuba (2/9).
Investigation of the completeness of extinct species last-occurrence dates, and the degree to which these can be temporally correlated with first regional human arrival or other major human-caused environmental impacts, should indicate whether past mammal extinction data can be of use for identifying causative threat factors in modern faunas. The available radiometric dates demonstrate that representatives of all of the major mammal clades present in the Late Quaternary fossil record of Madagascar, Cuba and Hispaniola persisted through the Pleistocene and into the climatically stable Holocene, the interval when humans first reached each of these island systems (figure 5). It is therefore probable that most or all Late Quaternary mammal species losses on all three islands were driven by some form of prehistoric or historical-era human activities, rather than by non-anthropogenic environmental processes [16,69,88]. However, it is only possible to demonstrate evidence for a temporal overlap between extinct species and humans for 62 per cent of the extinct Late Quaternary mammal fauna of Madagascar, 38 per cent of Hispaniola and 10 per cent of Cuba. Extinction events may vary considerably in duration from rapid disappearance to gradual population attrition and decline in response to different drivers and ecologies , complicating any retrospective attempt to identify cause and effect in such events. For the remainder of the extinct mammal faunas of each island, in the absence of a meaningful timeline of last-occurrence dates we cannot therefore attempt to identify possible extinction drivers with any confidence, if humans were indeed responsible for these events.
Furthermore, the resolution of available data indicating temporal overlap between humans and extinct mammals differs significantly between the three islands, demonstrating inter-regional biases in fossil data quality as well as overall incompleteness of fossil data (χ2 = 7.796, d.f. = 2, p < 0.05). This difference probably represents varying levels of past research effort in different regions, rather than intrinsic preservational differences in the regional fossil record for these large islands; the higher resolution of last-occurrence dates available for Madagascar matches the greater amount of historical research that has been carried out on this island's Holocene mammal fauna, due in large part to long-term anthropological interest in extinct lemurs , in contrast to the relatively limited radiometric investigations conducted to date for extinct Caribbean mammal faunas. Similar biases are likely to extend to other regional Holocene faunas, in particular for the Indonesia-New Guinea region, for which only 38 per cent (11/29) of the mammal species interpreted as having become extinct during the Holocene have even been formally described (see the electronic supplementary material, table S1). These inter-regional patterns of incompleteness in our understanding of Holocene extinctions therefore present further constraints to identifying causative processes in past extinction events, and further restrict the usefulness of the recent fossil record for informing current-day conservation of phylogenetically or ecologically similar taxa. They also emphasize caution when interpreting the conclusions we have drawn in the last two sections, as our analyses of Holocene mammal extinctions are largely based on fossil data and are therefore probably highly incomplete for many regions.
Our results demonstrate that mammal extinctions across the Holocene have been phylogenetically and spatially concentrated in specific taxa and geographical areas, and are generally consistent with the hypothesis that large mammals have been disproportionately vulnerable to extinction since the end of the last glaciation. An extinction filter operating for much of the Late Quaternary has successively removed susceptible large-bodied species from extinction-prone regions. These have tended to be either island systems with naive mammal faunas, such as Madagascar, the Caribbean or southeast Asia, or continental regions with long Holocene histories of high human population density and agricultural activity. The size selectivity of extinction risk throughout the Holocene has led to current mammalian threat hotspots being located in tropical regions with more recent histories of anthropogenic intensification, such as southern Asia and the Neotropics. However, phylogenetic and spatial patterns of past extinction and current risk also support the idea that extinction processes have changed through time, with the focus shifting from islands to continents and from introduced species and overexploitation to widespread habitat loss and fragmentation. Finally, data on the quality of the Holocene fossil record emphasize that there are still substantial limitations to describing the dynamics of mammal extinctions even in the recent past.
Given the complexity of these results, can the past provide useful new insights for understanding current mammal extinctions and developing predictive models for the future? Overall, our results seem to caution against using Holocene extinction data at face value to make oversimplified predictions; although there are consistent similarities in the body-size selectivity of past extinction and current risk, it may be difficult to predict risk for particular countries and taxa, with past extinctions only able to inform us about current risk at a relatively broad scale. Further research is required to investigate whether mammal faunas that have survived historical extinction filters are now threatened by different anthropogenic threats, which may therefore alter the phylogenetic, ecological or geographical pattern of extinction-proneness traits. Data biases and incompleteness pose an additional barrier to the usefulness of past extinction data for understanding long-term patterns or trends in anthropogenic impacts on regional mammal faunas, especially for hotspots of mammalian diversity such as southeast Asia, which is currently experiencing extremely high levels of mammalian extinction risk . Although the direct use of fossil data in future projections of extinction risk is currently not straightforward, additional insights into extinction processes gained from this unique source of data in more locally and taxonomically restricted, high-resolution studies will still be useful in understanding current species extinction risk.
This study was funded by a NERC Postdoctoral Fellowship (NE/D009456/1) and a Royal Society University Research Fellowship (UF080320) to STT, a Marie Curie Fellowship from the European Commission to SAF (FP6 Early-Stage Training Network ‘Understanding and Conserving Earth's Biodiversity Hotspots’), and funding from the Danish National Research Foundation to the Centre for Macroecology, Evolution and Climate. We thank Jennifer Crees for data on Holocene distributions of extinct European mammal species, Katie Marske and anonymous reviewers for helpful comments on the manuscript, and William Turvey for technical assistance.