How the temperate world was colonised by bindweeds: biogeography of the Convolvuleae (Convolvulaceae)
BMC Evolutionary Biology volume 16, Article number: 16 (2016)
At a global scale, the temperate zone is highly fragmented both between and within hemispheres. This paper aims to investigate how the world’s disjunct temperate zones have been colonised by the pan-temperate plant group Convolvuleae, sampling 148 of the c. 225 known species. We specifically determine the number and timing of amphitropical and transoceanic disjunctions, investigate the extent to which disjunctions in Convolvuleae are spatio-temporally congruent with those in other temperate plant groups and determine the impact of long-distance dispersal events on diversification rates.
Eight major disjunctions are observed in Convolvuleae: two Northern Hemisphere, two Southern Hemisphere and four amphitropical. Diversity in the Southern Hemisphere is largely the result of a single colonisation of Africa 3.1–6.4 Ma, and subsequent dispersals from Africa to both Australasia and South America. Speciation rates within this monophyletic, largely Southern Hemisphere group (1.38 species Myr−1) are found to be over twice those of the tribe as a whole (0.64 species Myr-1). Increased speciation rates are also observed in Calystegia (1.65 species Myr−1).
The Convolvuleae has colonised every continent of the world with a temperate biome in c. 18 Myr and eight major range disjunctions underlie this broad distribution. In keeping with other temperate lineages exhibiting disjunct distributions, long-distance dispersal is inferred as the main process explaining the patterns observed although for one American-Eurasian disjunction we cannot exclude vicariance. The colonisation of the temperate zones of the three southern continents within the last c. 4 Myr is likely to have stimulated high rates of diversification recovered in this group, with lineage accumulation rates comparable to those reported for adaptive radiations.
The successful colonisation of temperate biomes by tropical lineages has involved the crossing of a significant physiological barrier that has acted as an important filter . As a consequence, approximately half of all plant families remain restricted to the tropics . Lineages that have made the transition to temperate biomes have experienced different fates with some lineages expanding their ranges to occupy highly disjunct areas where a suitable climate occurs, with distributions spanning both different continents and different hemispheres. Whilst the processes responsible for such patterns are complex [3–6], long-distance dispersal (LDD) events have been proposed for many such disjunctions and they may have acted as triggers for diversification [7, 8].
Thorne  recognised fifteen temperate disjunction patterns, several of which have since been the focus of molecular phylogenetic studies to understand the extent to which vicariance and dispersal explain biogeographic patterns in temperate plant lineages (e.g., the eastern North American–East Asian disjunction ; the temperate North and South American disjunction ; the western North American–East Asian disjunction ). One of the patterns documented by Thorne  was the ‘North–South Temperate disjunction’ which describes groups that are widespread in the northern temperate region and that also occur in one or more of the southern temperate zones (i.e., those located in South America, Africa and Australasia).
The Convolvuleae (Choisy) Choisy is one of 12 tribes within the predominantly tropical plant family Convolvulaceae Juss. . It is one of only two tribes within the family to have widely colonised temperate regions (the other being the parasitic Cuscuteae) and it exhibits Thorne’s ‘North–south Temperate disjunction’ pattern. The Convolvuleae comprises three genera namely Convolvulus L., Calystegia R.Br. and Polymeria R.Br.. Convolvulus is the largest, comprising 190 species . It has a main centre of diversity in the Mediterranean and western Asia, with further centres of diversity in eastern Asia and in temperate South America, southern and eastern Africa and Australasia; i.e., the three temperate zones of the southern hemisphere. Species also occur in North America, although they are few in number. Calystegia is readily distinguished from Convolvulus based on morphological characters (namely polypantoporate pollen and stigma shape) but molecular analyses suggest it is nested within the larger Convolvulus clade [13, 15, 16]. Calystegia is taxonomically complex  with c. 26 species and more than 65 distinct taxa currently accepted . The centre of diversity for Calystegia is in California where nearly half of the described taxa occur . Other centres of diversity for Calystegia are found in eastern Asia and, to a lesser extent Europe and the Mediterranean. Calystegia also occurs in temperate regions of the Southern Hemisphere. Finally, the Australasian endemic Polymeria is the smallest of the three genera of Convolvuleae with eight species recognised . Molecular analyses place it as sister group to the remainder of Convolvuleae [13, 15, 16].
A recent study by Williams et al.  established a robust phylogenetic hypothesis of the Convolvuleae that sampled 62 % of species diversity in the tribe and was based on data from both the nuclear ITS region and the chloroplast matK and rbcL regions. The goal of this paper is to utilise that phylogenetic framework to determine how the North–South Temperate disjunction pattern displayed by Convolvuleae was generated. Specifically, we aim to (i) determine the number, timing and cause (dispersal versus vicariance) of amphitropical and transoceanic disjunctions in the pan-temperate Convolvuleae and (ii) determine how major disjunctions in the history of the group may have impacted on diversification rates.
An alignment comprising 153 species of Convolvulaceae (of which 109 were Convolvuleae; eight Polymeria, 11 Calystegia, 90 Convolvulus) and 343 species of Solanaceae (126 Solanoideae) and 1328 characters from the matK and rbcL regions (of which 538 were parsimony informative) was used to establish divergence times within Convolvulaceae. The rbcL region was coded with missing data for 241 taxa, of which four were Convolvuleae. A chronogram with major groups is summarised in Additional file 1. The Convolvulaceae are resolved to have arisen 44.1 (95 % HPD 33.9–51.2) Ma, in agreement with Särkinen et al. . Age estimates established for nodes within Solanaceae are also in agreement with Särkinen et al. . Within Convolvulaceae, the Convolvuloideae sensu Stefanović et al.  is resolved as 20.9 (14.3–27.5) Myr old, with the Convolvuleae crown group (corresponding to the split between Polymeria and the Convolvulus + Calystegia clade) resolved at 17.9 (11.8–23.7) Ma.
Five areas of endemism within Convolvuleae were delimited using UPGMA clustering of species by country distribution data (Fig. 1a). A Convolvuleae alignment consisted of 148 species of Convolvuleae (11 Polymeria, 18 Calystegia, 119 Convolvulus) and 2033 characters from the rbcL, matK and ITS regions (matrix deposited in TreeBASE, study 18623). Divergence times estimated in beast using calibration points derived from the Solanaceae-Convolvulaceae analysis above and ancestral area reconstructions estimated using lagrange are provided in Fig. 1b with Table 1 summarising the information for key nodes of interest.
The ancestral area of the Convolvulus + Calystegia clade is inferred to be the Mediterranean-and-Middle-East (area A in Fig. 1a; node 2). Dispersal between contiguous areas (i.e., A–B, A–C, C–D (Fig. 1a)) occurred frequently (Fig. 1b). Movement between disjunct (i.e., non-contiguous) areas has been much less common and eight such events are inferred. These are, in order of recency: (i) amphitropical disjunction between Australasia (Area G) and the Mediterranean-and-Middle-East (Area A) (posterior probability for node (PP) = 1) dated 17.61 Ma (95 % highest posterior density (HPD): 13.50–21.56 Ma) (node 1; Fig. 1b; Table 1), (ii) Northern hemisphere disjunction between the Mediterranean-and-Middle-East (A) and North America (F) (PP = 1) dated 5.92 Ma (3.53–8.6 Ma) (node 16), (iii) Northern hemisphere disjunction between Central-and-North-Eastern-Asia (B) and North America (F) (PP = 1), dated 5.56 Ma (3.45–8.26 Ma) (node 12), (iv) amphitropical disjunction between the Mediterranean-and-Middle-East (A) and Southern-and-Eastern-Africa (D) (PP = 1) dated 4.62 Ma (3.12–6.41 Ma) (node 5), (v) transoceanic southern hemisphere disjunction between Southern-and-Eastern-Africa (D) and Australasia (G) (PP = 1), dated 3.06 Ma (1.99–4.41 Ma) (node 6), (vi) weakly supported transoceanic southern hemisphere disjunction between Southern-and-Eastern-Africa (D) and South America (E) (PP = 0.52) dated as 1.15 Ma (0.65–1.78 Ma) (node 8), (vii) weakly supported (PP = 0.59) amphitropical disjunction between South America (E) and North America (F), dated 0.57 Ma (0.18–1.03 Ma) (node 9) and (viii) weakly supported (PP = 0.79) amphitropical disjunction between North America (F) and the Southern Hemisphere (Australasia (G), South America (E) or both, node 15), dated 0.82 Ma (0.3–1.53 Ma)
Diversification rate analysis
BAMM analysis found support for rate shifts within the Convolvuleae phylogeny compared with a single rate model (Bayes factors (BF) > 30 for 2–6 shifts). The maximum a posteriori (MAP) probability rate shift configuration, which alone explains 56 % of the data, infers two rate shifts: one on the stem branch of the southern hemisphere clade (Fig. 2; group A; BF = 755) and one on the stem branch of Calystegia (Fig. 2; group B; BF = 424) (Additional file 2a). Mean speciation rates (λ) within the circum-South Temperate (cST) clade (1.38, 90 % HPD: 0.71–2.07) and Calystegia (1.65, 90 % HPD: 0.61–2.73) are over twice those of the tribe as a whole (0.64, 90 % HPD: 0.5–0.83). Extinction rates (μ) however are also slightly higher in both the cST clade (0.47, 90 % HPD: 0.04–1.19) and Calystegia (0.86, 90 % HPD: 0.14–2.07) than in Convolvuleae as a whole (0.31, 90 % HPD: 0.11–0.58). Mean diversification rates (λ - μ) in the cST clade (0.91 species Myr−1) are therefore nearly three times greater than in the tree as a whole (0.34 species Myr−1), while those in Calystegia are over twice as fast (0.8 species Myr−1) as Convolvuleae in general. Extinction rates are inferred to have remained fairly constant over the history of the tribe, however speciation rates appear to have increased considerably in the last 2.5 Myr (Additional file 2b).
Both the world’s oceans and the equatorial tropics present barriers to dispersal of temperate lineages, potentially limiting exchange between the disjunct temperate zones of the world. In Convolvuleae, we observe four amphitropical and four transoceanic disjunctions in the history of the group with an increase in diversification rate associated with one amphitropical disjunction (the main group of Convolvulus in the southern hemisphere). A second increase in diversification rate is observed in Calystegia which also exhibits a transoceanic disjunction.
The four amphitropical disjunctions are spread throughout the history of the tribe (Fig. 1b, nodes 1, 5, 9 and 16). The earliest dates to the mid-Miocene (17.61 (13.5–21.56) Ma; node 1; Fig. 1b; Table 1) and separates the Australasian endemic Polymeria from the remainder of Convolvuleae (Convolvulus + Calystegia), for which the Mediterranean-and-Middle-East is resolved as the ancestral distribution area. This spatio-temporal pattern is consistent with the inferred timing of the disjunction between Australia and Eurasia in Carex subsect. Spirostachyae (c. 16–26 Ma ) and in Halosarcia (c. 15–20 Ma ). Escudero et al.  invoked LDD to explain this disjunction although the timing is coincident with the mid-Miocene Climatic Optimum 15–17 Ma, which saw the expansion of tropical forests, an event that is thought to have facilitated the dispersal of tropical plant and animal groups between Africa and Asia . It is plausible that a corresponding contraction of temperate areas may have led to the disjunction apparent in these groups. Särkinen et al.  resolved a similar sister group relationship in Solanum between the Western Mediterranean–Macaronesian endemic Normania clade and the Australasian endemic Archaesolanum clade and an Australian–Northern Hemisphere disjunction was also inferred for Atriplex . However, the timing of these was more recent (8.3 Ma and 9.8–7.8 Ma respectively) suggesting that the history of Australasian–Northern Hemisphere disjunctions is complex with multiple, temporally distinct events likely involved.
A second amphitropical disjunction in Convolvuleae is the result of the colonisation of the Southern Hemisphere by Convolvulus during the late Miocene to Pliocene c. 4.62 (3.12–6.41) Ma (node 5; Fig. 1b; Table 1). The most probable scenario involves dispersal from the Mediterranean-and-Middle-East into Southern-and-Eastern-Africa, followed by dispersal from there to Australasia c. 3.06 (1.99–4.41) Ma (node 6) and South America c. 1.15 (0.65–1.78) Ma (node 8) although the precise relationships of African and American taxa are not well supported. The mountains of the East African rift system, which link Southern Africa with the Horn of Africa are thought to have originated c. 12–40 Ma  and they provide a plausible trans-African dispersal corridor for Convolvuleae as has been suggested for other temperate taxa (e.g., Senecio ; Disa, Irideaeae p.p., Pentaschistis, Restionaceae ; Androcymbium ; Scabiosa ). It is notable that whilst some African amphitropical disjunctions are inferred to be pre-Pliocene (e.g., Juniperus, 30.5 (14.0–47.0) Ma ; Hyacinthoideae, 18.7 (18.8–18.7) Ma ; Thamnosma, 8.7 (5.3–12.1) Ma ), Plio– Pleistocene disjunctions, consistent with that observed in Convolvuleae have been reported in a number of groups. For example, colonisation of southern Africa from the north through this corridor has been inferred in Apium (4.1 (1.2–7.0) Ma ), Ranunculus (3.9 (2.6–5.3) Ma ) and Scabiosa (1.6 (0.7–2.6) Ma ) whilst South to North colonisation has been inferred for Androcymbium (4.0 (2.5–5.5) Ma and 3.0 (1.5–4.5) Ma ) and, very recently in Senecio (0.2 (0–0.4) Ma ).
The remaining amphitropical disjunctions are observed in the New World. Pleistocene dispersal (0.57 (0.18–1.03) Ma) from South America to North America is inferred in the cST group (node 9; Fig. 1b; Table 1). A southwards dispersal in Calystegia (node 15) is inferred to have occurred at a similar time c. 0.82 (0.3–1.53) Ma resulting in the colonisation of South America and Australasia. However, limited support (Convolvulus; PP = 0.59, Calystegia; PP = 0.79) or taxonomic uncertainty (Calystegia) means that these patterns should be interpreted with caution and the evolution and biogeography of the Calystegia clade in particular would benefit from further research.
These limitations notwithstanding, amphitropical American disjunctions of recent origin have been inferred in a range of groups with evidence for dispersal in both directions . Bird mediated dispersal has frequently been proposed as responsible for such disjunctions due to the seasonal migration of birds between the Northern and Southern hemispheres (e.g., [11, 36, 37]). Whilst evidence for this is largely anecdotal, epizoochoric bird-mediated LDD between California and Chile has been demonstrated in Lepidium . In the case of Convolvulus, viable seeds of Convolvulus arvensis have been recovered from the digestive tract of migratory killdeer (Charadrius vociferus) up to six days after ingestion . Importantly however, long-distance internal transport of seeds, even in generally larger waterbirds has been shown to be limited to around 300 km making extreme long-distance endozoochoric dispersal unlikely . Montane South American species of Convolvulus such as C. crenatifolius and C. montanus are frequently found above 1500 m  and the Andean high mountains, which are of late Miocene origin , may have provided a suitable route for the dispersal of temperate Convolvuleae lineages across the neotropics.
Remarkably few transoceanic dispersal events are necessary to explain the global distribution of the Convolvuleae, in contrast to groups such as Fabaeae . Between the major Southern Hemisphere landmasses, we infer only two such dispersal events which both occurred 0.65–4.41 Ma, long after the breakup of the Gondwanan landmass and too recent to involve an Antarctic corridor . Divergence time estimates for southern temperate plant groups indicate a wide range of ages, with Convolvuleae disjunctions among the most recent . Whilst detailed information on dispersal mechanisms within the tribe are lacking, the variability of seed characteristics in Convolvulaceae  and evidence of long-distance oceanic seed dispersal elsewhere in the family  anecdotally support an oceanic dispersal hypothesis for the Southern Hemisphere distribution of Convolvuleae.
Two independent events in the late Miocene/Pliocene are inferred to be responsible for the earliest colonisation of North America. The most likely scenario for the Calystegia clade is dispersal from the Mediterranean-and-Middle-East into Eastern Asia c. 7.85 Ma (node 11; Fig. 1b; Table 1) and from there into North America c. 5.56 (3.45–8.26) Ma (node 12). The East Asian–North American disjunction is one of the best studied disjunctions with numerous examples of movement between the two regions throughout the Cenozoic . The continents of the Northern Hemisphere were connected until 5.4–5.5 Ma when the Bering Land Bridge joining North America and Eastern Asia was finally severed . A circum-Arctic floral region spanning this landmass prior to the severing of the land connections is frequently hypothesised as responsible for both the similarity and diversity of the flora in these regions (e.g., [9, 31, 46–48]). Given our estimated age for the Northern Hemisphere disjunction in Convolvuleae we are unable to reject a vicariance hypothesis for the origin of Calystegia in North America, in contrast to all other disjunctions we have inferred.
In the case of Convolvulus simulans (node 16), dispersal directly from the Mediterranean-and-Middle-East to North America c. 5.92 (3.53–8.6) Ma is inferred. A disjunction between the Mediterranean regions of North America and Europe (Madrean–Tethyan) is well documented (see ) and long distance dispersal to North America from the Mediterranean during the late Miocene/Pliocene has been inferred in a number of lineages (e.g., Exaculum/Schenkia–Zeltnera c. 9 Ma , Crocanthemum–Hudsonia c. 5.2–9.2 Ma , Eobassia/Chenolea/Spirobassia–Neokochia c. 8.8–13.1 Ma ).
With regards widespread and naturalised Convolvulus that were excluded from our analyses, comparison between the phylogeny estimated in Williams et al.  and our biogeographic inference suggests that both excluded species (C. arvensis and C. lineatus) probably originated in the Mediterranean-and-Middle-East region.
The diversification rate analysis reveals mean diversification rates for Convolvuleae of 0.34 species Myr−1, above the estimated diversification rates of angiosperms as a whole (0.077–0.089 species Myr−1, . Furthermore, two shifts to increased diversification rates are supported within the tribe, with strong support for a rate shift in the southern hemisphere clade (Fig. 2; group A) leading to mean diversification rates (0.91 species Myr−1) which exceed those of adaptive radiations such as the Hawaiian radiation of Bidens (0.3–0.8 species Myr−1) . BAMM suggests the elevated diversification rates are linked to an increase in speciation rate as opposed to a decrease in extinction rate (Table 2). The southern hemisphere clade contains at least two long-distance oceanic dispersal events within the Southern Hemisphere (nodes 6 and 8; Fig. 1b; Table 1) and at least two amphitropical dispersal events (nodes 5 and 9), with the shift to elevated diversification rates associated with the initial dispersal into the southern hemisphere in the late Miocene or Pliocene. This is consistent with other studies demonstrating the impact of Miocene dispersal events important in promoting diversification [8, 53].
The second diversification rate shift is observed in Calystegia. BAMM finds support, albeit less strongly, for a shift to increased diversification rates on the stem branch of Calystegia (Fig. 2; group B), leading to mean diversification rates (0.8 species Myr−1), over twice as high as those found in Convolvuleae as a whole. Most of the diversity of Calystegia is in North America and specifically California  and dispersal from East Asia into North America, again in the Miocene/Pliocene could also have been an important trigger for diversification within the group.
In summary, our results indicate that the Convolvuleae has successfully colonised every continent of the world with a temperate biome in c. 18 Myr. The tropics and major oceans have been significant dispersal barriers for the group with only eight major disjunctions underlying this broad ‘North–south temperate’ distribution pattern. In keeping with many other disjunct temperate lineages, long-distance dispersal is inferred as the main process explaining the patterns observed although for one American-Eurasian disjunction we cannot exclude vicariance resulting from the severing of the Bering Land Bridge. Even though dispersal is the primary process generating the patterns observed, spatio-temperal congruence is observed with other temperate disjunct groups suggesting a common explanation for the patterns observed. The Convolvuleae exhibits high diversification rates overall when compared to other angiosperm groups and the colonisation of the temperate zones of the three southern continents within the last c. 4 Myr is associated with an increase in diversification rate with lineage accumulation rates in the clade comparable to those reported for adaptive radiations. The Calystegia clade also exhibits a high diversification rate which probably reflects rapid diversification following colonisation of western North America. The Convolvuleae thus provide a striking example of the ability of temperate lineages to rapidly colonise highly disjunct areas worldwide and to diversify.
Divergence time estimation
The fossil record of Convolvulaceae is poor with none of the fossils assigned to the family (Convolvulites orichitus , Tricolpites trioblatus , Calystegiapollis microechinatus (in ) able to be accurately placed within a phylogeny. We therefore adopted a two-step calibration procedure. We first utilised a recent phylogenetic study of Solanaceae, the sister group to Convolvulaceae  which reviewed all 50 of the known fossils assigned to the family, as the basis for calibration points for divergence time estimation within Convolvulaceae based on chloroplast data. Second, node age estimates from the chloroplast phylogeny were used to calibrate a combined nuclear ITS and plastid matK and rbcL phylogeny of Convolvuleae.
matK and rbcL sequences for Convolvuleae from Williams et al.  were manually aligned with sequences of the same regions for taxa across the remainder of the Convolvulaceae and Solanaceae, which were retrieved from GenBank. Details of all accessions sampled are included in Additional file 3. The datasets were concatenated, with taxa lacking matK sequences excluded and taxa lacking rbcL sequences coded with missing data for this region. Due to the lower levels of variation in the rbcL region  the missing data is unlikely to have any significant impact on the tree topology, as it will be overridden by the signal from the matK region .
Following Särkinen et al. , we used two calibration points reflecting the youngest age estimates of the oldest assignable fossils to constrain (i) the stem node of Solanoideae with a lognormal offset of 23.0 Ma, mean of 0.01, and standard deviation (SD) of 1.0 and (ii) the Solanaceae stem node with a lognormal offset of 46.0 Ma, mean of 0.01, and SD of 1.0. A gamma distribution (shape 0.001, scale 1000) was used as a prior for the mean mutation rate. Bayesian time estimation with an uncorrelated lognormal relaxed clock model was implemented in beast v1.8 . Two independent Markov Chain Monte Carlo (MCMC) runs of 200 million generations, sampling every 10,000 generations were conducted using a Speciation: Birth–Death process tree prior and the GTR + I + G model. A run sampling only from the prior probabilities was also performed to evaluate the performance of the priors. Mixing of the chains and convergence were assessed using Tracer v1.6  as was confirmation of an effective sample size (ESS) in the post burn-in samples exceeding 200 for all estimated parameters. The output tree files were combined using LogCombiner v1.8 (part of the beast software package) discarding the first 10 % of trees of each run as burn in. TreeAnnotator v1.8 (part of the beast software package) was used to combine post burn-in trees from the two runs, calculate the maximum clade credibility tree and the mean 95 % higher posterior density (HPD) intervals of node ages. Final trees were edited in FigTree v 1.4.0 .
A second divergence time analysis was performed on a concatenated ITS, matK and rbcL dataset modified from Williams et al.  since the ITS region included greater taxon sampling within Convolvuleae (see Additional file 3 for sampling details). Due to the separate modes of evolution, the manually aligned matrix was partitioned into nuclear and plastid regions and parameters estimated independently. Analysis in beast followed the protocol for the Convolvulaceae–Solanaceae analysis except that minimum age estimates from the aforementioned analysis were used to constrain the Convolvuleae root node and Polymeria, Convolvulus + Calystegia and Calystegia crown nodes with normally distributed prior at 17.89 Ma (SD = 3.0), 4.76 Ma (SD = 1.5), 15.18 Ma (SD = 2.5), and 2.7 Ma (SD = 0.8) respectively, with the distribution reflecting the 95 % HPD estimates, and MCMC runs were reduced to 20 million generations, sampling every 1000 generations.
The extant distributions of all accepted taxa at a country level were collated, largely from Wood et al. . Widespread taxa and those for which the natural distribution may have been obscured by frequent introductions/naturalisations (Convolvulus arvensis, Convolvulus lineatus, Calystegia soldanella, Calystegia pulchra, Calystegia sepium subsp. sepium, Calystegia sepium subsp. roseata and all Calystegia silvatica subspecies) were excluded as were countries with only a single taxon present. Areas of endemism were then delimited using Unweighted Pair-Group Method with Arithmetic Mean (UPGMA) clustering of a taxon × country distribution matrix using the Sørensen–Dice coefficient [61, 62] in DendroUPGMA . Seven areas were delimited (Fig. 1a) and each taxon was coded as belonging to one or more of these regions. Given our use of country borders as opposed to ecological boundaries to delimit areas, we only considered range shifts between non-contiguous regions as disjunctions.
The historical biogeography of the Convolvuleae was reconstructed using the dispersal–extinction–cladogenesis (DEC) model implemented in Lagrange v 2.0.1 , with taxon distributions coded as above. Three time slices were incorporated into the DEC model reflecting the presence or absence of Northern Hemisphere land connections between the Old and New Worlds. Subsequently, Northern America was isolated from the rest of the Northern Hemisphere 0–5.5 Ma, a connection existed between North America and Asia via the Bering Land Bridge (BLB) 5.5–15 Ma, and a connection also existed between North America and Europe via the North Atlantic Land Bridge (NALB) 15–25 Ma. All possible area combinations were permitted throughout. Dispersal probabilities followed Mao et al. , however given the relatively young age of Convolvuleae, and the subsequent reduction in major continental movements we simplified the model from five to three dispersal probabilities: 1.0 for connected areas, 0.1 for widely disjunct areas and 0.5 for three combinations of narrowly disjunct areas (area A–area D, area B–area F and area E–area F, Fig. 1b).
Diversification rate analysis
Bayesian Analysis of Macroevolutionary Mixtures (BAMM) v2.0  was used to model the dynamics of speciation and extinction on the time-calibrated Convolvuleae phylogenetic tree. Incomplete and non-random taxon sampling was incorporated directly into the likelihood calculations by utilising the recent monograph of Convolvulus  to place missing taxa into their respective clades. Two independent BAMM metropolis-coupled MCMC (MCMCMC) runs, with three heated and one cold chain, were run for 10 million generations and sampled every 1000 generations. Convergence of BAMM runs was assessed by computing ESS of log-likelihoods and numbers of shifts using the CODA library for R: both parameters had effective sample sizes > 1000. The first 10 % of samples were discarded as burn-in. Post-run analysis and visualisation was performed using the R package BAMMtools v2.0 .
Donoghue MJ. A phylogenetic perspective on the distribution of plant diversity. Proc Natl Acad Sci U S A. 2008;105:11549–55.
Ricklefs RE, Renner SS. Species richness within families of flowering plants. Evolution. 1994;48:1619–36.
Donoghue MJ, Smith SA. Patterns in the assembly of temperate forests around the Northern Hemisphere. Philos Trans R Soc London [Biol]. 2004;359(1450):1633–44.
Wen J, Ickert-Bond SM. Evolution of the Madrean–Tethyan disjunctions and the North and South American amphitropical disjunctions in plants. J Syst Evol. 2009;47:331–48.
Kadereit JW, Baldwin BG. Western Eurasian–western North American disjunct plant taxa: the dry-adapted ends of formerly widespread north temperate mesic lineages—and examples of long-distance dispersal. Taxon. 2012;61:3–17.
Manos PS, Meireles JE. Biogeographic analysis of the woody plants of the Southern Appalachians: implications for the origins of a regional flora. Am J Bot. 2015;102(5):780–804.
Knope ML, Morden CW, Funk VA, Fukami T. Area and the rapid radiation of Hawaiian Bidens (Asteraceae). J Biogeogr. 2012;39:1206–16.
Linder HP, Rabosky DL, Antonelli A, Wüest RO, Ohlemüller R. Disentangling the influence of climatic and geological changes on species radiations. J Biogeogr. 2014;41:1313–25.
Thorne RF. Major disjunctions in the geographic ranges of seed plants. Q Rev Biol. 1972;47:365–411.
Wen J, Ickert-Bond S, Nie Z-L, Li R. Timing and modes of evolution of eastern Asian–North American biogeographic disjunctions in seed plants. In: Long M, Gu H, Zhou Z, editors. Darwin’s Heritage Today: Proceedings of the Darwin 200 Beijing International Conference. Beijing: Higher Education Press; 2010. p. 252–69.
Popp M, Mirré V, Brochmann C. A single Mid-Pleistocene long-distance dispersal by a bird can explain the extreme bipolar disjunction in crowberries (Empetrum). Proc Natl Acad Sci U S A. 2011;108:6520–5.
Nie Z-L, Wen J, Sun H, Bartholomew B. Monophyly of Kelloggia Torrey ex Benth. (Rubiaceae) and evolution of its intercontinental disjunction between western North America and eastern Asia. Am J Bot. 2005;92:642–52.
Stefanović S, Krueger L, Olmstead RG. Monophyly of the Convolvulaceae and circumscription of their major lineages based on DNA sequences of multiple chloroplast loci. Am J Bot. 2002;89:1510–22.
Wood JRI, Williams BRM, Mitchell TC, Carine MA, Harris DJ, Scotland RW. A foundation monograph of Convolvulus L. (Convolvulaceae). PhytoKeys. 2015;51:1–282.
Carine MA, Russell SJ, Santos-Guerra A, Francisco-Ortega J. Relationships of the Macaronesian and Mediterranean floras: molecular evidence for multiple colonizations into Macaronesia and back-colonization of the continent in Convolvulus (Convolvulaceae). Am J Bot. 2004;91:1070–85.
Williams BRM, Mitchell TC, Wood JRI, Harris DJ, Scotland RW, Carine MA. Integrating DNA barcode data in a monographic study of Convolvulus L. Taxon. 2014. doi:10.12705/636.9.
Brown JM, Brummitt RK, Spencer M, Carine MA. Disentangling the bindweeds: hybridization and taxonomic diversity in British Calystegia (Convolvulaceae). Bot J Linn Soc. 2009;160:388–401.
Brummitt RK, Staples GW. Convolvulaceae. In: Heywood VH, Brummitt RK, Culham A, Seberg O, editors. Flowering plant families of the world. London: Royal Botanic Gardens, Kew; 2007. p. 108–10.
Brummitt RK. A remarkable new species of Calystegia (Convolvulaceae) from California. Kew Bull. 1974;29:499–502.
Särkinen T, Bohs L, Olmstead RG, Knapp S. A phylogenetic framework for evolutionary study of the nightshades (Solanaceae): a dated 1000-tip tree. BMC Evol Biol. 2013;13:214.
Stefanović S, Austin DF, Olmstead RG. Classification of Convolvulaceae: a phylogenetic approach. Syst Botany. 2003;28:791–806.
Escudero M, Valcárcel V, Vargas P, Luceño M. Significance of ecological vicariance and long-distance dispersal in the diversification of Carex sect. Spirostachyae (Cyperaceae). Am J Bot. 2009;96:2100–14.
Kadereit G, Mucina L, Freitag H. Phylogeny of Salicornioideae (Chenopodiaceae): diversification, biogeography, and evolutionary trends in leaf and flower morphology. Taxon. 2006;55:617–42.
Zhou Z, Wen J, Li G, Sun H. Phylogenetic assessment and biogeographic analyses of tribe Peracarpeae (Campanulaceae). Plant Syst Evol. 2012;298:323–36.
Kadereit G, Mavrodiev EV, Zacharias EH, Sukhorukov AP. Molecular phylogeny of Atripliceae (Chenopodioideae, Chenopodiaceae): implications for systematics, biogeography, flower and fruit evolution, and the origin of C4 photosynthesis. Am J Bot. 2010;97:1664–87.
Roberts EM, Stevens NJ, O’Connor PM, Dirks PHGM, Gottfried MD, Clyde WC, et al. Initiation of the western branch of the East African Rift coeval with the eastern branch. Nat Geosci. 2012;5:289–94.
Coleman M, Liston A, Kadereit JW, Abbott RJ. Repeat intercontinental dispersal and Pleistocene speciation in disjunct Mediterranean and desert Senecio (Asteraceae). Am J Bot. 2003;90:1446–54.
Galley C, Bytebier B, Bellstedt DU, Linder HP. The Cape element in the Afrotemperate flora: from Cape to Cairo? Proc R Soc B. 2007;274:535–43.
del Hoyo A, García-Marín JL, Pedrola-Monfort J. Temporal and spatial diversification of the African disjunct genus Androcymbium (Colchicaceae). Mol Phylogenet Evol. 2009;53:848–61.
Carlson SE, Linder HP, Donoghue MJ. The historical biogeography of Scabiosa (Dipsacaceae): implications for Old World plant disjunctions. J Biogeogr. 2012;39:1086–100.
Mao K, Hao G, Liu J, Adams RP, Milne RI. Diversification and biogeography of Juniperus (Cupressaceae): variable diversification rates and multiple intercontinental dispersals. New Phytol. 2010;188:254–72.
Ali SS, Yu Y, Pfosser M, Wetschnig W. Inferences of biogeographical histories within subfamily Hyacinthoideae using S-DIVA and Bayesian binary MCMC analysis implemented in RASP (Reconstruct Ancestral State in Phylogenies). Ann Bot. 2011;109:95–107.
Thiv M, van der Niet T, Rutschmann F, Thulin M, Brune T, Linder HP. Old–New World and trans-African disjunctions of Thamnosma (Rutaceae): Intercontinental long-distance dispersal and local differentiation in the succulent biome. Am J Bot. 2011;98:76–87.
Spalik K, Piwczyński M, Danderson CA, Kurzyna-Młynik R, Bone TS, Downie SR. Amphitropic amphiantarctic disjunctions in Apiaceae subfamily Apioideae. J Biogeogr. 2010;37:1977–94.
Emadzade K, Gehrke B, Linder HP, Hörandl E. The biogeographical history of the cosmopolitan genus Ranunculus L. (Ranunculaceae) in the temperate to meridional zones. Mol Phylogenet Evol. 2011;58:4–21.
Carlquist S. Intercontinental dispersal. In: Kubitzki K, editor. Dispersal and distribution. Hamburg: Parey; 1983. p. 37–47.
Wen J, Li PPL, Walck JL, Yoo K-O. Phylogenetic and biogeographic diversification in Osmorhiza (Apiaceae). Ann Missouri Bot Gard. 2002;89:414.
Proctor VW. Long-distance dispersal of seeds by retention in digestive tract of birds. Science. 1968;160:321–2.
Clausen P, Nolet BA, Fox AD, Klaassen M. Long-distance endozoochorous dispersal of submerged macrophyte seeds by migratory waterbirds in northern Europe—a critical review of possibilities and limitations. Acta Oecol. 2002;23:191–203.
Gregory-Wodzicki KM. Uplift history of the Central and Northern Andes: A review. Geol Soc Am Bull. 2000;112:1091–105.
Schaefer H, Hechenleitner P, Santos-Guerra A, de Sequeira MM, Pennington RT, Kenicer G, et al. Systematics, biogeography, and character evolution of the legume tribe Fabeae with special focus on the middle-Atlantic island lineages. BMC Evol Biol. 2012;12:250.
Winkworth RC, Hennion F, Prinzing A, Wagstaff SJ. Explaining the disjunct distributions of austral plants: the roles of Antarctic and direct dispersal routes. J Biogeogr. 2015;42(7):1197–209.
Abdel Khalik K, Osman AK. Seed morphology of some species of Convolvulaceae from Egypt (Identification of species and systematic significance). Feddes Repert. 2007;118:24–37.
Miryeganeh M, Takayama K, Tateishi Y, Kajita T. Long-distance dispersal by sea-drifted seeds has maintained the global distribution of Ipomoea pes-caprae subsp. brasiliensis (Convolvulaceae). PLoS ONE. 2014;9(4):e91836.
Gladenkov AY, Oleinik AE, Marincovich Jr L, Barinov KB. A refined age for the earliest opening of Bering Strait. Palaeogeogr Palaeoclimatol Palaeoecol. 2002;183:321–8.
Azuma H, García-Franco JG, Rico-Gray V, Thien LB. Molecular phylogeny of the Magnoliaceae: the biogeography of tropical and temperate disjunctions. Am J Bot. 2001;88:2275–85.
Davis CC, Fritsch PW, Li J, Donoghue MJ. Phylogeny and biogeography of Cercis (Fabaceae): evidence from nuclear ribosomal ITS and chloroplast ndhF sequence data. Syst Botany. 2002;27:289–302.
Milne RI. Northern Hemisphere plant disjunctions: a window on Tertiary land bridges and climate change? Ann Bot. 2006;98:465–72.
Mansion G, Zeltner L. Phylogenetic relationships within the New World endemic Zeltnera (Gentianaceae-Chironiinae) inferred from molecular and karyological data. Am J Bot. 2004;91:2069–86.
Guzmán B, Vargas P. Historical biogeography and character evolution of Cistaceae (Malvales) based on analysis of plastid rbcL and trnL-trnF sequences. Org Divers Evol. 2009;9:83–99.
Kadereit G, Freitag H. Molecular phylogeny of Camphorosmeae (Camphorosmoideae, Chenopodiaceae): Implications for biogeography, evolution of C4-photosynthesis and taxonomy. Taxon. 2011;60:51–78.
Magallón S, Sanderson MJ. Absolute diversification rates in Angiosperm clades. Evolution. 2001;55:1762–80.
Bacon CD, Baker WJ, Simmons MP. Miocene dispersal drives island radiations in the palm tribe Trachycarpeae (Arecaceae). Syst Biol. 2012;61:426–42.
MacGinitie HD. Fossil plants of the Florissant beds, Colorado. Carnegie Institution of Washington: Washington; 1953.
Martin HA. Re-assignment of the affinities of the fossil pollen type Tricolpites trioblatus Mildenhall and Pocknall to Wilsonia (Convolvulaceae) and a reassessment of the ecological interpretations. Rev Palaeobot Palynol. 2000;111:237–51.
Martin HA. The family Convolvulaceae in the Tertiary of Australia: evidence from pollen. Austral J Bot. 2001;49:221–34.
Wiens JJ, Moen DS. Missing data and the accuracy of Bayesian phylogenetics. J Syst Evol. 2008;46:307–14.
Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29:1969–73.
Rambaut A, Suchard MA, Xie D, Drummond AJ. Tracer v1.6. 2014. http://beast.bio.ed.ac.uk/Tracer.
Rambaut A. FigTree v1.4. 2012. http://tree.bio.ed.ac.uk/software/figtree.
Dice LR. Measures of the amount of ecologic association between species. Ecology. 1945;26:297–302.
Sørensen T. A method of establishing groups of equal amplitude in plant sociology based on similarity of species and its application to analyses of the vegetation on Danish commons. Kong Dansk Vidensk. 1948;5:1–34.
Garcia-Vallvé S, Puigbo P. DendroUPGMA: a dendrogram construction utility. http://genomes.urv.cat/UPGMA/.
Ree RH, Smith SA. Maximum Likelihood inference of geographic range evolution by dispersal, local extinction, and cladogenesis. Syst Biol. 2008;57:4–14.
Mao K, Milne RI, Zhang L, Peng Y, Liu J, Thomas P, et al. Distribution of living Cupressaceae reflects the breakup of Pangea. Proc Natl Acad Sci U S A. 2012;109:7793–8.
Rabosky DL. Automatic detection of key innovations, rate shifts, and diversity-dependence on phylogenetic trees. PLoS ONE. 2014;9:e89543.
Rabosky DL, Grundler M, Anderson C, Title P, Shi JJ, Brown JW, et al. BAMMtools: an R package for the analysis of evolutionary dynamics on phylogenetic trees. Methods Ecol Evol. 2014;5:701–7.
This work was supported by the German Research Foundation (DFG) and the Technische Universität München within the funding programme Open Access Publishing. This study was partially funded by a Syntax grant. We gratefully acknowledge the Natural History Museum (London), Royal Botanical Garden Edinburgh and Royal Botanical Gardens, Kew for providing accesses to herbarium specimens and tissue and the staff of the Molecular Biology Laboratories and the Wolfson-Wellcome Sequencing Facility at the Natural History Museum for technical support.
The authors declare that they have no competing interests.
MAC, RWS, JRIW and DJH conceived the ideas; TCM, BRMW and JRIW collected the data; TCM, BRMW and MAC analysed the data and led the writing. All authors read and approved the final manuscript.
The authors have broad interests in plant taxonomy, evolution and biogeography. This study formed part of a collaborative project to monograph Convolvulus.
Thomas C. Mitchell and Bethany R. M. Williams contributed equally to this work.
The data sets supporting the results of this article are available in the TreeBASE repository, study 18623; http://purl.org/phylo/treebase/phylows/study/TB2:S18623 and the following additional files.
Dated phylogeny of Convolvulaceae and Solanaceae inferred in beast from analysis of the concatenated mat K and rbcL dataset. Node bars represent 95 % HPD estimates. Bayesian Posterior Probabilities (BPP) ≥ 0.95 are given by their respective nodes. Scale bar represents millions of years before present. * indicate the location of fossil-calibrated nodes (Särkinen et al. ). + indicate the location of nodes used to calibrate the Convolvuleae phylogeny. (PDF 87 kb)
BAMM outputs. Phylorate plots (2a) and speciation rate through time curve for Convolvuleae (2b). Additional file 2a represents the distinct shift configurations that account for 95% of the probability of the data (f-values denote the posterior probability of each shift configuration). Branches are scale colour-coded to indicate rate variation from red (acceleration) to blue (deceleration). Circles indicate the location of core rate shifts and are similarly colour-coded, with circle size proportional to the marginal probability of a shift. Additional file 2b represents a speciation rate through time curve (red) for Convolvuleae. Blue shading represents the confidence on speciation rate at any point in time. (ZIP 474 kb)
List of accessions included in this study. GenBank numbers are provided in the respective columns. A dash indicates that no sequence data was included for that region. (DOCX 26 kb)
About this article
Cite this article
Mitchell, T.C., Williams, B.R.M., Wood, J.R.I. et al. How the temperate world was colonised by bindweeds: biogeography of the Convolvuleae (Convolvulaceae). BMC Evol Biol 16, 16 (2016). https://doi.org/10.1186/s12862-016-0591-6
- Diversification rates