Below is the uncorrected machine-read text of this chapter, intended to provide our own search engines and external engines with highly rich, chapter-representative searchable text of each book. Because it is UNCORRECTED material, please consider the following text as a useful but insufficient proxy for the authoritative book pages.
14 Phylogenetic Trees and the Future of Mammalian Biodiversity T. JONATHAN DAVIES,* SUSANNE A. FRITZ,â RICHARD GRENYER,â¡ C. DAVID L. ORME,â JON BIELBY,â Â§ OLAF R. P. BININDA-EMONDS,î© MARCEL CARDILLO,â # KATE E. JONES,Â§ JOHN L. GITTLEMAN,** GEORGINA M. MACE,â¡ and ANDY PURVISâ Phylogenies describe the origins and history of species. However, they can also help to predict speciesâ fates and so can be useful tools for managing the future of biodiversity. This chapter starts by sketching how phylogenetic, geographic, and trait information can be combined to elucidate present mammalian diversity pat- terns and how they arose. Recent diversification rates and stand- ing diversity show different geographic patterns, indicating that cradles of diversity have moved over time. Patterns in extinction risk reflect both biological differences among mammalian lineages and differences in threat intensity among regions. Phylogenetic comparative analyses indicate that for small-bodied mammals, extinction risk is governed mostly by where the species live and the intensity of the threats, whereas for large-bodied mammals, ecological differences also play an important role. This modeling approach identifies species whose intrinsic biology renders them particularly vulnerable to increased human pressure. We outline *National Center for Ecological Analysis and Synthesis, University of California, 735 State Street, Santa Barbara, CA 93101; â Division of Biology and â¡Natural Environment Research Council Centre for Population Biology, Department of Life Sciences, Imperial College Lon- don, Silwood Park Campus, Ascot SL5 7PY, United Kingdom; Â§Institute of Zoology, Zoologi- cal Society of London, Regentâs Park, London NW1 4RY, United Kingdom; î©AG Systematik und Evolutionsbiologie, Institut fÃ¼r Biologie und Umweltwissenschaften, Carl von Ossietz- ky University, 26111 Oldenburg, Germany; #Centre for Macroevolution and Macroecology, School of Botany and Zoology, Australian National University, Canberra 0200, Australia; and **Odum School of Ecology, University of Georgia, Athens, GA 30602. 263
264â /â T. Jonathan Davies et al. how the approach might be extended to consider future trends in anthropogenic drivers, to identify likely future battlegrounds of mammalian conservation, and the likely casualties. This frame- work could help to highlight consequences of choosing among different future climatic and socioeconomic scenarios. We end by discussing priority-setting, showing how alternative currencies for diversity can suggest very different priorities. We argue that aiming to maximize long-term evolutionary responses is inap- propriate, that conservation planning needs to consider costs as well as benefits, and that proactive conservation of largely intact systems should be part of a balanced strategy. T he Tree of Lifeâphylogenyâis a powerful metaphor for lifeâs diversity, showing all species, including our own, as part of an interrelated whole. But phylogeny is more than a metaphor. It is also a research toolâthe result of evolutionary processes integrated over the history of life, it can be analyzed for insights into how those processes have shaped todayâs biota (Harvey et al., 1996). This approach is becom- ing increasingly powerful as the trees become ever more inclusive, built from rapidly accumulating databases using methods that continue to be improved (Bininda-Emonds, 2004; Cracraft and Donoghue, 2004). Speciesâ histories are of interest, but their futures are of more pressing concern. The Tree of Life is currently under sustained attack, as people increasingly dominate landscapes (Millennium Ecosystem Assessment, 2005a). Comparisons of extinction rates between today and geological history are difficult for many reasons (Jablonski, 1994), but the Tree of Life is already being pruned more quickly than it is growing (Purvis and Hector, 2000), and extinction rates are projected to rise by at least another order of magnitude over the next centuries (Mace et al., 2005). This chapter describes how phylogeny has a role to play in understanding the pattern of survivors and casualties and how it can help us both to predict speciesâ futures and to estimate some of the biodiversity value that would be lost if they went extinct. We focus on mammals as a model system. They are much better known than almost any other group, with a mature species-level tax- onomy (Wilson and Reeder, 2005), a largely complete evaluation of spe- ciesâ extinction risk (Baillie et al., 2004), a large database of ecological, life history, and geographical information (K.E.J., J.B., A.P., C.D.L.O., S.A.F., Christina Connolly, Amber Teacher, J.L.G., R.G., Elizabeth Boakes, Michael Habib, Janna Rist, Chris Carbone, Christopher A. Plaster, O.R.P.B.-E., Janine K. Foster, Elisabeth A. Rigby, Michael J. Cutts, Samantha A. Price, Wes Sechrest, Justin OâDell, Kamran Safi, M.C., and G.M.M., unpub- lished data), and a comprehensive species-level estimate of phylogeny
Phylogenic Trees and the Future of Mammalian Biodiversityâ /â 265 (Bininda-Emonds et al., 2007). [Many node ages in Bininda-Emonds et al. (2007) were slightly affected by a software bug; all our analyses use corrected dates (Bininda-Emonds et al., 2008).] As always with a model system, advantages come at a price. Mammals are atypical (e.g., they are much larger than most other species), so results from them cannot neces- sarily be extrapolated more broadly. However, mammals are a charismatic group of special interest to many people, so results have value even if they cannot be generalized. We start with a snapshot of present mammalian diversity and the (overwhelmingly anthropogenic) pressures that species face, before going on to describe recent and ongoing attempts to under- stand the present and possible future consequences of those pressures. A SNAPSHOT OF MAMMALIAN BIODIVERSITY Mammalian species are distributed very unevenly among genera, families, and orders (Purvis and Hector, 2000; Wilson and Reeder, 2005). Differences in age among taxa of a given rank (Avise, Chapter 15, this volume) confound evolutionary interpretation of the pattern, but the phy- logeny permits a test of whether the chances of diversification have indeed varied among lineages. Under the equal-rates Markov model (ERM), in which chances are equal, phylogenies should have a weighted mean I [the degree to which species are partitioned unequally between sister clades (Purvis et al., 2002)] not significantly >0.5. The estimate of phylogeny (Bininda-Emonds et al., 2007) has a weighted mean I of 0.657 (SE = Â±0.0131), well above 0.5 (weighted t test vs. 0.5: t848 = 11.98, p << 0.001), indicating that lineages have had different propensities to diversify. Such inequality is common throughout the Tree of Life (Mooers and Heard, 1997; Purvis, 1996) and prompts the search for traits that might be responsible. Phy- logenetic analyses reveal that large litter size and high abundance are both linked with high richness in sister-clade comparisons pooled across four orders (primates, carnivores, marsupials, and bats), whereas small body size and short gestation period also predict high richness within carnivores (Isaac et al., 2005). In common with most comparable studies on other taxa (Coyne and Orr, 2004), however, the biological traits leave most of the variance in richness unexplained, suggesting a possible role for the environment. The geographic distribution of mammalian species is also very uneven (Fig. 14.1). Mammals follow global trends for higher tropical diversity, with a strong latitudinal diversity gradient (Fig. 14.1a). Within the tropics, richness seems to correlate with productivity and waterâenergy dynamics, peaking in Amazonia at the base of the Andes, in the Great Rift Valley in Africa (where richness exceeds 250 species per 10,000 km2), and in an arc running from the Himalayas into southeastern Asia (Fig. 14.1b). These
266â /â T. Jonathan Davies et al. FIGURE 1 Species Richness log10 Range Size km2 FIGURE 14.1â Geographic patterns in mammalian biodiversity. (a) The latitudinal gradient in species richness. The solid line shows the median; gray bands demar- cate the interquartile range; species are counted in every latitudinal band (100 km north to south) in which they occur. (b) A map of species richness within 100 Ã 100-km grid cells, ranging from light (minimum = 0 species) to dark (maximum = 258 species). (c) The latitudinal pattern of mean geographic range size; details as for Fig. 1a. peaks suggest that topographical heterogeneity has also shaped species- richness patterns (Kerr and Packer, 1997). At higher latitudes, the relation- ship between richness and productivity is much weaker, and ambient energy is a better predictor of richness (Kerr and Packer, 1997). The geo- graphic pattern of mammalian richness is highly congruent with those of birds and amphibians (Grenyer et al., 2006), which are also well predicted by the same environmental variables (Hawkins et al., 2003; Buckley and Jetz, 2007). All these groups have many narrow-range species on the Neo- tropical mainland, but other ââhotspots of endemicityââ differâcontinental shelf island systems for mammals, oceanic archipelagoes for birds, and mainland locations for amphibians (Grenyer et al., 2006). The latitudinal gradient in median geographic range size in mammals (Fig. 14.1c) cor- relates strongly with the land area available within the same latitudinal bands (r = 0.68, although spatial autocorrelation complicates significance testing). However, median range size remains high in the northernmost bands despite the rapid decline in land area toward the Pole. Bird range sizes show similar patterns (Orme et al., 2006). The tropics have been described as an evolutionary powerhouse, act- ing as a ââcradleââ for diversity (Stebbins, 1974). Might Fig. 14.1, therefore, reflect longstanding geographic differences in diversification rates? Cra- dles of diversity could be of particular interest if conservation actions are to be targeted toward conserving the processes that generate diversity, a point to which we return below. The geographic pattern of recent diversifi- cation can be examined in terms of taxonomy or phylogeny, although both have their problems. Boundaries 38 higher taxa can be arbitrary (Avise, to
Phylogenic Trees and the Future of Mammalian Biodiversityâ /â 267 Chapter 15, this volume), whereas both the relative and absolute diver- gence times in molecular phylogenies can be controversial (PulquÃ©rio and Nichols, 2007). Because of these reservations, we use both approaches. Taxonomically, the ratio of species to genera should indicate the regional diversification rate over the past several million years, if gen- era are approximately equal in age. This ratio correlates strongly with log(species richness) among equal-area grid cells (Pearsonâs r = 0.61; Fig. 14.2a). This correlation remains highly significant (corrected F = 78.49, corrected P < 0.001, n = 4,152: based on a subsample of cells and exclud- ing single species occurrences) when degrees of freedom are reduced to account for spatial autocorrelation (Clifford et al., 1989). Moving to phy- logeny, places where a high proportion of species are on short terminal branches in the tree are likely to have rapid diversification, turnover, or immigration in their recent history (Roy and Goldberg, 2007). However, analyses are complicated by the low resolution (uncertain relationships) at the tips of the phylogeny, which introduces overestimates of the respective branch lengths. We ameliorated this problem by reducing ages of terminal polytomies using the correction suggested by Nee in Mooers and Atkins (2003) and by assuming that the descendants from each polytomy diver- sified under a Yule process (Nee, 2006). Fig. 14.2b highlights the Andean and Himalayan diversity peaks, but not the African great lakes, as recent evolutionary crucibles. Much of the temperate north stands out more than much of the tropics in this map, and there is a negative overall correlation between the proportion of short branches and log(species richness) (r = -0.38, corrected F = 21.01, corrected P < 0.001, n = 4,210, analyses as above), although this depended on how we corrected for terminal polytomies. This result partially echoes recent findings of higher recent speciation and extinction rates in temperate than in tropical mammals (Weir and Schluter, 2007). These maps also imply that some regions have seen marked shifts in net diversification rate, whereas others may have remained steady. Areas of ââoldââ and ââyoungââ diversity can be identified from the resid- uals of a loess regression across cells of total evolutionary history (i.e., total branch length in the phylogeny of a cellâs species) on species number. The African diversity peak emerges as old, whereas much of Andean diversity is young (Fig. 14.2c). Character disparityâamong-species variation in morphologyâalso shows geographic pattern: Fig. 14.2d maps one index of disparity, the variance in log(body mass). Disparity tends to be high where diversity is old (r = 0.29, corrected F = 5.10, corrected P < 0.05, n = 4,210, analyses as above), although tropical regions drive this relationship. Mammalian biodiversity, then, shows complex geographic and phy- logenetic patterns of richness, recent diversification, and character varia- tion. The African diversity peak is old and disparate, that in Asia is young and disparate, and the Andean peak is young with low disparity. These
FIGURE 2 268â /â T. Jonathan Davies et al. High : 2 a Low : 1 High : 1 b Low : 0 High : 708 c Low : -1101 High : 6.82 d Low : 0 FIGURE 14.2â Maps showing four aspects of mammalian diversification and di- ` versity. (a) Mean numbers of species per genus within cells. (b) Proportion of species with shorter terminal branch lengths than the overall median [phylogeny from Bininda-Emonds et al. (2007), with branch lengths corrected and modified as Davies et al. figure 2 (images are separate) explained in text]. (c) Residuals from a loess regression of cell total evolutionary history against cell species number (phylogeny as above). (d) Variance in log(body mass) (unpublished data), with missing data inferred as the mean of all closest relatives with data.
Phylogenic Trees and the Future of Mammalian Biodiversityâ /â 269 patterns reflect a complex history of speciation, extinction, anagenesis, and dispersal, with each factor probably shaped by biological traits and both biotic and abiotic environmental features in ways that have changed through time. Environmental features could provide much of the missing explanation for phylogenetic asymmetry (Ricklefs, 2003), whereas traits may help to complete the explanation of geographic patterns. Other, less well known taxa doubtless have different patterns, but there is no particu- lar reason to expect the patterns to be any simpler. Natural diversity patterns increasingly bear the stamp of widespread anthropogenic system change. The biota we see today is already affected by anthropogenic extinction. Large, slowly reproducing mammals went extinct almost everywhere (with the notable exception of Africa, perhaps linking to the patterns in Fig. 14.2c and d) â7,000 to 50,000 years ago (Barnosky, Chapter 12, this volume). In the past few centuries, mammalian extinctions have mostly been on islands, notably the West Indies, with continental extinctions largely confined to Australia (Baillie et al., 2004). Our analyses are therefore of an already reduced fauna. The next section considers the main ways in which human actions continue to reduce and reshape mammalian biodiversity. THREATS FACING MAMMALIAN BIODIVERSITY The terrestrial environment is now dominated by peopleâ1/4 to 1/3 of the land area has been transformed for human use (Vitousek et al., 1997). Additionally, human population density tends to be higher in species-rich areas, probably because productivity shapes both (Luck, 2007). Only a few mammal species fare well in human-dominated environments; the vast majority are vulnerable to the widespread and rapid anthropogenic changes. The main direct human-induced drivers that impact biodiver- sity now are habitat loss and fragmentation (the most important present threat), alien invasive species, overutilization, disease, pollution, and cli- mate change (Baillie et al., 2004). The International Union for Conserva- tion of Nature (IUCN) assessed how these drivers are affecting mammals. Nearly all mammal species have been evaluated and, provided enough information was available, placed in one of the following extinction risk categories: least concern (LC), near-threatened (NT), vulnerable (VU), endangered (EN), critically endangered (CR), extinct in the wild (EW), and extinct (EX). The resulting IUCN Red List (International Union for Conservation of Nature, 2007) lists 74 mammal species as having gone extinct in their native range since A.D. 1500, 1,094 (22.5% of those evalu- ated) as being threatened (i.e., VU, EN, or CR), and only 2,652 (54.5%) as giving no cause for concern (i.e., as being LC).
270â /â T. Jonathan Davies et al. Does extinction risk show any patterns that might help us to under- stand the processes affecting biodiversity loss? The ââfield of bulletsââ sce- nario, in which extinction strikes completely at random, is a widely used null model for extinction (Raup et al., 1973; Van Valen, 1976). The metaphor comes from trench warfare, where soldiersâ survival may have depended more on luck than skill. This scenario predicts that threatened species should constitute a stochastically constant fraction of any sample. Mam- malian extinction risk is not a simple field of bullets but shows both geo- graphic and phylogenetic patterns (Russell et al., 1998; Mace and Balmford, 2000; Baillie et al., 2004). The prevalence of risk is higher in the Old World than in the New World and higher on islands than on continents (Mace and Balmford, 2000; Baillie et al., 2004). It varies among clades too, being higher than average in primates and perissodactyls and lower than aver- age in rodents (Baillie et al., 2004). Species with few close relatives are also more likely to be at risk (Russell et al., 1998; Purvis et al., 2000a). These patterns reject the original field-of-bullets model, but the model lacks a geographic dimension because there may have been fields that were near to the battle but that nevertheless had no bullets. Likewise, human pressures have changed some places beyond recognition but left others almost untouched. Because closely related species often live in the same region, geographic heterogeneity in threat intensity could, by itself, cause taxonomic selectivity in extinction risk (Russell et al., 1998). Alter- natively, the selectivity could arise because biological differences among clades affect speciesâ abilities to withstand threats (McKinney, 1997). How important for speciesâ survival is staying out of the firing line, and how important is being bulletproof? Human population density predicts proportions of threatened mam- mal species among continental countries (McKinney, 2001), supporting the ââfiring-lineââ model. However, an analysis of extinction risk within terrestrial World Wildlife Fund (WWF) ecoregions shows that phyloge- netic nonrandomness is common within single ecoregions, where pres- sures may be more even than across the globe. Within each ecoregion, we generated phylogenetically independent contrasts (essentially, differences between sister clades) (Harvey and Pagel, 1991) in extinction risk (0 for LC species, 1 for species having a higher status; data-deficient and unevalu- ated species were excluded) on the phylogeny (with polytomies resolved arbitrarily and branch lengths set to unity) and compared the sum of the absolute values of standardized contrasts with the sums obtained from 1,000 randomizations of the risk data among the ecoregionâs species. If high-risk species are strongly clumped in the phylogeny, the observed sum will be lower than in 95% of the simulations. Of 691 ecoregions with at least three higher-risk species and some variance in extinction risk, 386 (54%) showed significant clumping. Interestingly, the strength (rather
FIGURE 3 Phylogenic Trees and the Future of Mammalian Biodiversityâ /â 271 FIGURE 14.3â Strength and significance of clumping in extinction risk within WWF ecoregions. Scale bar below the map indicates clumping strength. A value of 1 indicates randomness, and clumping is stronger for lower values. Circle size indicates the p value [radius is proportional to -ln(p)]; circle size for P = 0.05 is shown at the lower left. than the significance) of the clumping is high in most realms apart from the Nearctic (Fig. 14.3). It also appears to be stronger in ecoregions with high disparity (Spearmanâs Ï = 0.316) and with relatively old diversity (Spearmanâs Ï = 0.195). These correlations should not be taken as evidence of a functional syndrome unless confirmed at more local scales: Some of the signal probably derives from differences among, rather than within, Davies et al. Figure 3 major biogeographic realms. The prevalence of clumping of risk implies that, faced with approximately equal pressures, species differ in their ability to persist because of lineage-specific characteristics. This finding invites a search for biological correlates of extinction risk. COMPARATIVE ANALYSES OF MAMMALIAN EXTINCTION RISK Perhaps the most obvious proposed risk factor for extinction is large body size. The end-Pleistocene mass extinction of mammals removed mostly large species (Barnosky, Chapter 12, this volume), and declining mammals are an order of magnitude heavier, on average, than are non- threatened species (Cardillo et al., 2005). There are several possible rea- sons: Large-bodied species are more tempting targets than small ones for hunters; they are, on average, less abundant; and they take longer to reach sexual maturity, have smaller litters of larger offspring, and have larger individual home ranges. Narrow ecological tolerances are also a plausible risk factorâhabitat specialists may be more at risk than generalists from
272â /â T. Jonathan Davies et al. habitat loss (Fisher and Owens, 2004). A small geographic range size may reflect narrow tolerances and increase the risk that the whole of the spe- ciesâ range is in the firing line. Which of these features matter most for extinction risk, and are any associations consistent across all mammals? Cardillo et al. (2005) carried out the most comprehensive investigation to date. Threatened species were included only if they were on the IUCN Red List because of observed decline, to avoid autocorrelation with predic- tor variables. Red List status, on a 0â5 scale, was used as the response vari- able. Many facets of geography (including human population density), ecology, and life history were tested as predictors of extinction risk, by using phylogenetically independent contrasts. A phylogenetic approach is needed because, although extinction risk and some of the possible pre- dictors listed above (e.g., geographic range size) do not evolve along the phylogenyâs branches like, say, body size does, they nonetheless tend to show phylogenetic signal [i.e., they tend to take more similar values in close relatives than in species chosen at random; (McKinney, 1997; Purvis et al., 2000b; Fisher and Owens, 2004)]. Minimum adequate models were derived from a large initial set of predictors. This approach helps exclude variables that correlate only indirectly with extinction risk, for example, because another variable shapes both them and risk. The predictors of risk were significantly different for smaller and larger species, with the importance of many predictors changing markedly at a body size of â3 kg. Species smaller than this fit the firing-line model: They are more likely to be threatened if they have small geographic ranges, live in temperate areas, face high human population densities, and live where a high proportion of the other mammal species are also threatened. Larger species, however, face multiple jeopardy: Biology matters as well as geography, with high-abundance, small neonates, and many litters per year all independently helping to bullet-proof species. High abundance is predicted to bullet-proof species if the field-of-bullets model operates at the level of individuals rather than species (Erwin, 2006a), but such a model also predicts that no other biological traits would independently predict risk. For both large and small mammals, the most important single risk factor is small geographic range size. The firing-line model predicts that small-ranged species will be most at risk because a single localized threat can impact their entire distribution. However, range size itself varies sys- tematically among clades [although it shows weaker phylogenetic signal than, e.g., body size (Gaston, 2003; KE Jones et al., 2005)], suggesting that it is shaped, at least in part, by organismal traits such as dispersal ability (BÃ¶hning-Gaese et al., 2006) or niche breadth as well as by circumstances of geography. For example, small-ranged species are more common at low latitudes and within climatically stable regions. Any traits, including
Phylogenic Trees and the Future of Mammalian Biodiversityâ /â 273 geographic location, that confer large species ranges also help make spe- cies bulletproof (although geographical variation in speciesâ range sizes again complicates separation of geographical variability in threat intensity from intrinsic biological vulnerability). Large-scale analyses can find general predictors of extinction risk but can miss interesting variation among regions or clades, which more narrowly focused models might pick up (Fisher and Owens, 2004). Order- specific models typically have higher explanatory power than the large- scale models. These models have some common features, such as the importance of geographic range size, but also differ considerably (Cardillo et al., 2008). For example, body size is a predictor in bats but not in rodents, whereas different life history traits predict risk in carnivores (gestation length) and ungulates (weaning age, interlitter interval). Likewise, differ- ent environmental factors and measures of human impact are implicated in different taxa. The models also vary regionally, with life history matter- ing less in north temperate regions than elsewhere (Cardillo et al., 2008). One likely source of variation in models is that different drivers may select against different characteristics and show spatial variation in inten- sity. Broad-scale analyses may therefore lump competing signals together (Owens and Bennett, 2000). Within artiodactyls, predictors of extinction risk differ between hunted and nonhunted species: Late weaning age was the sole risk factor among the former, whereas low income levels among local people and small range size predicted risk among the latter (Price and Gittleman, 2007). More generally, low reproductive rates and large size are likely risk factors for overexploitation, but a specialized habitat may matter more under habitat loss (Owens and Bennett, 2000). Analyses focused more tightly on driver-specific responses often tend to consider far fewer species, in which case far fewer predictor variables can be con- sidered simultaneously without overfitting, and statistical power may be lower. On the plus side, the tighter focus can reduce the chance of mixed signals [although interaction terms can also do this (Price and Gittleman, 2007)], and more precise measures of driver intensity and extinction risk might be available than can be had globally (Fisher et al., 2003; Isaac and Cowlishaw, 2004). Broad- and narrow-scale analyses each give part of an obviously very complex picture. Furthermore, we have focused on phy- logenetic nonindependence, but to fully consider the interaction between biology and geography, the development of methods that also deal with spatial nonindependence in comparative data will be critical. Analyses modeling risk as a function of intrinsic biology (i.e., not including driver information) can highlight species at lower risk than expected from their geography, ecology, and life history (Cardillo et al., 2006). Such species may be particularly likely to decline rapidly if drivers intensify, because their attributes are repeatedly found in rapidly declining
274â /â T. Jonathan Davies et al. taxa. Cardillo et al. (2006) termed this ââlatent riskââ and proposed 20 regions with largely intact, but intrinsically susceptible, mammalian faunas. These include the Nearctic boreal forests and the island arc of Southeast Asia, and are mostly not exceptionally high in numbers of total, endemic, or threat- ened species. Many have much less than 10% of their land within reserves, and some (especially in Southeast Asia) face very rapid human popula- tion growth. As such, latent-risk hotspots might represent cost-effective options for long-term conservation. However, these analyses do not yet consider realistic scenarios of future driver patterns; rather, they implicitly assume that places with low intensity will experience an increase to more typical levels (Cardillo et al., 2006). The next section discusses how more policy-relevant predictions could be obtained by projecting future driver patterns based on explicit scenarios. MODELING FUTURE DECLINES Predicting future declines is more complex than explaining present declines, because the future is not just a linear extrapolation from the past and present. Past extinctions were largely caused by invasive species and overexploitation; habitat alteration is now a more important driver (Baillie et al., 2004). Changes in land use have been mapped historically (www.mnp.nl/hyde) and are tracked in the present day (http://glcf. umiacs.umd.edu/data), but analogous spatial data for other main drivers are more problematic. Wild species might be most vulnerable to overex- ploitation where people live at high density and have few other protein sources, suggesting that predictive models can be developed at regional scales (Fa et al., 2003; Ling and Milner-Gulland, 2006). The patterns and driving processes behind invasive species have varied over time (Mack and Lonsdale, 2001), and, although there are clear associations with global movements from human migration and trade, identifying clear predictive methods for the intensity of invasives is a work in progress (Hastings et al., 2005). The same is true for disease (Pedersen et al., 2007). Given the difficulties in obtaining spatial data on the present intensity of direct drivers, let alone future projections, an alternative is to work with information on indirect driversâin particular, human population density and growth and patterns of land conversion. Projections of these drivers are available under a range of socioeconomic scenarios (Millen- nium Ecosystem Assessment, 2005b). Intensity data alone are not enough, however; the response curves linking intensity to decline are also needed, and responses will depend on how bulletproof the biota is. Thus, declines need to be modeled as a function of both driver intensity and relevant biological attributes. A first step (Cardillo et al., 2004) considered a single driver (human population density) under a single growth scenario, cou-
Phylogenic Trees and the Future of Mammalian Biodiversityâ /â 275 pling explanatory models of carnivore extinction risk from comparative analyses with human population projections to identify species whose conservation status was likely to worsen. Here, we enlarge this approach in a preliminary analysis of two driv- ers and all mammals. Across ecoregions, the proportion of species with risk status higher than LC was modeled (as a binomial denominator) as a function of two drivers and two summaries of biological vulnerability by using generalized additive models (Wood, 2006) to avoid forcing any particular form on the relationship. A smooth relationship was fitted to link risk level to mean human population density (Center for Inter- national Earth Science Information Network and Centro Internacional de Agricultura Tropical, 2005) and the proportion of land converted to urban or cropland (European Commission, Joint Research Centre, 2003; Center for International Earth Science Information Network et al., 2004). A second smooth relationship was fitted to control for two biotic variables [proportion of species weighing >3 kg, the size at which ecology and life history begin to influence risk strongly (Cardillo et al., 2005), and the proportion of species in the lowest quartile of global range size (K.E.J., J.B., A.P., C.D.L.O., Susanne A. Fritz, Christina Connolly, Amber Teacher, J.L.G., R.G., Elizabeth Boakes, Michael Habib, Janna Rist, Chris Carbone, Christopher A. Plaster, O.R.P.B.-E., Janine K. Foster, Elisabeth A. Rigby, Michael J. Cutts, Samantha A. Price, Wes Sechrest, Justin OâDell, Kamran Safi, M.C., and G.M.M., unpublished data)]. Fig. 14.4 shows the marginal effect of the drivers on extinction risk. The two drivers are strongly corre- lated across ecoregions and interact strongly. As expected, risk is low when both drivers are at the very lowest levels. However, risk rises rapidly as either driver increases. Medium levels of land conversion and low density are associated with high levels of risk, but risk falls as land conversion rises further. This suggests that land conversion is an extinction filter (Balmford, 1996), removing one set of species sufficiently thoroughly that highly converted regions can again have low levels of risk, with only the more bulletproof taxa remaining. Scenario analysis will obviously need to count projected extinctions as well as declines and may need to consider historical as well as present driver patterns. As human density reaches high levels, risk levels become uniformly higher. A more refined model, perhaps incorporating other drivers, could be combined with projected future driver intensity to predict where a high proportion of species will decline. Such an approach uses the spa- tial heterogeneity in present driver intensity as a surrogate time series, with high-intensity ecoregions suggesting what will happen elsewhere as conditions deteriorate. However, incorporating climate change into this modeling approach presents major challenges. Because it has not been a major driver of present risk patterns, we have not yet seen how intensity
276â /â T. Jonathan Davies et al. FIGURE 14.4â The surface relating the prevalence of extinction risk to urban and agricultural land use and human population density (HPD), controlling for two indices of biological susceptibility, among ecoregions. White regions in the top left and bottom right corners contain no ecoregions. Light = low risk; dark = high will relate to impact and cannot use spatial heterogeneity as a guide. Bio- climatic envelope models suggest that climate change is already affecting many species including mammals (Parmesan and Yohe, 2003) and may soon be the dominant driver, possibly exacerbated by interactions with invasive species and other threats (Thuiller, 2007). Climate change is likely to particularly impact species that face geographical or biological barriers to dispersal or that depend on environmental cues that may break down Davies et al. figure 4 as climate changes (Bradshaw and Holzapfel, 2006). SETTING PRIORITIES AMONG AREAS Conservation spending is not nearly enough to maintain even the current inadequate network of reserves and protected areas (James et al., 2001). The identification and prioritization of global conservation net- works are a major focus of conservation dollars (Halpern et al., 2006). However, rankings inherently depend on the currency used to evaluate regions. A simple demonstration can be provided by ranking ecoregions, in this case by applying a greedy complementarity algorithm, to maxi- mize the capture of seven possible currencies for mammalian conserva-
Phylogenic Trees and the Future of Mammalian Biodiversityâ /â 277 tion: species richness; total numbers of species extinctions predicted in the next 100 years from current Red List status [using the probabilities of extinction in Redding and Mooers (2006)]; two measures combining evolutionary uniqueness and present extinction risk [EDGE (Isaac et al., 2007) and ELEH (Redding and Mooers, 2006), both summed across spe- cies in the ecoregion]; evolutionary dynamism (indexed as the sum of the reciprocal of the terminal branch lengths of species in the ecoregion); evo- lutionary history (estimated as the total branch length in the phylogeny of the ecoregionâs species, adjusted for polytomies as described above); and total latent risk (Cardillo et al., 2006). Species-rich areas are likely to sum to higher values and, hence, rank highly across all currencies. Nonetheless, Table 14.1 shows that rank order differs considerably even for similar measures designed to optimize the same currency (e.g., EDGE and ELEH). The biggest difference is between latent risk and numbers of extinctions. This is perhaps unsurprising, because latent risk gives high weights to species that are less threatened (although not all species were weighted for latent risk, weakening the comparison). A greedy algorithm makes choices based on maximizing the immediate gain in currency at each step and so may give globally suboptimal complementarity networks (Pressey et al., 1996). In addition, ecoregions are, in any case, larger than most current conservation-management units (Grenyer et al., 2006). Our aim here is simply to show that surrogacy among currencies may be low. Furthermore, other more subjective or hard-to-quantify criteria might also be considered importantâfor example, the ecological function of a species or lineage, its cultural significance, or its charismatic appeal to humans (Avise, 2005)âfurther impeding comparisons between areas with similar numbers of species. Even if a currency can be agreed upon, a reserve network optimized for one clade is likely to be suboptimal for another (Grenyer et al., 2006). It is often tempting to use a single group as a surrogate for biodiversity TABLE 14.1â Numbers of Ecoregions Jointly Ranked in the Top 50 (of a Total of 791) by Each Pair of Currencies, According to Greedy Complementarity Searches Extinc- Diversi- Latent Evolutionary Richness EDGE ELEH tions fication Risk History Richness 50 41 22 21 34 15 36 EDGE 50 26 22 33 17 39 ELEH 50 34 18 10 26 Extinctions 50 20 9 21 Diversification 50 16 30 Latent risk 50 16 Evolutionary 50 history
278â /â T. Jonathan Davies et al. as a whole, but other clades are likely to have very differentâand equally complexâpatterns of diversity and extinction risk. The lack of strong surrogacy among groups introduces extra uncertainty into the measured biodiversity value of the regions being considered. In addition, we have focused on conservation benefits rather than costs, but costs vary spatially by several orders of magnitude (Balmford et al., 2003a; Halpern et al., 2006). Costâbenefit models can suggest very different priorities from allocations based solely on perceived biodiversity value (Possingham et al., 2002). So, if we assume that rational decision making must consider both benefits and costs, perhaps the most sensible investment would be in intact but susceptible regions (Balmford et al., 2003a; Cardillo et al., 2006). Public health-care systems may provide a useful analogy: A balanced health-care strategy includes money for preventative medicine as well as for hospital wards and life support. BEYOND THE DECLINES: THE EVOLUTIONARY FUTURE? The primary goal for most conservation management has been to maximize preservation of current diversity. However, by altering the envi- ronment, humans also influence future evolution (Smith and Bernatchez, 2008). A previous National Academy of Sciences Colloquium (Cowling and Pressey, 2001; Myers and Knoll, 2001) raised an important ques- tion: Should conservation goals be extended to consider the evolutionary future? A range of timescales might be considered. In the short term, species-recovery plans can address the requirements needed for continued adaptive evolution within populations (Mace and Purvis, 2008). But what of the longer term? We can identify clades that have recently diversified and the regions in which they are found (Fig. 14.2b). These lineages or areas might represent engines of speciation: Are they therefore conserva- tion priorities? The distinction between maximizing evolutionary his- tory versus centers of diversification is nontrivial. A network of reserves designed to capture maximal evolutionary history would look very dif- ferent from one designed to capture rapidly speciating lineages (Table 14.1), because rapid diversification results in low phylogenetic diversity per species (Mace et al., 2003; Roy and Goldberg, 2007). However, past ââsuccessââ may be a poor indicator of future performance because of the contingent nature of evolution (de Queiroz, 2002; Mace et al., 2003). The geographic pattern of mammalian diversification rates has changed mark- edly through time (Fig. 14.2aâc), and different lineages have also radiated at different times (Alroy, 2000; Bininda-Emonds et al., 2007). One secure prediction is that future environmental conditions will almost certainly differ from those in the past. Nonetheless, extrapolating current trends, key environmental changes are likely to include increas-
Phylogenic Trees and the Future of Mammalian Biodiversityâ /â 279 ing habitat loss and fragmentation, drastic shifts in species abundances and distributions, and climate change (see above) (Myers and Knoll, 2001; Barraclough and Davies, 2005). Changes in the biotic and abiotic compo- sition of the environment (including extinctions) may restructure niche space and moderate constraints on diversification imposed by niche satu- ration. Major turnovers in species composition may therefore be expected (Tilman and Lehman, 2001). Unfortunately, we lack detailed information on the past states of these attributes in geological history and so cannot eas- ily construct quantitative empirical models that can be projected forward. Predicting the evolutionary future is hampered by large uncertainty about the magnitude and form of environmental change and by lineage-specific responses. If we wished to safeguard the evolutionary future, a sensible strategy would be to maintain a set of species that is overdispersed with respect to their ecological adaptations and (as a simple proxy) their phy- logeny, maximizing the possibility of having the right set of features in an uncertain future (Crozier, 1997). Returning to our original question, should conservation goals con- sider the long-term evolutionary future? We sound two notes of caution. First, although we have restricted our focus to mammals, they are only a tiny branch on the Tree of Life, and many of the major limbs, from ciliates to Chlamydia, may be better insulated from anthropogenic disturbances (Nee, 2005). The evolutionary future of life on Earth is therefore unlikely to be in serious jeopardy. Second, anthropogenic environmental change and extinctions are occurring on the order of tens to hundreds of years, but times to speciation are frequently estimated in thousands to millions of years (Barraclough and Davies, 2005), and recovery times after previ- ous mass extinction events were perhaps 5â10 million years (Erwin, 2001). These timescales are too great for practical management. Diversity will almost certainly rebound after the current extinction event; however, it may be composed of species descended from a different, as yet unknown, subset of lineages from those that dominate now, and humans will likely not be included among them. Practical conservation should retain its focus on minimizing declines and extinctions in the present day. ACKNOWLEDGMENTS We thank Natalie Cooper, Lynsey McInnes, Jim Regetz, Gavin Thomas, and Nicola Toomey for discussion and John Avise and Doug Erwin for comments on the manuscript. This work was supported by the Natural Environment Research Council (U.K.), the National Science Foundation (U.S.), the European Union HOTSPOTS Early Stage Training Network, a Deutsche Forschungsgemeinschaft Heisenberg Scholarship, and the Lever- hulme Trust.