Research Article | Open | Published:
Why are some species older than others? A large-scale study of vertebrates
BMC Evolutionary Biologyvolume 16, Article number: 90 (2016)
Strong variations are observed between and within taxonomic groups in the age of extant species and these differences can clarify factors that render species more vulnerable to extinction. Understanding the factors that influence the resilience of species is thus a key component of evolutionary biology, but it is also of prime importance in a context of climate change and for conservation in general. We explored the effect of extrinsic and intrinsic factors on the timing of the oldest diversification event in over 600 vertebrate species distributed worldwide. We used phylogenetic comparative methods to show that color polymorphism, latitude and reproduction (the latter through its interaction with latitude) affected the timing of the oldest diversification event within a species.
Species from higher latitudes tended to be younger, and colour-polymorphic species were older than monomorphic species. Mode of reproduction was important also, in that the age of oviparous species decreased with latitude, whereas no pattern was apparent for viviparous species. Organisms which have already persisted for a long time may be more likely to deal with future modifications of their environment.
Species that are colour polymorphic, viviparous, and/or live at low latitudes have exhibited resilience to past environmental changes, and hence may be better able to deal with current climate change.
If we wish to predict the resilience of a species to current challenges, such as those caused by climate change, one important source of information may be the past history of that taxon. A species that has persisted through very long periods – and thus, survived across a wide range of environmental conditions – may be more resilient than a recently evolved taxon whose vulnerability has not yet been tested to the same degree. The age of a species can then be considered as a consequence of resilience, which, in turn, will be determined by specific phenotypic traits. All else being equal then, we might thus expect “older” species to be more resilient. An important question is then what are the factors that affect a species’ age?
Previous research shows that a species’ geographic location is linked to how long it has persisted since being formed (defined as the age of species in [1–3]). On average, terrestrial vertebrate species are older in the Southern than in the Northern Hemisphere (invertebrates, plants, fungi and marine vertebrates have never been considered). The severity of the Last Glacial Maximum was different in these two Hemispheres: extensive ice sheets were almost entirely absent from the Southern Hemisphere, whereas large portions of North America and Eurasia were covered by glaciers [4, 5]. These differences placed fewer constraints on organismal viability, enabling populations to persist over extensive areas in the Southern Hemisphere. As a result, the lack of glaciation reduced the loss of genetic diversity [1, 6]. For similar reasons, the latitudinal distribution of organisms might also impact the age of species. For example, Weir and Schluter  found that within New World birds and mammals, the age of species, as well as the divergence between sister species is older towards the tropics.
Intrinsic factors might also influence the resilience of species. In particular, it has been proposed that the reproductive strategies and neonatal behaviour, phenotypic characteristics, as well as the presence of intraspecific phenotypic variations, can have an effect on species ages (e.g. [8–12]). In this respect, recent studies have shown that intraspecific colour polymorphism allows species to exploit different habitat types or have broader distributions because of potential differences in behaviour, thermoregulation capacities, or prey–predator interactions (e.g. [13–15]). In addition, a study strictly focusing on snakes highlighted that the age of colour polymorphic species was older than those of monomorphic taxa . In combination, these results suggest that being colour polymorphic enhances the resilience of species. Concerning the reproduction strategies and precociality of neonates, studies on the evolution of viviparity in reptiles suggest that it is a key innovation compared to oviparity, buffering species against the negative impact of past climate fluctuations or enabling them to better exploit cold areas [10, 16, 17]. However, it may also be a dead-end in some circumstances as the transition to viviparity is often irreversible, meaning that cold adapted viviparous species might be at high risk in a context of global warming . Lastly, species with independent (precocial) neonates would be expected to be more resilient. Under changing conditions, independent neonates can move by themselves to find a better environment. In contrast, dependent (altricial) neonates need to be carried during moves, which burdens the parents and lessens their survival rate.
Any investigation of intrinsic factors associated with species age must also consider the indirect effect of body size. It is well known that body size is correlated with generation time and life history traits, as well as with population density. For example, small-bodied species have been shown to be more resilient than large ones (e.g. [9, 18]). Nonetheless, this factor has been neglected in previous studies looking at the age of species.
In the present study, we aim to understand how extrinsic and intrinsic factors influence the resilience of vertebrate taxa and hence their age. Based on a dataset of 600 species worldwide and using phylogenetic comparative methods, we tested for an effect of the geographic location, body size, reproduction mode (oviparity versus viviparity), newborn dependence behaviour (precocial, altricial), and presence of intraspecific colour polymorphism on the timing of the oldest diversification event within species. Our study differs from previous studies by including taxa that are distributed in all areas of the world and we consider in our model all the factors that have been hypothesized to affect the age of species. It is also the first to consider all classes of vertebrates simultaneously and consider for the first time the phylogenetic effects due to shared ancestry between species. It therefore provides a much more complete and clear picture of the factors shaping the resilience of species.
The phylogenetic tree for the 601 vertebrate species (including the shark species Squalus acanthias as an outgroup) was congruent with existing knowledge on the relationships of the taxonomic groups included (Fig. 1). Most of the nodes were well supported and the 100 trees sampled from the posterior distribution of trees that were used to account for phylogenetic uncertainty showed small differences (mainly occurring in nodes close to the tips of the tree). The divergence times obtained by maximum likelihood are within the 95 % interval of the posterior samples obtained with BEAST. The two analyses resulted in very similar divergence times and we therefore use the dates obtained by maximum likelihood for the remaining of the analyses. The dates obtained for the main lineages of vertebrates are also congruent with recent large scale phylogenetic reconstructions [19–22]. We used phylogenetic GLS to test the factors affecting the age of species, which was used as the response variable in our analyses. The backward model selection retained three single factors and one interaction in the final phylogenetic GLS model (Fig. 2a; Additional file 1). Colour polymorphism and latitude were highly significant (p-value < 0.001; Fig. 2a), while mode of reproduction was maintained due to its significant interaction with latitude. Polymorphic species had a much higher resilience than monomorphic ones (fig. 2b) and, on average, the age of the polymorphic species was (when transformed back into the original scale) 1.86 Myr older than the age of monomorphic ones (Fig. 2b; Additional file 2). Latitude was also significant, but its effect was influenced by the mode of reproduction. Oviparous species showed a strong latitudinal trend, with species living far from the Equator being younger than species living in tropical climates (Fig. 2c). In contrast, viviparous species did not show any latitudinal trend in species age (Fig. 2c).
Our study, based on a dataset of 600 vertebrate species, allowed us to identify several extrinsic and intrinsic factors that influence significantly the age of vertebrate species. It also clarified results from previous studies that were not based on comparative methods, focused on particular taxa or considered fewer factors (see e.g. [1–3, 7, 12]). Overall, absolute mean latitude of the species distribution was negatively correlated with species age and colour-polymorphic species were older than monomorphic species. Moreover, the age of oviparous species decreased with latitude, whereas no such pattern was apparent for viviparous ones.
The latitudinal gradient found in this study is consistent with the results of Weir and Schluter , who found that within New World birds and mammals, the age of species, as well as the divergence between sister species, is younger towards the poles. Using interspecific divergences, they estimated that the highest recent speciation and extinction rates occurred at high latitudes and declined toward the tropics. This pattern is likely explained by the presence of more stable climatic conditions in the tropics than in temperate areas during Pleistocene climatic fluctuations . According to Weir and Schluter , this climatic stability leads to a lower species turnover in the tropics compared to temperate areas, resulting in a lower diversity in the latter. Based on the similarity between our results and those of Weir and Schluter  in term of intraspecific latitudinal gradient, this pattern might be generally valid for vertebrate species from both Hemispheres. However, the situation is likely more complex. Indeed, the age of viviparous vertebrate species was not influenced by their mean latitudinal distribution, whereas the age of oviparous species decreased towards the poles. This result suggests a higher resilience of viviparous species to cold climate, and can be attributed to their buffering capacity against inadequate environmental conditions during the embryogenesis. This is particularly true for heliothermic vertebrates, for which a cold year might result in an absence of reproduction if favourable nesting sites are absent. In contrast, viviparous species might be able to better deal with poor environmental conditions by modifying their sun-basking activities (for a review see e.g. [10, 16, 23]). For instance, Pincheira-Donoso et al., , found that, in a widespread South American lizard genus (Liolaemus), the evolution of viviparity was associated with a radiation into cold climates. This argument is not relevant to birds or placental and marsupial mammals (as opposed to egg-laying monotremes) because their breeding type is constant and they are more widely distributed. However, the argument can be extended to other ectothermic vertebrates (fish and amphibians).
Body colour plays a major role in the evolution of organisms as it is involved, for instance, in thermoregulation, prey–predator interaction, behaviour, protection against UV-light or abrasive agents . Therefore, the presence of intraspecific colour polymorphism may promote the ecological success of species as previously shown for squamates [14, 15], amphibians  and owls , through their capacity to exploit a larger range of habitat types. In turn, polymorphic species might be more resilient to past and future climatic fluctuations. Our results are also in agreement with the study of Pizzatto & Dubey , which found that colour-polymorphic snake species were older than monomorphic species.
Recent studies focusing on temperate vertebrate taxa highlighted that species from the southern Hemisphere are older than those from the northern, and that differences in past climate are likely responsible for this pattern [1–3]. Interestingly, the Hemisphere of origin is not impacting the age of species in our study, which is likely due to the inclusion of additional intrinsic and extrinsic factors and to the use of phylogenetic comparative methods. Indeed, when the Hemisphere of origin is included in our analyses without accounting for the phylogenetic relationships between taxa, southern species remained significantly older than northern ones (see Additional file 3). The same is true when strictly temperate species of our dataset are included (the use of a comparative method leading to a marginally significant result; Additional file 3). These differences illustrate the importance of considering an array of factors in such analyses. In addition, the phylogenetic relationships of studied taxa might deeply influence such studies, keeping in mind that northern and southern Hemisphere temperate fauna are distantly related genetically, numerous families and genera being solely distributed in one of the two Hemispheres.
In conclusion, our findings suggest that intrinsic and extrinsic variables might contribute to the resilience of species to past environmental changes and hence are relevant in a context of climate change. Overall and independently of current threats, organisms which exhibit ancient intraspecific diversification events are more likely to deal with future modifications of their environment, those species being colour polymorphic, viviparous, or living at low latitudes.
Species sampling and traits
We gathered published phylogenetic and phylogeographic studies for vertebrate species using searches within Scopus and ISI Web of Knowledge. We defined the age of a species (in Myr) as the oldest intraspecific diversification event, which is an estimate of the amount of time since extant populations last shared a common ancestor [1–3, 7]. Species ages were obtained directly from the studies included in our dataset and hence were not estimated in the present study. In addition, only studies including a reasonable number of individuals and covering most of the distribution of the species have been included. We excluded non-monophyletic species and used the Catalogue of Life  to eliminate synonyms and ambiguities. Our final sampling contained species from 71 freshwater fish, 73 amphibians, 156 reptiles, 145 birds and 155 mammals (See Fig. 1 & Additional file 4).
We collected for each species their geographical distributions (Additional file 4), as represented by the minimum and maximum latitude taken from IUCN distributional polygon data. We calculated the mean latitude (in absolute values) and latitudinal range with ArcGIS® (ESRI 2006). We further categorized the mean latitude of each species to allocate them to Northern (between 10° and 90° N) or Southern (between 10° and 90° S) Hemispheres. Tropical species (between 10°N and 10°S) as well as species found on both Hemispheres were defined as intermediates. Latitude is considered here as a covariate that can represent several aspects of the life-history that cannot be easily associated with the more than 600 species sampled in our study.
We also collected data on life-history traits from published taxonomies and available databases. As our measure of mean adult size, we used snout-vent length for amphibian, reptile and mammal species, common size for fish species and length for bird species. Both sexes were pooled together and averaged. We also recorded the breeding type for each species and we combined ovoviviparous and viviparous species. The mean litter or clutch size was defined as the mean number of offspring for one female and we only considered the number of offspring in a single clutch/litter when multiple birthing took place in one year. We considered as well the behaviour of newborns and modified the scale of Starck & Ricklefs  into four levels to account for the larger taxonomic breadth of our study (Precocial - neonates eat by themselves and move freely; semi-precocial - neonates move freely by themselves but rely on the parents for food; semi-altricial - neonates unable to move freely by themselves, but still have the capacity to cling or act to survive; altricial - neonates are fully dependant on their parents). Finally, we defined a species as colour-polymorphic when it exhibited at least two different colour morphs within the same sex and age class. Thus, at least one of the sexes had to be polymorphic for the species to be considered polymorphic. No distinction was made between geographical colour variation and within-population colour-polymorphism. See Additional file 5 for factors retained in our final analyses.
The phylogenetic tree for the 601 species sampled was based on four mitochondrial genes (cytochrome b, COI, ND2, and ND4) and one nuclear gene (Rag1), but not all species in the dataset were represented for the five genes (TreeBase accession number). The shark Squalus acanthias was used as an outgroup. The sequences were downloaded from GenBank and aligned using MUSCLE . The best model for each gene defined using Akaike Information Criterion (AIC) as implemented in the phymltest function from the R package APE . We inferred the Maximum Likelihood tree using the BEST tree search algorithm as implemented in PhyML 2.3  and we performed 100 bootstrap replicates. We also inferred a Bayesian phylogenetic tree using MrBayes 3.2  and we allowed the model of substitution to take different parameter values for each gene present in the concatenated dataset. The MCMC was run for 50 millions generations and we sampled trees and parameter values every 1000 generations. We repeated the analyses twice to assess the convergence of the chain.
We estimated divergence times using 11 fossil calibration points obtained from the literature (Additional file 6). Given the size of the phylogenetic tree, a full estimation of the tree topology and the divergence times was not practical. We therefore decided to first use a Maximum Likelihood implementation of a relaxed molecular clock  as implemented in the function chronos from the R package APE . We also sampled 100 trees from the posterior distribution of trees obtained by MrBayes and repeated the divergence time analysis to obtain a set of calibrated trees. Second, we used BEAST v. 1.8  to estimate the divergence times using a log-normal relaxed molecular clock while fixing the tree topology to the one obtained by MrBayes. We applied the 11 fossil calibration points using log-normal priors. The means and standard deviations of these prior distributions were selected to include the fossil uncertainties in the 95 % of their density. We used a yule process for the prior on divergence times and ran the MCMC for 10 million generations. We assessed the convergence of the chain using Tracer v. 1.7.5 by looking at the trace of the log-likelihood and the ESS values for each parameter.
We used general least squares (GLS) to understand the effects of intrinsic and extrinsic factors on the age of species. We modified species ages with a Box-Cox transformation (lambda = 0.1) and we included the average latitude (in absolute value), Hemisphere, reproduction mode, newborn behavior, body size (Box-Cox transformed, lambda = −0.3) and the presence of colour polymorphism as covariables in the GLS analysis. We started with a full model that included the single factors and all pairwise interactions. We corrected the correlation structure of the GLS by considering the evolutionary relationships between species and we assumed a Brownian correlation structure based on the phylogenetic trees (lambda parameter = 1). We performed all analyses using the gls function from the R package NLME  and used backward model selection to test the factors. Model validation was performed using graphical tools (as proposed in ). We finally tested the effect of phylogenic tree construction on the best model output by repeating the GLS analysis 100 times sampling each time a new phylogenetic tree randomly from their posterior distribution.
Dubey S, Shine R. Geographic variation in the age of temperate-zone reptile and amphibian species: southern Hemisphere species are older. Biol Letters. 2011;7:96–7. doi:10.1098/rsbl.2010.0557.
Dubey S, Shine R. Are reptile and amphibian species younger in the Northern Hemisphere than the Southern Hemisphere? J Evol Biol. 2012;25:220–6. doi:10.1111/j.1420-9101.2011.02417.x.
Tingley R, Dubey S. Disparity in the timing of vertebrate diversification events between the northern and southern hemispheres. BMC Evol Biol. 2012;12:244. doi:10.1186/1471-2148-12-244.
Markgraf V, McGlone M, Hope G. Neogene paleoenvironmental and paleoclimatic change in southern temperate ecosystems – a southern perspective. Trends Ecol Evol. 1995;10:143–7. doi:10.1016/S0169-5347(00)89023-0.
Hewitt GM. The genetic legacy of the Quaternary ice ages. Nature. 2000;405:907–13. doi:10.1038/35016000.
Byrne M. Evidence for multiple refugia at different time scales during Pleistocene climatic oscillations in southern Australia inferred from phylogeography. Quaternary Science Rev. 2008;27:2576–85. doi:10.1016/j.quascirev.2008.08.032.
Weir JT, Schluter D. The latitudinal gradient in recent speciation and extinction rates of birds and mammals. Science. 2007;315:1574–6. doi:10.1126/science.1135590.
Balon EK. Dynamics of biodiversity and mechanisms of change: A plea for balanced attention to form creation and extinction. Biol Cons. 1993;66:5–8. doi:10.1016/0006-3207(93)90129-O.
Isaac JL, Vanderwal J, Johnson CN, Williams SE. Resistance and resilience: quantifying relative extinction risk in a diverse assemblage of Australian tropical rainforest vertebrates. Divers Distrib. 2009;15:280–8. doi:10.1111/j.1472-4642.2008.00531.x.
Lynch VJ. Live-birth in vipers (Viperidae) is a key innovation and adaptation to global cooling during the Cenozoic. Evolution. 2009;63:2457–65. doi:10.1111/j.1558-5646.2009.00733.x.
Gardner JL, Peters A, Kearney MR, Joseph L, Heinsohn R. Declining body size: a third universal response to warming? Trends Ecol Evol. 2011;26:285–91.
Pizzatto L, Dubey S. Colour-polymorphic snake species are older. Biol J Linn Soc. 2012;107:210–8. doi:10.1111/j.1095-8312.2012.01936.x.
Shreeve TG. Microhabitat use and hindwing phenotype in Hipparchia semele (Lepidoptera, Satyrinae): thermoregulation and background matching. Ecol Entom. 1990;15:201–13. doi:10.1111/j.1365-2311.1990.tb00801.x.
Forsman A, Aberg V. Associations of variable coloration with niche breadth and conservation status among Australian reptiles. Ecology. 2008;89:1201–7. doi:10.1890/07-1670.1.
Forsman A, Åberg V. Variable coloration is associated with more northerly geographic range limits and larger range sizes in North American lizards and snakes. Evol Ecol Res. 2008;10:1025–36.
Shine R. Does viviparity evolve in cold climate reptiles because pregnant females maintain stable (not high) body temperatures? Evolution. 2004;58:1809–18. ISSN: 0014–3820.
Pincheira-Donoso D, Tregenza T, Witt MJ, Hodgson DJ. The evolution of viviparity opens opportunities for lizard radiation but drives it into a climatic cul-de-sac. Global Ecol Biogeogr. 2013;22:857–67. doi:10.1111/geb.12052.
Damuth J. Population density and body size in mammals. Nature. 1981;290:699–700. doi:10.1038/290699a0.
O'leary MA et al. The Placental Mammal Ancestor and the Post-K-Pg Radiation of Placentals. Science. 2013;339:662–7. doi:10.1126/science.1229237.
Mulcahy DG et al. Estimating divergence dates and evaluating dating methods using phylogenomic and mitochondrial data in squamate reptiles. Mol Phylogenet Evol. 2012;65:974–91. doi:10.1016/j.ympev.2012.08.018.
Zhang P, Zhou H, Chen YQ, Liu YF, Qu LH. Mitogenomic Perspectives on the Origin and Phylogeny of Living Amphibians. Syst Biol. 2005;54:391–400. doi:10.1080/10635150590945278.
Zhang GJ et al. Comparative genomics reveals insights into avian genome evolution and adaptation. Science. 2014;346:1311–20. doi:10.1126/science.1251385.
Shine R. Life-history evolution in reptiles. Annu Rev Ecol Evol, S. 2005;36:23–46.
Dubey S, Roulin A. Evolutionary and biomedical consequences of internal melanins. Pigment Cell Melanoma R. 2014;27:327–38. doi:10.1111/pcmr.12231.
Forsman A, Hagman M. Association of coloration mode with population declines and endangerment in Australian frogs. Cons Biol. 2009;23:1535–43. doi:10.1111/j.1523-1739.2009.01244.x.
Roulin A, Burri R, Antoniazza S. Owl melanin-based plumage redness is more frequent near than away from the equator: implications on the effect of climate change on biodiversity. Biol J Linn Soc. 2011;102:573–82. doi:10.1111/j.1095-8312.2010.01614.x.
Roskov Y, Abucay L, Orrell T, Nicolson D, Flann C, Bailly N, Kirk P, Bourgoin T, DeWalt RE, Decock W, De Wever A. Species 2000 & ITIS Catalogue of Life, 2016 Annual Checklist. 2016. Digital resource at www.catalogueoflife.org/annual-checklist/2016. Species 2000: Naturalis, Leiden, the Netherlands. ISSN 2405-884X.
Starck JM, Ricjlefs RE. Avian Growth and Development: Evolution Within the Altricial-precocial Spectrum. Oxford University Press. 1998.
Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–7. doi:10.1093/nar/gkh340.
Paradis E, Claude J, Strimmer K. APE: Analyses of Phylogenetics and Evolution in R language. Bioinformatics. 2004;20:289–90. doi:10.1093/bioinformatics/btg412.
Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52:696–704. doi:10.1080/10635150390235520.
Ronquist F et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42. doi:10.1093/sysbio/sys029.
Paradis E. Molecular dating of phylogenies by likelihood methods: A comparison of models and a new information criterion. Mol Phylogenet Evol. 2013;67:436–44. doi:10.1016/j.ympev.2013.02.008.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73. doi:10.1093/molbev/mss075.
Pinheiro J, Bates D, DebRoy S, Sarkar D, the R Development Core Team. Nlme: Linear and Nonlinear Mixed Effects Models. R package version. 2013;3:1–111.
Zuur AF, Ieno EN, Walker N, Saveliev AA, Smith GM. In: Gail M, Krickeberg K, Samet JM, Tsiatis A, Wong W, editors. Mixed Effects Models and Extensions in Ecology with R. New York: Springer; 2009.
We thank Richard Shine for helpful comments. This work was financially supported by the Swiss National Science Foundation (SNSF; n° PZ00P3_136649 and PZ00P3_157244 to SD).
The authors declare that they have no competing interests.
SD, NS, LC, JS contributed with the conceptual development of the work and the writing of the manuscript and carried out the analyses. All authors read and approved the final version of the manuscript.
Cattin and Schuerch are co-first authors
Salamin and Dubey are co-last authors
Final phylogenetic GLS model. (DOCX 13 kb)
Age variation within factors. (DOCX 13 kb)
The effect of hemisphere on age of species, considering or not the phylogenic correction or/and species from the tropics (below 23° latitude in absolute value). (DOCX 18 kb)
Taxa included in our final dataset. (XLSX 72 kb)
Correlation between intrinsic and extrinsic fatcors relative to the species sample. Covariables retained in the model are represented in bold. (DOCX 15 kb)
Calibration points used in the phylogenetic analyses. (DOCX 23 kb)