Oligocene niche shift, Miocene diversification – cold tolerance and accelerated speciation rates in the St. John’s Worts (Hypericum, Hypericaceae)
© Nürk et al.; licensee BioMed Central. 2015
Received: 19 November 2014
Accepted: 22 April 2015
Published: 6 May 2015
Our aim is to understand the evolution of species-rich plant groups that shifted from tropical into cold/temperate biomes. It is well known that climate affects evolutionary processes, such as how fast species diversify, species range shifts, and species distributions. Many plant lineages may have gone extinct in the Northern Hemisphere due to Late Eocene climate cooling, while some tropical lineages may have adapted to temperate conditions and radiated; the hyper-diverse and geographically widespread genus Hypericum is one of these.
To investigate the effect of macroecological niche shifts on evolutionary success we combine historical biogeography with analyses of diversification dynamics and climatic niche shifts in a phylogenetic framework. Hypericum evolved cold tolerance c. 30 million years ago, and successfully colonized all ice-free continents, where today ~500 species exist. The other members of Hypericaceae stayed in their tropical habitats and evolved into ~120 species. We identified a 15–20 million year lag between the initial change in temperature preference in Hypericum and subsequent diversification rate shifts in the Miocene.
Contrary to the dramatic niche shift early in the evolution of Hypericum most extant species occur in temperate climates including high elevations in the tropics. These cold/temperate niches are a distinctive characteristic of Hypericum. We conclude that the initial release from an evolutionary constraint (from tropical to temperate climates) is an important novelty in Hypericum. However, the initial shift in the adaptive landscape into colder climates appears to be a precondition, and may not be directly related to increased diversification rates. Instead, subsequent events of mountain formation and further climate cooling may better explain distribution patterns and species-richness in Hypericum. These findings exemplify important macroevolutionary patterns of plant diversification during large-scale global climate change.
At the onset of evolutionary theory it was observed that “species of the same genus have usually […] some similarity in habits and constitution” (, p 76). That is, closely related species or lineages are expected to be more similar in ecology than to more distantly related taxa because of their more recent common ancestry [2,3]. Accordingly, it has been observed that under changing environmental conditions organisms tend to retain their ancestral ecological characteristics rather than evolving into a new ecological niche [4-6].
Large-scale global climate change dramatically alters the distribution of major biomes , and thus the ecological niches available to entire taxonomic groups [8,9]. During the Early Eocene Climatic Optimum, c. 55 million years (Ma) ago, tropical biomes dominated the Earth’s surface even at high latitudes . Within the last 50 Ma, however, the world has experienced a transition: a fluctuating but overall decrease in mean temperatures [7,11,12] that shifted the distribution of frost-intolerant plants towards the equatorial zones. The flowering plant families Araceae , Chloranthaceae , and Malpighiaceae  are well-studied examples of formerly more widely distributed lineages that are currently mostly restricted to the tropics.
Only some lineages of flowering plants have managed the transition from tropical to temperate climates [16,17], despite presumably having had ample opportunities to do so with the Neogene expansion of temperate habitats [18,19]. Adaptation to cold is supposed to involve complex reorganizations of the genome and physiology , implying that the evolution of tolerances to temperate climates with highly seasonal conditions may pose particular problems [21,22], especially for warm-adapted plant taxa confronted with cooling climate . In contrast, all extant plant lineages that do occur in temperate areas underwent adaptation to cooler climates at some point [23,24]. Also, in an area that undergoes environmental changes and that lacks suitable migration routes, only such resident lineages will survive that can develop the relevant traits necessary to persist [25-27]. In a recent study of angiosperms, Zanne et al. suggest that “many of these solutions were probably acquired before their foray into the cold” . On the other hand, it has been suggested that in the absence of obvious key innovations, climate cooling may have acted as an important driver of diversification, e.g., in the hyper-diverse genus Carex .
Hypericum has a nearly cosmopolitan distribution (Figure 1) with a primary center of species-richness in the temperate regions of Eurasia [35,36]. In temperate regions, Hypericum species are mostly native to low- and mid-elevation areas, while in the tropics they are almost always confined to high-elevation mountains, such as the South American Andes or mountains in topical Africa . The divergence time of the Hypericum crown node has been estimated to the Upper Eocene, c. 35 Ma ago , based on a single (seed-) fossil calibration. The same study reconstructed the Western Palearctic as the ancestral area for Hypericum. However, Sánchez Meseguer et al.  employed a Bayesian approach [39,40] that limits ancestral areas to only single states in the reconstruction (i.e. it is not possible to infer the occurrence of ancestral populations in multiple areas). For a genus like Hypericum that has several species that span very wide geographic ranges, this assumption is likely inadequate.
Within an otherwise pantropic clade  Hypericum is a cold-adapted but species-rich lineage. Because the clade is of worldwide cold/temperate distribution and evolved during the Paleaogene to Neogene (23.3 Ma ) climate transition, Hypericum is ideal for studying the effects of broad-scale climate change on plant distribution and diversification. Furthermore, the evolution of morphological characters in the genus has been studied at length [36-38,42], providing a foundation to further investigate evolutionary patters and processes in Hypericum.
In this study, we test the hypotheses that (i) Hypericum originated in the Western Palearctic prior to the Oligocene, as suggested by Meseguer et al. , that (ii) the occupation of temperate environments is derived in the Hypericaceae and a distinctive characteristic of Hypericum, and that (iii) it has stimulated the radiation of this lineage globally. To do so, we estimate a species phylogeny based on nuclear and chloroplast sequence variation. We calibrate our time-tree using six fossils and compare the effect on node ages of the specific assignment of the oldest known fossil in Hypericum. To evaluate historical biogeography, we employ parametric models that allow for the incorporation of paleogeographic information. Applying recently developed Bayesian approaches, we estimate the magnitude and placement of climatic niche shifts, and investigate the impact of the transition from the tropical into the cold/temperate climate niche by assessing the placement of shifts in diversification rates.
Maximum likelihood (ML) topologies were similar for both chloroplast (petD + trnL–trnF) and nuclear (ITS) data. Discordance between nuclear and chloroplast trees was present only in two places, although without strong support (Additional file 1: Phylogenetic inference). Thus, we concatenated the chloroplast and nuclear sequence data into a combined data set, which contained 24.1% missing data (0.1% missing in ITS, 21.2% in petD, and 59.4% in trnL–trnF).
Divergence times depend on fossil assignment
Results of divergence time estimation using different fossil assignments
Bayesian relaxed clock – crown ages
‘no missing’ data
‘no missing’ data
Brathys s.l. + Myriandra
core Hypericum + Ascyreia s.l.
Summary statistics detailing node support, age estimation, diversification rates, and niche shifts
Bioclimatic niche shifts*
core Hypericum – Brathys s.l.
core Hypericum + Ascyreia s.l.
Ancestral area estimation is equivocal at the Hypericum crown node
Ancestral areas were optimized over both divergence time estimations taking phylogenetic uncertainty into account by analyzing a posterior subset of 1000 trees. Additionally, two models were compared which differed in their dispersal/extinction probabilities both between areas and over time; a stratified M1 model that accounts for varying connectivity of areas over time (Additional file 1: Figure S3), and a unconstrained M2 model with equal probabilities of movement between areas at any time (i.e. equal dispersal/extinction probabilities). The ancestral areas estimated per node were highly congruent in all analyses, except for four nodes for which the ancestral states differed (Additional file 1: Table S2). Two of these affected nodes are located at basal dichotomies, one at the Hypericum crown node (Figure 2), and one at the ‘Ascyreia s.l.’ crown node. In all four cases of incongruence, and independent of the used age estimation, the evidence ratio  is generally higher in the stratified model when compared to unconstrained model (Additional file 1: Table S2).
Niche shifts during the evolution of the Hypericum stem lineage in the Oligocene
We used the first principal components (PC1) scores obtained from a phylogenetic PCA as an approximation of the climatic niche. The PC1 scores were dominated by the thermal bioclim variables (bio1–11; with the most influence from mean annual temperature [bio1]), and account for c. 34% of the variation within the data. Bayesian fitting of multi-optima OU models to the approximation of the climatic niche space estimated the placement and magnitude of two shifts. Both shifts were estimated to have occurred in the Oligocene, the first shift during the evolution of the stem lineage of Hypericum, and the second shift one divergence event latter at the stem lineage of the MRCA of ‘core Hypericum’ – ‘Triadenum (Table 2; for a comparison of results produced under the two age estimations see Additional file 1: Figure S5, Table S3). Both shift magnitudes were negative, indicating a niche shift into colder climates during the early evolution of Hypericum (Table 2 and Additional file 1: Table S3). The phylogenetic half-life (ln(2)/α = 0.009 Ma; based on a rate of adaptation of α = 1.525, and a total tree length of 52.31 Ma) signifies that the movement to the primary climatic optimum in Hypericum was rapid and resembles an OU process rather than Brownian motion .
Diversification rate shifts in Hypericum in the Miocene
In the Hypericaceae, two shifts in diversification dynamics were detected within Hypericum with strong support: (i) at the stem lineage of ‘core Hypericum + Ascyreia s.l.’, and (ii) with the divergence of Central- and South American ‘Brathys s.l.’ from its relatives (Figure 3, Table 2). In both dating approaches, the same clades are identified to show diversification rate shifts comparing the Bayesian credible sets of distinct shift configurations (Additional file 1: Figure S4). In the analysis using the dating approach A, the first increase in diversification rates is inferred to have occurred at about 13.93 Ma with a marginal probability of 0.93, and the second at about 12.54 Ma with a marginal probability of 0.64, resulting in a rate of 1.04 and 1.06 species per million year, respectively (for a comparison of results see Additional file 1: Figure S4, Table S3).
Historical biogeography of Hypericaceae
There is strong evidence that the MRCA of Hypericaceae occurred in Africa in the Eocene (Figure 2), a result also recovered in a previous study . The MRCAs of the tropical taxa within the Hypericaceae (Cratoxyleae and Vismieae) are revealed to have occurred in Africa + Southeastern Asia and Africa + South America, respectively. During the Early Eocene the higher thermal maximum allowed megathermal organisms, such as the MRCA of Hypericaceae, to disperse throughout the Northern Hemisphere, which explains the vicariance of many tropical lineages [10,19,48].
At the onset of worldwide climate cooling, c. 40 Ma the Hypericum stem lineage split from its tropical relatives. The ancestral area estimated for the MRCA of Hypericum is in the Northern Hemisphere, either solely in the Western Palaearctic, or more widely distributed between the Western Palaearctic and the Nearctic (Figure 2, Additional file 1: Table S2). The Nearctic has been connected to the Western Palearctic through islands that could have easily functioned as stepping-stones facilitating dispersal over the North Atlantic [49-51]. Seeds of Hypericum are tiny, ca. 1.5 to 2.5 mm long, and are easily dispersed by birds or strong winds , promoting long-distance dispersal. However, an ancestral distribution from the Western Palaearctic to the Nearctic is somewhat surprising given the fact that the oldest known fossil is from Siberia and SW China . Differential extinction may explain the partial absence of reconstructed ancestral populations in the areas of the oldest fossil record . Hence, the assumption of a widespread MRCA of Hypericum from the Palaearctic to the Nearctic connected via Beringia [19,53] would be plausible as well. On the other hand, it is likely that the fossil record does not accurately represent actual ancestral distributions, and thus, may be misleading for ancestral area estimations. Regardless, Hypericum is of Northern Hemisphere origin, likely around the late Tethys estuaries in the Palaearctic , perhaps as part of a deciduous mixed-mesophytic forest . The ecology of basal lineages within Hypericum, however, is diverse reaching from dry rocky Mediterranean to shallow aquatic habitats, and montane cloud forests. Thus, extensive differential extinction of intermediate ancestral populations (ecologically and geographically) and/or intercontinental long-distance colonization early in the evolutionary history is needed to explain biogeographic patterns in Hypericum.
Cold adaptation in the Hypericum stem lineage but later diversification
Initially, no speciation seems to have taken place in Hypericum – or extinction may have erased the evidence of earlier divergence in the genus leading to a stem lineage of c. 15 million years (c. 40–25 Ma). However, the genesis of modern species diversity in Hypericum traces to c. 25 Ma (Figures 2 and Additional file 1: Figure S2), after the Oligocene climate cooling caused a substantial decrease in mean temperatures . This worldwide cooling, which led to a period of rather constant cold temperatures , initiated the expansion of temperate habitats in the Northern Hemisphere, replacing the more tropical vegetation dominant during the Eocene . With the exception of Hypericum, the remaining lineages of Hypericaceae stayed in tropical climates, i.e. their distributions experienced a restriction with the retreat of tropical areas towards equatorial zones after the Eocene Thermal Maximum . In contrast, Hypericum adapted to colder climates c. 30 Ma (Figure 3) and consequently dispersed and diversified in temperate habitats. The adaptation towards the new climatic optimum (estimated by the phylogenetic half-life ) was rapid, meaning that Hypericum species were well adapted to the cold/temperate niche . Hence, following this initial niche shift, Hypericum never completely left the colder climates. Even in equatorial areas, Hypericum is only found at high elevations with a cool climate, such as the South American Andes and the high mountains of Africa.
Two shifts in diversification rates were detected with strong support in Hypericum (Figure 3). The first is an increase in speciation rates that coincides with a sudden decrease of temperature during the Middle Miocene Climate Transition at about 14 Ma . The second increase in speciation rates was likely the result of dispersal into the South American continent and may have been triggered by the orogeny of the Andes . That is, the adaptation of Hypericum to colder climates first evolved during the Oligocene but diversification rates increased c. 13–17 million years later during the Miocene, in nested clades within Hypericum. The evolution of a shrubby habit at the Hypericum stem lineage  might well be a life history trait that facilitated the initial adaptation to cold , but no obvious intrinsic (morphological or physiological) traits that could be interpreted as key innovations were identified for these rapidly diversifying clades [37,38]. However, our results suggest that cold tolerance is likely an important initial adaptation that was exploited when new biogeographic opportunities were presented (i.e. Andean orogeny). Only after a long lag phase when global temperatures dramatically decreased following the Mid-Miocene Thermal Optimum , did changing environments and expanding temperate regions (including tropical high mountains) allow Hypericum to spread and diversify. Because c. 42% of extant Hypericum species are found in montane biomes (c. 72% of this diversity in tropical high mountains), we postulate that the onset of extensive mountain formation in Eurasia and the Americas  is likely to have contributed to both detected rate shifts.
After divergence from its tropical relatives, adaptation towards colder climates c. 30 Ma ago enabled Hypericum to stay in the Northern Hemisphere, while its tropical relatives experienced habitat restriction towards equatorial zones. This niche shift offered dispersal and possible diversification opportunities in the expanding temperate areas in the Northern Hemisphere during the Oligocene. After the last thermal maximum in the Miocene, massive mountain formations and further climate cooling may have stimulated the radiation of this lineage globally. As a consequence, cold-adapted Hypericum contains 80% of the species present in the family. Despite its worldwide distribution and tropical ancestry, even the species growing in the tropics retain their temperate climate niche by growing exclusively in cool climates of higher elevation habitats.
Higher species richness in temperate climates is rather atypical and contrary to the usual pattern observed in plants where tropical groups generally tend to have more species than temperate relatives [6,23]. However, our findings mirror patterns described in sedges (Carex) , and to varying extents in buttercups (Ranunculaceae) , grasses (Poaceae) [59-62], and heaths (Ericales) .
We have demonstrated that a pronounced lag phase is present between the initial niche shift and diversification in new habitats. Therefore, we conclude that the disproportionate species numbers of Hypericum in comparison to its tropical relatives are not only a result of initial cold-climate adaptation. Our analyses indicate a relatively late increase in speciation rate, i.e. with the onset of further cooling and especially mountain formation during the Upper Miocene (and Pliocene). Likewise, Carex, the most diverse non-tropical sedge lineage (Cyperaceae)  has a stem lineage of c. 20 million years, with most extant lineages having diversified subsequent to the Oligocene in montane biomes. Thus, we emphasize that both (i) a niche shift following the Eocene Thermal Optimum as a precondition for the presence of cold-adapted lineages in temperate regions, and (ii) an extrinsic trait (perhaps in addition to lineage-specific intrinsic traits), i.e. the availability of emerging temperate and/or mountain habitats, are key events potentially triggering diversification rates in these cold-adapted, species-rich plant groups.
Taxon sampling and species richness
Plant material from herbarium specimens or silica dried samples were chosen for 100 species (and 3 subspecies) representing the distribution range of all major lineages present in Hypericaceae [30,37] (see Additional file 1: Voucher). Following Xi et al. , Garcinia xanthochymus Hook.f. ex T.Anderson (Clusiaceae Lindl.) was chosen as outgroup in the phylogenetic analyses.
The monograph of Hypericum [31,32,42,64-72] was used for data on species richness and distribution (for Hypericum sensu Robson 2012). Stevens  lists information for the remaining taxa of Hypericeae Choisy (Triadenum Raf., Thornea Breedlove & E.M.McClint., and Lianthus N.Robson; included in Hypericum in Ruhfel et al. ), as well as the tropical genera of the family Hypericaceae. Xi et al.  provides information on distributions for Clusiaceae and Calophyllaceae J.Agardh (Additional file 1: Voucher). Total species numbers with evidence from taxonomy , morphological cladistics , and molecular phylogenetic analyses [37,38] were used to assign species richnesses to the major clades defined in the diversification rate analyses.
Molecular marker and sequencing
We sequenced two fragments from the chloroplast genome, namely petD (including the petB–petD intergenic spacer, the petD-5′-exon, and the petD intron) and trnL–trnF (including the trnLUAA intron and the intergenic spacer between the trnLUAA 3′ exon and trnFGAA gene), and the nuclear rDNA internal transcribed spacer region (including ITS-1, 5.8S rDNA, and ITS-2). Extraction of DNA was done according to Nürk et al. . In the case of well preserved herbarium material or silica dried samples, entire regions were amplified using the primers ITS-A(F) and ITS-B(R) , PIpetB1411F and PIpetD738R , c(F) and f(R) designed by Taberlet  for the trnL–trnF region. In the case of degraded herbarium materials, ITS-1 and ITS-2 were amplified separately using in addition two internal primers, ITS-C(R) and ITS-D(F) binding in the conserved 5.8S rDNA . Similarly for petD, using the two internal primers SALpetD599F and OpetD897R designed by Korotkova et al. . PCR amplification of ITS was performed as described in Nürk et al. . PCR reaction mixes for petD and trnL–trnF were chosen according to Nürk et al. , but without adding MgCl2, and PCR profiles consisted of an initial denaturation at 96°C for 1.5 min, followed by 35 cycles of 95°C for 30s, 50°C for 60s, 73°C for 90s and a final step at 72°C for 10 min. Primer combinations as described above for poorly preserved samples were used for cycle sequencing. DNA sequencing was done by Eurofins MWG Operon (Ebersberg, Germany). All newly generated sequences have been submitted to the to the EMBL nucleotide database (Accession No. LK871650–LK871782).
Sequences were assembled and edited with Geneious v5.4 , aligned using the automatic selection of an appropriate strategy in Mafft v6.903b [78,79] and manually adjusted using PhyDE v0.996 (available online: http://www.phyde.de). In order to remove poorly aligned or length variable data partitions the alignments were subjected to Gblocks 0.91b sever  with the ‘less stringent’ options selected.
Phylogenetic analyses were performed under maximum likelihood (ML)  and Bayesian inference (BI)  to reveal confidence limits of the data. ML analyses were performed with the RAxML GUI v1.1 [83,84] and BI in MrBayes 3.2.2 . To test for discordance we analyzed the nuclear (ITS) and chloroplast data partitions (petD, trnL–trnF) separately by ML search under the GTRCAT model of sequence evolution. Clade support was evaluated with 1000 rapid bootstrap replicates .
The combined data set (ITS + petD + trnL–trnF) was analyzed under ML with the partitions defined and the settings chosen as described for the ML analysis above. For BI we started 4 simultaneous runs, each with 4 chains, set to run 108 cycles, with sampling every 104 cycle, setting temperature to 0.01, and with the appropriate model of sequence evolution specified per partition: GTR + I + Γ for ITS and petD and HKY + I + Γ for trnL–trnF; selected in MrModeltest  under the Akaike Information Criterion (AIC) . We used the ML tree as a starting tree, but introduced random perturbations into it to enable detection of possible convergence problems (using the command “mcmcp nperts = 5”). A ‘corrected’ exponential prior on a branch length of 1/λ = 0.1 [“prset brlenspr = Unconstrained:Exp(100)”] was specified . Convergence of parameter estimates was monitored using Tracer v1.5 . After discarding 25% of the sampled trees as burnin, posterior probabilities were calculated on the BI stationary sample. Trees and alignments are available at TreeBASE study number 16298.
Divergence time estimation and fossil assignment
The likelihood ratio test  conducted on the BI consensus tree in PAUP*  rejected a global molecular clock (P < 0.05) for the combined data set. Therefore, divergence times were estimated under a relaxed molecular clock employing the uncorrelated lognormal model  in BEAST v1.7 . Eight external time-constraints were imposed for calibration, comprising six fossils [52,95-97] and two secondary calibrations  (for details see Additional file 1: Age estimation, calibration). Fossil calibrations were constrained by hard minimum bounds and secondary calibrations by lognormal distributions to incorporate the uncertainty reported in the original study . Two approaches were designed differing only in the assignment of the seed fossil Hypericum antiguum Balueva & V.P. Nikitin : (i) to the stem node of Hypericum in analysis A, and (ii) to the crown node of Hypericum in analysis B (other calibrations remained unchanged in the two analyses; for a discussion of fossil assignment see Additional file 1: Age estimation, calibration).
A birth and death model of speciation considering incomplete species sampling  was set as tree prior. Both divergence time estimations (A and B) were started with two independent Monte Carlo Markov Chain (MCMC) runs, each set to run 108 cycles with sampling every 104 cycles. The substitution and clock models were not linked between the partitions. The ML tree was used as starting tree. To ensure that the prior branching times of the starting tree fulfilled the constraints imposed by the calibration priors we transformed branch length into absolute time using penalized likelihood  with the chronopl command in the R  package ape . Convergence of parameter estimates was monitored using Tracer . The resulting trees were combined in LogCombiner with a burnin of 50%. On the remaining 10,002 trees means and confidence intervals were calculated in TreeAnnotator  to obtain the final consensus tree, the ultrametric time calibrated maximum clade credibility (MCC) chronogram that has 95% of the highest posterior density (HPD). Because missing data can have deleterious effects on analyses that depend on branch length , we tested the effect on age estimates of missing sequences in our data set. We repeated both analyses (A and B) using a data set that did not contain missing data (‘no missing’ data) and that had therefore a reduced species sampling containing only 73 accessions (xml input files are available at the Dryad repository ).
Ancestral area estimation
Historical biogeography was analyzed by classifying the species to be distributed within six biogeographic regions, following Brummit et al.  for area subdivision. The region were defined to reflect general biogeographic entities, and to be meaningful for the study group: (A) Afrotropical [central Africa, the southern Arabian peninsula, Madagascar, and the West Indian Ocean islands], (WP) western Palearctic, (EP) eastern Palearctic, (IP) Indo-Pacific [SE tropical Asia, Australasia, and the Pacific], (NA) North America [Nearctic], (SA) South America [Neotropic].
We employed a parametric likelihood approach that uses the DEC model  implemented in Lagrange . Two models were designed, differing in dispersal probabilities between areas and time to take into account the impact of dispersal/extinction probabilities. The first model (M1) was stratified, i.e. incorporating paleogeographic information on area connectivity, e.g., the existence of land bridges, by varying the dispersal probabilities between areas and over time (Additional file 1: historical biogeography; Additional file 1: Figure S3). The second model (M2) was unconstrained, assuming equal dispersal probabilities between areas over time. Ancestral areas were optimized under both models over 1000 post-burnin posterior trees randomly chosen  and generated by divergence time estimation A and, in a second analysis, by estimation B. Composite Akaike weights  were used to summarize the likelihood of range optimizations. Evidence ratios (the ratio of Akaike weights w i /w j ) were used to evaluate which scenario is the most favored.
Bioclimatic niche analysis
We used the WorldClim  global climate layers considering all 19 Bioclim variables as an approximation of the climatic niche of the species. Collection sites were taken from the voucher specimens, or manually geo-referenced (if not documented on the specimen) while comparing voucher information with collection sites per species recorded in the Global Biodiversity Information Facility (GBIF). The bioclimatic variables for each collection locality were extracted using the R package raster . To identify significant changes in the adaptive landscape (i.e. the climatic niche space) within the study group we employed a reversible-jump Bayesian method of fitting multi-optima Ornstein-Uhlenbeck (OU) models to the bioclim variables . To do so, we first transformed the bioclim variables using a phylogenetic principal components analysis (PCA)  in the R package phytools  to obtain principal components (PC) scores per species while correcting for phylogenetic history (i.e. the correlation of independent contrasts ). Then, we analyzed the PC1 scores using in the R package bayou  with a standard error of 0.5 to estimate the placement and magnitude of shifts in the climatic niche directly from the data . We allowed only one shift per branch and assigned an equal probability of each branch having a shift. We placed a corrected Poisson distribution as prior on adaptive optima and a probability density for a half-Cauchy distribution on the number of shifts between adaptive regimes [using make.prior((tree, dists = list(dalpha = “dhalfcauchy”, dsig2 = “dhalfcauchy”, dsb = “dsb”, dk = “cdpois”, dtheta = “dnorm”), param = list(dalpha = list(scale = 1), dsig2 = list(scale = 1), dk = list(lambda = 15, kmax = 200), dsb = list(bmax = 1, prob = 1), dtheta = list(mean = mean(dat), sd = 2*sd(dat))))], and run MCMC for 106 cycles. To verify that MCMC analyses converged to the same posterior distribution, we applied the Gelman diagnostic  in the R package coda . After discarding the first 30% of samples per run as burnin parameter estimates of the two runs were combined and summarized using the Lposterior command in bayou  to obtain a posterior of shift locations and magnitudes, the rate of adaptation (α), and phylogenetic half-life (the amount of time it takes for the expected trait value to get halfway to the phenotypic optimum, defined as ln(2)/α units of time ).
Diversification rate analysis
We used the Bayesian approach for studying patterns of rate variation through time and among lineages implemented in BAMM (Bayesian analysis of macroevolutionary mixtures) [116,117]. The method aims at detecting and quantifying heterogeneity in evolutionary rates by using reversible-jump Markov Chain Monte Carlo to detect subclades that share common parameters of speciation and extinction. BAMM identifies sets of shifts (i.e. configurations) that are sampled together and allows one to compute relative and marginal probability of those configurations (given in the 95% credible set of distinct shift configurations). We analyzed both the MCC consensus tree produced by divergence time estimation A, and B, respectively (outgroup taxa were removed prior to the analyses to increase statistical power). To do so, we assigned species richness (i.e. sampling fraction) to ten well-supported clades to account for incomplete sampling, and ran MCMC over 108 cycles with default settings. Convergence of parameter estimates was evaluated by means of effective samples size (ESS) diagnostics using the R package coda  after discarding 10% of samples as burnin. Results were summarized and visualized using the R package BAMMtools .
We thank the herbaria ANDES, B, BM, GH, HEID, K, KYO, PE, TI, Z for allowing us to use their collections. We are grateful to Michael Pirie for discussions on fossil constrains, and Josef Uyeda, Matt Pennell and Luke Harmon for assistance regarding bayou. Tony Verboom and two anonymous reviewers are thanked for valuable comments on the manuscript. This work was supported by the Deutsche Forschungsgemeinschaft (DFG; grant number BL462/8, NU292/1).
- Darwin C. On the Origin of Species. London: J. Murray; 1859.Google Scholar
- Harvey PH, Pagel MD. The Comparative Method in Evolutionary Biology. Oxford: Oxford University Press; 1991.Google Scholar
- Moen D, Morlon H. Why does diversification slow down? Trends Ecol Evol. 2014;29:190–7.PubMedGoogle Scholar
- Peterson A, Soberón J, Sánchez-Cordero V. Conservatism of ecological niches in evolutionary time. Science. 1999;285:1265–7.PubMedGoogle Scholar
- Wiens JJ, Graham CH. Niche conservatism: integrating evolution, ecology, and conservation biology. Annu Rev Ecol Evol Syst. 2005;36:519–39.Google Scholar
- Jablonski D, Roy K, Valentine JW. Out of the tropics: evolutionary dynamics of the latitudinal diversity gradient. Science. 2006;314:102–6.PubMedGoogle Scholar
- Graham A. A Natural History of the New World. London: The University of Chicago Press; 2011.Google Scholar
- Donoghue MJ, Bell CD, Li J. Phylogenetic patterns in Northern Hemisphere plant geography. Int J Plant Sci. 2001;162:41–52.Google Scholar
- Ohlemüller R. Running out of climate space. Science. 2011;334:613–4.PubMedGoogle Scholar
- Morley RJ. Interplate dispersal paths for megathermal angiosperms. Perspect Plant Ecol Evol Syst. 2003;6:5–20.Google Scholar
- Zachos J, Dickens G, Zeebe R. An early Cenozoic perspective on greenhouse warming and carbon-cycle dynamics. Nature. 2008;451:279–83.PubMedGoogle Scholar
- Zachos J, Pagani M, Sloan L, Thomas E, Billups K. Trends, rhythms, and aberrations in global climate 65 Ma to present. Science. 2001;292:686–93.PubMedGoogle Scholar
- Nauheimer L, Metzler D, Renner SS. Global history of the ancient monocot family Araceae inferred with models accounting for past continental positions and previous ranges based on fossils. New Phytol. 2012;195:938–50.PubMedGoogle Scholar
- Antonelli AA, Sanmartín I. Mass extinction, gradual cooling, or rapid radiation? Reconstructing the spatiotemporal evolution of the ancient angiosperm genus Hedyosmum (Chloranthaceae) using empirical and simulated approaches. Syst Biol. 2011;60:596–615.PubMedGoogle Scholar
- Manchester SR, Wilde V, Collinson ME. Fossil cashew nuts from the Eocene of Europe: Biogeographic links between Africa and South America. Int J Plant Sci. 2007;168:1199–206.Google Scholar
- Judd W, Sanders R, Donoghue MJ. Angiosperm family pairs: preliminary phylogenetic analyses. Harvard Pap Bot. 1994;5:1–51.Google Scholar
- Crisp MD, Arroyo MTK, Cook LG, Gandolfo MA, Jordan GJ, McGlone MS, et al. Phylogenetic biome conservatism on a global scale. Nature. 2009;458:754–6.PubMedGoogle Scholar
- Ricklefs R, Renner SS. Species richness within families of flowering plants. Evolution. 1994;48:1619–36.Google Scholar
- Donoghue MJ. A phylogenetic perspective on the distribution of plant diversity. Proc Natl Acad Sci U S A. 2008;105:11549–55.PubMed CentralPubMedGoogle Scholar
- Sakai A, Larcher W. Frost Survival of Plants: Responses and Adaptation to Freezing Stress. Ecological Studies 62. Berlin: Springer Verlag; 1987.Google Scholar
- Wiens JJ, Donoghue MJ. Historical biogeography, ecology and species richness. Trends Ecol Evol. 2004;19:639–44.PubMedGoogle Scholar
- Wiens JJ, Kuczynski CA, Townsend T, Reeder TW, Mulcahy DG, Sites JW. Combining phylogenomics and fossils in higher-level Squamate reptile phylogeny: molecular data change the placement of fossil taxa. Syst Biol. 2010;59:674–88.PubMedGoogle Scholar
- Jansson R, Rodriguez-Castaneda G, Harding LE. What can multiple phylogenies say about the latitudinal diversity gradient? A new look at the tropical conservatism, out-of-the-tropics and diversification rate hypotheses. Evolution. 2013;67:1741–55.PubMedGoogle Scholar
- Zanne AE, Tank DC, Cornwell WK, Eastman JM, Smith SA, FitzJohn RG, et al. Three keys to the radiation of angiosperms into freezing environments. Nature. 2014;506:89–92.PubMedGoogle Scholar
- Engler HGA. Versuch einer Entwicklungsgeschichte der extratropischen Florengebiete der nördlichen Hemisphäre. Leipzig: Wilhelm Engelmann; 1879.Google Scholar
- Baldwin BG, Sanderson MJ. Age and rate of diversification of the Hawaiian silversword alliance (Compositae). Proc Natl Acad Sci U S A. 1998;95:9402–6.PubMed CentralPubMedGoogle Scholar
- Mannion PD, Upchurch P, Benson RBJ, Goswami A. The latitudinal biodiversity gradient through deep time. Trends Ecol Evol. 2013;29:42–50.PubMedGoogle Scholar
- Escudero M, Hipp AL, Waterway MJ, Valente LM. Diversification rates and chromosome evolution in the most diverse angiosperm genus of the temperate zone (Carex, Cyperaceae). Mol Phylogenet Evol. 2012;63:650–5.PubMedGoogle Scholar
- Wurdack KJ, Davis CC. Malpighiales phylogenetics: gaining ground on one of the most recalcitrant clades in the angiosperm tree of life. Am J Bot. 2009;96:1551–70.PubMedGoogle Scholar
- Xi Z, Ruhfel BR, Schaefer H, Amorim AM, Sugumaran M, Wurdack KJ, et al. Phylogenomics and a posteriori data partitioning resolve the Cretaceous angiosperm radiation Malpighiales. Proc Natl Acad Sci U S A. 2012;109:17519–24.PubMed CentralPubMedGoogle Scholar
- Robson NKB. Studies in the genus Hypericum L. (Guttiferae): 1. Infrageneric classification. Bull Brit Mus (Nat Hist) Bot. 1977;5:291–355.Google Scholar
- Robson NKB. Studies in the genus Hypericum L. (Hypericaceae) 9. Addenda, corrigenda, keys, lists and general discussion. Phytotaxa. 2012;72:1–111.Google Scholar
- Ruhfel B, Bittrich V, Bove C, Gustafsson M, Philbrick C, Rutishauser R, et al. Phylogeny of the clusioid clade (Malpighiales): evidence from the plastid and mitochondrial genomes. Am J Bot. 2011;98:306–25.PubMedGoogle Scholar
- Stevens PF. Hypericaceae. In: Kubitzki K, editor. The Families and Genera of Vascular Plants. Volume XI. Berlin, Heidelberg: Springer Verlag; 2007. p. 194–201.Google Scholar
- Robson NKB. Hypericum Botany. In: Ernst E, editor. Hypericum – The genus Hypericum. London, New York: Taylor and Francis; 2003. p. 1–22 [Hardman R (Series Editor): Medicinal and Aromatic Plants – Industrial Profiles].Google Scholar
- Nürk NM, Blattner FR. Cladistic analysis of morphological characters in Hypericum (Hypericaceae). Taxon. 2010;59:1495–507.Google Scholar
- Nürk NM, Madriñán S, Carine MA, Chase MW, Blattner FR. Molecular phylogenetics and morphological evolution of St. John’s wort (Hypericum; Hypericaceae). Mol Phylogenet Evol. 2013;66:1–16.PubMedGoogle Scholar
- Sanchez Meseguer A, Aldasoro JJ, Sanmartín I. Bayesian inference of phylogeny, morphology and range evolution reveals a complex evolutionary history in St. John’s wort (Hypericum). Mol Phylogenet Evol. 2013;67:379–403.Google Scholar
- Lewis P. A likelihood approach to estimating phylogeny from discrete morphological character data. Syst Biol. 2001;50:913–25.PubMedGoogle Scholar
- Lemey P, Rambaut A, Drummond AJ, Suchard MA. Bayesian phylogeography finds its roots. PLoS Comput Biol. 2009;5:e1000520.PubMed CentralPubMedGoogle Scholar
- Cohen KM, Finney SC, Gibbard PL, Fan JX. The ICS international chronostratigraphic chart. Episodes. 2013;36:199–204.Google Scholar
- Robson NKB. Studies in the genus Hypericum L. (Guttiferae): 2. Characters of the genus. Bull Brit Mus (Nat Hist) Bot. 1981;8:55–226.Google Scholar
- Nürk NM, Scheriau C, Madriñán S. Explosive radiation in high Andean Hypericum – Rates of diversification among New World lineages. Front Genet. 2013;4:1–14.Google Scholar
- Meseguer AS, Sanmartín I, Marcussen T, Pfeil BE. Utility of low-copy nuclear markers in phylogenetic reconstruction of Hypericum L. (Hypericaceae). Plant Syst Evol. 2014;300:1503–14.Google Scholar
- Magallón SS, Sanderson MJ. Absolute diversification rates in Angiosperm clades. Evolution. 2001;55:1762–80.PubMedGoogle Scholar
- Burnham KP, Anderson DR. Model Selection and Multimodel Inference: A Practical Information-Theoretic Approach. 2nd ed. Heidelberg: Springer Verlag; 2002.Google Scholar
- Uyeda JC, Harmon LJ. A novel Bayesian method for inferring and interpreting the dynamics of adaptive landscapes from phylogenetic comparative data. Syst Biol. 2014;63:902–18.PubMedGoogle Scholar
- Davis CC, Bell CD, Mathews S, Donoghue MJ. Laurasian migration explains Gondwanan disjunctions: evidence from Malpighiaceae. Proc Natl Acad Sci U S A. 2002;99:6833–7.PubMed CentralPubMedGoogle Scholar
- Tiffney BH. The Eocene North Atlantic land bridge – its importance in Tertiary and modern phytogeography of the Northern Hemisphere. J Arnold Arbor. 1985;66:243–73.Google Scholar
- Tiffney BH. Phylogeography, fossils, and northern hemisphere biogeography: the role of physiological uniformitarianism 1. Ann Mo Bot Gard. 2008;95:135–43.Google Scholar
- Denk T, Grímsson F, Zetter R. Episodic migration of oaks to Iceland: evidence for a North Atlantic “land bridge” in the latest Miocene. Am J Bot. 2010;97:276–87.PubMedGoogle Scholar
- Sanchez Meseguer A, Sanmartín I. Paleobiology of the genus Hypericum (Hypericaceae): a survey of the fossil record and its palaeogeographic implications. An Jard Bot Madr. 2012;69:97–106.Google Scholar
- Tiffney BH, Manchester S. The use of geological and paleontological evidence in evaluating plant phylogeographic hypotheses in the Northern Hemisphere Tertiary. Int J Plant Sci. 2001;162:3–17.Google Scholar
- Wolfe JA. Some aspects of plant geography of the Northern Hemisphere during the Late Cretaceous and Tertiary. Ann Mo Bot Gard. 1975;62:264–79.Google Scholar
- Mosbrugger V, Utescher T, Dilcher DL. Cenozoic continental climatic evolution of Central Europe. Proc Natl Acad Sci U S A. 2005;102:14964–9.PubMed CentralPubMedGoogle Scholar
- Hansen TF, Pienaar J, Orzack SH. A comparative method for studying adaptation to a randomly evolving environment. Evolution. 2008;62:1558–5646.Google Scholar
- Stanley SM. Earth System History. 3rd ed. London: Macmillan Publishers; 2008.Google Scholar
- Ziman SN, Keener CS. A geographical analysis of the family Ranunculaceae. Ann Mo Bot Gard. 1989;125:1198–205.Google Scholar
- Kellogg EA. Evolutionary history of the grasses. Plant Physiol. 2001;125:1198–205.PubMed CentralPubMedGoogle Scholar
- Bouchenak-Khelladi Y, Muasya AM, Linder HP. A revised evolutionary history of Poales: origins and diversification. Bot J Linn Soc. 2014;175:4–16.Google Scholar
- Bouchenak-Khelladi Y, Verboom GA, Savolainen V, Hodkinson TR. Biogeography of the grasses (Poaceae): a phylogenetic approach to reveal evolutionary history in geographical space and geological time. Bot J Linn Soc. 2010;162:543–57.Google Scholar
- Humphreys AM, Linder HP. Evidence for recent evolution of cold tolerance in grasses suggests current distribution is not limited by (low) temperature. New Phytol. 2013;198:1261–73.PubMedGoogle Scholar
- Sytsma KJ, Walker JB, Schonenberger J, Anderberg AA. Phylogenetics, Biogeography, and Radiation of Ericales. In: Botany 2006: 28 July-02 August 2006; Chico. 2006. p. 435.Google Scholar
- Robson NKB. Studies in the genus Hypericum L. (Guttiferae): 3. Sections 1. Campylosporus to 6a. Umbraculoides. Bull Brit Mus (Nat Hist) Bot. 1985;12:1–325.Google Scholar
- Robson NKB. Studies in the genus Hypericum L. (Guttiferae): 7. Section 29. Brathys (part 1). Bull Brit Mus (Nat Hist) Bot. 1987;16:1–106.Google Scholar
- Robson NKB. Studies in the genus Hypericum L. (Guttiferae): 8. Sections 29. Brathys (part 2) and 30. Trigynobrathys. Bull Nat Hist Mus London Bot. 1990;20:1–151.Google Scholar
- Robson NKB. Studies in the genus Hypericum L. (Guttiferae): 6. Sections 20. Myriandra to 28. Elodes. Bull Nat Hist Mus London Bot. 1996;26:75–217.Google Scholar
- Robson NKB. Studies in the genus Hypericum L. (Guttiferae) 4(1). Sections 7. Roscyna to 9. Hypericum sensu lato (part 1). Bull Nat Hist Mus London Bot. 2001;31:37–88.Google Scholar
- Robson NKB. Studies in the genus Hypericum L. (Guttiferae) 4(2). Section 9. Hypericum sensu lato (part 2): subsection 1. Hypericum series 1. Hypericum. Bull Nat Hist Mus London Bot. 2002;32:61–123.Google Scholar
- Robson NKB. Studies in the genus Hypericum L. (Clusiaceae). Section 9. Hypericum sensu lato (part 3): subsection 1. Hypericum series 2. Senanensia, subsection 2. Erecta and section 9b. Graveolentia. Syst Biodivers. 2006;4:19–98.Google Scholar
- Robson NKB. Studies in the genus Hypericum L. (Hypericaceae) 5(1). Sections 10. Olympia to 15/16. Crossophyllum. Phytotaxa. 2010;4:5–126.Google Scholar
- Robson NKB. Studies in the genus Hypericum L (Hypericaceae) 5(2). Sections 17. Hirtella to 19. Coridium. Phytotaxa. 2010;4:127–258.Google Scholar
- Blattner FR. Direct amplification of the entire ITS region from poorly preserved plant material using recombinant PCR. Biotechniques. 1999;27:1180–6.PubMedGoogle Scholar
- Löhne C, Borsch T. Molecular evolution and phylogenetic utility of the petD group II intron: a case study in basal angiosperms. Mol Biol Evol. 2005;22:317–32.PubMedGoogle Scholar
- Taberlet P, Gielly L, Pautou G, Bouvet J. Universal primers for amplification of three non-coding regions of chloroplast DNA. Plant Mol Biol. 1991;17:1105–9.PubMedGoogle Scholar
- Korotkova N, Schneider J, Quandt D, Worberg A, Zizka G, Borsch T. Phylogeny of the eudicot order Malpighiales: analysis of a recalcitrant clade with sequences of the petD group II intron. Plant Syst Evol. 2009;282:201–28.Google Scholar
- Drummond AJ, Ashton B, Buxton S, Cheung M, Cooper A, Heled J, et al.: Geneious v5.4 [http://www.geneious.com]. Accessed 17 Nov 2011.
- Katoh K. MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005;33:511–8.PubMed CentralPubMedGoogle Scholar
- Katoh K, Misawa K, Kuma K-i, Miyata T. MAFFT: a novel method for rapid multiple sequence alignment based on fast Fourier transform. Nucleic Acids Res. 2002;30:3059–66.PubMed CentralPubMedGoogle Scholar
- Castresana J. Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. Mol Biol Evol. 2000;17:540–52.PubMedGoogle Scholar
- Felsenstein J. Evolutionary trees from DNA sequences: a maximum likelihood approach. J Mol Evol. 1981;17:368–76.PubMedGoogle Scholar
- Mau B, Newton MA, Larget B. Bayesian phylogenetic inference via Markov chain Monte Carlo methods. Biometrics. 1999;55:1–12.PubMedGoogle Scholar
- Stamatakis A. RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models. Bioinformatics. 2006;22:2688–90.PubMedGoogle Scholar
- Silvestro D, Michalak I. raxmlGUI: a graphical front-end for RAxML. Org Divers Evol. 2011;12:335–7.Google Scholar
- Ronquist F, Teslenko M, van der Mark P, Ayres DL, Darling A, Höhna S, et al. MrBayes 3.2: efficient Bayesian phylogenetic inference and model choice across a large model space. Syst Biol. 2012;61:539–42.PubMed CentralPubMedGoogle Scholar
- Stamatakis A, Hoover P, Rougemont J. A rapid bootstrap algorithm for the RAxML web servers. Syst Biol. 2008;57:758–71.PubMedGoogle Scholar
- Nylander JAA: MrModeltest v2 [https://github.com/nylander/MrModeltest2]. Accessed 1 Jan 2013.
- Posada D, Buckley TR. Model selection and model averaging in phylogenetics: advantages of Akaike information criterion and Bayesian approaches over likelihood ratio tests. Syst Biol. 2004;53:793–808.PubMedGoogle Scholar
- Brown JM, Hedtke SM, Lemmon AR, Lemmon EM. When trees grow too long: investigating the causes of highly inaccurate bayesian branch-length estimates. Syst Biol. 2010;59:145–61.PubMedGoogle Scholar
- Rambaut A, Drummond AJ: Tracer v1.5 [http://beast.bio.ed.ac.uk/Tracer]. Accessed 10 Aug 2009.
- Huelsenbeck JP, Rannala B. Phylogenetic methods come of age: testing hypotheses in an evolutionary context. Science. 1997;276:227–32.PubMedGoogle Scholar
- Swofford DL: PAUP*. Phylogenetic analyses using parsimony (*and other methods) v4.0b10. Sinauer Associates; 2002.Google Scholar
- Drummond AJ, Ho SYW, Phillips MJ, Rambaut A. Relaxed phylogenetics and dating with confidence. PLoS Biol. 2006;4:e88.PubMed CentralPubMedGoogle Scholar
- Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214.PubMed CentralPubMedGoogle Scholar
- Miller NJ, Calkin PE. Paleoecological interpretation and age of an interstadial lake bed in western New York. Quat Res. 1992;37:75–88.Google Scholar
- Wijninga VM. Neogene ecology of the Salto de Tequendama site (2475 m altitude, Cordillera Oriental, Colombia): the paleobotanical record of montane and lowland forests. Rev Palaeobot Palynol. 1996;92:97–156.Google Scholar
- Arbuzova O: Hypericum. In Iskopaemye tsvetkovye rastenija Rossii i sopredel’nyh gosudarstv [Fossil flowering plants of Russia and adjacent countries]. Vol IV, Nyctaginaceae-Salicaceae. Edited by Budantsev L. Saint Petersburg; 2005: 48–52Google Scholar
- Stadler T. On incomplete sampling under birth-death models and connections to the sampling-based coalescent. J Theor Biol. 2009;261:58–66.PubMedGoogle Scholar
- Sanderson MJ. A nonparametric approach to estimating divergence times in the absence of rate constancy. Mol Biol Evol. 1997;14:1218–31.Google Scholar
- R Development Core Team: R v3.0.2 [http://www.R-project.org]. Accessed 26 Sep 2013.
- Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–90.PubMedGoogle Scholar
- Lemmon AR, Brown JM, Stanger-Hall K, Moriarty Lemmon E. The effect of ambiguous data on phylogenetic estimates obtained by maximum likelihood and Bayesian inference. Syst Biol. 2009;58:130–45.PubMedGoogle Scholar
- Nürk NM, Uribe-Convers S, Gehrke B, Tank DC, Blattner FR: Data from: Oligocene niche shift, Miocene diversification - Cold tolerance and accelerated speciation rates in the St. John’s Worts (Hypericum, Hypericaceae). Dryad Data Repository; 2015 [http://dx.doi.org/10.5061/dryad.4tm8j]
- Brummitt RK, Pando F, Hollis S, Brummitt NA. World Geographical Scheme for Recording Plant Distributions. Pittsburgh: Hunt Institute for Botanical Documentation; 2001.Google Scholar
- Ree RH, Moore BR, Webb CO, Donoghue MJ. A likelihood framework for inferring the evolution of geographic range on phylogenetic trees. Evolution. 2005;59:2299–311.PubMedGoogle Scholar
- Ree RH, Smith SA. Maximum likelihood inference of geographic range evolution by dispersal, local extinction, and cladogenesis. Syst Biol. 2008;57:4–14.PubMedGoogle Scholar
- Smith SA, Dunn CW. Phyutility: a phyloinformatics tool for trees, alignments and molecular data. Bioinformatics. 2008;24:715–6.PubMedGoogle Scholar
- Beaulieu JM, Tank DC, Donoghue MJ. A Southern Hemisphere origin for campanulid angiosperms, with traces of the break-up of Gondwana. BMC Evol Biol. 2013;13:80.PubMed CentralPubMedGoogle Scholar
- Hijmans RJ, Cameron SE, Parra JL, Jones PG, Jarvis A. Very high resolution interpolated climate surfaces for global land areas. Int J Climatol. 2005;25:1965–78.Google Scholar
- Hijmans RJ: raster: Geographic data analysis and modeling [http://cran.r-project.org/web/packages/raster/index.html]. Accessed 14 Aug 2014.
- Revell LJ. Size-correction and principal components for interspecific comparative studies. Evolution. 2009;63:3258–68.PubMedGoogle Scholar
- Revell LJ. phytools: an R package for phylogenetic comparative biology (and other things). Meth Ecol Evol. 2011;3:217–23.Google Scholar
- Uyeda JC, Eastman J, Harmon LJ: bayou: Bayesian fitting of Ornstein-Uhlenbeck models to phylogenies [http://CRAN.R-project.org/package=bayou]. Accessed 7 Jul 2014.
- Gelman A, Rubin DB. Inferences from iterative simulation using multiple sequences. Stat Sci. 1992;7:457–511.Google Scholar
- Plummer M, Best N, Cowles K, Vines K. CODA: convergence diagnosis and output analysis for MCMC. R News. 2006;6:7–11.Google Scholar
- Rabosky DL. BAMMtools: an R package for the analysis of evolutionary dynamics on phylogenetic trees. J Am Stat Assoc. 2014;90:773–95.Google Scholar
- Rabosky DL. Automatic detection of key innovations, rate shifts, and diversity-dependence on phylogenetic trees. PLoS One. 2014;9:e89543.PubMed CentralPubMedGoogle Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.