Skip to main content
  • Other Publications
    • Philosophical Transactions B
    • Proceedings B
    • Biology Letters
    • Open Biology
    • Philosophical Transactions A
    • Proceedings A
    • Royal Society Open Science
    • Interface
    • Interface Focus
    • Notes and Records
    • Biographical Memoirs

Advanced

  • Home
  • Content
    • Latest issue
    • All content
    • Subject collections
    • Special features
    • Videos
  • Information for
    • Authors
    • Reviewers
    • Readers
    • Institutions
  • About us
    • About the journal
    • Editorial board
    • Author benefits
    • Policies
    • Citation metrics
    • Publication times
    • Open access
  • Sign up
    • Subscribe
    • eTOC alerts
    • Keyword alerts
    • RSS feeds
    • Newsletters
    • Request a free trial
  • Submit
You have accessRestricted access

Primate extinction risk and historical patterns of speciation and extinction in relation to body mass

Luke J. Matthews, Christian Arnold, Zarin Machanda, Charles L. Nunn
Published 13 October 2010.DOI: 10.1098/rspb.2010.1489
Luke J. Matthews
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Christian Arnold
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Zarin Machanda
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
Charles L. Nunn
  • Find this author on Google Scholar
  • Find this author on PubMed
  • Search for this author on this site
  • Article
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF
Loading

Abstract

Body mass is thought to influence diversification rates, but previous studies have produced ambiguous results. We investigated patterns of diversification across 100 trees obtained from a new Bayesian inference of primate phylogeny that sampled trees in proportion to their posterior probabilities. First, we used simulations to assess the validity of previous studies that used linear models to investigate the links between IUCN Red List status and body mass. These analyses support the use of linear models for ordinal ranked data on threat status, and phylogenetic generalized linear models revealed a significant positive correlation between current extinction risk and body mass across our tree block. We then investigated historical patterns of speciation and extinction rates using a recently developed maximum-likelihood method. Specifically, we predicted that body mass correlates positively with extinction rate because larger bodied organisms reproduce more slowly, and body mass correlates negatively with speciation rate because smaller bodied organisms are better able to partition niche space. We failed to find evidence that extinction rates covary with body mass across primate phylogeny. Similarly, the speciation rate was generally unrelated to body mass, except in some tests that indicated an increase in the speciation rate with increasing body mass. Importantly, we discovered that our data violated a key assumption of sample randomness with respect to body mass. After correcting for this bias, we found no association between diversification rates and mass.

1. Introduction

Several studies have implicated body size as a risk factor for extinction in mammals [1–4]. Large body size is thought to correlate with higher extinction risk through a number of life-history covariates of increased body mass, such as longer generation times and smaller litters [2,5]. These factors should increase the time needed to recover from stochastic demographic reductions in population size, thus increasing the probability of extinction. To test this hypothesis, studies have quantified extinction risk among living species as a ranked variable derived from the IUCN Red List conservation status categories (www.iucnredlist.org; [1,2,4]). In mammals, Cardillo et al. [2] found that ‘intrinsic’ biological life-history traits influenced extinction risk only when species were greater than 3 kg in mass. Below 3 kg mass, the primary determinants of extinction risk were ‘extrinsic’ factors that are not genetically heritable, such as geographical range size and nearby human population density [2]. Studies also have failed to detect a consistent association between diversification rate and body mass across mammalian clades [5,6].

Understanding the predictors of present-day extinction risk, as reflected by the IUCN Red List, has clear relevance given the pressing need for primate conservation [7]. For example, the great apes are all large bodied and critically endangered (IUCN Red List). An important question is whether large-bodied primate species experience higher extinction rates irrespective of human activities. If so, then large-bodied primates may be especially vulnerable to anthropogenic drivers of extinction, such as reductions in habitat or climate change. It is possible to infer speciation and extinction rates from a dated phylogeny because higher rates of extinction relative to speciation should produce longer internal branches on a tree with an apparent burst of diversification close to the tips [8,9]. Recent advances provide a way to integrate the study of speciation and extinction that is dependent upon another trait [10], such as body mass. Previous studies have investigated patterns of diversification more generally in primates [6,11–15] using an older inference of primate phylogeny [16].

Using a newly inferred primate phylogeny that enabled us to incorporate phylogenetic uncertainty, we investigated the links between body mass, extinction risk and diversification rates. First, we tested for an association between body mass and present-day extinction risk categories obtained from the IUCN Red List. As with previous studies, we predicted that larger bodied primates are at greater risk of extinction [1,2,4]. In our analyses, we used methods that incorporate phylogenetic uncertainty, better model character evolution on the tree and test for spurious results arising from the ordinal (rather than integer) measurement of extinction threat categories.

Second, we predicted that primate lineages characterized by greater body mass have experienced higher extinction rates throughout evolutionary history. For this, we applied a new method to estimate extinction rates on a phylogenetic tree in relation to a biological characteristic [10], again running the test in a way that incorporated phylogenetic uncertainty.

Last, we turned our attention to speciation by predicting that lineages with smaller body mass have higher speciation rates. We made this prediction based on several studies showing a trend for smaller bodied mammalian clades to be more diverse, possibly because they experience vicariance events more frequently or they are able to partition the environment into more niches [1,2,5,6,17,18].

2. Methods

(a) Tree inference

As our hypotheses pertained to the order Primates as a whole, we needed a phylogeny that included as many species as possible. To this end, we used the trees available from the 10kTrees project, v. 1 (http://10ktrees.fas.harvard.edu/). The website provides extensive documentation on the tree inference, a graphical interface for downloading trees and a number of visualizations of the trees. Details regarding tree inference are available in Arnold et al. [19]. Our use of Bayesian tree inference enabled us to deal with phylogenetic uncertainty by running comparative tests on multiple trees saved from the Markov chain. Incorporating such topological and branch length uncertainty is important because the phylogeny used can affect the conclusions that are drawn from a comparative analysis [20,21].

We used two sets of trees from the Bayesian analysis: a sample of 100 trees distributed evenly along the post-burn-in Markov chain and a consensus tree of all nodes with clade credibility support greater than 0.5. We dated all trees prior to comparative analysis by using seven fossil calibration points employed by previous phylogenetic studies (table 1; [22–26]). We conducted molecular dating with the software r8s [27] using the penalized-likelihood algorithm with a smoothing parameter of 100, chosen because this value best recovered dates inferred from phylogenetic analyses of smaller taxonomic samples but with more extensive sequence data [23–25].

View this table:
  • View inline
  • View popup
Table 1.

Fossil calibration ranges used to date molecular phylogenies. (n.a., not applicable.)

(b) Body mass and IUCN Red List data

Body mass data were obtained from the work of Smith & Jungers [28]. We calculated the mean female body mass across study sites for our analysis. In this manner, we obtained body mass data on 160 species that could be matched to our phylogeny through translation via the taxonomy of Wilson & Reeder [29]. Analyses used the natural log of female body mass. We obtained IUCN conservation status for the species on our phylogeny from the IUCN Red List website (www.iucnredlist.org). From each IUCN category, we constructed an ordinal variable with higher ranks corresponding to greater extinction threat. Specifically, ranked from smallest to highest, we used the following categories of threat: least concern, near threatened, vulnerable, endangered and critically endangered. Six species labelled ‘data deficient’ in the IUCN Red List were removed from the analysis, reducing the sample size to 154 species. The highest rank was assigned to ‘critically endangered’ rather than ‘extinct in the wild’ because our study included only extant species as data points.

(c) Comparative analyses

We tested for an association between body size and IUCN extinction risk through a phylogenetic generalized least-squares (PGLS) model applied across our Bayesian tree block. This analysis of the primate data replicates previous independent contrasts analyses across mammals (e.g. [2,4]). However, PGLS implemented in the program BayesTraits [30] enabled us to incorporate Bayesian estimation of the branch-scaling parameter λ [31] across a tree block, which should improve upon contrasts-based approaches. In particular, we suspected that Brownian motion along the branch lengths in our block of 100 trees might not reflect evolutionary change in a character such as threat status. Brownian motion along branch lengths is an assumption of independent contrasts. The scaling parameter λ adjusts the internal branch lengths with a length multiplier such that the data meet the assumption of Brownian motion. By systematically searching through many possible values for λ, the program locates the value that makes the data most likely and thus best accommodates the assumption of Brownian motion. In primates, Purvis et al. [32] showed that IUCN rank alone exhibited λ = 0.77 that excluded both zero and one (95% CI 0.51–0.90). Thus, we considered it important to estimate λ in our PGLS, and especially to do so in a way that incorporates uncertainty in phylogenetic relationships and branch lengths.

The use of PGLS or independent contrasts treats IUCN extinction risk as a continuous variable [33,34]. Counter to this assumption, extinction risk codes in the IUCN Red List are not continuously varying. Instead, extinction risk is an ordinal variable in which ranks probably vary in the amount of difference in the actual underlying extinction risk [32]. For example, the true (quantitative) difference in extinction risk may differ between categories of near threatened and vulnerable, when compared with endangered and critically endangered. Treating ordinal variables as continuous can produce elevated type 1 error rates because such treatment applies arithmetic operations that do not preserve the variance structure of the original ordinal ranks [35,36]. The problem occurs specifically when the ordinal ranks are separated by unequal distances along the underlying continuous variable that they measure, a point acknowledged by Purvis et al. [32].

We used computer simulations to assess whether the treatment of IUCN threat categories as continuously varying may have introduced error into this study and previous studies [2,4,32]. Specifically, we tested whether treating threat status as continuously varying resulted in elevated type 1 error rates. We conducted 1000 simulations of two uncorrelated continuous characters on our consensus tree (function ‘sim.char’ in the R package ‘geiger’; [37,38]). Characters evolved randomly with a constant accumulation of variance (set at 1.0 per unit; branch length) and an initial state of zero. We then rescored one character from each pair into a set of ordinal ranks with the same number of species in each state as was observed for each corresponding IUCN threat status in our dataset. The lowest continuous values of the simulated character thus became rank ‘0’, while the highest became rank ‘4’, with other ranks derived from intermediate values and all in matching proportion to the observed frequency of each rank in the IUCN data. This rendered one character a true continuous distribution, such as body size, but the other character had been rescored into a set of ordinal ranks with different real distances between the means of each rank. We then conducted 1000 PGLS tests to test for a significant association of the two characters, one of which was ordinal and one continuous. Significant associations were counted as type 1 errors because the characters evolved independently during the simulations. An alternative approach to deal with ordinal IUCN data would be to include additional parameters that model the ordinal nature of the dependent variable. While this procedure in principle can be implemented [39], given the prominent prior research that treated IUCN as continuous, we preferred to use simulations to test this approach more generally. Additionally, ordinal data do not necessarily violate the assumptions of linear regression if the states are sufficiently evenly spaced [35].

To test for effects of body size on the actual speciation and extinction rates experienced by different lineages on the primate tree, we employed the ‘binary-state speciation and extinction’ (BiSSE) test [10], as implemented in the R package diversitree [40,41]. This procedure uses likelihood methods to test a six-parameter model of speciation, extinction and trait evolution. Following the procedure recommended in Maddison et al. [10], we constructed five- and six-parameter models to test body-size-dependent speciation and extinction rates, which we investigated separately. For each state of a binary character, BiSSE can model a rate of speciation, extinction and character state transitions. Thus, six rates are possible for the most complex BiSSE model. The five-parameter models constrained the speciation (or extinction) rates to be equal in the two body size categories, while the six-parameter models allowed speciation (or extinction) rates to vary for different body sizes. Because these models were nested, we assessed statistical significance using likelihood ratio tests.

To use BiSSE with our measures of body mass, we binned species into categories of ‘small’ and ‘large’ body mass. As it is not immediately clear what cut-off should be used for these categories, we initially tested our hypotheses with three different cut-offs derived from the primatological literature. The first break point, 500 g, reflects the point above which primate species are not strictly faunivorous (i.e. Kay's threshold) [42]. We thought this widely accepted energetic constraint might reflect the life-history variables that underlie previously observed associations between body mass and extinction risk. The second break point, 984 g, was the phylogenetic mean body mass at the root of the dated consensus tree, as inferred through squared change parsimony reconstruction [43] of log body mass values in Mesquite v. 2.6 [44]. The third break point of 3000 g came from a study of extinction risk across mammals [2], in which the authors found that biologically intrinsic life-history traits only influenced extinction risk above a body mass of 3000 g. We ran these analyses of the three a priori break points across the 100 dated trees from the Bayesian tree search.

To assess the sensitivity of our results to the break point used for binning species, we conducted an analysis with a sliding break point that divided the body mass data into small and large categories along 20 intervals of 0.2 log units of body mass. We quantified the precision of the sliding break point estimates by sampling the six-parameter model with a Markov chain Monte Carlo (MCMC) that we initialized with the maximum-likelihood estimates for parameter values (R package diversitree; [40,41]). We used only the dated consensus tree for the sliding break point analysis. We constructed 95 per cent credible intervals for parameter values by calculating the range of each parameter across the MCMC samples after excluding the largest and smallest 2.5 per cent of sampled values.

The implementation of BiSSE in diversitree allows one to specify the degree of species sampling employed in the study. Doing so is important because methods to test trait-dependent extinction are biased by incomplete sampling of the study group [8,41]. The diversitree package includes modified-likelihood equations to account for this bias [41]. Under the Wilson & Reeder [29] taxonomy, which served as the taxonomy for the present study, we sampled 42 per cent of extant primates. We used this value as our sampling percentage in all BiSSE analyses.

The BiSSE correction for incomplete sampling is statistically valid only when the species have been sampled randomly. Two issues are relevant here: random sampling of species from the phylogeny, and random sampling with respect to the character hypothesized to affect speciation and extinction rates.

To assess the randomness of our species sample with respect to phylogeny, we conducted a G-test of proportions of species from each genera within our observed sample ([45] R script by P. Hurd, http://www.psych.ualberta.ca/∼phurd/cruft/). The G-test uses maximum-likelihood techniques to assess whether an observed proportion of species within genera can be viewed as a random sample of species given their known frequencies from the complete taxonomy [29]. We conducted this test at the generic level because, based on our primate phylogeny [19], generic classifications reflect monophyletic clades.

To test whether sampling was random with respect to body mass, we calculated the deviation of the observed number of species sampled from a genus compared with the number expected to be sampled given the number of species in each genus in the complete taxonomy. Thus, a positive deviation indicated that more species were sampled from a genus than expected under random sampling, whereas a negative number indicated fewer were sampled than expected. We then conducted a PGLS test for an association between body mass and sampling deviation using our dated consensus tree to control for phylogenetic non-independence (R packages ape [46] and nlme [47]).

We then retested our hypothesis regarding speciation and body size using the diversification test of Freckleton et al. [15]. This test is simply a PGLS of a character trait as a dependent variable with the number of nodes from the tree root to each tip as the independent variable. If a particular character state, such as body size, is associated positively with diversification rates, then a positive association should exist between the character value and lineages that exhibit more nodes. The advantage of this test is that, because it is conducted through a PGLS framework, we were able to include the deviations from the randomly expected sampling as an independent variable; in other words, we could investigate the effects of body mass on diversification while controlling for body-mass-related biased sampling if it exists, even if not statistically significant in the above tests. Another advantage of this test is it treats body size as a continuous variable, whereas the use of BiSSE required us to bin body size into a set of binary categories. The disadvantage of the Freckleton et al. [15] test is that it investigates the association between net diversification rate and body size, rather than specifically investigating speciation and extinction separately (as in BiSSE). We used the dated consensus tree for this analysis.

3. Results

(a) Extinction risk and body mass

In our dataset, 68 species were categorized as ‘least concern’, nine as ‘near threatened’, 30 as ‘vulnerable’, 40 as ‘endangered’ and seven as ‘critically endangered’ (table S1, electronic supplementary material). We found a significant and positive association between IUCN extinction risk and body mass in a PGLS analysis across our dated Bayesian tree block (n = 154, β = 0.40, p = 0.004). Additionally, we found that λ differed from zero and one (λ = 0.78, 95% credible interval λ = 0.63–0.89). This indicates significant phylogenetic signal and thus highlights the importance of controlling for phylogeny, but further shows that a non-Brownian component may contribute to variation in residual values. In our 1000 simulations that tested the effect of ordinal ranking for IUCN, only 5 per cent of statistical tests were deemed significant at the 0.05 nominal level. Thus, the type 1 error rate was as expected.

(b) Extinction rate and body mass

In our BiSSE analyses of 160 primate species, likelihood ratio tests supported no significant relationship between log body mass and historical extinction rate either at our a priori break points or in our sliding break point analysis (figure 1a). Similarly, across the 100 trees, extinction rates at the a priori break points were significantly different relative to body mass for only 1 of the 100 trees sampled from our tree block. This was true even though extinction rates were estimated to be different from zero in many cases (figure 1a). Several papers have indicated that methods for inferring trait-dependent extinction may suffer from low statistical power [10,11,18,48,49], a point to which we return in §4.

Figure 1.
  • Download figure
  • Open in new tab
  • Download powerpoint
Figure 1.

Sliding break point analysis of body-mass-dependent speciation and extinction rates on the consensus tree. Solid lines connecting points show rates for species that have larger body mass than the given point on the body mass in the x-axis (black 95% credible interval bars), while dashed lines are for species with smaller body mass (grey 95% credible interval bars). Triangles in (a) show inferred extinction rates and squares in (b) show inferred speciation rates. Thick vertical arrows correspond to our three a priori hypotheses.

(c) Speciation rate and body mass

In contrast to the results from the extinction analysis, we found significant speciation rate differences at Kay's threshold a priori break point on our dated consensus phylogeny and across 95 of the 100 trees (for the consensus tree, p = 0.016). Differences in the speciation rate at the phylogenetic mean and 3 kg break points were not significant on the consensus tree or on any of the trees from the tree block. Importantly, the result for Kay's threshold was opposite to our prediction, with higher estimated speciation rates in larger bodied species (figure 1b). This pattern seemed to apply across most of the variation in primate body mass, but not at the largest body mass categories (figure 1b). However, MCMC sampling of the speciation estimates showed that 95 per cent credible intervals overlapped across most of the body mass sliding break points (figure 1b).

In testing for speciation and extinction rates, BiSSE also estimated two rates of transition (small to large and large to small) in all analyses. The transition rate estimates showed overlapping credible intervals across almost the entire range of sliding break points (figure S1, electronic supplementary material). This indicates rates of body size increase and decrease are generally equivalent regardless of the break point used to identify large and small species.

The BiSSE results are not compromised by non-random sampling of species in our dataset with respect to phylogeny. The G-test for independence showed that the observed sampling of species within genera was within the variation expected from random sampling of taxa (p = 0.147). However, the observed sample of species was biased towards greater sampling of species with larger body mass in a PGLS test (β = 0.80, p = 0.037). Thus, biased sampling towards species with larger body mass may have influenced the BiSSE results; specifically, the denser sampling of the larger species may have biased the test to estimate a higher rate of speciation because these lineages would have occurred more commonly on the tree representing the sample of species.

When we controlled for the observed sampling deviations from the random sampling expectation, our PGLS test of diversification effects on body size did not find any association between these variables (β = 0.00, p = 0.71). When we removed the sampling variable from the analysis, however, the PGLS obtained a trend more consistent with the results of the BiSSE analysis (β = 0.08, p = 0.081). Although this PGLS test cannot technically distinguish between body-mass-related speciation and extinction effects, we would expect to find a positive association between diversification and body mass if speciation rate correlates positively with body mass and extinction rates do not (as suggested by the above tests). Additionally, simulations by Freckleton et al. [15] showed that their test was more sensitive to speciation effects than to extinction effects. Thus, given the lower power of the test to detect an effect of extinction, any positive result would more probably result from speciation effects.

4. Discussion

We tested a series of predictions involving the links among body mass, phylogenetic diversification and current extinction risk in primates. With regard to current patterns of extinction, we found that larger bodied species experience higher extinction risk. This result replicates previous findings, but for the first time, to our knowledge, incorporates phylogenetic uncertainty and uses the scaling parameter λ to control for the non-Brownian distribution of variation in IUCN threat status categories [2–4]. We also tested, for the first time, whether the ordinal nature of IUCN threat status categories impacts the statistical performance of independent contrasts and PGLS analyses that use it as a dependent variable. Our 1000 phylogenetic simulations show the ordinal coding itself does not produce elevated Type 1 error in our primate dataset.

Despite finding a strong association between threat categories and body mass in primates, we failed to find evidence for an historical association between body mass and extinction rates. Instead, in a limited number of tests, we found that speciation rates may increase with body mass, which was opposite to our predictions. These analyses point to an association between larger body mass and higher speciation rate, which contrasts with most predictions, including ours, that speciation rates covary negatively with body mass (see also [50]). After controlling for biased sampling of species in a PGLS model, however, we found no association between the diversification rate—a function of speciation and extinction—and body mass.

Previous studies found that intrinsic biological variables can explain a substantial proportion of variation in risk status [1,2]. These previous statistical models found that high risk variables, including large body mass, high trophic level, small geographical range and slow life history, explain threat status to a substantial degree and independently of anthropogenic effects [1,2,51]. Thus, the present-day extinction risk status should be related to the historically experienced extinction rates of lineages.

We found it surprising that current patterns of extinction risk covaried with body mass in our study and previous studies, while estimates of extinction rates across primate phylogeny failed to show this effect. It could be that higher threat categories among large-bodied primates have resulted only from anthropogenic effects in the present. Alternatively, the tests may have failed to detect differences because available methods to estimate historical extinction rates suffer from low statistical power ([10], see also [11,18,48,49,52,53]). Using simulations, for example, Maddison et al. [10] showed that the power to detect variation in speciation rates is relatively high (power for speciation in simulations of 500 species was approx. 58%), when compared with the power of detecting variation in extinction rates (approx. 21%). Rabosky [53] recently showed that, in addition to low power, current methods cannot accurately estimate extinction rates when these rates vary across the tree. His analyses were conducted with complete sampling and phylogenetic information of every extant species. Given the fossil evidence for a plethora of extinct primate species, we agree with Rabosky [53] that future analyses should incorporate fossil information as a means to increase power and to estimate extinction rates more accurately (see also [51]).

Our study also offers a cautionary tale about the importance of testing the assumptions of methods for studying speciation and extinction rates. When using a new method to account for incomplete sampling of taxa, we found some support for an association between speciation and body mass. However, our data violated another assumption of these methods, namely that sampling of species is random with respect to body mass. When we controlled for non-random sampling, the association between diversification rates and body mass became non-significant. Thus, we offer a valuable empirical example of shortcomings to these methods that complement previous simulation studies and provide a procedure for investigating biased sampling for others to use in future studies. The potential effects of non-random sampling are a serious concern for previous studies of trait-dependent diversification [11,14,15,54].

We began this study with the prediction that populations of smaller bodied animals may speciate more quickly, as they experience vicariance events more frequently or they are able to better partition the environment into more niches [18]. Previous studies provided hints that the speciation rate could increase with body mass in primates, rather than showing a strict negative correlation with body mass. For example, Purvis et al. [11] found that the family Cercopithecidae (Old World monkeys) experienced larger speciation rates than did other primate lineages. While their study did not assess how speciation covaried with characteristics such as body mass, our results are consistent with their findings given that the cercopithecids are generally larger bodied than extant platyrrhine and strepsirrhine primates that more frequently fall below Kay's threshold. The cercopithecids are also typically smaller than the (relatively) species depauperate hominoids. These clade-specific patterns could be tested more rigorously with the new statistical software MEDUSA [55]. Our results are consistent with Freckleton et al.'s [15] findings of a significant correlation between body mass and diversification rate. The results imply that some degree of species selection on body mass may have occurred in primates, in that the influence of body mass on cladogenesis explains the distribution of body mass even after allowing for anagenic change as it is incorporated into the model (figure S1, electronic supplementary material; [10,11,56]).

In summary, we used new methods to incorporate phylogenetic uncertainty and to assess how body mass influences primate diversification. Our results add to a growing body of evidence that larger bodied animals are more susceptible to extinction. The lack of significance for historical patterns of extinction may indicate that methods to detect extinction rates from extant species are severely compromised in statistical power [10,11,18,48,49,53]. Intriguingly, we found some evidence that speciation rates appear to increase with body mass in primates (see also [50]). It will be interesting to see if similar patterns occur in other groups of vertebrates and, if so, how and why the effects of body size on speciation vary among clades.

Acknowledgements

This research was supported by Harvard University and by the National Science Foundation (BCS-0923791). We thank Daniel Lieberman for his helpful discussions regarding diversification mechanisms, Richard Wrangham's laboratory group for helpful comments, and two anonymous referees and Andy Purvis for valuable suggestions during the review process.

  • Received July 13, 2010.
  • Accepted September 22, 2010.
  • This Journal is © 2010 The Royal Society

References

  1. ↵
    1. Purvis A.,
    2. Gittleman J. L.,
    3. Cowlishaw G.,
    4. Mace G. M.
    2000 Predicting extinction risk in declining species. Proc. R. Soc. Lond. B 267, 1947–1952. doi:10.1098/rspb.2000.1234 (doi:10.1098/rspb.2000.1234)
    OpenUrlAbstract/FREE Full Text
  2. ↵
    1. Cardillo M.,
    2. Mace G. M.,
    3. Jones K. E.,
    4. Bielby J.,
    5. Bininda-Emonds O. R. P.,
    6. Sechrest W.,
    7. Orme C. D. L.,
    8. Purvis A.
    2005 Multiple causes of high extinction risk in large mammal species. Science 309, 1239–1241. doi:10.1126/science.1116030 (doi:10.1126/science.1116030)
    OpenUrlAbstract/FREE Full Text
    1. Clauset A.,
    2. Erwin D. H.
    2008 The evolution and distribution of species body size. Science 321, 399–401. doi:10.1126/science.1157534 (doi:10.1126/science.1157534)
    OpenUrlAbstract/FREE Full Text
  3. ↵
    1. Davies T. J.,
    2. et al.
    2008 Phylogenetic trees and the future of mammalian biodiversity. Proc. Natl Acad. Sci. USA 105, 11 556–11 563. doi:10.1073/pnas.0801917105 (doi:10.1073/pnas.0801917105)
    OpenUrlAbstract/FREE Full Text
  4. ↵
    1. Isaac N. J. B.,
    2. Jones K. E.,
    3. Gittleman J. L.,
    4. Purvis A.
    2005 Correlates of species richness in mammals: body size, life history, and ecology. Am. Nat. 165, 600–607. doi:10.1086/429148 (doi:10.1086/429148)
    OpenUrlCrossRefPubMedWeb of Science
  5. ↵
    1. Gittleman J. L.,
    2. Purvis A.
    1998 Body size and species-richness in carnivores and primates. Proc. R. Soc. Lond. B 265, 113–119. doi:10.1098/rspb.1998.0271 (doi:10.1098/rspb.1998.0271)
    OpenUrlAbstract/FREE Full Text
  6. ↵
    1. Cowlishaw G.,
    2. Dunbar R. I. M.
    2000 Primate conservation biology. Chicago, IL: University of Chicago Press.
  7. ↵
    1. Nee S.,
    2. Holmes E. C.,
    3. May R. M.,
    4. Harvey P. H.
    1994 Extinction rates can be estimated from molecular phylogenies. Phil. Trans. R. Soc. Lond. B 344, 77–82. doi:10.1098/rstb.1994.0054 (doi:10.1098/rstb.1994.0054)
    OpenUrlAbstract/FREE Full Text
  8. ↵
    1. Nee S.
    2001 Inferring speciation rates from phylogenies. Evolution 55, 661–668. doi:10.1554/0014-3820(2001)055[0661:ISRFP]2.0.CO;2 (doi:10.1554/0014-3820(2001)055[0661:ISRFP]2.0.CO;2)
    OpenUrlCrossRefPubMedWeb of Science
  9. ↵
    1. Maddison W. P.,
    2. Midford P. E.,
    3. Otto S. P.
    2007 Estimating a binary character's effect on speciation and extinction. Syst. Biol. 56, 701–710. doi:10.1080/10635150701607033 (doi:10.1080/10635150701607033)
    OpenUrlAbstract/FREE Full Text
  10. ↵
    1. Purvis A.,
    2. Nee S.,
    3. Harvey P. H.
    1995 Macroevolutionary inferences from primate phylogeny. Proc. R. Soc. Lond. B 260, 329–333. doi:10.1098/rspb.1995.0100 (doi:10.1098/rspb.1995.0100)
    OpenUrlAbstract/FREE Full Text
    1. Chan K. M. A.,
    2. Moore B. R.
    2002 Whole-tree methods for detecting differential diversification rates. Syst. Biol. 51, 855–865. doi:10.1080/10635150290102555 (doi:10.1080/10635150290102555)
    OpenUrlAbstract/FREE Full Text
    1. Bininda-Emonds O.
    1. Moore B. R.,
    2. Chan K. M. A.,
    3. Donoghue M. J.
    2004 Detecting diversification rate variation in supertrees. In Phylogenetic supertrees: combining information to reveal the tree of life (ed. Bininda-Emonds O.), pp. 487–533. Dordrecht, The Netherlands: Kluwer Academic Publishers.
  11. ↵
    1. Paradis E.
    2005 Statistical analysis of diversification with species traits. Evolution 59, 1–12.
    OpenUrlPubMedWeb of Science
  12. ↵
    1. Freckleton R. P.,
    2. Phillimore A. B.,
    3. Pagel M.
    2008 Relating traits to diversification: a simple test. Am. Nat. 172, 102–115. doi:10.1086/588076 (doi:10.1086/588076)
    OpenUrlCrossRefPubMedWeb of Science
  13. ↵
    1. Purvis A.
    1995 A composite estimate of primate phylogeny. Phil. Trans. R. Soc. Lond. B 348, 405–421. doi:10.1098/rstb.1995.0078 (doi:10.1098/rstb.1995.0078)
    OpenUrlAbstract/FREE Full Text
  14. ↵
    1. Maurer B. A.,
    2. Brown J. H.,
    3. Rusler R. D.
    1992 The micro and macro in body size evolution. Evolution 46, 939–953. doi:10.2307/2409748 (doi:10.2307/2409748)
    OpenUrlCrossRefWeb of Science
  15. ↵
    1. Gardezi T.,
    2. da Silva J.
    1999 Diversity in relation to body size in mammals: a comparative study. Am. Nat. 153, 110–123. doi:10.1086/303150 (doi:10.1086/303150)
    OpenUrlCrossRefWeb of Science
  16. ↵
    1. Arnold C.,
    2. Matthews L. J.,
    3. Nunn C. L.
    2010 The 10kTrees website: a new online resource for primate phylogeny. Evol. Anthropol. 19, 114–118. doi:10.1002/evan.20251 (doi:10.1002/evan.20251)
    OpenUrlCrossRefWeb of Science
  17. ↵
    1. Lutzoni F.,
    2. Pagel M.,
    3. Reeb V.
    2001 Major fungal lineages are derived from lichen symbiotic ancestors. Nature 411, 937–940. doi:10.1038/35082053 (doi:10.1038/35082053)
    OpenUrlCrossRefPubMedWeb of Science
  18. ↵
    1. Lässig M.,
    2. Valleriani A.
    1. Pagel M.,
    2. Lutzoni F.
    2002 Accounting for phylogenetic uncertainty in comparative studies of evolution and adaptation. In Biological evolution and statistical physics (eds Lässig M., Valleriani A.), pp. 148–161. Berlin, Germany: Springer-Verlag.
  19. ↵
    1. Seiffert E. R.,
    2. Simons E. L.,
    3. Attia Y.
    2003 Fossil evidence for an ancient divergence of lorises and galagos. Nature 422, 421–424. doi:10.1038/nature01489 (doi:10.1038/nature01489)
    OpenUrlCrossRefGeoRefWeb of Science
  20. ↵
    1. Yang Z. H.,
    2. Yoder A. D.
    2003 Comparison of likelihood and Bayesian methods for estimating divergence times using multiple gene loci and calibration points, with application to a radiation of cute-looking mouse lemur species. Syst. Biol. 52, 705–716. doi:10.1080/10635150390235557 (doi:10.1080/10635150390235557)
    OpenUrlAbstract/FREE Full Text
    1. Yoder A. D.,
    2. Yang Z. H.
    2004 Divergence dates for Malagasy lemurs estimated from multiple gene loci: geological and evolutionary context. Mol. Ecol. 13, 757–773. doi:10.1046/j.1365-294X.2004.02106.x (doi:10.1046/j.1365-294X.2004.02106.x)
    OpenUrlCrossRefPubMedWeb of Science
  21. ↵
    1. Hodgson J. A.,
    2. Sterner K. N.,
    3. Matthews L. J.,
    4. Burrell A. S.,
    5. Jani R. A.,
    6. Raaum R. L.,
    7. Stewart C. B.,
    8. Disotell T. R.
    2009 Successive radiations, not stasis, in the South American primate fauna. Proc. Natl Acad. Sci. USA 106, 5534–5539. doi:10.1073/pnas.0810346106 (doi:10.1073/pnas.0810346106)
    OpenUrlAbstract/FREE Full Text
  22. ↵
    1. Godinot M.
    2006 Lemuriform origins as viewed from the fossil record. Folia Primatologica 77, 446–464. doi:10.1159/000095391 (doi:10.1159/000095391)
    OpenUrlCrossRefPubMedWeb of Science
  23. ↵
    1. Sanderson D. S.
    2002 Estimating absolute rates of molecular evolution and divergence times: a penalized likelihood approach. Mol. Biol. Evol. 19, 101–109.
    OpenUrlAbstract/FREE Full Text
  24. ↵
    1. Smith R. J.,
    2. Jungers W. L.
    1997 Body mass in comparative primatology. J. Hum. Evol. 32, 523–559. doi:10.1006/jhev.1996.0122 (doi:10.1006/jhev.1996.0122)
    OpenUrlCrossRefPubMedWeb of Science
  25. ↵
    1. Wilson D. E.,
    2. Reeder D. M.
    2005 Mammal species of the world. Baltimore, MD: Johns Hopkins University Press.
  26. ↵
    1. Pagel M.,
    2. Meade A.
    2007 BayesTraits. Reading, UK, version 1.0. See www.evolution.rdg.ac.uk.
  27. ↵
    1. Pagel M.
    1999 Inferring the historical patterns of biological evolution. Nature 401, 877–884. doi:10.1038/44766 (doi:10.1038/44766)
    OpenUrlCrossRefPubMed
  28. ↵
    1. Purvis A.,
    2. Gittleman J. L.,
    3. Brooks T.
    1. Purvis A.,
    2. Cardillo M.,
    3. Grenyer R.,
    4. Collen B.
    2005 Correlates of extinction risk: phylogeny, biology, threat and scale. In Conservation biology 8: phylogeny and conservation (eds Purvis A., Gittleman J. L., Brooks T.), pp. 295–316. New York, NY: The Zoological Society of London.
  29. ↵
    1. Felsenstein J.
    1985 Phylogenies and the comparative method. Am. Nat. 125, 1–15. doi:10.1086/284325 (doi:10.1086/284325)
    OpenUrlCrossRefWeb of Science
  30. ↵
    1. Garland T.,
    2. Midford P. E.,
    3. Ives A. R.
    1999 An introduction to phylogenetically based statistical methods, with a new method for confidence intervals on ancestral values. Am. Zoologist 39, 374–388.
    OpenUrlWeb of Science
  31. ↵
    1. Stevens S. S.
    1946 On the theory of scales of measurement. Science 103, 677–680. doi:10.1126/science.103.2684.677 (doi:10.1126/science.103.2684.677)
    OpenUrlFREE Full Text
  32. ↵
    1. Siegel S.
    1956 Nonparametric statistics for the behavioral sciences, pp. 23–26. New York, NY: McGraw-Hill Book Company, Inc.
  33. ↵
    1. Harmon L.,
    2. Weir J.,
    3. Brock C.,
    4. Glor R.,
    5. Challenger W.,
    6. Hunt G.
    2009. geiger: analysis of evolutionary diversification. R package version 1.3-1. See http://CRAN.R-project.org/package=geiger.
  34. ↵
    R Development Core Team. 2009 R: a language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing. ISBN 3-900051-07-0. See URL http://www.R-project.org.
  35. ↵
    1. Hadfield J. D.,
    2. Nakagawa S.
    2010 General quantitative genetic methods for comparative biology: phylogenies, taxonomies and multi-trait models for continuous and categorical characters. J. Evol. Biol. 23, 494–508. doi:10.1111/j.1420-9101.2009.01915.x (doi:10.1111/j.1420-9101.2009.01915.x)
    OpenUrlCrossRefPubMedWeb of Science
  36. ↵
    1. FitzJohn R. G.
    2009. Diversitree: comparative phylogenetic tests of diversification. R package version 0.2-1. See http://www.zoology.ubc.ca/prog/diversitree.
  37. ↵
    1. FitzJohn R. G.,
    2. Maddison W. P.,
    3. Otto S. P.
    2009 Estimating trait-dependent speciation and extinction rates from incompletely resolved phylogenies. Syst. Biol. 58, 595–611. doi:10.1093/sysbio/syp067 (doi:10.1093/sysbio/syp067)
    OpenUrlAbstract/FREE Full Text
  38. ↵
    1. Rodman P.,
    2. Cant J.
    1. Kay R. F.
    1984 On the use of anatomical features to infer foraging behavior in extinct primates. In Adaptations for foraging in nonhuman primates (eds Rodman P., Cant J.), pp. 21–53. New York, NY: Columbia University Press.
  39. ↵
    1. Maddison W.
    1991 Squared-change parsimony reconstructions of ancestral states for continuous-valued characters on a phylogenetic tree. Syst. Zool. 40, 304–314. doi:10.2307/2992324 (doi:10.2307/2992324)
    OpenUrlAbstract/FREE Full Text
  40. ↵
    1. Maddison W. P.,
    2. Maddison D. R.
    2009. Mesquite: a modular system for evolutionary analysis, version 2.6. See http://mesquiteproject.org.
  41. ↵
    1. Sokal R. R.,
    2. Rohlf F. J.
    1995 Biometry. New York, NY: W. H. Freeman and Company.
  42. ↵
    1. Paradis E.,
    2. Claude J.,
    3. Strimmer K.
    2004 APE: analyses of phylogenetics and evolution in R language. Bioinformatics 20, 289–290. doi:10.1093/bioinformatics/btg412 (doi:10.1093/bioinformatics/btg412)
    OpenUrlAbstract/FREE Full Text
  43. ↵
    1. Pinheiro J.,
    2. Bates D.,
    3. DebRoy S.,
    4. Sarkar S.
    & the R Core team 2009 nlme: linear and nonlinear mixed effects models. R package, v. 3.1–96. Vienna, Austria: R Foundation for Statistical Computing.
  44. ↵
    1. Kubo T.,
    2. Iwasa Y.
    1995 Inferring the rates of branching and extinction from molecular phylogenies. Evolution 49, 694–704. doi:10.2307/2410323 (doi:10.2307/2410323)
    OpenUrlCrossRefWeb of Science
  45. ↵
    1. Paradis E.
    2004 Can extinction rates be estimated without fossils? J. Theor. Biol. 229, 19–30. doi:10.1016/j.jtbi.2004.02.018 (doi:10.1016/j.jtbi.2004.02.018)
    OpenUrlCrossRefPubMedWeb of Science
  46. ↵
    1. Liow L. J.,
    2. Fortelius M.,
    3. Bingham E.,
    4. Lintulaakso K.,
    5. Mannila H.,
    6. Flynn L.,
    7. Stenseth N. C.
    2008 Higher origination and extinction rates in larger mammals. Proc. Natl Acad. Sci. USA 105, 6097–6102. doi:10.1073/pnas.0709763105 (doi:10.1073/pnas.0709763105)
    OpenUrlAbstract/FREE Full Text
  47. ↵
    1. Purvis A.
    2008 Phylogenetic approaches to the study of extinction. Annu. Rev. Ecol. Evol. Syst. 39, 301–319. doi:10.1146/annurev.ecolsys.063008.102010 (doi:10.1146/annurev.ecolsys.063008.102010)
    OpenUrlCrossRefWeb of Science
  48. ↵
    1. Nee S.,
    2. Mooers A. O.,
    3. Harvey P. H.
    1992 Tempo and mode of evolution revealed from molecular phylogenies. Proc. Natl Acad. Sci. USA 89, 8322–8326. doi:10.1073/pnas.89.17.8322 (doi:10.1073/pnas.89.17.8322)
    OpenUrlAbstract/FREE Full Text
  49. ↵
    1. Rabosky D. L.
    2010 Extinction rates should not be estimated from molecular phylogenies. Evolution 64, 1816–1824. doi:10.1111/j.1558-5646.2009.00926.x (doi:10.1111/j.1558-5646.2009.00926.x)
    OpenUrlCrossRefPubMedWeb of Science
  50. ↵
    1. Nakagawa S.,
    2. Freckleton R. P.
    2008 Missing inaction: the dangers of ignoring missing data. Trends Ecol. Evol. 23, 592–596. doi:10.1016/j.tree.2008.06.014 (doi:10.1016/j.tree.2008.06.014)
    OpenUrlCrossRefPubMedWeb of Science
  51. ↵
    1. Alfaro M. E.,
    2. Santini F.,
    3. Brock C.,
    4. Alamillo H.,
    5. Dornburg A.,
    6. Rabosky D. L.,
    7. Carnevale G.,
    8. Harmon L. J.
    2009 Nine exceptional radiations plus high turnover explain species diversity in jawed vertebrates. Proc. Natl Acad. Sci. USA 106, 13 410–13 414. doi:10.1073/pnas.0811087106 (doi:10.1073/pnas.0811087106)
    OpenUrlAbstract/FREE Full Text
  52. ↵
    1. Vrba E. S.
    1984 What is species selection. Syst. Zool. 33, 318–328. doi:10.2307/2413077 (doi:10.2307/2413077)
    OpenUrlFREE Full Text
  53. ↵
    1. Haile-Selassie Y.
    2001 Late Miocene hominids from the Middle Awash, Ethiopia. Nature 412, 178–181. doi:10.1038/35084063 (doi:10.1038/35084063)
    OpenUrlCrossRefPubMedWeb of Science
  54. ↵
    1. Senut B.,
    2. Pickford M.,
    3. Gommery D.,
    4. Mein P.,
    5. Cheboi K.,
    6. Coppens Y.
    2001 First hominid from the Miocene (Lukeino Formation, Kenya). Comptes Rendus De L Academie Des Sciences Serie Ii Fascicule a-Sciences De La Terre Et Des Planetes 332, 137–144.
    OpenUrlGeoRef
  55. ↵
    1. Vignaud P.,
    2. et al.
    2002 Geology and palaeontology of the Upper Miocene Toros-Menalla hominid locality, Chad. Nature 418, 152–155. doi:10.1038/nature00880 (doi:10.1038/nature00880)
    OpenUrlCrossRefPubMed
  56. ↵
    1. Brunet M.,
    2. et al.
    2002 A new hominid from the Upper Miocene of Chad, central Africa. Nature 418, 145–151. doi:10.1038/nature00879 (doi:10.1038/nature00879)
    OpenUrlCrossRefPubMedWeb of Science
  57. ↵
    1. Hartwig W. C.
    1. Kelley J.
    2002 The hominoid radiation in Asia. In The primate fossil record (ed. Hartwig W. C.), pp. 369–384. Cambridge, UK: Cambridge University Press.
  58. ↵
    1. Jablonski N. G.
    1. Leakey M. G.
    1993 Evolution of Theropithecus in the Turkana Basin. In Theropithecus: the rise and fall of a primate genus (ed. Jablonski N. G.), pp. 85–123. Cambridge, UK: Cambridge University Press.
  59. ↵
    1. Young N. M.,
    2. MacLatchy L.
    2004 The phylogenetic position of Morotopithecus. J. Hum. Evol. 46, 163–184. doi:10.1016/j.jhevol.2003.11.002 (doi:10.1016/j.jhevol.2003.11.002)
    OpenUrlCrossRefPubMedWeb of Science
  60. ↵
    1. Hartwig W. C.
    1. Benefit B. R.,
    2. McCrossin M. L.
    2002 The Victoriapithecidae, Cercopithecoidea. In The primate fossil record (ed. Hartwig W. C.), pp. 241–253. Cambridge, UK: Cambridge University Press.
  61. ↵
    1. Hartwig W. C.
    1. Hartwig W. C.,
    2. Meldrum D. J.
    2002 Miocene platyrrhines of the northern Neotropics. In The primate fossil record (ed. Hartwig W. C.), pp. 175–188. Cambridge, UK: Cambridge University Press.
View Abstract
PreviousNext
Back to top
PreviousNext
22 April 2011
Volume 278, issue 1709
Proceedings of the Royal Society B: Biological Sciences: 278 (1709)
  • Table of Contents
  • Table of Contents (PDF)
  • About the Cover
  • Index by author
  • Back Matter (PDF)
  • Ed Board (PDF)
  • Front Matter (PDF)
Share
Primate extinction risk and historical patterns of speciation and extinction in relation to body mass
Luke J. Matthews, Christian Arnold, Zarin Machanda, Charles L. Nunn
Proc. R. Soc. B 2011 278 1256-1263; DOI: 10.1098/rspb.2010.1489. Published 11 March 2011
del.icio.us logo Digg logo Reddit logo Twitter logo CiteULike logo Connotea logo Facebook logo Google logo Mendeley logo
Email

Thank you for your interest in spreading the word on Proceedings of the Royal Society of London B: Biological Sciences.

NOTE: We only request your email address so that the person you are recommending the page to knows that you wanted them to see it, and that it is not junk mail. We do not capture any email address.

Enter multiple addresses on separate lines or separate them with commas.
Primate extinction risk and historical patterns of speciation and extinction in relation to body mass
(Your Name) has sent you a message from Proceedings of the Royal Society of London B: Biological Sciences
(Your Name) thought you would like to see the Proceedings of the Royal Society of London B: Biological Sciences web site.
Print
Manage alerts

Please log in to add an alert for this article.

Sign In to Email Alerts with your Email Address
Citation tools

Primate extinction risk and historical patterns of speciation and extinction in relation to body mass

Luke J. Matthews, Christian Arnold, Zarin Machanda, Charles L. Nunn
Proc. R. Soc. B 2011 278 1256-1263; DOI: 10.1098/rspb.2010.1489. Published 11 March 2011

Citation Manager Formats

  • BibTeX
  • Bookends
  • EasyBib
  • EndNote (tagged)
  • EndNote 8 (xml)
  • Medlars
  • Mendeley
  • Papers
  • RefWorks Tagged
  • Ref Manager
  • RIS
  • Zotero
Download

Article reuse

  • Article
    • Abstract
    • 1. Introduction
    • 2. Methods
    • 3. Results
    • 4. Discussion
    • Acknowledgements
    • References
  • Figures & Data
  • Info & Metrics
  • eLetters
  • PDF

See related subject areas:

  • taxonomy and systematics
  • evolution
  • ecology

Related articles

Cited by

Large datasets are available through Proceedings B's partnership with Dryad

Open biology

  • PROCEEDINGS B
    • About this journal
    • Contact information
    • Purchasing information
    • Submit
    • Author benefits
    • Open access membership
    • Recommend to your library
    • FAQ
    • Help

Royal society publishing

  • ROYAL SOCIETY PUBLISHING
    • Our journals
    • Open access
    • Publishing policies
    • Conferences
    • Podcasts
    • News
    • Blog
    • Manage your account
    • Terms & conditions
    • Cookies

The royal society

  • THE ROYAL SOCIETY
    • About us
    • Contact us
    • Fellows
    • Events
    • Grants, schemes & awards
    • Topics & policy
    • Collections
    • Venue hire

Copyright © 2018 The Royal Society