Skip to main content

Phylogenetic and paleobotanical evidence for late Miocene diversification of the Tertiary subtropical lineage of ivies (Hedera L., Araliaceae)



Hedera (ivies) is one of the few temperate genera of the primarily tropical Asian Palmate group of the Araliaceae, which extends its range out of Asia to Europe and the Mediterranean basin. Phylogenetic and phylogeographic results suggested Asia as the center of origin and the western Mediterranean region as one of the secondary centers of diversification. The bird-dispersed fleshy fruits of ivies suggest frequent dispersal over long distances (e.g. Macaronesian archipelagos), although reducing the impact of geographic barriers to gene flow in mainland species. Genetic isolation associated with geographic barriers and independent polyploidization events have been postulated as the main driving forces of diversification. In this study we aim to evaluate past and present diversification patterns in Hedera within a geographic and temporal framework to clarify the biogeographic history of the genus.


Phylogenetic (biogeographic, time divergence and diversification) and phylogeographic (coalescence) analyses using four DNA regions (nrITS, trnH-psbA, trnT-trnL, rpl32) revealed a complex spatial pattern of lineage divergence. Scarce geographic limitation to gene flow and limited diversification are observed during the early-mid Miocene, followed by a diversification rate increase related to geographic divergence from the Tortonian/Messinian. Genetic and palaeobotanical evidence points the origin of the Hedera clade in Asia, followed by a gradual E-W Asian extinction and the progressive E-W Mediterranean colonization. The temporal framework for the E Asia - W Mediterranean westward colonization herein reported is congruent with the fossil record. Subsequent range expansion in Europe and back colonization to Asia is also inferred. Uneven diversification among geographic areas occurred from the Tortonian/Messinian onwards with limited diversification in the newly colonized European and Asian regions. Eastern and western Mediterranean regions acted as refugia for Miocene and post-Miocene lineages, with a similar role as consecutive centers of centrifugal dispersal (including islands) and speciation.


The Miocene Asian extinction and European survival of Hedera question the general pattern of Tertiary regional extinction of temperate angiosperms in Europe while they survived in Asia. The Tortonian/Messinian diversification increase of ivies in the Mediterranean challenges the idea that this aridity period was responsible for the extinction of the Mediterranean subtropical Tertiary flora. Differential responses of Hedera to geographic barriers throughout its evolutionary history, linked to spatial isolation related to historical geologic and climatic constraints may have shaped diversification of ivies in concert with recurrent polyploidy.


Hedera (ivies) is an Old World root-climber genus that extends from N Africa to Europe and S Asia [14]. The main diagnostic characters for species identification and recognition are morphological features from foliar trichomes and vegetative leaves [5]. However, ploidy level and geographic distribution provide fundamental information for species delimitation [5, 6] (Fig. 1). For example, ploidy level was essential for the identification of two morphologically similar species that were traditionally considered as the same species (H. helix: 2×, H. hibernica: 4×) [7], or for the segregation of two N African endemics (H. algeriensis: 4×, H. maroccana: 2×) [810]. In addition, geographic isolation helped distinguish two closely related species (H. iberica, SW Iberian Peninsula; H. maderensis, Madeira) [11], or disclose incipient speciation processes (H. nepalensis var. nepalensis, Himalaya; var. sinensis, E & SE China) [12, 13]. The combination of morphological and cytogenetic variation together with geographic information, help characterize 12 species (14 taxa): 6 diploid species (3 island endemics), 2 tetraploid species, 4 hexaploid taxa (2 island endemics), and 1 octoploid species (Fig. 1). The numerous island endemics (five) and the strong geographic structure detected in the DNA sequence variation [4, 11, 1416] are interpreted as an imprint of the geographic barriers in the diversification process of Hedera. However, the endozoochorus dispersal syndrome of ivies, mainly mediated by birds [17, 18], together with the winter ripping of their fleshy fruits when food is scarce for animals, suggests that small geographic obstacles might not be such effective barriers to gene flow.

Fig. 1
figure 1

Geographic distribution of the samples of Hedera included in the phylogeographic study. In parenthesis level of ploidy. Coloured envelopes delimit the five biogeographic regions recognised in Hedera

Ivies usually occupy shaded and humid understories of temperate and subtropical woodlands and riparian vegetation. However, they can also live in extremely dry environments such as fully sun-exposed rocks [19, 20]. Not all species are equally tolerant to both deep shaded and drought environments. Indeed, whereas H. helix occurs under the above-described contrasted environments, other species show very strict environmental requirements. For example, H. iberica is restricted to warm humid places (“Canuto”) in southern Iberia [21], where remnants of the Tertiary flora also occur [22, 23]. Another example is H. canariensis that only occurs in humid, warm and shaded understories of the Macaronesian subtropical woodland ‘Laurisilva’ [24]. The fact that some species display the subtropical affinity that characterized the Asian Palmate group of the Araliaceae while others show a strict temperate tendency, may be suggesting that climate might have also contributed to the speciation in Hedera.

Different phylogenetic studies have described a very complex evolutionary history for this genus [4, 6, 1316, 25]. The origin of the Hedera clade has been estimated in the Oligocene in Asia during the deep radiation of the Asian Palmate group of the Araliaceae [25, 26]. In previous phylogeographic studies conducted within Hedera, Asia was also suggested as the ancestral area of the extant species of ivies [16, 25]. It has been hypothesized that extinctions in Asia together with Mediterranean colonization and diversification have led to the present distribution. In this scenario, current Asian species of the genus would be the result of a re-colonization from the Mediterranean [15].

In this study we aim to reconstruct ancient and recent diversification patterns in Hedera under the working hypothesis that geographic barriers have determined the main patterns of diversification by promoting speciation. Under this working hypothesis, divergence events in Hedera would be expected to occur preferentially between areas after a colonization event. The resulting geographic isolation would have led to a strong geographic structure in the genetic variation. To evaluate this hypothesis we analyzed three plastid DNA spacer/intron regions (trnH-psbA, trnT-trnL, and rpL12) and the nuclear ribosomal ITS (nrDNA) region. First, a nrITS dated phylogeny was reconstructed and used as a starting point to conduct biogeographic and diversification analyses. Phylogenetic results based on nrITS data were examined together with those obtained from the phylogeographic analysis of the plastid dataset. To achieve our ultimate goal of clarifying the biogeographic and phylogeographic history of Hedera, the following specific objectives were addressed, the: (1) study of the Mediterranean and Asian areas using a geographically-balanced targeted sampling, (2) reconstruction of past and present diversification patterns, and (3) evaluation of the importance of geographic barriers in promoting isolation.


Taxon sampling and sequencing

Phylogenetic sampling

A phylogenetic-based study was performed to provide a temporal context to conduct the biogeographic and diversification analyses for reconstructing the evolution of the Hedera clade. Wide sampling of outgroup is needed for biogeographic inferences when the nodes of interest approximate the root of the ingroup tree because node’s estimates partly rely on the optimization of their stems [27]. Therefore, all the generic-lineages of the Asian Palmate group have been included, as well as the putative sister-group of the Asian Palmate group (the Aralia-Panax group). The phylogenetic-based analyses (biogeographic, diversification and divergence age analyses) used the nrITS region because: (1) there is a large number of available sequences of Araliaceae, (2) it is more variable than the fastest evolving plastid regions and provides more resolved tree topologies, (3) it better complements the evolutionary history of the genus where nuclear and plastid incongruence has been previously reported due to hybridization [15, 25], and (4) main diagnostic characters in the taxonomy of Hedera come from foliar trichomes and trichomes are genetically controlled by nuclear genes [28, 29]. In any case, all the analyses were also done with the plastid dataset used for the phylogeographic study (see below) but not included in the study because the lack of branch support may have resulted in inaccurate interpretations of the biogeographic, divergence age and diversification results. The nrITS dataset included 34 samples representing 12 species of Hedera, 44 of the other 20 generic-lineages of the Asian Palmate group, and 12 other Araliaceae genera (Additional file 1). Harmsiopanax ingens was used as the outgroup. All the 90 nrITS sequences were obtained from previous studies [6, 15, 25, 26] and downloaded from GenBank (, Additional file 1).

Phylogeographic sampling

A phylogeographic study was conducted to reconstruct the geographic pattern of genetic diversity within Hedera, including 153 samples representing the 12 species (14 taxa) recognised (Fig. 1 and Additional file 2). The number of samples per species varied between 5 and 40, except for H. algeriensis (endemic to N Algeria and N Tunisia) for which only two samples were available. Sampling effort was more intensive on the two most widespread species, leading to the inclusion of 40 samples of the European H. helix and 32 of the Asian H. nepalensis. Samples were selected to represent the whole geographic range of each species with an emphasis on the areas considered as Tertiary refugia both in the Mediterranean [30] and in China [31]. To investigate the geographic origin of Hedera, Kalopanax septemlobus was included as the outgroup (Additional file 2) [25]. Additionally, to evaluate the potential impact of the uncertainty on the sister group of Hedera [25], six other Asian Palmate genera were also included (Additional file 2).

Three plastid DNA regions were analyzed (rpL32, trnH-psbA and trnT-trnL) for this part of the study. The primers used for the amplifications were as follows: (1) trn a and trn b for trnT-L spacer [32], (2) rpL32F and trnL(UAG) for the rpL32 intron [33], and (3) trnHR and psbAF for the trnH-psbA spacer [34]. Amplifications and sequencing protocols followed Valcárcel et al. [15] for the trnT-trnL region, Mitchell et al. [26] for the trnH-psbA spacer and Shaw et al. [33] for the rpL32 intron. As a result 270 sequences were newly generated in this study (76 for rpL32, 89 for trnT-trnL, and 111 for trnH-psbA). The three plastid DNA regions of K. septemlobus were taken from Li et al. [35] and downloaded from GenBank (, as well as for the other six genera of Araliaceae included. The sampling for the trnT-trnL spacer was completed by the addition of 44 sequences from our previous phylogeographic studies [15, 16]. Three DNA matrices were compiled using only Kalopanax as the outroup: trnT-trnL (134 samples, 89 new sequences), trnH-psbA (112 samples, 111 new sequences), and rpL32 (77 samples, 76 new sequences). Alignments were carried out with MUSCLE [36] followed by manual revision in Geneious v9.0.5 ( Sequences were concatenated into a fourth matrix with the program Sequence Matrix [37], only including samples with the three DNA regions sequenced (66 samples). A fifth matrix was additionally built to check for the impact of the different rootings of the Hedera network on the geographic interpretations, including 65 samples of Hedera plus 7 different genera of Araliaceae.

Phylogenetic-based analyses

Divergence age estimates

Divergence age estimates were inferred from the nrITS matrix through a relaxed molecular clock implemented in Beast v.1.7.5 [38]. The substitution rate variation was modeled using an uncorrelated lognormal distribution and a Birth-Death process was applied to model speciation. The best evolutionary model for each of the DNA regions was selected by jModeltest setting a threshold of 3 ∆AIC (Additional file 3) [39]. The analyses were run in the absence of topological constraints, except for the calibration nodes. Two MCMC analyses were run for 100 million generations sampled every 10,000 generations. Convergence, mixing and effective sample size (ESS) of model parameters were assessed using Tracer 1.5 [40]. Samples from the two independent runs were pooled after removing a 25% burn-in using Log Combiner 1.7.5 [38]. Trees were summarized in a maximum clade credibility (MCC) tree obtained in TreeAnotator 1.7.5. Seven leaf macrofossils and two pollen grain fossils have been recorded in Hedera (Table 1). The taxonomy of ivies is mainly based on foliar trichomes that are not well preserved in fossils. Also, pollen grains do not show any morphological variation between the extant species of ivies [3]. Therefore, certainty on the phylogenetic placements of these fossils is limited. Only the oldest fossil found (Oligocene, Table 1) can be placed with certainty as a calibration point at the stem of Hedera, as inferred from the age recovered for the lineage of Hedera in previous studies [between 36.6 Mya and 51.55 Mya; 25]. However, we decided not to use this node as a calibration point in the final analyses because the stem of Hedera represents an uncertain node in the phylogeny [25]. Instead, two previous divergence time estimates obtained from plastid DNA were employed as secondary calibration points [25]. This secondary calibration approach is more conservative since the age estimates in which it is based were obtained from a fossil-based calibration with a certain placement of fossils on robust nodes [25]. Accordingly, the crown groups of the Asian Palmate and Hedera clades were set as normal distributions of 72.55 ± 9.0 Myr and 7.65 ± 3.5 Myr, respectively. Calibration accuracy was tested by comparing divergence times herein estimated to the ages of reliable Hedera fossils [4150] (Table 1). A second Beast analysis was performed using the oldest fossil of Hedera as the minimum age for the stem of Hedera to double check the posterior age recovered for the crown of Hedera in the secondary calibration estimate.

Table 1 Fossil records of Hedera (entries arranged in chronological order)

Biogeographic range estimation

Estimation of spatial patterns of geographic diversification in Hedera was conducted using a model-based likelihood method (Lagrange) and the nrITS dataset. A Dispersal-Extinction-Cladogenesis analysis was performed over multiple trees using a script provided by Richard Ree (pers. com.). A multi-tree approach was essential to account for the impact of phylogenetic uncertainty on the ancestral range estimate of the stem group of Hedera, due to its ambiguous sister-group relationship [25]. For this purpose, 1000 post-burnin trees from the nrITS Beast analysis were randomly selected with the R-package ape [51] and used as input trees for Lagrange. The geographic range of the Araliaceae was divided into eight regions based on floristic endemicity with special emphasis on Hedera distribution: (A) tropical Africa, (B) Neotropics, (C) Australia, (D) western Mediterranean region (including Macaronesia, hereafter W Mediterranean), (E) eastern Mediterranean region (hereafter E Mediterranean), (F) Europe, (G) western Asia (hereafter W Asia), and (H) eastern Asia (hereafter E Asia). The codification of areas for each sample is provided in Additional file 1. In range constraints, adjacency of areas was allowed only between areas that share the edge (i.e., between W Mediterranean and Europe, and between E Mediterranean and W Asia). Maximum range size was set to two areas. Ranges allowed in the analysis included all possible combinations within those imposed by adjacency and maximum range size. Because we were only interested in recovering the early evolution of Hedera, results were only computed for the most internal nodes of Hedera with support (posterior probability (PP) >0.95; Fig. 2). Particularly, eight nodes were analyzed: the stem and crown groups of the Hedera clade (nodes 0 and 1; Fig. 2), two main clades (nodes 2 and 3; Fig. 2) and two main subclades (nodes 4 and 5; Fig. 2). The output file obtained from the Lagrange analysis of 1000 trees was read and parsed with a new R script herein designed for parsing multi-tree Lagrange results (Additional file 4). This script makes automatic parsing Lagrange multi-tree results, which simplifies the process of summarizing results saving time. The results were summarized as the mean of probabilities estimated by Lagrange for the posterior trees analyzed. The specific ancestral areas of the descendants of a given node are provided when a congruent biogeographic pattern is consistently recovered over the multi-trees analyzed (e.g., nodes 0 and 1; Fig. 2). If the estimated ancestral areas of the descendants of a given node i resulted in incongruent biogeographic patterns when the multi-trees results were analyzed together, only the ancestral areas of node i are provided with no specification to the descendant lineages. The ancestral area for a given node i is assumed to be the combination of the ancestral areas estimated for its descendants. For example, if the most probable biogeographic patterns for a given node are (1) E|D at a mean probability of 0.39, (2) D|E at a mean probability of 0.16, and (3) ED|D at a mean probability of 0.14; then, a simplification of the results are shown by providing the ancestral area for the node as ED with a mean probability of 0.69.

Fig. 2
figure 2

Beast Maximum Clade Credibility chronogram of the Hedera clade from the nrITS dataset. Mean ages and 95% CI are only represented for clades with >0.95 Posterior Probability support. Nodes of interest are labelled. Ancestral areas and estimated probabilities obtained from the Lagrange multi-tree analysis are provided only for the nodes of interest. The complete Asian Palmate group MCC tree including the 90 sequences of Araliaceae is shown in Additional file 5. Branch lengths were modified from the original tree (Additional file 5) to better fit the biogeographic results

Diversification analyses

The 6500 post-burnin Beast trees were pruned to only consider the clade of Hedera and one tip per species, except for non-monophyletic species for which one tip per species-lineages was kept (18 tips, 14 taxa). We decided to only analyze the clade of Hedera because the independent analysis of a particular clade is recommended to isolate its diversification pattern from the heterogeneous diversification patterns of other clades in the phylogeny [52]. The stem of Hedera was included in the analysis because ignoring long branches before crown nodes may result in inaccurate interpretations of the diversification pattern within the crown group [53, 54]. The resulting 6500 post-burnin pruned trees were used as the inputs for the diversification analyses. The log-transformed number of extant taxa was plotted against time (LTT plot) for the 6500 post-burnin pruned trees using the R-package ape [51]. Fitness to speciation models with one, two or three diversification rates was also tested as implemented in the R-package LASER [55]. The best evolutionary model was selected for each of the 6500 post-burnin pruned trees based on the AIC using a ∆AIC of 4. Two contingency table tests with one dimension and two levels were performed over the results of the 6500 trees using chi-square goodness of fit [56] in R [57]: (1) number of trees with constant vs. variable rate models and (2) number of trees with Yule two rates vs. Yule three rates variable models.

The Phylogenetic Diversity (PD) index of Faith [58] measures the length of evolutionary pathways that connect a given set of taxa as the sum of branch lengths connecting taxa in a given area. In this study, Faith’s index was used to account for the amount of PD of Hedera represented in each of the five major endemicity areas delimited for ivies in the biogeographic analyses (W Mediterranean, Europe, E Mediterranean, W Asia, E Asia). The aim of this analysis was to evaluate whether the most species-rich endemicity areas also hold the greatest evolutionary diversity. Faith’s PD was estimated as implemented in the R-package Picante [59]. To assess if regions have significantly higher or lower PDs than at random expectations for a given number of species, PDs were calculated over the 1000 randomized posterior trees used for the Lagrange analysis but pruned to only consider the clade of Hedera and one tip per species-lineages. Two-tailed test was used to compare the observed PDs to the null distribution of the 1000 random replicates (significance level of 0.05).

Network reconstructions based on the coalescence [58] analysis of plastid haplotypes were performed on the five-plastid DNA matrices. Statistical Parsimony (SP) was applied as implemented in TCS 1.13 [60]. The 95% probability limit of parsimonious connections was applied and gaps were coded as missing data. Predictions from coalescent theory were applied to deals with homoplasy [61, 62]. To test the hypothesis that interior and tip haplotypes are equally frequent we used a contingency table test with one dimension and two levels using chi-square goodness of fit [56] in R [57].


Estimates of divergence times

The early divergence age estimates of the Hedera clade (Fig. 2, Additional file 5) fit the timing set by the fossil record (Table 1). The two main clades of Hedera (diploid and polyploid) diverged in parallel during the late Miocene (9.47 / 9.15 Mya, 5.02–14.91 / 5.48–13.31 Mya 95% CI; nodes 2 / 3, Fig. 2), five million years after the crown age divergence of the Hedera clade in the early-mid Miocene (14.23, 9.96–19.1 Mya 95% CI; node 1; Fig. 2). The divergence time estimate for the eastern polyploid subclade is 4.79 Mya (1.66–9.04 Mya 95% CI; node 4), whereas 3.37 Mya (1.05–7.08 Mya 95% CI; node 5) is recovered for the western polyploid subclade. The diploid clade displays a large basal polytomy that prevents from any further divergence time or biogeographic estimates (Fig. 2). The posterior estimates on the calibration node 1 lay outside the bound of the calibration, which is interpreted as low influence of this calibration prior on the posterior. The MCC Beast tree obtained after removing this secondary calibration point and including the oldest fossil of Hedera as the calibration point at the stem of Hedera revealed similar posterior node estimates (Additional file 6).

Ancestral range inference

Results from the multi-tree Lagrange analysis (1000 posterior random trees) plotted on the nrITS MCC Beast tree of Hedera posit Asia as the most likely ancestral area for the stem of the genus (|H or |G, P = 0.66; node 0; Fig. 2). The most probable range inferred for the crown group of Hedera is the combination of the E Mediterranean and W Asia (P = 0.77; node 1). This result suggests that an extinction event occurred in E Asia between Hedera’s stem and crown along with a dispersal event to the E Mediterranean. Subsequently, two equally plausible biogeographic scenarios are inferred for the diploid clade (from node 1 to 2; Fig. 2): (1) persistence in the ancestral area (EG, P = 0.43; node 2) or (2) W Asian extinction and dispersal to the W Mediterranean (ED, P = 0.34; node 2). The most probable biogeographic scenario for the polyploid clade (from node 1 to 3; Fig. 2) is the E Mediterranean persistence (E, P = 0.93; node 3), W Asian extinction (≠ G, P = 0.71), and the W Mediterranean colonization (D, P = 0.49). The W Mediterranean is recovered as part of the ancestral area for the divergence of the western polyploid subclade (node 5; D, P = 0.85; Fig. 2), whereas divergence of the eastern polyploid subclade may have occurred in the E Mediterranean and W Asia (EG, P = 0.66; node 4) or in the E Mediterranean (EE, P = 0.33).

Diversification analyses

Rate variable models (Yule two rates and Yule three rates) are more frequently selected as the best fitting speciation process than constant ones (pureBirth and birth-death) among the 6500 posterior pruned trees (4003 vs. 2497, respectively; U = 359.5319, p-value <0.001; Additional file 7). The Yule two rates model estimates an initial diversification of 0.04 per lineage per unit of time (−0.14–0.23 95% CI) that increased to 0.29 (0.024–0.55 95% CI) at around 7.82 Mya (5.94–9.71 95% CI). The equally probable Yule three rates (U = 0.89, p-value = 0.4) recovered an initial diversification of 0.04 (−0.29–0.37 95% CI) that increased to 0.94 (−2.96–4.45 95% CI) at around 10.19 Mya (7.42–12.97 95% CI) and decreased to 0.36 (−1.11–1.84 95% CI) at 3.20 Mya (0.35–6.05 95% CI). The LTT plot of the 6500 posterior pruned trees describes a flat initial diversification rate period followed by an ever-increasing diversification pattern from the Miocene onwards for the Hedera clade (Fig. 3). This graphical representation is congruent with the recovery of Yule two rates and Yule three rates as the best fitting speciation models. The E Mediterranean and W Asia are the only regions that reveal greater PDs than expected at random for the given number of species (Table 2). The largest value of Faith’s PD index is estimated for W Mediterranean (72.05, Table 2), while the smallest value is computed for E Asia (32.54). The remaining three main biogeographic areas display similar low values of PD index (Table 2).

Fig. 3
figure 3

Lineage through time plots of the 6500 posterior pruned trees of the Hedera clade from the nrITS dataset

Table 2 Summary of different diversity indicators in Hedera according to the five biogeographic regions recognized
Fig. 4
figure 4

Phylogeographic network of the Hedera plastid matrix (rpL32, trnH-psbA, trnT-trnL) including Kalopanax septemlobus as outgroup. Haplotype numbering is according to Additional file 8. Circle dimensions are proportional to the number of samples displaying each haplotype as indicated at the bottom. Lines indicate a single nucleotide substitution and dots () represent extinct or not-detected haplotypes. Dashed lines indicate ambiguities resolved under predictions of the Coalescent Theory [65]. Abbreviation names of the taxa displaying each haplotype are as follows: ALG, H. algeriensis; AZO, H. azorica; CAN, H. canariensis; COL, H. colchica; CYP, H. pastuchowii subsp. cypria; HEL, H. helix; HIB, H. hibernica; IBE, H. iberica; MAD, H. maderensis; MAR, H. maroccana; NEP, H. nepalensis var. nepalensis; PAS, H. pastuchowii subsp. pastuchowii; RHO, H. rhombea; SIN, H. nepalensis var. sinensis

Phylogeographic networks

Twenty-four haplotypes are detected within Hedera when the three-plastid DNA regions are concatenated (Fig. 4, Additional file 8). The highest number of haplotypes is detected in the Mediterranean region with 15 haplotypes, 57% of them exclusive to E Asia. Seven haplotypes are detected in Asia, 67% of them exclusive. Finally, five haplotypes occur in Europe, with only the 2 of them exclusive. The most frequent haplotype (Hp 8) is detected in 14 Asian samples (9 from E Asia and 5 from W Asia) of two species (H. colchica, H. nepalensis), followed by two Euro-Mediterranean haplotypes (Hp 22: 8 Mediterranean samples, 3 European; Hp 16: 4 European samples, 2 Mediterranean). Comparable geographic patterns of haplotype diversity are obtained when individual plastid matrices are analyzed (Additional files 9, 10 and 11). The W Mediterranean is consistently recovered as the geographic region with the highest number of haplotypes (rpL32: 6 Hps, trnH-psbA: 5 Hps, trnT-trnL: 7 Hps; Additional files 9, 10 and 11). The area with the second highest number of haplotypes is E Asia according to trnH-psbA (4 Hps; Additional file 10) and trnT-trnL (8 Hps; Additional file 11) while Europe according to rpL32 (4 Hps; Additional file 9).

The number of interior haplotypes is 10 (6 unambiguous, 29 samples; 4 ambiguous, 4 samples) while 14 are tip haplotypes (23 samples). Five of the six unambiguous interior haplotypes are from the W Mediterranean (5) and Europe (2) while only one from Asia (Table 2, Fig. 4). The mean number of sequences per haplotype varies from 4.8 for unambiguous interior haplotypes and 1.64 for tip haplotypes. Also the ratio of single-sample haplotyes (singleton) vs. non-singleton haplotypes varies between 0.2 for interior haplotypes and 0.6 for tip haplotypes. Interior haplotypes are more frequent than tip ones (U = 14.087, p-value <0.001). This result together with the application of Templeton’s rules allows us to resolve the network uncertainties (Fig. 4).

Limited taxonomic congruence is detected in the haplotype network, with the widespread Euro-Mediterranean H. helix and H. hibernica and the E Asian H. rhombea displaying between two and eight unrelated haplotypes (Fig. 4 and Additional file 8). However, strong geographic structure and environmental affinity is detected (Fig. 4). The W Mediterranean haplotypes are scattered through the network, whereas the majority of the European and Asian haplotypes are organized in two star-like groups. The Euro-Mediterranean and Asian start-like groups are connected to each other through their most widespread haplotypes (Hp 16 and Hp 8, respectively; Fig. 4). The haplotypes of all the samples included from the southwest Mediterranean are connected with no missing haplotype needed (hereafter “relict SW Mediterranean haplotypes”). The independent analysis of each plastid region consistently recovers the group of relict SW Mediterranean haplotypes and the Asian star-like group connected to the widespread Euro-Mediterranean star-like group (Additional files 9, 10 and 11).

The Hedera network is connected to Kalopanax (outgroup) through the relict SW Mediterranean Hp 1 with 16 missing haplotypes needed for connection. Similarly, a relict SW Mediterranean haplotype is connected to an Asian outgroup when the seven genera of Araliaceae are included (Additional file 12). Connection to outgroup slightly differs when plastid datasets are analyzed independently (Additional files 9, 10 and 11). The biogeographic connection to the outgroup is inferred through a widespread Euro-Mediterranean haplotype in rpL32 (Additional file 9) or through E Asian haplotypes in trnH-psbA (Additional file 10) and trnT-trnL (Additional file 11). The Taiwanese sample that connects to the outgroup in the trnH-psbA network (Additional file 10) is recovered as a tip haplotype when the three DNA regions are analyzed together (Hp 7; Fig. 4).


Uneven geographic diversification of Hedera: E Asia as both ancestral and sink area

Eastern Asia is inferred as the most likely ancestral area for the Hedera clade (Fig. 2), which is congruent with the E Asian location of the oldest fossil record of Hedera (Table 1). This geographic context is also consistent with the fact that all the remaining 20 generic-lineages of the Asian Palmate group occur in E Asia. Indeed, tropical and subtropical SE Asian environments are inferred to be the center of diversification of the Asian Palmate group and some of its genera [26, 6369]. This E Asian placement of the ancestor of the temperate Hedera clade is congruent with the Asian origin of numerous lineages of the flora of Europe [7072]. Westward dispersal and E Asia extinction appear to have occurred early in the evolution of the clade of Hedera leading to an initial divergence during the Lower-Middle Miocene in W Asia and the E Mediterranean (Fig. 2). The fossil record also supports the timing of this westward migration, since reliable fossils of Hedera are found for the first time in W Asia (Georgia) [47] and the E Mediterranean (Greece) [48] during the Lower-Middle Miocene (Table 1). The two early descendant lineages (polyploid and diploid clades) persisted in the E Mediterranean with a probable and independent W Mediterranean colonization coupled with W Asian extinction in the polyploid clade. The low diversification rate found for this period in the Hedera clade (r0 = 0.03–0.04, >10.19 Mya - >7.82 Mya; Fig. 3, Additional file 7) is better explained by a high extinction rate in Asia. This Asian extinction scenario supports the previous biogeographic hypothesis of Asian extinction and Mediterranean differentiation of Hedera [15, 16], but challenges the general pattern of Tertiary regional extinction of temperate genera in Europe and survival in Asia or North America [7375]. The E-W gradual extinction in Asia is not only supported by the biogeographic analysis (Fig. 2), but also by the high number of missing haplotypes that are needed to connect the relict SW Mediterranean group to the Asian outgroup in the network (Fig. 4, Additional file 12). Soon after the westward migration and Asian extinction, a diversification rate increase is detected for Hedera starting in the Tortonian/Messinian (10.19 or 7.82 Mya; Fig. 3, Additional file 7), when ivies were already established at both sides of the Mediterranean as inferred from biogeographic reconstruction (Fig. 2) and fossil record (Table 1). The increase in aridity that synchronously started across the whole basin by that time [7678] has been proposed to be responsible for the extinction of most of the sub-tropical Tertiary plants that previously inhabited this area [7981]. The diversification rate increase of Hedera is intriguing since it places the intensification of diversification under an arid climate, which is an unpredicted scenario for a group of plants that need humid environments. Confinement to wet habitats in locally isolated environments could have promoted isolation and thus speciation for those elements of the sub-tropical Tertiary flora that were also tolerant to arid conditions. This could be the case of ivies that are considered part of the dry sub-tropical Tertiary flora [69, 82, 83]. Indeed, certain Mediterranean traits have been observed in Hedera species: low stomatal density primarily on leaf underside (Virginia Valcárcel pers. obs.), high density of indumentum on young shoots (Hugh A. McAllister pers. obs.), and flowering only in sun-exposed branches (Pablo Vargas pers. obs.). Independent colonization events from the Mediterranean to other areas of Europe and independent back colonizations to Asia are inferred coinciding with the highest diversification periods (10.19–3.20 Mya or 7.82 Mya - present, Additional file 7). Interestingly, the intense diversification that accompanied this post-Miocene geographic range expansion did not occur evenly in all areas. As a result, an unbalanced diversity pattern is currently observed with the ancestral area (E Asia) operating as a sink region for the diversity of extant ivies (Table 2) and the distant Mediterranean basin acting as its source of re-colonization.

E and W Mediterranean as refugia in the Miocene and Pliocene

Two main regions are considered as Tertiary and Quaternary refugia for temperate and subtropical plant lineages in Eurasia, the Mediterranean basin [30] and SW China [73, 84]. The fact that Pliocene and Pleistocene extinctions were more intense in Europe than in Asia [75] has been explained by several factors, such as climatic and topographic changes, complex geological history, or the antiquity of the flora in Asia ([84], but see [85]). As a whole, there is a greater diversity of temperate genera in Asia than in Europe [73]. After a thorough full sampling of Hedera in the Mediterranean and Asian Tertiary refugia (Fig. 1), the Mediterranean region is not only still inferred as the main center of diversification for ivies [16] but also as a cumulative lineage refugium. The E and W Mediterranean are interpreted as consecutive refugia with differential roles in the evolution of the genus. Whereas the E Mediterranean acted as the refugium for Hedera during the dramatic changes of the Miocene, the W Mediterranean is identified as a more recent refugium for ivies during the climatic changes that shaped the Mediterranean flora from the Pliocene onwards. Our biogeographic analysis suggests that the E Mediterranean was part of the ancestral area back during the early differentiation of Hedera in the Miocene (nodes 2–3; Fig. 2), which is also supported by the time of the first fossils of Hedera in the Mediterranean that are from the eastern side [48] (Table 1). This pattern reveals that the E Mediterranean was a phylogenetic center of diversification in the past that served as a Miocene refugium for the ancestors of the lineages of extant ivies. This scenario is also supported by the fact that the E Mediterranean has a greater Faith’s PD than expected by chance (Table 2). The reason why the E Mediterranean currently shows low diversity might be due to higher extinction rates in this area and higher levels of differentiation in the W Mediterranean. Currently, the E Mediterranean only has three species, three ploidy levels, two-plastid tip haplotypes and a relatively low PD index (45, Table 2). The limited diversity of extant ivies in the E Mediterranean contrasts the results of some studies showing this side of the basin as a center of diversification for different groups of organisms [8692]. From the ancestral refugium in the E Mediterranean a centrifugal dispersal resulted in the colonization of the W Mediterranean in the first place and the back colonization of W Asia in more recent times. A similar centrifugal dispersal from the E Mediterranean has been described for Anthemis also in the Miocene [93]. Despite the relevance of the E Mediterranean in the early evolution of the genus, the W Mediterranean is the region that accounts for the highest diversity in ivies. Indeed, the W Mediterranean displays the greatest number of extant species (eight), ploidy levels (four), plastid haplotypes (13, Additional files 7), and Faith’s PD (72, Table 2). This is consistent with the pattern observed in the phylogeographic study where the W Mediterranean samples are scattered through the network displaying internal and tip haplotypes both in relict and derived groups (Fig. 4). Also, seven out of the 10 internal haplotypes and six out of the 13 tip haplotypes are widespread or unique to western Mediterranean samples (Fig. 4). All these results support the W Mediterranean as a secondary center of diversification [6], which might also be the reason for this area to have the greatest PD index but with this index lower than the one expected by chance (Table 2). All sources of evidence agree with the W Mediterranean as a source area of dispersal in the post-Miocene colonization of Europe and E Asia. Indeed, the star-like organization of the Asian and Euro-Mediterranean haplotypes (Fig. 4) suggest a recent and fast colonization from the W Mediterranean.

The relative importance of the E and W Mediterranean regions as refugia for the Mediterranean flora varies among groups [see 30 for details] with a general trend towards considering the W Mediterranean as the main Tertiary refugium (33 cases vs. 19 in E Mediterranean) [30]. It seems reasonable to assume that the role of the E Mediterranean as an ancestral refugium might have been overlooked because the more recent diversification processes might have erased footprints of earlier events, as in the case of Hedera.

Geographic isolation related to new habitats as the driver for Hedera evolution

The pattern of endemicity in Hedera, with 10 endemic taxa (five island endemics) and three widespread species (H. helix, H. hibernica and H. nepalensis), can only be explained by a key role of geographic isolation in speciation. However, geographic isolation seems less likely for a bird-dispersed fleshy-fruited plant group, like Hedera, due to the expected long-distance dispersability of endozoochorous dispersal syndromes mediated by birds [94]. Successful dispersal is supported by the occurrence of three species of Hedera in three archipelagoes of Macaronesia. Nevertheless, complex and different geographic patterns of genetic variation have been suggested for fleshy-fruited plants [14, 95, 96]. Evidence points to additional external factors that may alter the efficiency of dispersal mediated frugivory, such as site availability (spatial limitation) [97]. The geographic structure inferred from the biogeographic and phylogeographic reconstructions of Hedera is not spatially or temporarily constant. Geographic isolation is interpreted for the recent past of the genus from the geographic congruence of main clades diverged during the Pliocene (nodes 4, 5; Fig. 2). In contrast, back to the early evolution of Hedera, diversification occurred within the same geographic area (nodes 1–3, Fig. 2). This may be explained by a combination of historical and contemporary geological and ecological factors. The early divergence of the genus that led to the diploid and polyploid clades most likely took place in the E Mediterranean and W Asia during the Lower-Middle Miocene (Fig. 2, Table 1). By that time the Arabian microplate collided with Eurasia resulting in a major change to a colder climate and the uplift of E Mediterranean basin (~16 Mya) [98]. This may have promoted regional or habitat-dependent isolation that cannot be observed at the scale of biogeographic analyses. The independent range expansion to the W Mediterranean that occurred in the diploid and polyploid clades before the end of the Miocene (Fig. 2) does not seem to result in the observation of any divergence pattern either. During that period the Mediterranean region was under seasonal sub-tropical climate, which may have been crucial for Hedera establishment and dispersal across the basin, hindering geographic isolation. Given that the diploid and polyploid clades persisted in the E Mediterranean, the similarity in the diversification pattern detected could be related to the short compressional event that took place in Eastern Aegean 9–8 Mya [99] that may have promoted a secondary contact between the two clades. Such contact provided a reliable context for the inter-lineage hybridization already proposed in Hedera [6, 15] that might have eventually resulted in genome duplication [100]. Also, it provides a likely explanation for the incongruence detected between plastid and nuclear results [15] (cf. Figs. 2 and 4). The evolution of the diploid clade remains unclear, whereas a vicariant event is found for the divergence of the polyploid clade leading to the eastern and western polyploid subclades. Interestingly the time of lineage differentiation for eastern polyploids appears to be older than that of the western one, which was unexpected given that the W Mediterranean basin is much older [98]. This may be the result of a geographic filter due to the closer proximity of the E Mediterranean to the ancestral range of ivies (E Asia).


The consequences of the geological changes occurred in the Mediterranean area during Lower-Middle Miocene were dramatic for the Mediterranean Tertiary plant lineages [79, 81]. However, the subtropical Hedera clade survived and diversified in an arid and relatively constant climate. Indeed, the diversification rate increase of ivies since the Tortonian/Messinian under the increasing aridity in the Mediterranean, suggests that a climate-driven spatial limitation (i.e. habitat availability) may have enabled geographic speciation. Cumulative Miocene and post-Miocene refugia are detected in E and W Mediterranean, respectively, that also acted as consecutive dispersal centers.



Million years


Million years ago


Phylogenetic index


Posterior probability


  1. Meusel H, Jäger E, Weinert E. Vergleichende Chorologie Der Zentraleuropäischen Flora. Jena: Veb Gustav Fischer Verlag; 1965.

    Google Scholar 

  2. Mabberley DJ. The plant-book. Cambridge: Cambridge University Press; 1997.

    Google Scholar 

  3. Valcárcel V: Taxonomy, systematics and evolution of Hedera L. (Araliaceae). Universidad Pablo de Olavide; 2008.

  4. Green AF, Ramsey TS, Ramsey J. Phylogeny and biogeography of ivies (Hedera spp., Araliaceae), a polyploid complex of woody vines. Syst Bot. 2011;36:1114–27.

    Article  Google Scholar 

  5. Rutherford A, McAllister HA, Mill RR. New ivies from the Mediterranean area and Macaronesia. The Plantsman. 1993;15:115–28.

    Google Scholar 

  6. Vargas P, Mcallister HA, Morton C, Jury SL, Wilkinson MJ. Polyploid speciation in Hedera (Araliaceae): phylogenetic and biogeographic insights based on chromosome counts and ITS sequences. Plant Syst Evol. 1999;219:165–79.

    Article  CAS  Google Scholar 

  7. McAllister HA, Rutherford A. Hedera helix L. and H. hibernica (Kirchner) bean (Araliaceae) in the British isles. Watsonia. 1990;18:7–15.

    Google Scholar 

  8. Rutherford A. The history of the Canary Islands ivy and its relatives. Ivy J. 1984;10:13–8.

    Google Scholar 

  9. Rutherford A. The ivies of Andalusia (southern Spain). Ivy J. 1989;15:7–17.

    Google Scholar 

  10. McAllister HA. Canary and Algerian ivies. The Plantsman. 1988;10:27–9.

    Google Scholar 

  11. Ackerfield J, Wen J. Evolution of Hedera (the ivy genus, Araliaceae): insights from chloroplast DNA data. Int J Plant Sci. 2003;164:593–602.

    Article  CAS  Google Scholar 

  12. Ackerfield J. Trichome morphology in Hedera (Araliaceae). Edinb J Bot. 2001;58:259–67.

    Article  Google Scholar 

  13. Ackerfield J, Wen J. A morphometric analysis of Hedera L. (the ivy genus, Araliaceae). Adansonia Sér. 2002;324:197–212.

    Google Scholar 

  14. Grivet D, Petit RJ. Phylogeography of the common ivy (Hedera sp.) in Europe: genetic differentiation through space and time. Mol Ecol. 2002;11:1351–62.

    Article  CAS  PubMed  Google Scholar 

  15. Valcárcel V, Fiz O, Vargas P. Chloroplast and nuclear evidence for multiple origins of polyploids and diploids of Hedera (Araliaceae) in the Mediterranean basin. Mol Phylogenet Evol. 2003;27:1–20.

    Article  PubMed  Google Scholar 

  16. Valcárcel V, Vargas P. Phylogenetic reconstruction of key traits in the evolution of ivies (Hedera L.). Plant Syst Evol. 2013;229:447–58.

    Article  Google Scholar 

  17. Ridley HN. The dispersal of plants throughout the world. L. Reeve & Co: Kent; 1930.

    Google Scholar 

  18. Guitián J. Dispersal of ivy Hedera helix seeds by birds: time spent in the plant and seed removal efficiency. Ardeola. 1987;34:25–35.

    Google Scholar 

  19. Kollmann J, Grubb PJ. Recruitment of fleshy-fruited species under different shrub species: control by under-canopy environment. Ecol Res. 1999;14:9–21.

    Article  Google Scholar 

  20. Sack L, Grubb PJ, Maranon T. The functional morphology of juvenile plants tolerant of strong summer drought in shaded forest understories in southern Spain. Plant Ecol. 2003;168:139–63.

    Article  Google Scholar 

  21. Valcárcel V, Rutherford A, Miller R, McAllister HA. Hedera L. in Flora iberica. Vol. X. Araliaceae-Umbelliferae. Edited by Nieto Feliner G. Madrid: Departamento de publicaciones del CSIC; 2003.

    Google Scholar 

  22. McAllister HA. New work on ivies. Int Dendrol Soc Yearb. 1981;1981:106–9.

    Google Scholar 

  23. Rodríguez-Sánchez F, Hampe A, Jordano P, Arroyo J: Past tree range dynamics in the Iberian Peninsula inferred through phylogeography and palaeodistribution modelling: A review. Rev Palaeobot Palynol 2010, 162:507–521. [Iberian Floras through Time: Land of Diversity and Survival].

  24. Bramwell D, Bramwell Z. Flores Silvestres de Las Islas Canarias. Editorial Rueda: Madrid, Spain; 1990.

    Google Scholar 

  25. Valcárcel V, Fiz O, Wen J. The origin of the early differentiation of ivies (Hedera L.) and the radiation of the Asian Palmate group (Araliaceae). Mol Phylogenet Evol. 2014;70:492–503.

    Article  PubMed  Google Scholar 

  26. Mitchell A, Li R, Brown JW, Schönberger I, Wen J. Ancient divergence and biogeography of Raukaua (Araliaceae) and close relatives in the southern hemisphere. Aust Syst Bot. 2012;25:432–46.

    Article  Google Scholar 

  27. Ronquist F. Dispersal-vicariance analysis: a new approach to the quantification of historical biogeography. Syst Biol. 1997;45:195–203.

    Article  Google Scholar 

  28. Marks MD. Molecular genetic analysis of trichome development in Arabidopsis. Annu Rev Plant Physiol Plant Mol Biol. 1997;48:137–63.

    Article  CAS  PubMed  Google Scholar 

  29. Ishida T, Kurata T, Okada K. Wada: a genetic regulatory network in the development of trichomes and root hairs. Annu Rev Plant Biol. 2008;59:365–86.

    Article  CAS  PubMed  Google Scholar 

  30. Médail F, Diadema K. Glacial refugia influence plant diversity patterns in the Mediterranean Basin. J Biogeogr. 2009;36:1333–45.

    Article  Google Scholar 

  31. López-Pujol J, Zhang F-M, Sun H-Q, Ying T-S, Ge S. Centres of plant endemism in China: places for survival or for speciation? J Biogeogr. 2011;38:1267–80.

    Article  Google Scholar 

  32. 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.

    Article  CAS  PubMed  Google Scholar 

  33. Shaw J, Lickey EB, Schilling EE, Small RL. Comparison of whole chloroplast genome sequences to choose noncoding regions for phylogenetic studies in angiosperms: the tortoise and the hare III. Am J Bot. 2007;94:275–88.

    Article  CAS  PubMed  Google Scholar 

  34. Sang T, Crawford D, Stuessy T. Chloroplast DNA phylogeny, reticulate evolution, and biogeography of Paeonia (Paeoniaceae). Am J Bot. 1997;84:1120.

    Article  CAS  PubMed  Google Scholar 

  35. Li R, Ma P-F, Wen J, Yi T-S. Complete sequencing of five Araliaceae chloroplast genomes and the phylogenetic implications. PLoS One. 2013;8:e78568.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Edgar RC. MUSCLE: multiple sequence alignment with high accuracy and high throughput. Nucleic Acids Res. 2004;32:1792–979.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  37. Vaidya G, Lohman DJ, Meier R. SequenceMatrix: concatenation software for the fast assembly of multi-gene datasets with character set and codon information. Cladistics. 2010;27:171–80.

    Article  Google Scholar 

  38. Drummond AJ, Rambaut A. BEAST: Bayesian evolutionary analysis by sampling trees. BMC Evol Biol. 2007;7:214.

    Article  PubMed  PubMed Central  Google Scholar 

  39. Posada D. jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008;25:1253–6.

    Article  CAS  PubMed  Google Scholar 

  40. Rambaut A, Drummond AJ: Tracer v1.4. 2007.

  41. Szafer W. Miocenska Flora ze Starych Gliwic na Slasku (Miocene Flora of stare Gliwice in upper Silesia). Pr Inst Geol. 1961;33:1–205.

    Google Scholar 

  42. Barrón E, Postigo-Mijarra JM, Diéguez C. The late Miocene macroflora of the Cerdanya Basin (eastern Pyrenees, Spain): towards a synthesis. Paleontogr Abt B Paleobotany - Paleophytology. 2014;291:85–129.

    Google Scholar 

  43. Mai DH, Walther H. Die pliozaenen Floren von Thueringen, Deutsche Demokratische Republik. Quartaerpalaeontologie. 1988;7:55–297.

  44. Rim KH. Fossils of North Korea. Pyongyang: Science and Technology Press; 1994.

    Google Scholar 

  45. Müller J. Fossil pollen records of extant angiosperms. Bot Rev. 1981;47:1–142.

    Article  Google Scholar 

  46. Kong WS. The vegetational and environmental history of the pre-Holocene period in the Korean Peninsula. Korean J Quat Res. 1992;6:12.

    Google Scholar 

  47. Kolakovskii AA, Shakryl AK. Kimmeriyskaya flora Gul’ripsha (Bagazhishta). Tr Sukum Bot Sada. 1978;24:134–56.

    Google Scholar 

  48. Kvaček Z, Velitzelos D, Velitzelos E. Late Miocene Flora of Vegora Macedonia N. Greece: University of Athens, Greece; 2002.

    Google Scholar 

  49. Kovar-Eder J, Kvaček Z, Martinetto E, Roiron P: Late Miocene to Early Pliocene vegetation of southern Europe (7–4 Ma) as reflected in the megafossil plant record. Palaeogeogr Palaeoclimatol Palaeoecol 2006, 238:321–339. [Late Miocene to Early Pliocene Environment and Climate Change in the Mediterranean Area].

  50. Vieira M, Poças E, Pais J, Pereira D. Pliocene flora from S. Pedro da Torre deposits (Minho, NW Portugal). Geodiversitas. 2011;33:71–85.

    Article  Google Scholar 

  51. Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–90.

  52. Ricklefs RE. History and diversity: explorations at the intersection of ecology and evolution. Am Nat. 2007;170:S56–70.

    Article  PubMed  Google Scholar 

  53. Crisp MD, Cook LG. Explosive radiation or cryptic mass extinction? Interpreting signatures in molecular mhylogenies. Evolution. 2009;63:2257–65.

    Article  PubMed  Google Scholar 

  54. Fiz-Palacios O, Valcárcel V. From Messinian crisis to Mediterranean climate: a temporal gap of diversification recovered from multiple plant phylogenies. Perspect Plant Ecol Evol Syst. 2013;15:130–7.

    Article  Google Scholar 

  55. Rabosky DL. LASER: a maximum likelihood toolkit for detecting temporal shifts in diversification rates from molecular phylogenies. Evol Bioinforma. 2006;2:257–60.

    CAS  Google Scholar 

  56. Crawley MJ. The R book. Chichester: Wiley; 2007.

    Book  Google Scholar 

  57. R Development Core Team: R: A Language and Environment for Statistical Computing. Vienna, Austria; 2011.

  58. Faith DP. Conservation evaluation and phylogenetic diversity. Biol Con. 1992;61:1–10.

    Article  Google Scholar 

  59. Kembel SW, Cowan PD, Helmus MR, Cornwell WK, Morlon H, Ackerly DD, et al. Picante: R tools for integrating phylogenies and ecology. Bioinformatics. 2010;26:1463–4.

    Article  CAS  PubMed  Google Scholar 

  60. Clement M, Posada D, Crandall KA. TCS: A computer program to estimate gene genealogies. Mol Ecol. 2000;9:1657–9.

    Article  CAS  PubMed  Google Scholar 

  61. Crandall KA, Templeton AR. Empirical tests of some predictions from coalescent theory with applications to intraspecific phylogeny reconstruction. Genetics. 1993;134:959–69.

    CAS  PubMed  PubMed Central  Google Scholar 

  62. Posada D, Crandall KA. Intraspecific gene genealogies: trees grafting into networks. Trends Ecol Evol. 2001;16:37–45.

    Article  CAS  PubMed  Google Scholar 

  63. Li R, Wen J. Phylogeny and biogeography of Asian Schefflera (Araliaceae) based on nuclear and plastid DNA sequences data. J Syst Evol. 2014;52:431–49.

  64. Li R, Wen J. Phylogeny and diversification of Chinese Araliaceae based on nuclear and plastid DNA sequence data. J Syst Evol. 2016;54:453–67.

    Article  Google Scholar 

  65. Li R, Wen J. Phylogeny and biogeography of Dendropanax (Araliaceae), an Amphi-Pacific Disjunct genus between tropical/subtropical Asia and the Neotropics. Syst Bot. 2013;38:536–51.

  66. Mitchell A, Wen J. Phylogeny of Brassaiopsis (Araliaceae) in Asia based on nuclear ITS and 5S–NTS DNA sequences. Syst Bot. 2005;30:872–86.

  67. Plunkett GM, Wen J, Ii PPL. Infrafamilial classifications and characters in Araliaceae: insights from the phylogenetic analysis of nuclear (ITS) and plastid (trnL-trnF) sequence data. Plant Syst Evol. 2004;245:1–39.

  68. Wen J, Plunkett GM, Mitchell AD, Wagstaff SJ. The evolution of Araliaceae: a phylogenetic analysis based on ITS sequences of nuclear ribosomal DNA. Syst Bot. 2001;26:144–67.

    Google Scholar 

  69. Wen J, Nie Z-L, Soejima A, Meng Y. Phylogeny of Vitaceae based on the nuclear GAI1 gene sequences. Can J Bot. 2007;85:731–45.

  70. Axelrod D. Evolution and biogeography of Madrean–Tethyan sclerophyll vegetation. Ann Mo Bot Gard. 1975;62:280–334.

    Article  Google Scholar 

  71. Harris AJ, Xiang QY, Thomas DT: Phylogeny, origin, and biogeographic history of Aesculus L. (Sapindales) an update from combined analysis of DNA sequences, morphology and fossils. Taxon 2009, 58:1–19.

  72. Harris A, Wen J. Xiang Q-Y (Jenny): inferring the biogeographic origins of inter-continental disjunct endemics using a Bayes-DIVA approach. J Syst Evol. 2013;51:117–33.

    Article  Google Scholar 

  73. Wen J. Evolution of eastern Asian and eastern North American disjunct pattern in flowering plants. Annu Rev Ecol Syst. 1999;30:421–55.

    Article  Google Scholar 

  74. Latham RE, Ricklefs RE. Global patterns of tree species richness in moist forests: energy-diversity theory does not account for variation in species richness. Oikos. 1993;67:325–33.

    Article  Google Scholar 

  75. Qian H, Ricklefs RE. Large-scale processes and the Asian bias in temperate plant species diversity. Nature. 2000;407:180–2.

    Article  CAS  PubMed  Google Scholar 

  76. Krijgsman W. The onset of the Messinian salinity crisis in the eastern Mediterranean (Pissouri Basin, Cyprus). Earth Planet Sci Lett. 2002;194:299–310.

    Article  CAS  Google Scholar 

  77. Ivanov D, Ashraf AR, Mosbrugger V, Palamarev E. Palynological evidence for Miocene climate change in the Forecarpathian Basin (central Paratethys, NW Bulgaria). Palaeogeogr Palaeoclimatol Palaeoecol. 2002;178:19–37.

    Article  Google Scholar 

  78. Van Dam JA. Geographic and temporal patterns in the late Neogene (12.3 ma) aridification of Europe: the use of small mammals as paleoprecipitation proxies. Palaeogeogr Palaeoclimatol Palaeoecol. 2006;238:190–218.

    Article  Google Scholar 

  79. Thompson JD. Plant evolution in the Mediterranean. Oxford: Oxford University Press; 2005.

    Book  Google Scholar 

  80. Rodríguez-Sánchez F, Pérez-Barrales R, Ojeda F, Vargas P, Arroyo J: The Strait of Gibraltar as a melting pot for plant biodiversity. Quat Sci Rev 2008, 27:2100–2117. [The Coastal Shelf of the Mediterranean and Beyond: Corridor and Refugium for Human Populations in the Pleistocene].

  81. Jiménez-Moreno G, Fauquette S, Suc JP. Miocene to Pliocene vegetation reconstruction and climate estimates in the Iberian Peninsula from pollen data. Rev Paleobot Palynol. 2010;162:410–5.

    Article  Google Scholar 

  82. Herrera C. Historical effects and sorting processes as explanations for contemporary ecological patterns: character syndromes in Mediterranean woody plants. Am Nat. 1992;140:421–46.

    Article  Google Scholar 

  83. Verdú M, Dávila P, García-Fayos P, Flores-Hernández N, Valiente-Banuet A. “convergent” traits of Mediterranean woody plants belong to pre-Mediterranean lineages. Biol J Linn Soc. 2003;78:415–27.

    Article  Google Scholar 

  84. Qiu YX, Fu CX, Comes HP. Plant molecular phylogeography in China and adjacent regions: tracing the genetic imprints of quaternary climate and environmental change in the world’s most diverse temperate flora. Mol Phylogenet Evol. 2011;59:225–44.

    Article  PubMed  Google Scholar 

  85. Qian H. A comparison of the taxonomic richness of temperate plants in East Asia and North America. Am J Bot. 2002;89:1818–25.

    Article  PubMed  Google Scholar 

  86. Torrecilla P, López-Rodríguez JA, Catalán P. Phylogenetic relationships of Vulpia and related genera (Poeaae; Poaceae) based on analysis of ITS and trnL-F sequences. Ann Mo Bot Gard. 2004;91:124–58.

    Google Scholar 

  87. Oberprieler C. Temporal and spatial diversification of Circum-Mediterranean Compositae-anthemideae. Taxon. 2005;54:951–66.

    Article  Google Scholar 

  88. Inda LA, Torrecilla P, Catalán P, Ruiz-Zapata T. Phylogeny of Cleome L. and ITS close relatives Podandrogyne Ducke and Polanisia Raf. (Cleomoideae, Cleomaceae) based on analysis of nuclear ITS sequences and morphology. Plant Syst Evol. 2008;274:111–26.

    Article  CAS  Google Scholar 

  89. Mansion G, Rosenbaum G, Schoenenberger N, Bacchetta G, Rosselló JA, Conti E. Phylogenetic analysis informed by geological history supports multiple, sequential invasions of the mediterranean basin by the angiosperm family Araceae. Syst Biol. 2008;57:269–85.

    Article  CAS  PubMed  Google Scholar 

  90. Micó E, Sanmartín I, Galante E. Mediterranean diversification of the grass-feeding Anisopliina beetles (Scarabaeidae, Rutelinae, Anomalini) as inferred by bootstrap-averaged dispersal-vicariance analysis. J Biogeogr. 2013;36:546–60.

    Article  Google Scholar 

  91. Ree RH, Sanmartín I. Prospects and challenges for parametric models in historical biogeographical inference. J Biogeogr. 2009;36:1211–20.

    Article  Google Scholar 

  92. Sanmartín I. Dispersal vs. vicariance in the Mediterranean: historical biogeography of the Palearctic Pachydeminae (Coleoptera, Scarabaeoidea). J Biogeogr. 2003;30:1883–97.

    Article  Google Scholar 

  93. Lo Presti RM, Oberprieler C. Evolutionary history, biogeography and eco-climatological differentiation of the genus Anthemis L. (Compositae, Anthemideae) in the circum-Mediterranean area. J Biogeogr. 2009;36:1313–32.

  94. Costa JM, Ramos JA, da Silva LP, Timoteo S, Araújo PM, Felgueiras MS, et al. Endozoochory largely outweighs epizoochory in migrating passerines. J Avian Biol. 2014;45:59–64.

    Article  Google Scholar 

  95. Hampe A, Arroyo J, Jordano P, Petit R. Rangewide phylogeography of a bird-dispersed Eurasian shrub: contrasting Mediterranean and temperate glacial refugia. Mol Ecol. 2003;12:3415–26.

    Article  CAS  PubMed  Google Scholar 

  96. Dubreuil M, Riba M. González–Martínez SC, Vendramin GG, Sebastiani F, Mayol M. Genetic effects of chronic habitat fragmentation revisited: strong genetic structure in a temperate tree, Taxus baccata (Taxaceae), with great dispersal capability. Am J Bot. 2010;97:303–10.

  97. García C, Jordano P, Arroyo JM, Godoy JA. Maternal genetic correlations in the seed rain: effects of frugivore activity in heterogeneous landscapes. J Ecol. 2009;97:1424–35.

    Article  Google Scholar 

  98. Krijgsman W. The Mediterranean: Mare nostrum of earth sciences. Earth Planet Sci Lett. 2002;205:1–12.

    Article  CAS  Google Scholar 

  99. Ring U, Laws S, Bennet M. Structural analysis of a complex nappe sequence and late-orogenic basins from the Aegean Island of Samos, Greece. J Struct Geol. 1999;21:1575–601.

    Article  Google Scholar 

  100. Escudero M, Martín-Bravo S, Mayrose I, Fernández-Mazuecos M, Fiz-Palacios O, Hipp AL, et al. Karyotypic changes through dysploidy persist longer over evolutionary time than polyploid changes. PLoS One. 2014;9:e85266.

    Article  PubMed  PubMed Central  Google Scholar 

Download references


Part of the data was obtained with the financial support of a grant to V. Valcárcel from Junta de Andalucía.

Availability of data and materials

The sequence datasets generated during the current study are available in the GenBank database repository,

The Lagrange splits generated during the current study and needed to run the new R script herein provided is available in ResearchGate database repository,

Authors’ contributions

VV obtained half of the data, performed the biogeographic, diversification and phylogeographic analyses, helped developing the R-script, interpreted all results and lead the writing process. BG obtained half of the data, performed the divergence age analysis, wrote the methodology for this analysis in the manuscript and contributed to the results interpretation and writing of the manuscript. NGM developed the R-script for summarizing the multi-tree results from biogeographic analysis and contributed to the writing of the manuscript. PV and JW led the elaboration of the hypothesis, provided funding for obtaining the data and were both major contributors to help interprete the results and improve the writing of the manuscript. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Consent for publication

Not applicable.

Ethics approval and consent to participate

Not applicable.

Author information

Authors and Affiliations


Corresponding author

Correspondence to V. Valcárcel.

Additional files

Additional file 1:

List of the studied material included in the phylogenetic-based analyses. Localities, geographic area codification and GenGank accession numbers of the nuclear Internal Transcribed Spacer are provided. Areas abbreviation are as follows: A, Neotropics; B, Tropical Africa; C, Asutralia: D, W Mediterranean; E, E Mediterranean; F, Europe; G, W Asia; H, E Asia. Papers of reference are provided as superscript as follows: (1) Li et al. [35], (2) Mitchell et al. [26], (3) Mitchell et al. [66], (4) Valcárcel et al. [25], (5) Vargas et al. [7], and (6) Valcárcel et al. [16]. (DOCX 141 kb)

Additional file 2:

GenBank accession numbers of the studied material included in the phylogeographic study. (DOCX 171 kb)

Additional file 3:

Summary of DNA sequences variation and evolutionary models best fitting the nuclear and plastid matrices. (DOCX 41 kb)

Additional file 4:

R script for parsing Lagrange results from multi-tree analyses. (R 22 kb)

Additional file 5:

Beast Maximum Clade Credibility chronogram of the nrITS dataset of Araliaceae with the secondary calibration approach. Legend: Mean ages and 95% CI are only represented for clades with >0.5 Posterior Probability support. (PDF 255 kb)

Additional file 6:

Beast Maximum Clade Credibility chronogram of the nrITS dataset of Araliaceae with the fossil calibration approach. Mean ages and 95% CI are only represented for clades with >0.5 Posterior Probability support. (PDF 236 kb)

Additional file 7:

Fitness to speciation models. Summary of the results of the fitness to speciation models estimated from 6500 posterior pruned trees of the Hedera clade. N indicates the number of trees that recover a given evolutionary model. Mean values and 95% CI are provided for AIC and each of the parameters of the models. (DOCX 49 kb)

Additional file 8:

List of the haplotypes detected in the samples included in the phylogeographic study. Taxa name and general distribution is provided. Locality and voucher is specified for each sample as well as the number of haplotype detected individually for the three regions (rpL32, trnH-psbA, trnT-trnL). Last column indicates the number of haplotype when combining the three-plastid regions rpL32, trnH-psbA and trnT-trnL (HP3). (DOCX 106 kb)

Additional file 9:

Hedera phylogeographic network obtained from the Statistical Parsimony analysis of the rpL32 plastid region. Haplotype numbering is according to Additional file 8. Circle dimensions are proportional to the number of samples displaying each haplotype. Kalopanax septembolus is used as outgroup. (PDF 96 kb)

Additional file 10:

Hedera phylogeographic network obtained from the Statistical Parsimony analysis of the trnH-psbA plastid region. Haplotype numbering is according to Additional file 8. Circle dimensions are proportional to the number of samples displaying each haplotype. Kalopanax septembolus is used as outgroup. (PDF 93 kb)

Additional file 11:

Hedera phylogeographic network obtained from the Statistical Parsimony analysis of the trnT-trnL plastid region. Haplotype numbering is according to Additional file 8. Circle dimensions are proportional to the number of samples displaying each haplotype. Kalopanax septembolus is used as outgroup. (PDF 108 kb)

Additional file 12:

Hedera phylogeographic network obtained from the Statistical Parsimony analysis of the three-plastid regions (trnH-psbA, trnT-trnL, rpL32) and using seven Araliaceae genera as outgroup. Circle dimensions are proportional to the number of samples displaying each haplotype. Abbreviation names of the taxa displaying each haplotype are as follows: ALG, H. algeriensis; AZO, H. azorica; CAN, H. canariensis; COL, H. colchica; CYP, H. pastuchowii subsp. cypria; HEL, H. helix; HIB, H. hibernica; IBE, H. iberica; MAD, H. maderensis; MAR, H. maroccana; NEP, H. nepalensis var. nepalensis; PAS, H. pastuchowii subsp. pastuchowii; RHO, H. rhombea; SIN, H. nepalensis var. sinensis. (PDF 188 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (, which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Valcárcel, V., Guzmán, B., Medina, N.G. et al. Phylogenetic and paleobotanical evidence for late Miocene diversification of the Tertiary subtropical lineage of ivies (Hedera L., Araliaceae). BMC Evol Biol 17, 146 (2017).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: